*DECK DPPGQ8 SUBROUTINE DPPGQ8 (FUN, LDC, C, XI, LXI, KK, ID, A, B, INPPV, ERR, + ANS, IERR) C***BEGIN PROLOGUE DPPGQ8 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