Repatch (from linux) of OSX IRAF

1 files changed, 29 insertions, 0 deletions

diff --git a/math/deboor/setdat2.f b/math/deboor/setdat2.f
new file mode 100644
index 00000000..cf767a59
--- /dev/null
+++ b/math/deboor/setdat2.f
@@ -0,0 +1,29 @@
+      subroutine setdt2(icount)
+c  from  * a practical guide to splines *  by c. de boor
+c  to be called in main program  l 2 m a i n .
+c     this routine is set up to provide the specific data for example 3
+c     in chapter xiv.
+c
+      integer icount,  i,k,l,ntau
+      real break,coef,gtau,step,tau,totalw,weight,x,round
+c     parameter lpkmax=100,ntmax=200,ltkmax=2000
+c     common / data / ntau, tau(ntmax),gtau(ntmax),weight(ntmax),totalw
+c     common /approx/ break(lpkmax),coef(ltkmax),l,k
+      common / data / ntau, tau(200),gtau(200),weight(200),totalw
+      common /approx/ break(100),coef(2000),l,k
+      round(x) = float(ifix(x*100.))/100.
+      if (icount .gt. 0)		stop
+      icount = icount + 1
+      ntau = 65
+      step = 3./float(ntau-1)
+      do 10 i=1,ntau
+	 tau(i) = i*step
+	 gtau(i) = round(exp(tau(i)))
+   10	 weight(i) = 1.
+      totalw = ntau
+      l = 1
+      break(1) = tau(1)
+      break(2) = tau(ntau)
+      k = 3
+					return
+      end