Changeset 25509 in project

Timestamp:

11/16/11 12:10:14 (9 years ago)

Author:

Ivan Raikov

Message:

sundials: added M-L example with IDA solver

File:

• release/4/sundials/trunk/tests/ml.scm (modified) (4 diffs)

release/4/sundials/trunk/tests/ml.scm

@@ -66 +66 @@
 
 
-
-(define (main)
+(define (ida-main)
+
+  (define (ressc t yy yp)
+    (let ((dd (rhs t yy)))
+
+      (let ((v (- (f64vector-ref dd 0) (f64vector-ref yp 0)))
+            (w (- (f64vector-ref dd 1) (f64vector-ref yp 1))))
+        (f64vector v w))))
   
-  (let ((yy (f64vector -60.899 0.0149)) ;; v w
+  (let* ((yy (f64vector -60.899 0.0149)) ;; v w
+         (yp (rhs 0.0 yy))
+         
+         ;; Integration limits 
+         (t0  0.0)
+         (tf  TEND)
+         (dt  1e-2))
+    
+
+    ;; IDA initialization 
+    (let ((solver (ida-create-solver t0 tf yy yp ressc  
+                                     abstol: 1e-14
+                                     reltol: 1e-14)))
+
+      ;; In loop, call IDASolve, print results, and test for error. 
+      
+      (let recur ((tnext (+ t0 dt)) (iout 1))
+
+        (let ((flag  (ida-solve solver tnext)))
+          (if (negative? flag) (error 'main "IDA solver error" flag))
+
+          
+          (let ((yy (ida-yy solver)))
+            (printf "~A ~A ~A~%" 
+                    tnext
+                    (f64vector-ref yy 0)
+                    (f64vector-ref yy 1)
+                    ))
+
+
+          (if (< tnext tf)
+              (recur (+ tnext dt) (+ 1 iout)))
+          ))
+      
+      (ida-destroy-solver solver)
+      
+      )))
+
+
+
+(define (cvode-main)
+  
+  (let ((yy (f64vector -60.899 0.0149));; v w
 
         ;; Integration limits 
@@ -89 +137 @@
           (if (negative? flag) (error 'main "CVODE solver error" flag))
 
+
           (if (< tnext tf)
               (recur (+ tnext dt) (+ 1 iout)))
@@ -103 +152 @@
       )))
 
-(define (main/unsafe)
+(define (cvode-main/unsafe)
   
   (let ((yy (f64vector -60.899 0.0149)) ;; v w
@@ -151 +200 @@
             )))
       
-      
-(main)
-(main/unsafe)
+;(ida-main)      
+(cvode-main)
+(cvode-main/unsafe)