Changeset 25539 in project


Ignore:
Timestamp:
11/21/11 18:20:42 (9 years ago)
Author:
Ivan Raikov
Message:

sundials: fixes to unit tests

Location:
release/4/sundials/trunk
Files:
5 edited

Legend:

Unmodified
Added
Removed
  • release/4/sundials/trunk/sundials.scm

    r25538 r25539  
    373373    }
    374374
    375     flag = IDASetStopTime(solver_handle->ida_mem, time_stop);
     375    if (time_stop > 0.0)
     376    {
     377      flag = IDASetStopTime(solver_handle->ida_mem, time_stop);
     378    }
    376379    if (flag != IDA_SUCCESS)
    377380         chicken_error("could not set stop time with IDASetStopTime", C_fix(flag)) ;   
     
    699702
    700703(define (ida-create-solver
    701          tstart tstop variables derivatives
     704         tstart variables derivatives
    702705         residual-main #!key
     706         (tstop 0.0)
    703707         (residual-init #f)
    704708         (residual-event #f)
     
    798802
    799803(define (ida-create-solver/unsafe
    800          tstart tstop variables derivatives
     804         tstart variables derivatives
    801805         residual-main #!key
     806         (tstop 0.0)
    802807         (residual-init #f)
    803808         (residual-event #f)
     
    936941
    937942(define (cvode-create-solver
    938          tstart tstop variables 
     943         tstart variables 
    939944         rhs-fn #!key
     945         (tstop 0.0)
    940946         (lmm cvode-lmm/adams)
    941947         (iter cvode-iter/functional)
     
    10141020
    10151021(define (cvode-create-solver/unsafe
    1016          tstart tstop variables 
     1022         tstart variables 
    10171023         rhs-fn #!key
     1024         (tstop 0.0)
    10181025         (lmm cvode-lmm/adams)
    10191026         (iter cvode-iter/functional)
  • release/4/sundials/trunk/tests/e.scm

    r25082 r25539  
    3333   
    3434    ;; IDA initialization
    35     (let ((solver (ida-create-solver t0 tf yy yp ressc 
     35    (let ((solver (ida-create-solver t0 yy yp ressc 
     36                                     tstop: tf
    3637                                     abstol: 1e-14
    3738                                     reltol: 1e-14)))
     
    6768   
    6869    ;; IDA initialization
    69     (let ((solver (ida-create-solver/unsafe t0 tf yy yp ressc/unsafe
     70    (let ((solver (ida-create-solver/unsafe t0 yy yp ressc/unsafe
     71                                            tstop: tf
    7072                                            abstol: 1e-14
    7173                                            reltol: 1e-14)))
  • release/4/sundials/trunk/tests/hh.scm

    r25082 r25539  
    133133    ;; CVODE initialization
    134134    (let ((solver (cvode-create-solver
    135                    t0 tf yy rhs 
     135                   t0 yy rhs 
     136                   tstop: tf
    136137                   abstol: 1e-4
    137138                   reltol: 1e-4)))
     
    143144        (let ((flag  (cvode-solve solver tnext)))
    144145          (if (negative? flag) (error 'main "CVODE solver error" flag))
     146
     147          ;;(print-results solver tnext)
    145148
    146149          (if (< tnext tf)
     
    154157              (h (f64vector-ref yy 2))
    155158              (n (f64vector-ref yy 3)))
    156           (assert (< (abs (- v 12.8334007736321)) 1e-12))
    157           (assert (< (abs (- m 0.988277249273334)) 1e-12))
    158           (assert (< (abs (- h 0.151908759594835)) 1e-12))
    159           (assert (< (abs (- n 0.681407916363)) 1e-12))
     159          (assert (< (abs (- v 13.1818250584672)) 1e-12))
     160          (assert (< (abs (- m 0.988216867260496)) 1e-12))
     161          (assert (< (abs (- h 0.152871810132604)) 1e-12))
     162          (assert (< (abs (- n 0.680166108403409)) 1e-12))
    160163          ))
    161164     
     
    175178    ;; CVODE initialization
    176179    (let ((solver (cvode-create-solver/unsafe
    177                    t0 tf yy rhs/unsafe 
    178                    abstol: 1e-4
     180                   t0 yy rhs/unsafe 
     181                   tstop: tf
     182                   abstol: 1e-4
    179183                   reltol: 1e-4)))
    180184
     
    196200              (h (f64vector-ref yy 2))
    197201              (n (f64vector-ref yy 3)))
    198           (assert (< (abs (- v 12.8334007736321)) 1e-12))
    199           (assert (< (abs (- m 0.988277249273334)) 1e-12))
    200           (assert (< (abs (- h 0.151908759594835)) 1e-12))
    201           (assert (< (abs (- n 0.681407916363)) 1e-12))
     202          (assert (< (abs (- v 13.1818250584672)) 1e-12))
     203          (assert (< (abs (- m 0.988216867260496)) 1e-12))
     204          (assert (< (abs (- h 0.152871810132604)) 1e-12))
     205          (assert (< (abs (- n 0.680166108403409)) 1e-12))
    202206          ))
    203207     
  • release/4/sundials/trunk/tests/iaf.scm

    r25245 r25539  
    4343    (let ((dt 1e-2)
    4444          (solver (cvode-create-solver
    45                    t0 tf yy subthreshold 
     45                   t0 yy subthreshold 
     46                   tstop: tf
    4647                   abstol: abstol
    4748                   reltol: reltol
     
    6364              (cvode-reinit-solver solver (+ tnext trefractory) (f64vector vreset)) )
    6465         
    65           (print-results solver tnext)
     66          ;;(print-results solver tnext)
    6667
    6768          (if (< tnext tf)
  • release/4/sundials/trunk/tests/ml.scm

    r25509 r25539  
    8585
    8686    ;; IDA initialization
    87     (let ((solver (ida-create-solver t0 tf yy yp ressc 
     87    (let ((solver (ida-create-solver t0 yy yp ressc 
     88                                     tstop: tf
    8889                                     abstol: 1e-14
    8990                                     reltol: 1e-14)))
     
    9697          (if (negative? flag) (error 'main "IDA solver error" flag))
    9798
     99          ;;(print-results solver tnext)
    98100         
    99101          (let ((yy (ida-yy solver)))
     
    126128    ;; CVODE initialization
    127129    (let ((solver (cvode-create-solver
    128                    t0 tf yy rhs
     130                   t0 yy rhs
     131                   tstop: tf
    129132                   abstol: 1e-4
    130133                   reltol: 1e-4)))
     
    137140          (if (negative? flag) (error 'main "CVODE solver error" flag))
    138141
     142          ;;(print-results solver tnext)
    139143
    140144          (if (< tnext tf)
     
    145149        (let ((v (f64vector-ref yy 0))
    146150              (w (f64vector-ref yy 1)))
    147           (assert (< (abs (- v -30.0536210523504)) 1e-12))
    148           (assert (< (abs (- w 0.0402640532124851 )) 1e-12))))
     151          (assert (< (abs (- v -30.0333833471119 )) 1e-12))
     152          (assert (< (abs (- w 0.04023343585461 )) 1e-12))
     153          ))
    149154
    150155      (cvode-destroy-solver solver)
     
    163168    ;; CVODE initialization
    164169    (let ((solver (cvode-create-solver/unsafe
    165                    t0 tf yy rhs/unsafe
     170                   t0 yy rhs/unsafe
     171                   tstop: tf
    166172                   abstol: 1e-4
    167173                   reltol: 1e-4)))
     
    173179        (let ((flag  (cvode-solve solver tnext)))
    174180          (if (negative? flag) (error 'main "CVODE solver error" flag))
     181
     182          ;;(print-results solver tnext)
    175183
    176184          (if (< tnext tf)
     
    181189        (let ((v (f64vector-ref yy 0))
    182190              (w (f64vector-ref yy 1)))
    183           (assert (< (abs (- v -30.0536210523504)) 1e-12))
    184           (assert (< (abs (- w 0.0402640532124851 )) 1e-12))))
     191          (assert (< (abs (- v -30.0333833471119 )) 1e-12))
     192          (assert (< (abs (- w 0.04023343585461 )) 1e-12))
     193          ))
    185194
    186195      (cvode-destroy-solver solver)
Note: See TracChangeset for help on using the changeset viewer.