Viewing file:     bigdecimal_ja.html (28.92 KB)      -rw-r--r--

BigDecimal(可変長浮動少数点演算用拡張ライブラリ)

• http://www.ruby-lang.org/ja/：Ruby公式ページ
• http://kahori.com/ruby/ring/：Rubyに関するページを辿れます

目次

• はじめに
• 使用方法とメソッドの一覧
• 無限、非数、ゼロの扱い
• 内部構造
• 2進と10進
• 計算精度について

インストールについて

require 'bigdecimal'
a=BigDecimal::new("0.123456789123456789")
b=BigDecimal("123456.78912345678",40)
c=a+b

require "bigdecimal/math.rb"

require "bigdecimal/util.rb"

• +

加算（c = a + b）
c の精度については「計算精度について」を参照してください。

• -

減算（c = a - b）、または符号反転（c = -a）
c の精度については「計算精度について」を参照してください。

• *

乗算(c = a * b)
cの精度は(aの精度)+(bの精度)程度です。
詳しくは「計算精度について」を参照してください。

• /

除算(c = a / b)
c の精度については「計算精度について」を参照してください。

• add(b,n)

以下のように使用します。
c = a.add(b,n)
c = a + b を最大で n 桁まで計算します。
a + b の精度が n より大きいときは BigDecimal.mode で指定された方法で丸められます。
n がゼロなら + と同じです。

• sub(b,n)

以下のように使用します。
c = a.sub(b,n)
c = a - b を最大で n 桁まで計算します。
a - b の精度が n より大きいときは BigDecimal.mode で指定された方法で丸められます。
n がゼロなら - と同じです。

• mult(b,n)

以下のように使用します。
c = a.mult(b,n)
c = a * b を最大で n 桁まで計算します。
a * b の精度が n より大きいときは BigDecimal.mode で指定された方法で丸められます。
n がゼロなら * と同じです。

• div(b[,n])

以下のように使用します。
c = a.div(b,n)
c = a / b を最大で n 桁まで計算します。 a / b の精度が n より大きいときは BigDecimal.mode で指定された方法で丸められます。
n がゼロなら / と同じです。
n が省略されたときは Float#div と同様に結果が整数(BigDecimal)になります。

• fix

a の小数点以下の切り捨て。
c = a.fix

• frac

a の整数部分の切り捨て。
c = a.frac

• floor[(n)]

c = a.floor
a 以下の最大整数（BigDecimal 値）を返します。

c = BigDecimal("1.23456").floor  #  ==> 1
c = BigDecimal("-1.23456").floor #  ==> -2

以下のように引数 n を与えることもできます。
n>=0 なら、小数点以下 n+1 位の数字を操作します(少数点以下を、最大 n 桁にします)。
n が負のときは小数点以上 n 桁目を操作します(小数点位置から左に少なくとも n 個の 0 が並びます)。


 c = BigDecimal("1.23456").floor(4)   #  ==> 1.2345
 c = BigDecimal("15.23456").floor(-1) #  ==> 10.0

• ceil[(n)]

c = a.ceil
a 以上の整数のうち、最も小さい整数を計算し、その値（BigDecimal 値）を返します。

c = BigDecimal("1.23456").ceil  #  ==> 2
c = BigDecimal("-1.23456").ceil #  ==> -1

以下のように引数を与えて、小数点以下 n+1 位の数字を操作することもできます。
n>=0 なら、小数点以下 n+1 位の数字を操作します(少数点以下を、最大 n 桁にします)。
n が負のときは小数点以上 n 桁目をを操作します(小数点位置から左に少なくとも n 個の 0 が並びます)。


 c = BigDecimal("1.23456").ceil(4)   # ==> 1.2346
 c = BigDecimal("15.23456").ceil(-1) # ==> 20.0

• round[(n[,b])]

c = a.round
クラスメソッド BigDecimal::mode(BigDecimal::ROUND_MODE,flag) で指定した ROUND_MODE に従って丸め操作を実行します。 BigDecimal::mode(BigDecimal::ROUND_MODE,flag) で何も指定せず、かつ、引数 を指定しない場合は「小数点以下第一位の数を四捨五入して整数（BigDecimal 値）」にします。


 c = BigDecimal("1.23456").round  #  ==> 1
 c = BigDecimal("-1.23456").round #  ==> -1

以下のように引数を与えて、小数点以下 n+1 位の数字を操作することもできます。
n が正の時は、小数点以下 n+1 位の数字を丸めます(少数点以下を、最大 n 桁にします)。
n が負のときは小数点以上 n 桁目を丸めます(小数点位置から左に少なくとも n 個の 0 が並びます)。

