🚀 go-pugleaf

I would like to find the code source of two Akima spline procedures in
Delphi. These procedures are in the IMSLM, a fortran library. I couldn't find
the source code, even in Fortran, and what's more I would like to incorporate
them in my Delphi program, WITHOUT using dll.

* The 1st one computes the Akima cubic spline interpolant. It's name is :
              CSAKM(NDATA, XDATA, FDATA, BREAK, CSCOEF)
Arguments:
   NDATA  - Number of data points.  (Input)
   XDATA  - Array of length NDATA containing the data point
            abscissas.  (Input)
            The data point abscissas must be distinct.
   FDATA  - Array of length NDATA containing the data point
            ordinates.  (Input)
   BREAK  - Array of length NDATA containing the breakpoints for
            the piecewise cubic representation.  (Output)
   CSCOEF - Matrix of size 4 by NDATA containing the local
            coefficients of the cubic pieces.  (Output)


* The 2nd one evaluates a cubic spline. It's name is :
                CSVAL(X, NDATA, BREAK, CSCOEF, Y)
Arguments:
   X      - Point at which the spline is to be evaluated.  (Input)
   NDATA  - Number of polynomial pieces.  (Input)
   BREAK  - Array of length NINTV+1 containing the breakpoints for
            the piecewise cubic representation.  (Input)
            BREAK must be strictly increasing.
   CSCOEF - Matrix of size 4 by NINTV+1 containing the local
            coefficients of the cubic pieces.  (Input)
   CSVAL  - Value of the polynomial at X.  (Output)


I'm not used to spline thory, but somebody showed me how ti use these
procedures and I find them very convinient. Could somebody help me ?

Catherine

-----------== Posted via Deja News, The Discussion Network ==----------
http://www.dejanews.com/       Search, Read, Discuss, or Start Your Own