DEVPCK(3)

SYNOPSIS
SUBROUTINE DEVPCK (LOUT, KPRINT, X, Y, F, FX, FY, XE, YE, FE, DE, + FE2, FAIL)

NAME
devpck - Test usage of increment argument in DPCHFD and DPCHFE for DPCHQ1.

LIBRARY
SLATEC (PCHIP)

TYPE
DOUBLE PRECISION (EVPCCK-S, DEVPCK-D)

KEYWORDS
PCHIP EVALUATOR QUICK CHECK

AUTHOR
Fritsch, F. N., (LLNL)

DESCRIPTION
---- CODE TO TEST USAGE OF INCREMENT ARGUMENT IN DPCHFD AND DPCHFE --- EVALUATES A BICUBIC FUNCTION AND ITS FIRST PARTIAL DERIVATIVES ON A 4X6 MESH CONTAINED IN A 10X10 ARRAY. INTERPOLATION OF THESE DATA ALONG MESH LINES IN EITHER DIMENSION SHOULD AGREE WITH CORRECT FUNCTION WITHIN ROUNDOFF ERROR. ARRAYS ARE ARGUMENTS ONLY TO ALLOW SHARING STORAGE WITH OTHER TEST ROUTINES. NOTE: RUN WITH KPRINT=4 FOR FULL GORY DETAILS (10 PAGES WORTH). FORTRAN INTRINSICS USED: ABS. FORTRAN LIBRARY ROUTINES USED: (WRITE). SLATEC LIBRARY ROUTINES USED: DPCHFD, DPCHFE, D1MACH.

ROUTINES CALLED
D1MACH, DPCHFD, DPCHFE

REVISION HISTORY (YYMMDD)
820601 DATE WRITTEN
820714 CONVERTED TO QUICK CHECK FOR SLATEC LIBRARY.
820715 1. CORRECTED SOME FORMATS. 2. ADDED CALL TO D1MACH TO SET MACHEP.
890406 1. Modified to make sure final elements of X and XE agree, to avoid possible failure due to roundoff error. 2. Added printout of TOL in case of failure. 3. Removed unnecessary IMPLICIT declaration. 4. Corrected a few S.P. constants to D.P. 5. Minor cosmetic changes.
890706 Cosmetic changes to prologue. (WRB)
890911 Removed unnecessary intrinsics. (WRB)
891004 Cosmetic changes to prologue. (WRB)
891214 Prologue converted to Version 4.0 format. (BAB)
900315 Revised prologue and improved some output formats. (FNF)
900316 Additional minor cosmetic changes. (FNF)
900321 Made miscellaneous cosmetic changes. (FNF)
901130 Made many changes to output: (FNF) 1. Reduced amount of output for KPRINT=3. (Now need to use KPRINT=4 for full output.) 2. Added 1P's to formats and revised some to reduce maximum line length.
910708 Minor modifications in use of KPRINT. (WRB)
930317 Improved output formats. (FNF)