############################################################ # Demonstration of use of 'tlsoda' and 'chk-tlsoda' ############################################################ # The following disables the 'FORTRAN STOP' messages from # PGI-compiled Fortran code. % setenv NO_STOP_MESSAGE on % pwd; ls /home/phys410/ode/tlsoda Makefile Tlsoda* chk-tlsoda.f dvvfrom.f gpin gpine tlsoda.f % make pgf77 -g -Msecond_underscore -c tlsoda.f pgf77 -g -Msecond_underscore -L/usr/local/PGI/lib tlsoda.o -lp410f \ -lodepack -llinpack -lblas -o tlsoda Linking: pgf77 -g -Msecond_underscore -c chk-tlsoda.f pgf77 -g -Msecond_underscore -c dvvfrom.f pgf77 -g -Msecond_underscore -L/usr/local/PGI/lib chk-tlsoda.o \ dvvfrom.o -lp410f -o chk-tlsoda Linking: % 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.4948095973303711 -1.6788213251945621E-006 0.3750000000000000 0.7325475414412526 -2.4832691574360410E-006 0.5000000000000000 0.9588543029073683 -3.2256989623937526E-006 0.6250000000000000 1.170198397015519 -3.8511345951106923E-006 0.7500000000000000 1.363281824729051 -4.3046823822395854E-006 0.8750000000000000 1.535091537001702 -4.5325296476765367E-006 1.000000000000000 1.682946613665435 -4.6440496417296835E-006 rmserr: 3.2594244952664999E-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: 6.6996249710096993E-005 Tolerance: 1.0e-8 rmserr: 1.0073954017993227E-006 Tolerance: 1.0e-10 rmserr: 1.2902520732298939E-008 Tolerance: 1.0e-12 rmserr: 1.2675330006783407E-010 --------------------------------------------- Applying O(dt^2) approximation of ODE to tol=1.0e-12 results --------------------------------------------- dt rms(residual) 0.3125000000000000 5.6697977263394215E-003 0.1562500000000000 1.4121956022936352E-003 7.8125000000000000E-002 3.5224712324578329E-004 3.9062500000000000E-002 8.7952460543792181E-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: 6.6996249710096993E-005 65536 output times rmserr: 2.8837402677146864E-004 Tolerance: 1.0e-8 No additional output times rmserr: 5.4600515843744528E-009 256 output times rmserr: 1.0073954017993227E-006 65536 output times rmserr: 1.1280220551945170E-004 Tolerance: 1.0e-10 No additional output times rmserr: 9.8365551836988008E-010 256 output times rmserr: 1.2902520732298939E-008 65536 output times rmserr: 1.1271883090994448E-006 Tolerance: 1.0e-12 No additional output times rmserr: 3.9105783431699081E-012 256 output times rmserr: 1.2675330006783407E-010 65536 output times rmserr: 1.5508993178092734E-008