Changeset 30930 in project


Ignore:
Timestamp:
05/27/14 07:19:01 (6 years ago)
Author:
Ivan Raikov
Message:

sundials: small tweaks to example model parameters

Location:
release/4/sundials/trunk/tests
Files:
3 edited

Legend:

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

    r30926 r30930  
    1212
    1313(define  t0 0.0)
    14 (define  tf 1000.0)
     14(define  tf 500.0)
    1515
    1616
    17 (define C        281.0 )
    18 (define gL       30.0 )
    19 (define EL      -70.6 )
    20 (define VT      -50.4 )
    21 (define Isyn     1000.0 )
    22 (define Delta    2.0 )
    23 (define theta   (+ VT 5.0))
     17(define C       200.0 )
     18(define gL      10.0 )
     19(define EL      -58.0 )
     20(define VT      -50.0 )
     21(define Delta   2.0 )
     22(define theta   0.0 )
    2423(define trefractory  0.25)
    2524
    26 (define tau_w  144.0 )
    27 (define a      4.0 )
    28 (define b      0.0805 )
    29 (define Vr      EL )
     25(define a       2.0 )
     26(define tau_w  120.0 )
     27(define b       100.0 )
     28(define Vr     -46.0 )
    3029
    31 
     30(define Isyn     210.0 )
    3231
    3332(define (subthreshold t yy)
    3433  (let ((V (f64vector-ref yy 0))
    3534        (W (f64vector-ref yy 1)))
    36     (f64vector (/ (+ (* (- gL) (- V EL))
    37                      (* gL Delta (exp (/ (- V VT) Delta))) 
    38                      (- W) 
    39                      Isyn)
     35    (f64vector (/ (+ (* (- gL ) (- V EL))
     36                     (* gL Delta (exp (/ (- V VT) Delta)))
     37                     (- W) Isyn)
    4038                  C)
    4139               (/ (- (* a (- V EL)) W) tau_w)
     
    5048(define (main)
    5149 
    52   (let ((yy (f64vector -65.0 0.0)))
     50  (let ((yy (f64vector -70.0 0.0)))
    5351
    5452    ;; CVODE initialization
  • release/4/sundials/trunk/tests/hh.scm

    r29477 r30930  
    232232          (if (negative? flag) (error 'main "IDA solver error" flag))
    233233
    234 ;;        (print-results/ida solver)
     234          (print-results/ida solver)
    235235
    236236          (if (< tnext tf)
  • release/4/sundials/trunk/tests/izhfs.scm

    r28013 r30930  
    1212
    1313(define  t0 0.0)
    14 (define  tf 80.0)
     14(define  tf 800.0)
    1515
    1616(define k     1.0)
     
    3232(define  FS_U (* FS_b V))
    3333
    34 (define Iext 400.)
     34(define Iext 100.)
    3535
    3636(define (UU v)  (if (< v Vb) 0. (* FS_b (- v Vb) (- v Vb) (- v Vb))))
     
    6868
    6969    ;; CVODE initialization
    70     (let ((h 1e-3)
     70    (let ((h 1e-1)
    7171          (solver (cvode-create-solver
    7272                   t0 yy subthreshold 
     
    8989              (begin
    9090                ;; discrete event: re-initialize solver to integrate over the discontinuity
    91 ;;              (print-results solver)
     91                (print-results solver)
    9292                (cvode-reinit-solver solver (+ tnext h) (vreset (cvode-t solver) h (cvode-yy solver)) )
    9393                (recur tnext (+ iout 1)))
     
    9696                  (recur tnext (+ 1 iout))
    9797                  (begin
    98 ;;                  (print-results solver)
     98
     99                    (print-results solver)
    99100                    (if (< tnext tf)
    100101                        (recur (if (= 1 flag) (+ (+ tnext h) h) (+ tnext h)) (+ 1 iout)))
Note: See TracChangeset for help on using the changeset viewer.