The routine and its subsidiaries will now return an
error status as supplied by PDA_XERMSG.
SUBROUTINE PDA_DPLINT (N, X, Y, C, STATUS)
***BEGIN PROLOGUE PDA_DPLINT
***PURPOSE Produce the polynomial which interpolates a set of discrete
***TYPE DOUBLE PRECISION (POLINT-S, PDA_DPLINT-D)
***KEYWORDS POLYNOMIAL INTERPOLATION
***AUTHOR Huddleston, R. E., (SNLL)
Subroutine PDA_DPLINT is designed to produce the polynomial which
interpolates the data (X(I),Y(I)), I=1,...,N. PDA_DPLINT sets up
information in the array C which can be used by subroutine PDA_DPOLVL
to evaluate the polynomial and its derivatives and by subroutine
PDA_DPOLCF to produce the coefficients.
*** All TYPE REAL variables are DOUBLE PRECISION ***
N - the number of data points (N .GE. 1)
X - the array of abscissas (all of which must be distinct)
Y - the array of ordinates
C - an array of information used by subroutines
STATUS - Returned error status.
The status must be zero on entry. This
routine does not check the status on entry.
******* Dimensioning Information *******
Arrays X,Y, and C must be dimensioned at least N in the calling
***REFERENCES L. F. Shampine, S. M. Davenport and R. E. Huddleston,
Curve fitting by polynomials in one variable, Report
SLA-74-0270, Sandia Laboratories, June 1974.
***ROUTINES CALLED PDA_XERMSG
***REVISION HISTORY (YYMMDD)
740601 DATE WRITTEN
891006 Cosmetic changes to prologue. (WRB)
891006 REVISION DATE from Version 3.2
891214 Prologue converted to Version 4.0 format. (BAB)
900315 CALLs to XERROR changed to CALLs to PDA_XERMSG. (THJ)
920501 Reformatted the REFERENCES section. (WRB)
950403 Implement status. (HME)
***END PROLOGUE PDA_DPLINT