c = BigDecimal("1.23456").round(4)   #  ==> 1.2346
c = BigDecimal("15.23456").round(-1) #  ==> 20.0

２番目の引数を指定すると、BigDecimal#mode の指定を無視して、指定された方法で 丸め操作を実行します。

c = BigDecimal("1.23456").round(3,BigDecimal::ROUND_HALF_EVEN)   #  ==> 1.234
c = BigDecimal("1.23356").round(3,BigDecimal::ROUND_HALF_EVEN)   #  ==> 1.234

• truncate

c = a.truncate
小数点以下の数を切り捨てて整数（BigDecimal 値）にします。
以下のように引数を与えて、小数点以下 n+1 位の数字を操作することもできます。
n が正の時は、小数点以下 n+1 位の数字を切り捨てます(少数点以下を、最大 n 桁にします)。 n が負のときは小数点以上 n 桁目をを操作します(小数点位置から左に少なくとも n 個の 0 が並びます)。


c = BigDecimal("1.23456").truncate(4)   #  ==> 1.2345
c = BigDecimal("15.23456").truncate(-1) #  ==> 10.0

• abs

ａの絶対値
c = a.abs


• to_i

少数点以下を切り捨てて整数に変換します。
i = a.to_i
i は値に応じて Fixnum か Bignum になります。 a が Infinity や NaN のとき、i は nil になります。

• to_f

Float オブジェクトに変換します。 よりきめ細かい値が必要ならば split メソッドを利用して ください。

• to_s[(n)]

