* ext/bigdecimal/sample/linear.rb: resolve LoadError.
* ext/bigdecimal/sample/nlsolve.rb: ditto. * ext/bigdecimal/lib/bigdecimal/nlsolve.rb: removed because this file is sample script and same file exists in ext/bigdecimal/sample. git-svn-id: svn+ssh://ci.ruby-lang.org/ruby/trunk@7996 b2dd03c8-39d4-4d8f-98ff-823fe69b080e
This commit is contained in:
parent
64fbfe54fb
commit
c1ed63a4c4
10
ChangeLog
10
ChangeLog
@ -1,8 +1,14 @@
|
|||||||
Sat Feb 19 01:28:56 2005 Hirokazu Yamamoto <ocean@m2.ccsnet.ne.jp>
|
Sat Feb 19 01:28:56 2005 Hirokazu Yamamoto <ocean@m2.ccsnet.ne.jp>
|
||||||
|
|
||||||
* bigdecimal/newton.rb: resolve LoadError. [ruby-dev:25685]
|
* ext/bigdecimal/lib/bigdecimal/newton.rb: resolve LoadError.
|
||||||
|
[ruby-dev:25685]
|
||||||
|
|
||||||
* bigdecimal/nlsolve.rb: ditto.
|
* ext/bigdecimal/sample/linear.rb: ditto.
|
||||||
|
|
||||||
|
* ext/bigdecimal/sample/nlsolve.rb: ditto.
|
||||||
|
|
||||||
|
* ext/bigdecimal/lib/bigdecimal/nlsolve.rb: removed because this file
|
||||||
|
is sample script and same file exists in ext/bigdecimal/sample.
|
||||||
|
|
||||||
Thu Feb 17 22:15:34 2005 K.Kosako <sndgk@ybb.ne.jp>
|
Thu Feb 17 22:15:34 2005 K.Kosako <sndgk@ybb.ne.jp>
|
||||||
|
|
||||||
|
@ -1,38 +0,0 @@
|
|||||||
#!/usr/local/bin/ruby
|
|
||||||
|
|
||||||
#
|
|
||||||
# nlsolve.rb
|
|
||||||
# An example for solving nonlinear algebraic equation system.
|
|
||||||
#
|
|
||||||
|
|
||||||
require "bigdecimal"
|
|
||||||
require "bigdecimal/newton"
|
|
||||||
include Newton
|
|
||||||
|
|
||||||
class Function
|
|
||||||
def initialize()
|
|
||||||
@zero = BigDecimal::new("0.0")
|
|
||||||
@one = BigDecimal::new("1.0")
|
|
||||||
@two = BigDecimal::new("2.0")
|
|
||||||
@ten = BigDecimal::new("10.0")
|
|
||||||
@eps = BigDecimal::new("1.0e-16")
|
|
||||||
end
|
|
||||||
def zero;@zero;end
|
|
||||||
def one ;@one ;end
|
|
||||||
def two ;@two ;end
|
|
||||||
def ten ;@ten ;end
|
|
||||||
def eps ;@eps ;end
|
|
||||||
def values(x) # <= defines functions solved
|
|
||||||
f = []
|
|
||||||
f1 = x[0]*x[0] + x[1]*x[1] - @two # f1 = x**2 + y**2 - 2 => 0
|
|
||||||
f2 = x[0] - x[1] # f2 = x - y => 0
|
|
||||||
f <<= f1
|
|
||||||
f <<= f2
|
|
||||||
f
|
|
||||||
end
|
|
||||||
end
|
|
||||||
f = BigDecimal::limit(100)
|
|
||||||
f = Function.new
|
|
||||||
x = [f.zero,f.zero] # Initial values
|
|
||||||
n = nlsolve(f,x)
|
|
||||||
p x
|
|
@ -11,7 +11,7 @@
|
|||||||
#
|
#
|
||||||
|
|
||||||
require "bigdecimal"
|
require "bigdecimal"
|
||||||
require "ludcmp"
|
require "bigdecimal/ludcmp"
|
||||||
|
|
||||||
#
|
#
|
||||||
# NOTE:
|
# NOTE:
|
||||||
|
@ -6,7 +6,7 @@
|
|||||||
#
|
#
|
||||||
|
|
||||||
require "bigdecimal"
|
require "bigdecimal"
|
||||||
require "newton"
|
require "bigdecimal/newton"
|
||||||
include Newton
|
include Newton
|
||||||
|
|
||||||
class Function
|
class Function
|
||||||
|
Loading…
x
Reference in New Issue
Block a user