diff options
author | Karl Berry <karl@freefriends.org> | 2017-06-14 21:14:48 +0000 |
---|---|---|
committer | Karl Berry <karl@freefriends.org> | 2017-06-14 21:14:48 +0000 |
commit | a994c5902423e8e50bb47ac78e1a555cad25a06e (patch) | |
tree | d65488fe981855e9e75922871af3c6c103f6ef1e /Master/texmf-dist/dvips | |
parent | 4642b26ce44a10c1fc7be7663d08ad15fc604ef7 (diff) |
pst-ode (14jun17)
git-svn-id: svn://tug.org/texlive/trunk@44593 c570f23f-e606-0410-a88d-b1316a301751
Diffstat (limited to 'Master/texmf-dist/dvips')
-rw-r--r-- | Master/texmf-dist/dvips/pst-ode/pst-ode.pro | 9 |
1 files changed, 5 insertions, 4 deletions
diff --git a/Master/texmf-dist/dvips/pst-ode/pst-ode.pro b/Master/texmf-dist/dvips/pst-ode/pst-ode.pro index 848ac7d728c..df404a2de86 100644 --- a/Master/texmf-dist/dvips/pst-ode/pst-ode.pro +++ b/Master/texmf-dist/dvips/pst-ode/pst-ode.pro @@ -1,7 +1,7 @@ %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % % PostScript prologue for pst-ode.tex. -% Version 0.8, 2017/06/12 +% Version 0.9, 2017/06/14 % % Alexander Grahn (C) 2012--today % @@ -117,10 +117,11 @@ end %scaling vector for step size adjustment (Numerical Recipies) eabsvect k1 eabsvect addvect {1e-30 add} forall xlength array astore % => [x(t)] [x(t+ddt) by RKF4] [err] [xscale] - %maximum rel. error - /errmax 0 def + /errmax 0 def %maximum rel. error edivvect {abs dup errmax gt {/errmax exch def}{pop} ifelse} forall errmax % => [x(t)] [x(t+ddt) by RKF4] errmax + %Sometimes, the rel. error is zero (Why is that???). So we add 1e-30 + 1e-30 add end } bind def /ODEINT { % performs integration over output step [t,t+dt] @@ -137,7 +138,7 @@ end exch pop pshrink exp 0.1 max sfty mul ddt mul /ddt exch def ode@dict tcur ddt add tcur end eq { % error: step size underflow in ODEINT - (!) odeprint + (! t=) odeprint tcur 256 string cvs odeprint true }{ (-) odeprint |