diff --git a/examples/cvode/F2003_serial/cv_diurnal_kry_bp_f2003.out b/examples/cvode/F2003_serial/cv_diurnal_kry_bp_f2003.out new file mode 100644 index 0000000000..4eb0e6fa74 --- /dev/null +++ b/examples/cvode/F2003_serial/cv_diurnal_kry_bp_f2003.out @@ -0,0 +1,48 @@ + + Finished initialization, starting time steps + + t c1 (bottom left middle top right) | lnst lh + t c2 (bottom left middle top right) | lnst lh + ----------------------------------------------------------------------------------- + 7.200000E+03 1.839583E+04 3.009343E+04 1.582817E+04 191 1.632687E+02 + 4.440237E+11 7.263870E+11 3.820457E+11 + 1.440000E+04 1.149769E+07 1.249023E+07 1.329996E+07 223 3.726109E+02 + 4.473889E+11 4.861995E+11 5.178613E+11 + 2.160000E+04 2.525181E+07 7.837545E+07 3.019017E+07 242 3.726109E+02 + 2.824148E+11 9.233155E+11 3.419900E+11 + 2.880000E+04 1.189262E+07 2.395873E+07 1.042564E+07 262 3.726109E+02 + 4.627736E+11 9.345413E+11 4.054191E+11 + 3.600000E+04 2.342995E+04 2.291401E+04 2.522350E+04 285 1.992953E+02 + 5.651797E+11 5.527335E+11 6.084463E+11 + 4.320000E+04 2.608063E-04 4.048821E-05 -4.802744E-04 336 5.653068E+02 + 3.689229E+11 8.238195E+11 4.623800E+11 + 5.040000E+04 1.817289E-04 -1.516970E-06 3.755188E-04 349 5.653068E+02 + 3.862522E+11 1.017917E+12 3.536680E+11 + 5.760000E+04 3.523733E-10 1.237556E-11 -3.779706E-10 362 5.653068E+02 + 5.833332E+11 6.187025E+11 5.789597E+11 + 6.480000E+04 1.065643E-13 -1.378092E-13 -2.363348E-12 374 5.653068E+02 + 4.273617E+11 6.789731E+11 5.386164E+11 + 7.200000E+04 -2.565418E-14 -7.919368E-14 -2.889844E-13 387 5.653068E+02 + 3.453809E+11 1.030257E+12 3.448164E+11 + 7.920000E+04 -7.212118E-14 1.770733E-13 1.512967E-12 400 5.653068E+02 + 5.451911E+11 7.328821E+11 5.110635E+11 + 8.640000E+04 5.298797E-14 -2.802708E-14 2.918795E-13 413 5.653068E+02 + 5.089882E+11 5.756321E+11 5.984473E+11 + ----------------------------------------------------------------------------------- + + General Solver Stats: + Total internal steps taken = 413 + Total rhs function call = 502 + Total num preconditioner evals = 7 + Total num preconditioner solves = 937 + Num error test failures = 16 + Num nonlinear solver iters = 499 + Num linear solver iters = 491 + Avg Krylov subspace dim = 9.839679E-01 + Num nonlinear solver fails = 0 + Num linear solver fails = 0 + main solver real/int workspace sizes = 2689 53 + linear solver real/int workspace sizes = 2454 42 + CVBandPre real/int workspace sizes = 2800 622 + CVBandPre number of f evaluations = 35 + diff --git a/examples/cvode/F2003_serial/cv_diurnal_kry_f2003.out b/examples/cvode/F2003_serial/cv_diurnal_kry_f2003.out new file mode 100644 index 0000000000..c795a3f956 --- /dev/null +++ b/examples/cvode/F2003_serial/cv_diurnal_kry_f2003.out @@ -0,0 +1,46 @@ + + Finished initialization, starting time steps + + t c1 (bottom left middle top right) | lnst lh + t c2 (bottom left middle top right) | lnst lh + ----------------------------------------------------------------------------------- + 7.200000E+03 1.839591E+04 3.009355E+04 1.582807E+04 220 1.471813E+02 + 4.440256E+11 7.263896E+11 3.820433E+11 + 1.440000E+04 1.149781E+07 1.248992E+07 1.330010E+07 251 3.393888E+02 + 4.473937E+11 4.861874E+11 5.178666E+11 + 2.160000E+04 2.525117E+07 7.837531E+07 3.019003E+07 273 3.393888E+02 + 2.824072E+11 9.233129E+11 3.419886E+11 + 2.880000E+04 1.189255E+07 2.395900E+07 1.042540E+07 294 3.393888E+02 + 4.627711E+11 9.345518E+11 4.054103E+11 + 3.600000E+04 2.342993E+04 2.291370E+04 2.522348E+04 318 1.824177E+02 + 5.651802E+11 5.527262E+11 6.084465E+11 + 4.320000E+04 1.056286E-06 1.419210E-06 1.303552E-06 374 1.880511E+02 + 3.689281E+11 8.238008E+11 4.623838E+11 + 5.040000E+04 -4.908995E-09 -2.678148E-09 -2.964936E-09 387 6.366726E+02 + 3.862586E+11 1.017917E+12 3.536566E+11 + 5.760000E+04 -2.802667E-12 -4.033809E-13 -7.939800E-13 398 6.366726E+02 + 5.833440E+11 6.186405E+11 5.790094E+11 + 6.480000E+04 -3.960979E-16 -4.688990E-17 -1.039471E-16 409 6.366726E+02 + 4.273158E+11 6.790339E+11 5.385536E+11 + 7.200000E+04 -3.517091E-19 -1.092965E-19 -1.465617E-19 421 6.366726E+02 + 3.453981E+11 1.030299E+12 3.448146E+11 + 7.920000E+04 -2.182503E-22 4.319341E-23 -1.981237E-24 432 6.366726E+02 + 5.452334E+11 7.327463E+11 5.111093E+11 + 8.640000E+04 -3.382282E-26 4.678058E-26 3.177913E-26 443 6.366726E+02 + 5.089344E+11 5.756561E+11 5.984629E+11 + ----------------------------------------------------------------------------------- + + General Solver Stats: + Total internal steps taken = 443 + Total rhs function call = 617 + Total num preconditioner evals = 67 + Total num preconditioner solves = 1458 + Num error test failures = 21 + Num nonlinear solver iters = 614 + Num linear solver iters = 903 + Avg Krylov subspace dim = 1.470684E+00 + Num nonlinear solver fails = 66 + Num linear solver fails = 67 + main solver real/int workspace sizes = 2689 53 + linear solver real/int workspace sizes = 2454 42 + diff --git a/examples/cvode/F2003_serial/cv_roberts_dnsL_f2003.out b/examples/cvode/F2003_serial/cv_roberts_dnsL_f2003.out new file mode 100644 index 0000000000..1accffe897 --- /dev/null +++ b/examples/cvode/F2003_serial/cv_roberts_dnsL_f2003.out @@ -0,0 +1,43 @@ + + cv_roberts_dnsL_f2003.f90: Robertson CV ODE serial example problem for CVODE + Three equation chemical kinetics problem. + + Linear solver: LAPACK DENSE, with user-supplied Jacobian. +Tolerance parameters: rtol = 0.0001 atol = 1.E-06 1.E-11 1.E-05 +Initial conditions y0 = ( 1.00 0.00 0.00 ) + + --------------------------------------------------- + t y1 y2 y3 + --------------------------------------------------- + 2.6391E-01 9.8997E-01 3.4706E-05 1.0000E-02 + rootsfound[] = 0 1 + 4.0000E-01 9.8517E-01 3.3863E-05 1.4801E-02 + 4.0000E+00 9.0553E-01 2.2406E-05 9.4449E-02 + 4.0000E+01 7.1587E-01 9.1872E-06 2.8412E-01 + 4.0000E+02 4.5055E-01 3.2230E-06 5.4945E-01 + 4.0000E+03 1.8323E-01 8.9439E-07 8.1677E-01 + 4.0000E+04 3.8986E-02 1.6219E-07 9.6101E-01 + 4.0000E+05 4.9340E-03 1.9833E-08 9.9507E-01 + 4.0000E+06 5.1653E-04 2.0672E-09 9.9948E-01 + 1.9818E+07 1.0000E-04 4.0007E-10 9.9990E-01 + rootsfound[] = -1 0 + 4.0000E+07 4.4844E-05 1.7938E-10 9.9996E-01 + 4.0000E+08 6.8491E-06 2.7396E-11 9.9999E-01 + 4.0000E+09 -6.8850E+05 -4.0000E-06 6.8850E+05 + 4.0000E+10 -1.7968E+07 -4.0000E-06 1.7968E+07 + + ------------------------------------------------------ + Final y1' y2' y3' + ------------------------------------------------------ + -4.7999E-04 1.4590E-21 4.7999E-04 + + General Solver Stats: + Total internal steps taken = 745 + Total rhs function calls = 1206 + Total Jacobian function calls = 45 + Total root function calls = 781 + Total LU function calls = 190 + Num error test failures = 46 + Num nonlinear solver iters = 1203 + Num nonlinear solver fails = 38 + diff --git a/examples/cvode/F2003_serial/cv_roberts_dns_constraints_f2003.out b/examples/cvode/F2003_serial/cv_roberts_dns_constraints_f2003.out new file mode 100644 index 0000000000..36da8d04ff --- /dev/null +++ b/examples/cvode/F2003_serial/cv_roberts_dns_constraints_f2003.out @@ -0,0 +1,44 @@ + + cv_roberts_dns_constraints_f2003.f90: Robertson CV ODE serial example problem for CVODE + Three equation chemical kinetics problem. + + Linear solver: DENSE, with user-supplied Jacobian. +Tolerance parameters: rtol = 0.0001 atol = 1.E-06 1.E-11 1.E-05 +Initial conditions y0 = ( 1.00 0.00 0.00 ) +Constraints cval = ( 1.00 1.00 1.00 ) + + --------------------------------------------------- + t y1 y2 y3 + --------------------------------------------------- + 2.6391E-01 9.8997E-01 3.4706E-05 1.0000E-02 + rootsfound[] = 0 1 + 4.0000E-01 9.8517E-01 3.3863E-05 1.4801E-02 + 4.0000E+00 9.0553E-01 2.2406E-05 9.4449E-02 + 4.0000E+01 7.1587E-01 9.1872E-06 2.8412E-01 + 4.0000E+02 4.5055E-01 3.2230E-06 5.4945E-01 + 4.0000E+03 1.8323E-01 8.9439E-07 8.1677E-01 + 4.0000E+04 3.8986E-02 1.6219E-07 9.6101E-01 + 4.0000E+05 4.9340E-03 1.9833E-08 9.9507E-01 + 4.0000E+06 5.1659E-04 2.0675E-09 9.9948E-01 + 1.9607E+07 1.0000E-04 4.0004E-10 9.9990E-01 + rootsfound[] = -1 0 + 4.0000E+07 5.1479E-05 2.0593E-10 9.9995E-01 + 4.0000E+08 6.0238E-06 2.4095E-11 9.9999E-01 + 4.0000E+09 4.6711E-07 1.8685E-12 1.0000E+00 + 4.0000E+10 3.4054E-09 1.3622E-14 1.0000E+00 + + --------------------------------------------------- + Final y1' y2' y3' + --------------------------------------------------- + 7.4513E-20 2.9805E-25 -7.4513E-20 + + General Solver Stats: + Total internal steps taken = 387 + Total rhs function calls = 561 + Total Jacobian function calls = 12 + Total root function calls = 425 + Total LU function calls = 95 + Num error test failures = 15 + Num nonlinear solver iters = 558 + Num nonlinear solver fails = 6 + diff --git a/examples/cvode/F2003_serial/cv_roberts_klu_f2003.out b/examples/cvode/F2003_serial/cv_roberts_klu_f2003.out new file mode 100644 index 0000000000..5ec1165cb9 --- /dev/null +++ b/examples/cvode/F2003_serial/cv_roberts_klu_f2003.out @@ -0,0 +1,44 @@ + + cv_roberts_klu_f2003.f90: Robertson CV ODE serial example problem for CVODE + Three equation chemical kinetics problem. + + Linear solver: DENSE, with user-supplied Jacobian. +Tolerance parameters: rtol = 0.0001 atol = 1.E-06 1.E-12 1.E-04 +Initial conditions y0 = ( 1.00 0.00 0.00 ) + Constraints not used. + + --------------------------------------------------- + t y1 y2 y3 + --------------------------------------------------- + 2.6357E-01 9.8997E-01 3.4705E-05 1.0000E-02 + rootsfound[] = 0 1 + 4.0000E-01 9.8514E-01 3.3859E-05 1.4822E-02 + 4.0000E+00 9.0543E-01 2.2394E-05 9.4551E-02 + 4.0000E+01 7.1583E-01 9.1857E-06 2.8416E-01 + 4.0000E+02 4.5056E-01 3.2234E-06 5.4944E-01 + 4.0000E+03 1.8324E-01 8.9446E-07 8.1676E-01 + 4.0000E+04 3.8981E-02 1.6217E-07 9.6102E-01 + 4.0000E+05 4.9378E-03 1.9848E-08 9.9506E-01 + 4.0000E+06 5.1595E-04 2.0648E-09 9.9948E-01 + 2.0821E+07 1.0000E-04 4.0004E-10 9.9990E-01 + rootsfound[] = -1 0 + 4.0000E+07 5.1479E-05 2.0592E-10 9.9995E-01 + 4.0000E+08 5.4266E-06 2.1707E-11 9.9999E-01 + 4.0000E+09 4.6488E-07 1.8595E-12 1.0000E+00 + 4.0000E+10 4.8655E-08 1.9462E-13 1.0000E+00 + + --------------------------------------------------- + Final y1' y2' y3' + --------------------------------------------------- + -8.8151E-19 -3.5260E-24 8.8151E-19 + + General Solver Stats: + Total internal steps taken = 317 + Total rhs function calls = 941 + Total Jacobian function calls = 40 + Total root function calls = 357 + Total LU function calls = 117 + Num error test failures = 11 + Num nonlinear solver iters = 940 + Num nonlinear solver fails = 40 +