Demonstration program for DLSODE package Problem 1: Van der Pol oscillator: xdotdot - 3*(1 - x**2)*xdot + x = 0, x(0) = 2, xdot(0) = 0 neq = 2 itol = 1 rtol = 0.0E+00 atol = 0.1E-05 Solution with mf = 10 t x xdot nq h 0.13928E+01 0.16801E+01 -0.291E+00 3 0.123E+00 0.36076E+01 -0.77986E-04 -0.317E+01 5 0.217E-01 0.58224E+01 -0.16801E+01 0.291E+00 3 0.475E-01 0.80372E+01 0.11669E-03 0.317E+01 5 0.234E-01 Final statistics for this run: rwork size = 52 iwork size = 20 number of steps = 297 number of f-s = 352 (excluding J-s) = 352 number of J-s = 0 error overrun = 0.12E+03 Solution with mf = 11 t x xdot nq h 0.13928E+01 0.16801E+01 -0.291E+00 5 0.121E+00 0.36076E+01 -0.17732E-04 -0.317E+01 5 0.187E-01 0.58224E+01 -0.16801E+01 0.291E+00 6 0.963E-01 0.80372E+01 0.25894E-04 0.317E+01 5 0.190E-01 Final statistics for this run: rwork size = 58 iwork size = 22 number of steps = 203 number of f-s = 281 (excluding J-s) = 281 number of J-s = 29 error overrun = 0.26E+02 Solution with mf = 12 t x xdot nq h 0.13928E+01 0.16801E+01 -0.291E+00 5 0.121E+00 0.36076E+01 -0.17732E-04 -0.317E+01 5 0.187E-01 0.58224E+01 -0.16801E+01 0.291E+00 6 0.963E-01 0.80372E+01 0.25894E-04 0.317E+01 5 0.190E-01 Final statistics for this run: rwork size = 58 iwork size = 22 number of steps = 203 number of f-s = 339 (excluding J-s) = 281 number of J-s = 29 error overrun = 0.26E+02 Solution with mf = 13 t x xdot nq h 0.13928E+01 0.16801E+01 -0.291E+00 5 0.739E-01 0.36076E+01 0.34401E-04 -0.317E+01 6 0.260E-01 0.58224E+01 -0.16801E+01 0.291E+00 4 0.133E+00 0.80372E+01 -0.59053E-04 0.317E+01 5 0.205E-01 Final statistics for this run: rwork size = 56 iwork size = 20 number of steps = 198 number of f-s = 315 (excluding J-s) = 289 number of J-s = 26 error overrun = 0.59E+02 Solution with mf = 20 t x xdot nq h 0.13928E+01 0.16801E+01 -0.291E+00 5 0.549E-01 0.36076E+01 -0.56579E-04 -0.317E+01 5 0.143E-01 0.58224E+01 -0.16801E+01 0.291E+00 4 0.583E-01 0.80372E+01 0.10387E-03 0.317E+01 5 0.149E-01 Final statistics for this run: rwork size = 38 iwork size = 20 number of steps = 289 number of f-s = 321 (excluding J-s) = 321 number of J-s = 0 error overrun = 0.10E+03 Solution with mf = 21 t x xdot nq h 0.13928E+01 0.16801E+01 -0.291E+00 5 0.676E-01 0.36076E+01 -0.48977E-04 -0.317E+01 5 0.141E-01 0.58224E+01 -0.16801E+01 0.291E+00 5 0.126E+00 0.80372E+01 0.96867E-04 0.317E+01 5 0.142E-01 Final statistics for this run: rwork size = 44 iwork size = 22 number of steps = 262 number of f-s = 345 (excluding J-s) = 345 number of J-s = 30 error overrun = 0.97E+02 Solution with mf = 22 t x xdot nq h 0.13928E+01 0.16801E+01 -0.291E+00 5 0.676E-01 0.36076E+01 -0.48977E-04 -0.317E+01 5 0.141E-01 0.58224E+01 -0.16801E+01 0.291E+00 5 0.126E+00 0.80372E+01 0.96867E-04 0.317E+01 5 0.142E-01 Final statistics for this run: rwork size = 44 iwork size = 22 number of steps = 262 number of f-s = 405 (excluding J-s) = 345 number of J-s = 30 error overrun = 0.97E+02 Solution with mf = 23 t x xdot nq h 0.13928E+01 0.16801E+01 -0.291E+00 5 0.709E-01 0.36076E+01 -0.46705E-04 -0.317E+01 5 0.139E-01 0.58224E+01 -0.16801E+01 0.291E+00 3 0.719E-01 0.80372E+01 0.54700E-04 0.317E+01 5 0.154E-01 Final statistics for this run: rwork size = 42 iwork size = 20 number of steps = 271 number of f-s = 414 (excluding J-s) = 383 number of J-s = 31 error overrun = 0.55E+02 ---------------------------------------------------------------------- Problem 2: ydot = A * y , where A is a banded lower triangular matrix derived from 2-D advection PDE neq = 25 ml = 5 mu = 0 itol = 1 rtol = 0.0E+00 atol = 0.1E-05 Solution with mf = 10 t max.err. nq h 0.10000E-01 0.556E-06 2 0.766E-02 0.10000E+00 0.655E-05 3 0.249E-01 0.10000E+01 0.274E-05 4 0.520E-01 0.10000E+02 0.114E-05 3 0.117E+00 0.10000E+03 0.221E-05 2 0.262E+00 Final statistics for this run: rwork size = 420 iwork size = 20 number of steps = 524 number of f-s = 552 (excluding J-s) = 552 number of J-s = 0 error overrun = 0.65E+01 Solution with mf = 13 t max.err. nq h 0.10000E-01 0.839E-06 2 0.949E-02 0.10000E+00 0.208E-05 3 0.250E-01 0.10000E+01 0.127E-03 3 0.168E-01 0.10000E+02 0.113E-04 3 0.385E+00 0.10000E+03 0.145E-05 2 0.149E+02 Final statistics for this run: rwork size = 447 iwork size = 20 number of steps = 129 number of f-s = 235 (excluding J-s) = 201 number of J-s = 34 error overrun = 0.13E+03 Solution with mf = 14 t max.err. nq h 0.10000E-01 0.877E-06 2 0.965E-02 0.10000E+00 0.206E-05 3 0.250E-01 0.10000E+01 0.126E-05 5 0.935E-01 0.10000E+02 0.311E-06 6 0.442E+00 0.10000E+03 0.159E-07 2 0.291E+02 Final statistics for this run: rwork size = 697 iwork size = 45 number of steps = 92 number of f-s = 113 (excluding J-s) = 113 number of J-s = 18 error overrun = 0.21E+01 Solution with mf = 15 t max.err. nq h 0.10000E-01 0.877E-06 2 0.965E-02 0.10000E+00 0.206E-05 3 0.250E-01 0.10000E+01 0.126E-05 5 0.935E-01 0.10000E+02 0.311E-06 6 0.442E+00 0.10000E+03 0.160E-07 2 0.291E+02 Final statistics for this run: rwork size = 697 iwork size = 45 number of steps = 92 number of f-s = 221 (excluding J-s) = 113 number of J-s = 18 error overrun = 0.21E+01 Solution with mf = 20 t max.err. nq h 0.10000E-01 0.465E-06 2 0.483E-02 0.10000E+00 0.131E-05 3 0.148E-01 0.10000E+01 0.427E-05 5 0.635E-01 0.10000E+02 0.192E-05 4 0.351E+00 0.10000E+03 0.929E-07 1 0.455E+00 Final statistics for this run: rwork size = 245 iwork size = 20 number of steps = 330 number of f-s = 530 (excluding J-s) = 530 number of J-s = 0 error overrun = 0.43E+01 Solution with mf = 23 t max.err. nq h 0.10000E-01 0.101E-05 2 0.598E-02 0.10000E+00 0.446E-06 3 0.146E-01 0.10000E+01 0.153E-05 5 0.738E-01 0.10000E+02 0.578E-06 4 0.324E+00 0.10000E+03 0.908E-08 1 0.992E+02 Final statistics for this run: rwork size = 272 iwork size = 20 number of steps = 180 number of f-s = 325 (excluding J-s) = 274 number of J-s = 51 error overrun = 0.15E+01 Solution with mf = 24 t max.err. nq h 0.10000E-01 0.104E-05 2 0.608E-02 0.10000E+00 0.463E-06 3 0.146E-01 0.10000E+01 0.247E-05 5 0.666E-01 0.10000E+02 0.828E-06 5 0.391E+00 0.10000E+03 0.384E-09 1 0.108E+03 Final statistics for this run: rwork size = 522 iwork size = 45 number of steps = 118 number of f-s = 136 (excluding J-s) = 136 number of J-s = 18 error overrun = 0.25E+01 Solution with mf = 25 t max.err. nq h 0.10000E-01 0.104E-05 2 0.608E-02 0.10000E+00 0.463E-06 3 0.146E-01 0.10000E+01 0.247E-05 5 0.666E-01 0.10000E+02 0.828E-06 5 0.391E+00 0.10000E+03 0.384E-09 1 0.108E+03 Final statistics for this run: rwork size = 522 iwork size = 45 number of steps = 118 number of f-s = 244 (excluding J-s) = 136 number of J-s = 18 error overrun = 0.25E+01 Number of errors encountered = 0