From a994c5902423e8e50bb47ac78e1a555cad25a06e Mon Sep 17 00:00:00 2001 From: Karl Berry Date: Wed, 14 Jun 2017 21:14:48 +0000 Subject: pst-ode (14jun17) git-svn-id: svn://tug.org/texlive/trunk@44593 c570f23f-e606-0410-a88d-b1316a301751 --- Master/texmf-dist/dvips/pst-ode/pst-ode.pro | 9 +++++---- 1 file changed, 5 insertions(+), 4 deletions(-) (limited to 'Master/texmf-dist/dvips') 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 -- cgit v1.2.3