Initial import of new "odepack" package.

Systematized Collection of ODE Solvers

12 files changed, 160 insertions, 0 deletions

diff --git a/math/odepack/patches/patch-aa b/math/odepack/patches/patch-aa
new file mode 100644
index 00000000000..60544019dba
--- /dev/null
+++ b/math/odepack/patches/patch-aa
@@ -0,0 +1,13 @@
+$NetBSD: patch-aa,v 1.1.1.1 2000/11/30 00:51:17 jtb Exp $
+
+--- pjibt.f.orig	Fri Nov 17 17:05:55 2000
++++ pjibt.f
+@@ -13,7 +13,7 @@
+      1   ccmax, el0, h, hmin, hmxi, hu, rc, tn, uround
+       double precision con, fac, hl0, r, srur
+       dimension neq(1), y(1), yh(nyh,1), ewt(1), rtem(1),
+-     1   s(1), savr(1), wm(1), iwm(1)
++     1   s(1), savr(1), wm(*), iwm(*)
+       common /ls0001/ rowns(209),
+      2   ccmax, el0, h, hmin, hmxi, hu, rc, tn, uround,
+      3   iownd(14), iowns(6),
diff --git a/math/odepack/patches/patch-ab b/math/odepack/patches/patch-ab
new file mode 100644
index 00000000000..f55b4452154
--- /dev/null
+++ b/math/odepack/patches/patch-ab
@@ -0,0 +1,15 @@
+$NetBSD: patch-ab,v 1.1.1.1 2000/11/30 00:51:17 jtb Exp $
+
+--- prepj.f.orig	Fri Nov 17 17:05:55 2000
++++ prepj.f
+@@ -13,8 +13,8 @@
+      1   ccmax, el0, h, hmin, hmxi, hu, rc, tn, uround
+       double precision con, di, fac, hl0, r, r0, srur, yi, yj, yjj,
+      1   vnorm
+-      dimension neq(1), y(1), yh(nyh,1), ewt(1), ftem(1), savf(1),
+-     1   wm(1), iwm(1)
++      dimension neq(1), y(1), yh(nyh,*), ewt(1), ftem(1), savf(1),
++     1   wm(*), iwm(*)
+       common /ls0001/ rowns(209),
+      2   ccmax, el0, h, hmin, hmxi, hu, rc, tn, uround,
+      3   iownd(14), iowns(6),
diff --git a/math/odepack/patches/patch-ac b/math/odepack/patches/patch-ac
new file mode 100644
index 00000000000..bdd7c5a4566
--- /dev/null
+++ b/math/odepack/patches/patch-ac
@@ -0,0 +1,13 @@
+$NetBSD: patch-ac,v 1.1.1.1 2000/11/30 00:51:17 jtb Exp $
+
+--- prepji.f.orig	Fri Nov 17 17:05:55 2000
++++ prepji.f
+@@ -13,7 +13,7 @@
+      1   ccmax, el0, h, hmin, hmxi, hu, rc, tn, uround
+       double precision con, fac, hl0, r, srur, yi, yj, yjj
+       dimension neq(1), y(1), yh(nyh,1), ewt(1), rtem(1),
+-     1   s(1), savr(1), wm(1), iwm(1)
++     1   s(1), savr(1), wm(*), iwm(*)
+       common /ls0001/ rowns(209),
+      2   ccmax, el0, h, hmin, hmxi, hu, rc, tn, uround,
+      3   iownd(14), iowns(6),
diff --git a/math/odepack/patches/patch-ad b/math/odepack/patches/patch-ad
new file mode 100644
index 00000000000..68925ddb3dd
--- /dev/null
+++ b/math/odepack/patches/patch-ad
@@ -0,0 +1,15 @@
+$NetBSD: patch-ad,v 1.1.1.1 2000/11/30 00:51:17 jtb Exp $
+
+--- prjs.f.orig	Fri Nov 17 17:05:55 2000
++++ prjs.f
+@@ -16,8 +16,8 @@
+       double precision con0, conmin, ccmxj, psmall, rbig, seth
+       double precision con, di, fac, hl0, pij, r, r0, rcon, rcont,
+      1   srur, vnorm
+-      dimension neq(1), y(1), yh(nyh,1), ewt(1), ftem(1), savf(1),
+-     1   wk(1), iwk(1)
++      dimension neq(1), y(1), yh(nyh,*), ewt(1), ftem(1), savf(1),
++     1   wk(*), iwk(1)
+       common /ls0001/ rowns(209),
+      2   ccmax, el0, h, hmin, hmxi, hu, rc, tn, uround,
+      3   iownd(14), iowns(6),
diff --git a/math/odepack/patches/patch-ae b/math/odepack/patches/patch-ae
new file mode 100644
index 00000000000..e5e4f6cf5e8
--- /dev/null
+++ b/math/odepack/patches/patch-ae
@@ -0,0 +1,13 @@
+$NetBSD: patch-ae,v 1.1.1.1 2000/11/30 00:51:17 jtb Exp $
+
+--- rchek.f.orig	Fri Nov 17 17:05:55 2000
++++ rchek.f
+@@ -3,7 +3,7 @@
+       external g
+       integer job, neq, nyh, jroot, irt
+       double precision y, yh, g0, g1, gx
+-      dimension neq(1), y(1), yh(nyh,1), g0(1), g1(1), gx(1), jroot(1)
++      dimension neq(1), y(1), yh(nyh,*), g0(1), g1(1), gx(1), jroot(1)
+       integer iownd, iowns,
+      1   icf, ierpj, iersl, jcur, jstart, kflag, l, meth, miter,
+      2   maxord, maxcor, msbp, mxncf, n, nq, nst, nfe, nje, nqu
diff --git a/math/odepack/patches/patch-af b/math/odepack/patches/patch-af
new file mode 100644
index 00000000000..7feb7c024af
--- /dev/null
+++ b/math/odepack/patches/patch-af
@@ -0,0 +1,13 @@
+$NetBSD: patch-af,v 1.1.1.1 2000/11/30 00:51:17 jtb Exp $
+
+--- slsbt.f.orig	Fri Nov 17 17:05:55 2000
++++ slsbt.f
+@@ -3,7 +3,7 @@
+       integer iwm
+       integer lblox, lpb, lpc, mb, nb
+       double precision wm, x, tem
+-      dimension wm(1), iwm(1), x(1), tem(1)
++      dimension wm(*), iwm(*), x(1), tem(1)
+ c-----------------------------------------------------------------------
+ c this routine acts as an interface between the core integrator
+ c routine and the solbt routine for the solution of the linear system
diff --git a/math/odepack/patches/patch-ag b/math/odepack/patches/patch-ag
new file mode 100644
index 00000000000..0d43a13cdfd
--- /dev/null
+++ b/math/odepack/patches/patch-ag
@@ -0,0 +1,13 @@
+$NetBSD: patch-ag,v 1.1.1.1 2000/11/30 00:51:18 jtb Exp $
+
+--- slss.f.orig	Fri Nov 17 17:05:55 2000
++++ slss.f
+@@ -14,7 +14,7 @@
+      1   ccmax, el0, h, hmin, hmxi, hu, rc, tn, uround
+       double precision rlss
+       double precision di, hl0, phl0, r
+-      dimension wk(1), iwk(1), x(1), tem(1)
++      dimension wk(*), iwk(1), x(1), tem(1)
+       common /ls0001/ rowns(209),
+      2   ccmax, el0, h, hmin, hmxi, hu, rc, tn, uround,
+      3   iownd(14), iowns(6),
diff --git a/math/odepack/patches/patch-ah b/math/odepack/patches/patch-ah
new file mode 100644
index 00000000000..56caf015c0c
--- /dev/null
+++ b/math/odepack/patches/patch-ah
@@ -0,0 +1,13 @@
+$NetBSD: patch-ah,v 1.1.1.1 2000/11/30 00:51:17 jtb Exp $
+
+--- solsy.f.orig	Fri Nov 17 17:05:55 2000
++++ solsy.f
+@@ -9,7 +9,7 @@
+       double precision rowns,
+      1   ccmax, el0, h, hmin, hmxi, hu, rc, tn, uround
+       double precision di, hl0, phl0, r
+-      dimension wm(1), iwm(1), x(1), tem(1)
++      dimension wm(*), iwm(*), x(1), tem(1)
+       common /ls0001/ rowns(209),
+      2   ccmax, el0, h, hmin, hmxi, hu, rc, tn, uround,
+      3   iownd(14), iowns(6),
diff --git a/math/odepack/patches/patch-ai b/math/odepack/patches/patch-ai
new file mode 100644
index 00000000000..f4cf3d2a3fd
--- /dev/null
+++ b/math/odepack/patches/patch-ai
@@ -0,0 +1,13 @@
+$NetBSD: patch-ai,v 1.1.1.1 2000/11/30 00:51:18 jtb Exp $
+
+--- stoda.f.orig	Fri Nov 17 17:05:55 2000
++++ stoda.f
+@@ -18,7 +18,7 @@
+      1   r, rh, rhdn, rhsm, rhup, told, vmnorm
+       double precision alpha, dm1, dm2, exm1, exm2, pdh, pnorm, rate,
+      1   rh1, rh1it, rh2, rm, sm1
+-      dimension neq(1), y(1), yh(nyh,1), yh1(1), ewt(1), savf(1),
++      dimension neq(1), y(1), yh(nyh,*), yh1(1), ewt(1), savf(1),
+      1   acor(1), wm(1), iwm(1)
+       dimension sm1(12)
+       common /ls0001/ conit, crate, el(13), elco(13,12),
diff --git a/math/odepack/patches/patch-aj b/math/odepack/patches/patch-aj
new file mode 100644
index 00000000000..4f7df817b77
--- /dev/null
+++ b/math/odepack/patches/patch-aj
@@ -0,0 +1,13 @@
+$NetBSD: patch-aj,v 1.1.1.1 2000/11/30 00:51:17 jtb Exp $
+
+--- stode.f.orig	Fri Nov 17 17:05:55 2000
++++ stode.f
+@@ -12,7 +12,7 @@
+      2   ccmax, el0, h, hmin, hmxi, hu, rc, tn, uround
+       double precision dcon, ddn, del, delp, dsm, dup, exdn, exsm, exup,
+      1   r, rh, rhdn, rhsm, rhup, told, vnorm
+-      dimension neq(1), y(1), yh(nyh,1), yh1(1), ewt(1), savf(1),
++      dimension neq(1), y(1), yh(nyh,*), yh1(1), ewt(1), savf(1),
+      1   acor(1), wm(1), iwm(1)
+       common /ls0001/ conit, crate, el(13), elco(13,12),
+      1   hold, rmax, tesco(3,12),
diff --git a/math/odepack/patches/patch-ak b/math/odepack/patches/patch-ak
new file mode 100644
index 00000000000..172e77f70b4
--- /dev/null
+++ b/math/odepack/patches/patch-ak
@@ -0,0 +1,13 @@
+$NetBSD: patch-ak,v 1.1.1.1 2000/11/30 00:51:18 jtb Exp $
+
+--- stodi.f.orig	Fri Nov 17 17:05:55 2000
++++ stodi.f
+@@ -13,7 +13,7 @@
+       double precision dcon, ddn, del, delp, dsm, dup,
+      1   eljh, el1h, exdn, exsm, exup,
+      2   r, rh, rhdn, rhsm, rhup, told, vnorm
+-      dimension neq(1), y(1), yh(nyh,1), yh1(1), ewt(1), savf(1),
++      dimension neq(1), y(1), yh(nyh,*), yh1(1), ewt(1), savf(1),
+      1   savr(1), acor(1), wm(1), iwm(1)
+       common /ls0001/ conit, crate, el(13), elco(13,12),
+      1   hold, rmax, tesco(3,12),
diff --git a/math/odepack/patches/patch-al b/math/odepack/patches/patch-al
new file mode 100644
index 00000000000..09f3e34d36a
--- /dev/null
+++ b/math/odepack/patches/patch-al
@@ -0,0 +1,13 @@
+$NetBSD: patch-al,v 1.1.1.1 2000/11/30 00:51:16 jtb Exp $
+
+--- prja.f.orig	Fri Nov 17 17:22:20 2000
++++ prja.f
+@@ -16,7 +16,7 @@
+       double precision con, fac, hl0, r, r0, srur, yi, yj, yjj,
+      1   vmnorm, fnorm, bnorm
+       dimension neq(1), y(1), yh(nyh,1), ewt(1), ftem(1), savf(1),
+-     1   wm(1), iwm(1)
++     1   wm(*), iwm(*)
+       common /ls0001/ rowns(209),
+      2   ccmax, el0, h, hmin, hmxi, hu, rc, tn, uround,
+      3   iownd(14), iowns(6),