1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119
|
Robertson ODE test problem:
initial conditions: u0 = 1, v0 = 0, w0 = 0
t u v w
--------------------------------------------------
0.000e+00 1.00000e+00 0.00000e+00 0.00000e+00
1.000e+09 2.07857e-06 -1.09385e-10 9.99998e-01
2.000e+09 1.04202e-06 4.21741e-12 9.99999e-01
3.000e+09 6.94597e-07 2.66732e-12 9.99999e-01
4.000e+09 5.20918e-07 2.08504e-12 9.99999e-01
5.000e+09 4.16721e-07 1.50113e-12 1.00000e+00
6.000e+09 3.47261e-07 3.09313e-12 1.00000e+00
7.000e+09 2.97649e-07 1.12013e-12 1.00000e+00
8.000e+09 2.60438e-07 1.30410e-12 1.00000e+00
9.000e+09 2.31501e-07 6.01213e-15 1.00000e+00
1.000e+10 2.08348e-07 8.07940e-13 1.00000e+00
1.100e+10 1.89407e-07 5.27220e-13 1.00000e+00
1.200e+10 1.73622e-07 5.67013e-13 1.00000e+00
1.300e+10 1.60265e-07 6.62046e-13 1.00000e+00
1.400e+10 1.48820e-07 -1.96317e-12 1.00000e+00
1.500e+10 1.38895e-07 1.17103e-12 1.00000e+00
1.600e+10 1.30214e-07 3.98260e-13 1.00000e+00
1.700e+10 1.22554e-07 6.36369e-13 1.00000e+00
1.800e+10 1.15747e-07 -6.66047e-13 1.00000e+00
1.900e+10 1.09653e-07 4.39347e-13 1.00000e+00
2.000e+10 1.04171e-07 3.61443e-13 1.00000e+00
2.100e+10 9.92098e-08 5.00251e-13 1.00000e+00
2.200e+10 9.47001e-08 4.91884e-13 1.00000e+00
2.300e+10 9.05822e-08 8.71476e-13 1.00000e+00
2.400e+10 8.68083e-08 3.63347e-13 1.00000e+00
2.500e+10 8.33363e-08 -4.58969e-15 1.00000e+00
2.600e+10 8.01306e-08 3.78826e-13 1.00000e+00
2.700e+10 7.71624e-08 6.09154e-13 1.00000e+00
2.800e+10 7.44070e-08 2.04226e-13 1.00000e+00
2.900e+10 7.18410e-08 3.71843e-13 1.00000e+00
3.000e+10 6.94465e-08 4.86116e-14 1.00000e+00
3.100e+10 6.72060e-08 2.63403e-13 1.00000e+00
3.200e+10 6.51057e-08 3.32928e-13 1.00000e+00
3.300e+10 6.31329e-08 2.79447e-13 1.00000e+00
3.400e+10 6.12760e-08 2.69539e-13 1.00000e+00
3.500e+10 5.95252e-08 2.37903e-13 1.00000e+00
3.600e+10 5.78717e-08 2.27549e-13 1.00000e+00
3.700e+10 5.63075e-08 2.16412e-13 1.00000e+00
3.800e+10 5.48258e-08 2.20845e-13 1.00000e+00
3.900e+10 5.34200e-08 1.80732e-13 1.00000e+00
4.000e+10 5.20844e-08 2.07741e-13 1.00000e+00
4.100e+10 5.08141e-08 2.02952e-13 1.00000e+00
4.200e+10 4.96042e-08 1.74303e-13 1.00000e+00
4.300e+10 4.84505e-08 2.18218e-13 1.00000e+00
4.400e+10 4.73494e-08 1.66929e-13 1.00000e+00
4.500e+10 4.62972e-08 1.38049e-13 1.00000e+00
4.600e+10 4.52908e-08 1.08118e-13 1.00000e+00
4.700e+10 4.43271e-08 1.76051e-13 1.00000e+00
4.800e+10 4.34036e-08 1.70607e-13 1.00000e+00
4.900e+10 4.25178e-08 1.81842e-13 1.00000e+00
5.000e+10 4.16676e-08 -1.13757e-14 1.00000e+00
5.100e+10 4.08507e-08 -2.05065e-13 1.00000e+00
5.200e+10 4.00651e-08 -1.78624e-13 1.00000e+00
5.300e+10 3.93089e-08 1.41307e-13 1.00000e+00
5.400e+10 3.85809e-08 1.37592e-13 1.00000e+00
5.500e+10 3.78794e-08 1.22696e-13 1.00000e+00
5.600e+10 3.72032e-08 -1.10452e-13 1.00000e+00
5.700e+10 3.65501e-08 3.36533e-13 1.00000e+00
5.800e+10 3.59205e-08 -2.61161e-13 1.00000e+00
5.900e+10 3.53190e-08 -7.59722e-12 1.00000e+00
6.000e+10 3.47227e-08 1.85374e-13 1.00000e+00
6.100e+10 3.41536e-08 9.39629e-14 1.00000e+00
6.200e+10 3.36021e-08 7.25389e-13 1.00000e+00
6.300e+10 3.30696e-08 -1.63630e-13 1.00000e+00
6.400e+10 3.25554e-08 -2.65856e-12 1.00000e+00
6.500e+10 3.20517e-08 1.30711e-13 1.00000e+00
6.600e+10 3.15661e-08 1.26194e-13 1.00000e+00
6.700e+10 3.10950e-08 1.40543e-13 1.00000e+00
6.800e+10 3.06377e-08 1.09899e-13 1.00000e+00
6.900e+10 3.01937e-08 1.22297e-13 1.00000e+00
7.000e+10 2.97623e-08 1.43754e-13 1.00000e+00
7.100e+10 2.93432e-08 9.00254e-14 1.00000e+00
7.200e+10 2.89356e-08 1.09983e-13 1.00000e+00
7.300e+10 2.85392e-08 1.07976e-13 1.00000e+00
7.400e+10 2.81535e-08 1.10593e-13 1.00000e+00
7.500e+10 2.77782e-08 1.05388e-13 1.00000e+00
7.600e+10 2.74134e-08 -6.45069e-13 1.00000e+00
7.700e+10 2.70583e-08 -1.58650e-12 1.00000e+00
7.800e+10 2.67111e-08 -1.22796e-12 1.00000e+00
7.900e+10 2.63717e-08 8.74686e-14 1.00000e+00
8.000e+10 2.60420e-08 9.65943e-14 1.00000e+00
8.100e+10 2.57204e-08 1.41528e-13 1.00000e+00
8.200e+10 2.54068e-08 1.61138e-13 1.00000e+00
8.300e+10 2.51008e-08 -4.32908e-15 1.00000e+00
8.400e+10 2.48019e-08 6.26448e-14 1.00000e+00
8.500e+10 2.45101e-08 1.00567e-13 1.00000e+00
8.600e+10 2.42251e-08 9.76438e-14 1.00000e+00
8.700e+10 2.39466e-08 1.01689e-13 1.00000e+00
8.800e+10 2.36745e-08 1.26632e-13 1.00000e+00
8.900e+10 2.34085e-08 1.35620e-13 1.00000e+00
9.000e+10 2.31484e-08 8.73805e-14 1.00000e+00
9.100e+10 2.28938e-08 3.36553e-13 1.00000e+00
9.200e+10 2.26446e-08 7.21253e-13 1.00000e+00
9.300e+10 2.24013e-08 5.28310e-13 1.00000e+00
9.400e+10 2.21668e-08 -3.34767e-12 1.00000e+00
9.500e+10 2.19300e-08 1.72230e-13 1.00000e+00
9.600e+10 2.17016e-08 7.81589e-14 1.00000e+00
9.700e+10 2.14779e-08 1.07435e-13 1.00000e+00
9.800e+10 2.12602e-08 -1.41747e-12 1.00000e+00
9.900e+10 2.10435e-08 6.16285e-13 1.00000e+00
1.000e+11 2.08336e-08 8.08491e-14 1.00000e+00
--------------------------------------------------
Final Solver Statistics:
Internal solver steps = 372 (attempted = 563)
Total RHS evals: Fe = 0, Fi = 14565
Total linear solver setups = 581
Total RHS evals for setting up the linear system = 0
Total number of Jacobian evaluations = 391
Total number of Newton iterations = 12497
Total number of nonlinear solver convergence failures = 376
Total number of error test failures = 1
Total number of constraint test failures = 14
Total number of failed steps from solver failure = 176
|