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/util.rb"

• +

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

• -

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

• *

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

• /

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

• add

以下のように使用します。
c = a.add(b,n)
c = a + b を最大で n 桁まで計算します。 a + b の精度が n より大きいときは丸められます。

• sub

以下のように使用します。
c = a.sub(b,n)
c = a - b を最大で n 桁まで計算します。 a - b の精度が n より大きいときは丸められます。

• mult

以下のように使用します。
c = a.mult(b,n)
c = a * b を最大で n 桁まで計算します。 a * b の精度が n より大きいときは丸められます。

• div

以下のように使用します。
c,r = a.div(b,n)
c=a/b の計算をします。 r には剰余が代入されます。a/bは 必要ならn 桁まで計算されます。divmod メソッド と異なり、c は整数とは限りません。 また、 c は丸められることはありません。 a = c*b + r の関係は成立します。

• %

r = a%b
a/b の余りを計算します。以下の計算と同じものです。
r = a-((a/b).floor)*b


• 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::new("1.23456").floor(4)   #  ==> 1.2345
 c = BigDecimal::new("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::new("1.23456").ceil(4)   # ==> 1.2346
 c = BigDecimal::new("15.23456").ceil(-1) # ==> 20.0

• round[(n[,b])]

c = a.round
小数点以下第一位の数を四捨五入して整数（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::new("1.23456").round(4)   #  ==> 1.2346
c = BigDecimal::new("15.23456").round(-1) #  ==> 20.0

２番目の引数（デフォルトは 0）にゼロ以外を指定すると、いわゆる Banker's rounding になります。
Banker's rounding とは、四捨五入する数字を p として、p < 5 なら切り捨て p > 5 なら切り上げ、 p がちょうど５のときだけは切り上げ先の数字＋１が偶数になるときだけ切り上げます。

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

• truncate

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


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

• divmod

商と剰余の配列を返します。
c,r = a.divmod(b) # a = c*b + r
divmodメソッドは a = c * b + r となる a / b の浮動小数点型の商 c と剰余 r を 計算します。ここで c は整数（少数部分のない実数）になります。
c = (a/b).floor
r = a - c*b
で計算されます。

• remainder

r=a.remainder(b)
a/b の剰余 r を計算します。
c = (a/b).fix
r = a - c*b
で計算されます。

• 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"の形になります）。
s = a.to_s
n が指定されたときは、仮数部分を n 桁毎に空白で区切ります。
s = a.to_s(n)

• exponent

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

• prec

n,m = a.prec
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)


• <=>

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

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

(評価段階の)クラスメソッド

• E

自然対数の底e(=2.718281828....)を計算します（正直にテイラー展開で）。
e = BigDecimal::E(n)
nは必要な有効桁数を整数で指定します。

• PI

円周率(=3.14159265358979....)を計算します（J．Machinの公式を用います）。
e = BigDecimal::PI(n)
n は必要な有効桁数を整数で指定します。

(評価段階の)インスタンスメソッド

• sincos

a の有効桁 n 桁の sin と cos を同時に（テイラー展開で）計算して、 sin と cos の配列を返します。 n は必要な有効桁数です（ n の sin や cos を計算するわけではありません）。
sin,cos = a.sincos(n)
|a| < 2*3.1415....でないと正しい答えを計算できないこともあります。

• exp

自然対数の底e(=2.718281828....)の a 乗を計算します。
c = a.exp(n)
n は必要な有効桁数です。

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

#
# pi.rb
#  USAGE: ruby pi.rb n
#   where n is the number of digits required.
#  EX.: ruby pi.rb 1000
#

require "bigdecimal"
#
# Calculates 3.1415.... 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.exponent >= exp) 
    t   = t*m25
    u,r = 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.exponent >= exp )
    t,r = t.div(m57121,sig)
    u,r = t.div(k,sig)
    pi  = pi + u
    k   = k+two
  end
  pi
end

if $0 == __FILE__
  print "PI("+ARGV[0]+"):\n"
  p pi(ARGV[0].to_i)
end