NAG Library Routine Document
E02JFF
Note: before using this routine, please read the Users' Note for your implementation to check the interpretation of bold italicised terms and other implementation-dependent details.
1 Purpose
E02JFF calculates a mesh of values of a spline computed by
E02JDF.
2 Specification
INTEGER
NXEVAL, NYEVAL, IOPTS(*), IFAIL
REAL (KIND=nag_wp)
XEVALM(NXEVAL), YEVALM(NYEVAL), COEFS(*), FEVALM(NXEVAL,NYEVAL), OPTS(*)
3 Description
E02JFF calculates values on a rectangular mesh of a bivariate spline computed by
E02JDF. The points in the mesh are defined by
coordinates (
), for
, and
coordinates (
), for
. This routine is derived from the TSFIT package of O. Davydov and F. Zeilfelder.
4 References
Davydov O and Zeilfelder F (2004) Scattered data fitting by direct extension of local polynomials to bivariate splines Advances in Comp. Math. 21 223–271
Farin G and Hansford D (2000) The Essentials of CAGD Natic, MA: A K Peters, Ltd.
5 Parameters
- 1: NXEVAL – INTEGERInput
On entry: , the number of values in the direction forming the mesh on which the spline is to be evaluated.
Constraint:
.
- 2: NYEVAL – INTEGERInput
On entry: , the number of values in the direction forming the mesh on which the spline is to be evaluated.
Constraint:
.
- 3: XEVALM(NXEVAL) – REAL (KIND=nag_wp) arrayInput
On entry: the values forming the mesh on which the spline is to be evaluated.
Constraint:
for all
,
must lie inside, or on the boundary of, the spline's bounding box as determined by
E02JDF.
- 4: YEVALM(NYEVAL) – REAL (KIND=nag_wp) arrayInput
On entry: the values forming the mesh on which the spline is to be evaluated.
Constraint:
for all
,
must lie inside, or on the boundary of, the spline's bounding box as determined by
E02JDF.
- 5: COEFS() – REAL (KIND=nag_wp) arrayInput
On entry: the computed spline coefficients
COEFS as output from
E02JDF.
- 6: FEVALM(NXEVAL,NYEVAL) – REAL (KIND=nag_wp) arrayOutput
On exit: if on exit contains the computed spline value at .
- 7: IOPTS() – INTEGER arrayCommunication Array
-
On entry: must be the same array
IOPTS supplied in a previous call to
E02JDF. The contents of the array
must not have been modified either directly or indirectly, by a call to
E02ZKF, between calls to
E02JDF and E02JFF.
- 8: OPTS() – REAL (KIND=nag_wp) arrayCommunication Array
-
On entry: must be the same array
OPTS supplied in a previous call to
E02JDF. The contents of the array
must not have been modified either directly or indirectly, by a call to
E02ZKF, between calls to
E02JDF and E02JFF.
- 9: IFAIL – INTEGERInput/Output
-
On entry:
IFAIL must be set to
,
. If you are unfamiliar with this parameter you should refer to
Section 3.3 in the Essential Introduction for details.
For environments where it might be inappropriate to halt program execution when an error is detected, the value
is recommended. If the output of error messages is undesirable, then the value
is recommended. Otherwise, if you are not familiar with this parameter, the recommended value is
.
When the value is used it is essential to test the value of IFAIL on exit.
On exit:
unless the routine detects an error or a warning has been flagged (see
Section 6).
6 Error Indicators and Warnings
If on entry
or
, explanatory error messages are output on the current error message unit (as defined by
X04AAF).
Errors or warnings detected by the routine:
-
On entry, .
Constraint: .
-
On entry, .
Constraint: .
-
Option arrays are not initialized or are corrupted.
-
The fitting routine has not been called, or the array of coefficients has been corrupted.
-
On entry, was outside the bounding box.
Constraint: for all .
-
On entry, was outside the bounding box.
Constraint: for all .
-
Dynamic memory allocation failed.
7 Accuracy
A
real array of length is dynamically allocated by each invocation of E02JFF.
9 Example