Ошибка в bigdecimal / newton для реализации XIRR - PullRequest
2 голосов
/ 30 июля 2011

Я хотел бы использовать модуль bigdecimal/newton в Ruby для реализации XIRR .Я написал скрипт, чтобы попробовать его, следуя этому примеру .

Когда я запускаю код (Ruby 1.9.2 в Mac OS X 10.6), я получаю следующую ошибку:

/opt/local/lib/ruby1.9/1.9.1/bigdecimal/ludcmp.rb:84:in `div': wrong number of arguments(2 for 1) (ArgumentError)
   from /opt/local/lib/ruby1.9/1.9.1/bigdecimal/ludcmp.rb:84:in `block in lusolve'
   from /opt/local/lib/ruby1.9/1.9.1/bigdecimal/ludcmp.rb:78:in `downto'
   from /opt/local/lib/ruby1.9/1.9.1/bigdecimal/ludcmp.rb:78:in `lusolve'
   from /opt/local/lib/ruby1.9/1.9.1/bigdecimal/newton.rb:60:in `nlsolve'
   from gistfile1.rb:52:in `<main>'

Это из-за ошибки в моей программе или из-за чего-то в библиотеке BigDecimal?

#!/opt/local/bin/ruby1.9
require 'bigdecimal'
require 'bigdecimal/newton'
require 'time'
include Newton

transactions = [
  {amount:BigDecimal::new("-1000"), date:Time.parse("2011-01-01")},
  {amount:BigDecimal::new("100"), date:Time.parse("2011-02-11")},
  {amount:BigDecimal::new("100"), date:Time.parse("2011-03-11")},
  {amount:BigDecimal::new("100"), date:Time.parse("2011-04-11")},
  {amount:BigDecimal::new("100"), date:Time.parse("2011-05-11")},
  {amount:BigDecimal::new("100"), date:Time.parse("2011-06-11")},
  {amount:BigDecimal::new("100"), date:Time.parse("2011-07-11")},
  {amount:BigDecimal::new("100"), date:Time.parse("2011-08-11")},
  {amount:BigDecimal::new("100"), date:Time.parse("2011-09-11")},
  {amount:BigDecimal::new("100"), date:Time.parse("2011-10-11")},
  {amount:BigDecimal::new("100"), date:Time.parse("2011-11-11")},
  {amount:BigDecimal::new("100"), date:Time.parse("2011-12-11")},
  {amount:BigDecimal::new("100"), date:Time.parse("2012-01-11")},
]

class Function
  values ={ eps: "1.0e-16", one: "1.0", two: "2.0", ten: "10.0", zero: "0.0" }

  values.each do |key, value|
    define_method key do
      BigDecimal::new value
    end
  end

  def initialize(transactions)
    @transactions = transactions
  end

  def values(x)
    start = @transactions[0][:date]
    value = []

    value << @transactions.reduce(0) do |sum, t|
      pwr = (t[:date] - start) / 365
      r = t[:amount] / (1.0 + x[0]) ** pwr
      sum + r
    end

    value
  end
end

f = Function.new(transactions)
x = [f.zero]
n = nlsolve(f,x)
puts x

Любая помощь при отладке будет принята с благодарностью.

Спасибо!

1 Ответ

0 голосов
/ 30 июля 2011

Проблема в том, что вы смешиваете BigDecimals с "нормальными" числами в ваших вычислениях.Попробуйте заменить values на это:

values = { 
  eps: BigDecimal.new("1.0e-16"), 
  one: BigDecimal.new("1.0"), 
  two: BigDecimal.new("2.0"), 
  ten: BigDecimal.new("10.0"), 
  zero: BigDecimal.new("0.0") 
}

и блок в def values(x) на

value << @transactions.reduce(0) do |sum, t|
  pwr = (t[:date] - start) / BigDecimal.new("365.0")
  r = t[:amount] / (BigDecimal.new("1.0") + x[0]) ** pwr
  sum + r
end
...