rb-daspk 0.0.3-x86-darwin-10 → 0.0.4-x86-darwin-10

data/examples/springmass.rb

@@ -0,0 +1,130 @@
+# Copyright (c) 2009 Eric Todd Meyers
+#
+# Permission is hereby granted, free of charge, to any person
+# obtaining a copy of this software and associated documentation
+# files (the "Software"), to deal in the Software without
+# restriction, including without limitation the rights to use,
+# copy, modify, merge, publish, distribute, sublicense, and/or sell
+# copies of the Software, and to permit persons to whom the
+# Software is furnished to do so, subject to the following
+# conditions:
+#
+# The above copyright notice and this permission notice shall be
+# included in all copies or substantial portions of the Software.
+#
+# THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
+# EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES
+# OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
+# NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT
+# HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY,
+# WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
+# FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR
+# OTHER DEALINGS IN THE SOFTWARE.
+#
+#
+# More information, source code and gems @ http://rubyforge.org/projects/rb-daspk/
+#
+
+require 'rubygems'
+require 'rb-daspk'
+include DASPK
+
+#
+#  classic spring mass problem
+#
+#     d2x/dt + 2*l*w0*dxdt + w0*w0*x = 0
+#
+#     let v = dx/dt, then
+#
+#  1. dxdt - v = 0  
+#  2. dvdt + 2*l*w0*v + w0*w0*x = 0
+#
+#
+#  let's start with x at 1, dxdt at 0 (v=0)
+#  we'll ask DASPK for the state derivatives
+#
+y     = [1.0,0.0]
+yprime= [0.0,0.0] # these are gusses in this case. we will call solveForInitialValues to get them
+
+
+w0 = 1.0
+l  = 0.0 # undamped
+
+#
+#  in case we want to see details during the integration. in this case were just dumping values
+#  not required (can pass nil to Solver.new if you want)
+#
+inter = Proc.new do |time,y,yprime|
+   print " time = #{time}  x = #{y[0]}\n"
+end
+
+jac = Proc.new do |time,y,yprime,jacm,cj|
+
+   jacm[0][0] = cj
+   jacm[0][1] = -1.0
+   jacm[1][0] = w0**2.0
+   jacm[1][1] = (2*l*w0 + cj)
+
+end
+
+s=Solver.new(2,y,yprime,inter,jac) do |time,y,yprime,delta|
+#s=Solver.new(2,y,yprime) do |time,y,yprime,delta|
+             delta[0]  = yprime[0] - y[1] 
+             delta[1]  = yprime[1] + 2.0*l*w0*y[1] + w0*w0*y[0] 
+end
+
+
+#
+#  solve the initial value problem
+#  note, this isnt required if the y and yprimes are already consistent
+# 
+#  return code from solveFoeInitialValues
+#    4 success (see DDASPK for all possible values)
+#
+idid = s.solveForInitialValues(0,y,yprime,[0,0])
+print " ----------------------------------------\n"
+print " initial conditions\n"
+print " x            = #{y[0]}\n"
+print " vel          = #{y[1]}\n"
+print " accelaration = #{yprime[1]}\n"
+print " ----------------------------------------\n"
+
+#
+# integrate to PI seconds
+#
+#
+ans = s.solve(0.0,Math::PI) 
+#
+# answer is an array of arrays
+# ans[0] = end time
+# ans[1] = end y array
+# ans[2] = end yprime array
+# ans[3] = DASPK return code (see DDASPK.F for meaning)
+#
+print " ----------------------------------------\n"
+print " final conditions\n"
+print " time         = #{ans[0]}\n"
+print " x            = #{ans[1][0]}\n"
+print " vel          = #{ans[1][1]}\n"
+print " accelaration = #{ans[2][1]}\n"
+print " ----------------------------------------\n"
+
+#
+# continue to integrate to 2PI seconds
+#
+#
+ans = s.solve(0.0,2*Math::PI) 
+#
+# answer is an array of arrays
+# ans[0] = end time
+# ans[1] = end y array
+# ans[2] = end yprime array
+# ans[3] = DASPK return code (see DDASPK.F for meaning)
+#
+print " ----------------------------------------\n"
+print " final conditions\n"
+print " time         = #{ans[0]}\n"
+print " x            = #{ans[1][0]}\n"
+print " vel          = #{ans[1][1]}\n"
+print " accelaration = #{ans[2][1]}\n"
+print " ----------------------------------------\n"

data/lib/rb-daspk.rb

@@ -24,7 +24,6 @@
 #
 # More information, source code and gems @ http://rubyforge.org/projects/rb-daspk/
 #
-#:include:test.rb
 
 require 'rubygems'
 require 'ffi'
@@ -74,7 +73,7 @@ module DASPK
 
     include DASPK
 
