texlive[44593] Master/texmf-dist: pst-ode (14jun17)

[commits+karl at tug.org]

Revision: 44593
          http://tug.org/svn/texlive?view=revision&revision=44593
Author:   karl
Date:     2017-06-14 23:14:48 +0200 (Wed, 14 Jun 2017)
Log Message:
-----------
pst-ode (14jun17)

Modified Paths:
--------------
    trunk/Master/texmf-dist/doc/generic/pst-ode/ChangeLog
    trunk/Master/texmf-dist/doc/generic/pst-ode/pst-ode-doc.pdf
    trunk/Master/texmf-dist/dvips/pst-ode/pst-ode.pro
    trunk/Master/texmf-dist/tex/generic/pst-ode/pst-ode.tex

Modified: trunk/Master/texmf-dist/doc/generic/pst-ode/ChangeLog
===================================================================
--- trunk/Master/texmf-dist/doc/generic/pst-ode/ChangeLog	2017-06-14 00:17:45 UTC (rev 44592)
+++ trunk/Master/texmf-dist/doc/generic/pst-ode/ChangeLog	2017-06-14 21:14:48 UTC (rev 44593)
@@ -1,3 +1,7 @@
+2017-06-14
+	* version 0.9
+	* fix: occasional failure of the stepping routine (zero relative error)
+
 2017-06-12
 	* version 0.8
 	* fix: Integration step size underflow doesn't lead to ps2pdf crashing

Modified: trunk/Master/texmf-dist/doc/generic/pst-ode/pst-ode-doc.pdf
===================================================================
(Binary files differ)

Modified: trunk/Master/texmf-dist/dvips/pst-ode/pst-ode.pro
===================================================================
--- trunk/Master/texmf-dist/dvips/pst-ode/pst-ode.pro	2017-06-14 00:17:45 UTC (rev 44592)
+++ trunk/Master/texmf-dist/dvips/pst-ode/pst-ode.pro	2017-06-14 21:14:48 UTC (rev 44593)
@@ -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 @@
     %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 @@
       exch pop pshrink exp 0.1 max sfty mul ddt mul /ddt exch def
       ode at dict tcur ddt add tcur end eq {
         % error: step size underflow in ODEINT
-        (!) odeprint
+        (! t=) odeprint tcur 256 string cvs odeprint
         true
       }{
         (-) odeprint

Modified: trunk/Master/texmf-dist/tex/generic/pst-ode/pst-ode.tex
===================================================================
--- trunk/Master/texmf-dist/tex/generic/pst-ode/pst-ode.tex	2017-06-14 00:17:45 UTC (rev 44592)
+++ trunk/Master/texmf-dist/tex/generic/pst-ode/pst-ode.tex	2017-06-14 21:14:48 UTC (rev 44593)
@@ -11,8 +11,8 @@
 %% ODEs using the Runge-Kutta-Fehlberg (RKF45) method with automatic
 %% step size adjustment
 %%
-\def\fileversion{0.8}
-\def\filedate{2017/06/12}
+\def\fileversion{0.9}
+\def\filedate{2017/06/14}
 
 \csname PSTODELoaded\endcsname
 \let\PSTODELoaded\endinput
@@ -191,7 +191,7 @@
       ]
       globaldict exch /#1 exch cvx put
       false setglobal
-    } repeat
+    } repeat (\string\n) odeprint
     \ifPstODE at saveData statefile closefile \fi
     end % tx at odeDict
   }%