Initial import of new "odepack" package.

Systematized Collection of ODE Solvers

19 files changed, 248 insertions, 0 deletions

diff --git a/math/odepack/Makefile b/math/odepack/Makefile
new file mode 100644
index 00000000000..07325e18854
--- /dev/null
+++ b/math/odepack/Makefile
@@ -0,0 +1,28 @@
+# $NetBSD: Makefile,v 1.1.1.1 2000/11/30 00:51:16 jtb Exp $
+
+DISTNAME=		odepack-20001129
+CATEGORIES=		math
+MASTER_SITES=		${MASTER_SITE_LOCAL}
+
+MAINTAINER=		jtb@netbsd.org
+
+WRKSRC=			${WRKDIR}/odepack
+USE_LIBTOOL=		YES
+USE_FORTRAN=		YES
+
+FOPTS?=			=-O3
+
+.include "../../mk/bsd.prefs.mk"
+
+pre-build:
+	@${SED} -e 's:%%FORTRAN%%:'${FC}':g;s:%%FOPTS%%:'${FOPTS}':g' \
+	< ${FILESDIR}/Makefile > ${WRKSRC}/Makefile
+
+post-install:
+	${INSTALL_DATA_DIR} ${PREFIX}/share/doc/odepack
+	${INSTALL_DATA} ${WRKSRC}/doc ${PREFIX}/share/doc/odepack
+	${INSTALL_DATA} ${WRKSRC}/readme ${PREFIX}/share/doc/odepack
+	${INSTALL_DATA_DIR} ${PREFIX}/share/examples/odepack
+	${INSTALL_DATA} ${WRKSRC}/demo.f ${PREFIX}/share/examples/odepack
+
+.include "../../mk/bsd.pkg.mk"
diff --git a/math/odepack/files/Makefile b/math/odepack/files/Makefile
new file mode 100644
index 00000000000..58b9bf794d9
--- /dev/null
+++ b/math/odepack/files/Makefile
@@ -0,0 +1,28 @@
+LIB	= odepack
+
+LIBDIR	= ${PREFIX}/lib
+FORTRAN = %%FORTRAN%%
+FFLAGS	= %%FOPTS%%
+
+all:	lib$(LIB)
+
+OBJ = adjlr.o aigbt.o ainvg.o blkdta000.o bnorm.o cdrv.o cfode.o cntnzu.o \
+	decbt.o ewset.o fnorm.o intdy.o iprep.o jgroup.o lsoda.o	  \
+	lsodar.o lsode.o lsodes.o lsodi.o lsoibt.o md.o mdi.o mdm.o mdp.o \
+	mdu.o nnfc.o nnsc.o nntc.o nroc.o nsfc.o odrv.o pjibt.o prep.o	  \
+	prepj.o	prepji.o prja.o prjs.o rchek.o roots.o slsbt.o slss.o	  \
+	solbt.o solsy.o srcar.o srcma.o srcms.o srcom.o sro.o stoda.o	  \
+	stode.o stodi.o vmnorm.o vnorm.o xerrwv.o xsetf.o xsetun.o
+
+lib$(LIB):	$(OBJ)
+	@$(LIBTOOL) --mode=link $(FORTRAN) $(FFLAGS) $(FOPTS) \
+	-O -o lib$(LIB).la $(OBJ:.o=.lo) -rpath $(LIBDIR)
+
+.f.o:
+	@$(LIBTOOL) --mode=compile $(FORTRAN) $(FFLAGS) $(FOPTS) -c $*.f
+
+install:	lib$(LIB)
+	@$(LIBTOOL) --mode=install install -c lib$(LIB).la $(LIBDIR)
+
+clean:
+	rm -rf *.o *.lo .libs lib$(LIB).la
\ No newline at end of file
diff --git a/math/odepack/files/md5 b/math/odepack/files/md5
new file mode 100644
index 00000000000..a0fa99c1127
--- /dev/null
+++ b/math/odepack/files/md5
@@ -0,0 +1,3 @@
+$NetBSD: md5,v 1.1.1.1 2000/11/30 00:51:16 jtb Exp $
+
+MD5 (odepack-20001129.tar.gz) = 9ba3c480e911387b515e91bcfc69806e
diff --git a/math/odepack/files/patch-sum b/math/odepack/files/patch-sum
new file mode 100644
index 00000000000..163f3c59532
--- /dev/null
+++ b/math/odepack/files/patch-sum
@@ -0,0 +1,14 @@
+$NetBSD: patch-sum,v 1.1.1.1 2000/11/30 00:51:16 jtb Exp $
+
+MD5 (patch-aa) = 00a4a408ee7be7059ec065877029d1e5
+MD5 (patch-ab) = 5b378ad2a8788baba6b9ac7cdc7a1536
+MD5 (patch-ac) = 54d3cd855f8491e32d89121b4d90e236
+MD5 (patch-ad) = 3fcc0c8c51208cf5c1f6754428e72f63
+MD5 (patch-ae) = 44689456c551ccca51a44557682745ad
+MD5 (patch-af) = baae500a91c90f639825a68d6bc2a71d
+MD5 (patch-ag) = e23f1bc46916a7fa18642441dacb49dd
+MD5 (patch-ah) = 163dbe847e35f5ea1b8b834485c30841
+MD5 (patch-ai) = d946f2eb3bd1a362e391f20170abbbe3
+MD5 (patch-aj) = 7515896d8fc06cedda78b16955ceba89
+MD5 (patch-ak) = c046b34eed189f8fba26dde58314ffcf
+MD5 (patch-al) = 0f84b969ebe41ba523f81bdca2f62573
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),
diff --git a/math/odepack/pkg/COMMENT b/math/odepack/pkg/COMMENT
new file mode 100644
index 00000000000..3637ca57d70
--- /dev/null
+++ b/math/odepack/pkg/COMMENT
@@ -0,0 +1 @@
+Systematized Collection of ODE Solvers
diff --git a/math/odepack/pkg/DESCR b/math/odepack/pkg/DESCR
new file mode 100644
index 00000000000..ea7d4823b64
--- /dev/null
+++ b/math/odepack/pkg/DESCR
@@ -0,0 +1,5 @@
+ODEPACK is a collection of Fortran solvers for the initial value
+problem for ordinary differential equation (ODE) systems.  It
+currently includes six solvers, suitable for both stiff and nonstiff
+systems, and includes solvers for systems given in linearly implicit
+form as well as solvers for systems given in explicit form.
diff --git a/math/odepack/pkg/PLIST b/math/odepack/pkg/PLIST
new file mode 100644
index 00000000000..d5e441bf4ed
--- /dev/null
+++ b/math/odepack/pkg/PLIST
@@ -0,0 +1,9 @@
+@comment $NetBSD: PLIST,v 1.1.1.1 2000/11/30 00:51:16 jtb Exp $
+lib/libodepack.a
+lib/libodepack.la
+lib/libodepack.so.0.0
+share/doc/odepack/doc
+share/doc/odepack/readme
+share/examples/odepack/demo.f
+@dirrm share/examples/odepack
+@dirrm share/doc/odepack