############################################################ # Demonstration of use of 'tlsoda' and 'chk-tlsoda' ############################################################ % pwd; ls /home/phys410/ode/tlsoda Makefile Tlsoda* chk-tlsoda.f dvvfrom.f gpin gpine tlsoda.f % make pgf77 -g -c tlsoda.f pgf77 -g -L/usr/local/PGI/lib tlsoda.o -lp410f \ -lodepack -llinpack -lblas -o tlsoda pgf77 -g -c chk-tlsoda.f pgf77 -g -c dvvfrom.f pgf77 -g -L/usr/local/PGI/lib chk-tlsoda.o \ dvvfrom.o -lp410f -o chk-tlsoda % tlsoda usage: tlsoda % tlsoda 1.0 0.0 2.0 1.0d-6 3 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 0.1250000000000000 0.2493503338815517 -8.6711109632017607E-007 0.2500000000000000 0.4948090295080188 -1.1109989729458360E-006 0.3750000000000000 0.7325463645140672 -1.3063419721714098E-006 0.5000000000000000 0.9588526174206525 -1.5402122465304256E-006 0.6250000000000000 1.170196263557266 -1.7176763417157437E-006 0.7500000000000000 1.363279629007695 -2.1089610266150722E-006 0.8750000000000000 1.535088083213516 -1.0787414617844999E-006 1.000000000000000 1.682942631747096 -6.6213130311607298E-007 rmserr: 1.2937977308201230E-006 ############################################################ # Invoke 'Tlsoda' script to put 'tlsoda' through its paces ############################################################ % Tlsoda --------------------------------------------- Running tlsoda with the following tolerances: 1.0e-6 1.0e-8 1.0e-10 1.0e-12 --------------------------------------------- Tolerance: 1.0e-6 rmserr: 5.4901404035008383E-006 Tolerance: 1.0e-8 rmserr: 1.1801577890437745E-008 Tolerance: 1.0e-10 rmserr: 5.6420426412897609E-010 Tolerance: 1.0e-12 rmserr: 7.4219223237344245E-012 --------------------------------------------- Applying O(dt^2) approximation of ODE to tol=1.0e-12 results --------------------------------------------- dt rms(residual) 0.3125000000000000 5.6697977613005648E-003 0.1562500000000000 1.4121956372318181E-003 7.8125000000000000E-002 3.5224715816071033E-004 3.9062500000000000E-002 8.7952495474958567E-005 --------------------------------------------- Demonstrating dependence of results on number of requested output times --------------------------------------------- Tolerance: 1.0e-6 No additional output times rmserr: 1.1209185226999998E-005 256 output times rmserr: 5.4901404035008383E-006 65536 output times rmserr: 5.0924090121097052E-006 Tolerance: 1.0e-8 No additional output times rmserr: 5.4600515843744528E-009 256 output times rmserr: 1.1801577890437745E-008 65536 output times rmserr: 2.7422522181745856E-008 Tolerance: 1.0e-10 No additional output times rmserr: 9.8276284230249562E-010 256 output times rmserr: 5.6420426412897609E-010 65536 output times rmserr: 6.0607617739567266E-010 Tolerance: 1.0e-12 No additional output times rmserr: 3.8023204681436628E-012 256 output times rmserr: 7.4219223237344245E-012 65536 output times rmserr: 6.3700996654287074E-012