-    attr_accessor :rtol,:atol,:init,:jac,:maxstep,:mxnit,:mxnj,:mxnh,:lsoff,:stptol,:epinit
+    attr_accessor :rtol,:atol,:inter,:jac,:maxstep,:mxnit,:mxnj,:mxnh,:lsoff,:stptol,:epinit
 #
 #
 # * neq    = the number of equations. this is always the sum of the number of states and algebraics

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification 
 name: rb-daspk
 version: !ruby/object:Gem::Version 
-  version: 0.0.3
+  version: 0.0.4
 platform: x86-darwin-10
 authors: 
 - Eric Meyers
@@ -33,7 +33,7 @@ extra_rdoc_files: []
 files: 
 - lib/libdaspk.dylib
 - lib/rb-daspk.rb
-- lib/test.rb
+- examples/springmass.rb
 has_rdoc: true
 homepage: http://rb-daspk.rubyforge.org
 licenses: []

data/lib/test.rb

@@ -1,132 +0,0 @@
-# Copyright (c) 2009 Eric Todd Meyers
-#
-# Permission is hereby granted, free of charge, to any person
-# obtaining a copy of this software and associated documentation
-# files (the "Software"), to deal in the Software without
-# restriction, including without limitation the rights to use,
-# copy, modify, merge, publish, distribute, sublicense, and/or sell
-# copies of the Software, and to permit persons to whom the
-# Software is furnished to do so, subject to the following
-# conditions:
-#
-# The above copyright notice and this permission notice shall be
-# included in all copies or substantial portions of the Software.
-#
-# THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
-# EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES
-# OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
-# NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT
-# HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY,
-# WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
-# FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR
-# OTHER DEALINGS IN THE SOFTWARE.
-#
-#
-# More information, source code and gems @ http://rubyforge.org/projects/rb-daspk/
-#
-
-require 'rubygems'
-require 'rb-daspk'
-include DASPK
-
-y      = [1.0,0.5]
-yprime = [0.0,0.0]
-
-#
-#  in case the user wants to see details during the integration. in this case were just dumping values
-#  not required
-#
-inter = Proc.new do |time,y,yprime|
-    print "time = #{time}   y0=#{y[0]}  y1=#{y[1]}\n"
-end
-
-
-#
-#  let's solve and easy one
-#
-#     dy1/dy = y1
-#     dy2/dt = 2*y2
-#
-#  we all know the answer is (for the starting y's above)
-#     y1 = e**t
-#     y2 = e**(2t+ln(0.5))
-#
-#  at 2 seconds we should get y1 = 7.389056098930649, y2 = 27.29907501657209
-#
-#
-
-#
-#  explicit modified jacobian
-#  not required, but used if passed to solver
-#
-jac = Proc.new do |time,y,yprime,pd,cj|
-   pd[0][0] = 1.0 - cj   # first equation with respect to first var
-   pd[0][1] = 0.0        # first equation with respect to second var
-   pd[1][0] = 0.0        # second equation with respect to first var
-   pd[1][1] = 2.0 - cj   # second equation with respect to second var
-end
-
-s=Solver.new(2,y,yprime,inter, jac) do |time,yy,yyprime,delta| 
-#s=Solver.new(2,y,yprime,inter) do |time,yy,yyprime,delta|
-             delta[0]  =     yy[0] - yyprime[0] 
-             delta[1]  = 2.0*yy[1] - yyprime[1] 
-end
-
-
-#
-#  solve the initial value problem
-#  this is only required if the Y's and Y'primes are not
-#  already consistent
-#
-#  y and yprime will be modified
-#
-#
-#
-#
-#
-#
-ans = s.solveForInitialValues(0,y,yprime,[0,0],1)
-print " y0=#{y[0]}  y1=#{y[1]}\n"
-print " yp0=#{yprime[0]}  yp1=#{yprime[1]}\n"
-print "idid = #{ans}\n"
-
-#
-# integrate to 1 second
-#
-ans = s.solve(0.0,1.0) 
-
-
-#print final values
-time = ans[0]
-ya = ans[1]
-ypa = ans[2]
-print "time = #{time}   y0=#{ya[0]}  y1=#{ya[1]}\n"
-print "time = #{time}   yp0=#{ypa[0]}  yp1=#{ypa[1]}\n"
-
-#
-# lets keep going. should start at the last time we left off regardless of the sart time
-# integrate to 2 seconds
-#
-ans = s.solve(0.0,2.0) 
-#print final values
-time = ans[0]
-ya = ans[1]
-ypa = ans[2]
-print "time = #{time}   y0=#{ya[0]}  y1=#{ya[1]}\n"
-print "time = #{time}   yp0=#{ypa[0]}  yp1=#{ypa[1]}\n"
-
-#
-# reset integration. should go back to tzero
-# integrate to 2 seconds
-#
-s.reset
-
-ans = s.solve(0.0,2.0) 
-#print final values
-time = ans[0]
-ypa = ans[2]
-ya = ans[1]
-print "time = #{time}   y0=#{ya[0]}  y1=#{ya[1]}\n"
-print "time = #{time}   yp0=#{ypa[0]}  yp1=#{ypa[1]}\n"
-
-