文字列に変換します(デフォルトは "0.xxxxxEn" の形になります）。

BigDecimal("1.23456").to_s  #  ==> "0.123456E1"

引数 n に正の整数が指定されたときは、少数点で分けられる左右部分を、それぞれ n 桁毎 に空白で区切ります。

BigDecimal("0.1234567890123456789").to_s(10)   #  ==> "0.1234567890 123456789E0"

引数 n に正の整数を表す文字列を指定することもできます。

BigDecimal("0.1234567890123456789").to_s("10") #  ==> "0.1234567890 123456789E0"

文字列の最初に '+'（または ' '）を付けると、値が正の場合、先頭に '+'（または ' '）が付きます （負の場合は、常に '-' が付きます。）。

BigDecimal("0.1234567890123456789").to_s(" 10") #  ==> " 0.1234567890 123456789E0"
BigDecimal("0.1234567890123456789").to_s("+10") #  ==> "+0.1234567890 123456789E0"
BigDecimal("-0.1234567890123456789").to_s("10") #  ==> "-0.1234567890 123456789E0"

さらに文字列の最後に E(または e) か F(または f) を指定することで、以下のように 表示形式を変更することができます。

BigDecimal("1234567890.123456789").to_s("E")  #  ==> "0.1234567890123456789E10"
BigDecimal("1234567890.123456789").to_s("F")  #  ==> "1234567890.123456789"
BigDecimal("1234567890.123456789").to_s("5E") #  ==> "0.12345 67890 12345 6789E10"
BigDecimal("1234567890.123456789").to_s("5F") #  ==> "12345 67890.12345 6789"

• exponent

指数部を整数値で返します。 n = a.exponent
は a の値が 0.xxxxxxx*10**n を意味します。

• precs

n,m = a.precs
a の有効数字 (n) と最大有効数字 (m) の配列を返します。

• sign

値が正(sign > 0)、負(sign < 0)、その他(sigh==0)であるかの情報を返します。 n = a.sign
としたとき n の値は a が以下のときを意味します。
() の中の数字は、実際の値です(「内部構造」を参照)。
n = BigDecimal::SIGN_NaN(0) : a は NaN
n = BigDecimal::SIGN_POSITIVE_ZERO(1) : a は +0
n = BigDecimal::SIGN_NEGATIVE_ZERO(-1) : a は -0
n = BigDecimal::SIGN_POSITIVE_FINITE(2) : a は正の値
n = BigDecimal::SIGN_NEGATIVE_FINITE(-2) : a は負の値
n = BigDecimal::SIGN_POSITIVE_INFINITE(3) : a は+Infinity
n = BigDecimal::SIGN_NEGATIVE_INFINITE(-3) : a は-Infinity


• nan?

a.nan? は a がNaNのとき真を返します。

• infinite?

a.infinite? は a が+∞のとき 1 、-∞のときは -1、それ以外のときは nil を返します。

• finite?

a.finite? は a が∞または NaN でないとき真を返します。

• zero?

a が 0 なら true になります。
c = a.zero?

• nonzero?

a が 0 なら nil、0 以外なら a そのものが返ります。
c = a.nonzero?

• split

BigDecimal 値を 0.xxxxxxx*10**n と表現したときに、符号（NaNのときは 0、それ以外は+1か-1になります）、 仮数部分の文字列（"xxxxxxx"）と、基数（10）、更に指数 n を配列で 返します。
a=BigDecimal::new("3.14159265")
f,x,y,z = a.split
とすると、f=+1、x="314159265"、y=10、z=1になります。
従って、
s = "0."+x
b = f*(s.to_f)*(y**z)
で Float に変換することができます。

• inspect

デバッグ出力に使用されます。
p a=BigDecimal::new("3.14",10)
とすると、[0x112344:'0.314E1',4(12)]のように出力されます。 最初の16進数はオブジェクトのアドレス、次の '0.314E1' は値、 次の4は現在の有効桁数(表示より若干大きいことがあります)、 最後はオブジェクトが取り得る最大桁数になります。

• **

a の n 乗を計算します。ｎは整数。
c = a ** n
結果として c の有効桁は a の n 倍以上になるので注意。

• power

** と同じで、a の n 乗を計算します。ｎは整数。
c = a.power(n)
結果として c の有効桁は a の n 倍以上になるので注意。

• sqrt

aの有効桁 n 桁の平方根（n の平方根ではありません）を ニュートン法で計算します。
c = a.sqrt(n)


• divmod,quo,modulo,%,remainder

詳細は対応する Float の各メソッドを参照して下さい。

• <=>

a==b なら 0、a > b なら 1、a < b なら -1 になります。
c = a <=> b

• ==
• ===
「==」と同じですが case 文で使用されます。
• !=
• <
• <=
• >
• >=

coerceについて

   a = BigDecimal.E(20)
   c = a * "0.123456789123456789123456789" # 文字を BigDecimal に変換してから計算

  a = BigDecimal.E(20)
  c = "0.123456789123456789123456789" * a # エラー

require "bigdecimal"

aa  = %w(1 -1 +0.0 -0.0 +Infinity -Infinity NaN)
ba  = %w(1 -1 +0.0 -0.0 +Infinity -Infinity NaN)
opa = %w(+ - * / <=> > >=  < == != <=)

for a in aa
  for b in ba
    for op in opa
      x = BigDecimal::new(a)
      y = BigDecimal::new(b)
      eval("ans= x #{op} y;print a,' ',op,' ',b,' ==> ',ans.to_s,\"\n\"")
    end
  end
end

  typedef struct {
     unsigned long MaxPrec; // 最大精度(frac[]の配列サイズ)
     unsigned long Prec;    // 精度(frac[]の使用サイズ)
     short    sign;         // 以下のように符号等の状態を定義します。
                            //  ==0 : NaN
                            //    1 : +0
                            //   -1 : -0
                            //    2 : 正の値
                            //   -2 : 負の値
                            //    3 : +Infinity
                            //   -3 : -Infinity
     unsigned short flag;   // 各種の制御フラッグ
     int      exponent;     // 指数部の値(仮数部*BASE**exponent)
     unsigned long frac[1]; // 仮数部の配列(可変)
  } Real;

    0.1234 5678 4321*(10000)**1

#!/usr/local/bin/ruby

require "bigdecimal"
#
# Calculates 3.1415.... (the number of times that a circle's diameter
# will fit around the circle) using J. Machin's formula.
#
def big_pi(sig) # sig: Number of significant figures
  exp    = -sig
  pi     = BigDecimal::new("0")
  two    = BigDecimal::new("2")
  m25    = BigDecimal::new("-0.04")
  m57121 = BigDecimal::new("-57121")

  u = BigDecimal::new("1")
  k = BigDecimal::new("1")
  w = BigDecimal::new("1")
  t = BigDecimal::new("-80")
  while (u.nonzero? && u.exponent >= exp) 
    t   = t*m25
    u   = t.div(k,sig)
    pi  = pi + u
    k   = k+two
  end

  u = BigDecimal::new("1")
  k = BigDecimal::new("1")
  w = BigDecimal::new("1")
  t = BigDecimal::new("956")
  while (u.nonzero? && u.exponent >= exp )
    t   = t.div(m57121,sig)
    u   = t.div(k,sig)
    pi  = pi + u
    k   = k+two
  end
  pi
end

if $0 == __FILE__
  if ARGV.size == 1
    print "PI("+ARGV[0]+"):\n"
    p big_pi(ARGV[0].to_i)
  else
    print "TRY: ruby pi.rb 1000 \n"
  end
end