chem/KPP/kpp/kpp-2.1/test/stochastic_dimer.def

   1 {Decaying-dimerizing reaction set
   2  (Gillespie, J. Chem. Phys. 115(4), p. 1716, 2001)}
   3
   4 #DEFVAR
   5 S1  = ignore; {}
   6 S2  = ignore; {}
   7 S3  = ignore; {}
   8
   9 #EQUATIONS
  10 S1        = PROD        : (1);    {c1}
  11 S1 + S1   = S2          : (0.004);{2*c2}
  12 S2        = S1 + S1     : (0.5);  {c3}
  13 S2        = S3          : (0.04); {c4}
  14
  15 #LOOKATALL
  16
  17 #INITVALUES
  18   CFACTOR = 1.0;
  19   S1   = 1.0e+5;
  20   S2   = 0.0;
  21   S3   = 0.0;
  22
  23 #INLINE F77_INIT
  24       TSTART = 0.d0
  25       TEND   = 30.0d0
  26       DT     = 0.5d0;
  27       DO i=1,NVAR
  28         RTOL(i) = 1.0e-4
  29         ATOL(i) = 1.0e-8
  30       END DO
  31       Volume = 1.0d0
  32 #ENDINLINE
  33
  34 #INLINE F90_INIT
  35   TSTART = 0.d0
  36   TEND   = 30.0d0
  37   DT     = 0.5d0
  38   RTOL(1:NVAR) = 1.0e-4
  39   ATOL(1:NVAR) = 1.0e-8
  40   Volume = 1.0d0
  41 #ENDINLINE
  42
  43 #INLINE MATLAB_INIT
  44   global TSTART TEND DT TEMP
  45   TSTART = 0;
  46   TEND   = 30;
  47   DT     = 0.5;
  48   RTOL(1:NVAR) = 1.0e-4;
  49   ATOL(1:NVAR) = 1.0e-8;
  50   Volume = 1;
  51 #ENDINLINE
  52
  53 #INLINE C_INIT
  54   TSTART = 0.0;
  55   TEND   = 30.0;
  56   DT     = 0.5;
  57   for(i=0; i<NVAR; i++) {
  58       RTOL[i] = 1.0e-4;
  59       ATOL[i] = 1.0e-8;
  60   }
  61   Volume = 1.0;
  62 #ENDINLINE