dppgq8.f
SUBROUTINE DPPGQ8 (FUN, LDC, C, XI, LXI, KK, ID, A, B, INPPV, ERR,
+ ANS, IERR)
C***BEGIN PROLOGUE DPPGQ8
....
....
Warning: this routine is not intended to be user-callable.
....
....
C***SUBSIDIARY
C***PURPOSE Subsidiary to DPFQAD
C***LIBRARY SLATEC
C***TYPE DOUBLE PRECISION (PPGQ8-S, DPPGQ8-D)
C***AUTHOR Jones, R. E., (SNLA)
C***DESCRIPTION
C
C Abstract **** A DOUBLE PRECISION routine ****
C
C DPPGQ8, a modification of GAUS8, integrates the
C product of FUN(X) by the ID-th derivative of a spline
C DPPVAL(LDC,C,XI,LXI,KK,ID,X,INPPV) between limits A and B.
C
C Description of Arguments
C
C Input-- FUN,C,XI,A,B,ERR are DOUBLE PRECISION
C FUN - Name of external function of one argument which
C multiplies DPPVAL.
C LDC - Leading dimension of matrix C, LDC .GE. KK
C C - Matrix of Taylor derivatives of dimension at least
C (K,LXI)
C XI - Breakpoint vector of length LXI+1
C LXI - Number of polynomial pieces
C KK - Order of the spline, KK .GE. 1
C ID - Order of the spline derivative, 0 .LE. ID .LE. KK-1
C A - Lower limit of integral
C B - Upper limit of integral (may be less than A)
C INPPV- Initialization parameter for DPPVAL
C ERR - Is a requested pseudorelative error tolerance. Normally
C pick a value of ABS(ERR) .LT. 1D-3. ANS will normally
C have no more error than ABS(ERR) times the integral of
C the absolute value of FUN(X)*DPPVAL(LDC,C,XI,LXI,KK,ID,X,
C INPPV).
C
C
C Output-- ERR,ANS are DOUBLE PRECISION
C ERR - Will be an estimate of the absolute error in ANS if the
C input value of ERR was negative. (ERR Is unchanged if
C the input value of ERR was nonnegative.) The estimated
C error is solely for information to the user and should
C not be used as a correction to the computed integral.
C ANS - Computed value of integral
C IERR- A status code
C --Normal Codes
C 1 ANS most likely meets requested error tolerance,
C or A=B.
C -1 A and B are too nearly equal to allow normal
C integration. ANS is set to zero.
C --Abnormal Code
C 2 ANS probably does not meet requested error tolerance.
C
C***SEE ALSO DPFQAD
C***ROUTINES CALLED D1MACH, DPPVAL, I1MACH, XERMSG
C***REVISION HISTORY (YYMMDD)
C 800901 DATE WRITTEN
C 890531 Changed all specific intrinsics to generic. (WRB)
C 890911 Removed unnecessary intrinsics. (WRB)
C 891214 Prologue converted to Version 4.0 format. (BAB)
C 900315 CALLs to XERROR changed to CALLs to XERMSG. (THJ)
C 900326 Removed duplicate information from DESCRIPTION section.
C (WRB)
C 900328 Added TYPE section. (WRB)
C 910408 Updated the AUTHOR section. (WRB)
C***END PROLOGUE DPPGQ8