E02JFF (PDF version)
E02 Chapter Contents
E02 Chapter Introduction
NAG Library Manual

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.

+ Contents

    1  Purpose
    7  Accuracy
    9  Example

1  Purpose

E02JFF calculates a mesh of values of a spline computed by E02JDF.

2  Specification

SUBROUTINE E02JFF ( NXEVAL, NYEVAL, XEVALM, YEVALM, COEFS, FEVALM, IOPTS, OPTS, IFAIL)
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 x coordinates (xi), for i=1,2,,nx, and y coordinates (yj), for j=1,2,,ny. 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: nx, the number of values in the x direction forming the mesh on which the spline is to be evaluated.
Constraint: NXEVAL1.
2:     NYEVAL – INTEGERInput
On entry: ny, the number of values in the y direction forming the mesh on which the spline is to be evaluated.
Constraint: NYEVAL1.
3:     XEVALM(NXEVAL) – REAL (KIND=nag_wp) arrayInput
On entry: the xi values forming the mesh on which the spline is to be evaluated.
Constraint: for all i, XEVALMi 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 yj values forming the mesh on which the spline is to be evaluated.
Constraint: for all j, YEVALMj 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 IFAIL=0 on exit FEVALMij contains the computed spline value at xi,yj.
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 0, -1​ or ​1. 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 -1​ or ​1 is recommended. If the output of error messages is undesirable, then the value 1 is recommended. Otherwise, if you are not familiar with this parameter, the recommended value is 0. When the value -1​ or ​1 is used it is essential to test the value of IFAIL on exit.
On exit: IFAIL=0 unless the routine detects an error or a warning has been flagged (see Section 6).

6  Error Indicators and Warnings

If on entry IFAIL=0 or -1, explanatory error messages are output on the current error message unit (as defined by X04AAF).
Errors or warnings detected by the routine:
IFAIL=6
On entry, NXEVAL=value.
Constraint: NXEVAL1.
IFAIL=7
On entry, NYEVAL=value.
Constraint: NYEVAL1.
IFAIL=9
Option arrays are not initialized or are corrupted.
IFAIL=10
The fitting routine has not been called, or the array of coefficients has been corrupted.
IFAIL=13
On entry, XEVALMvalue=value was outside the bounding box.
Constraint: valueXEVALMivalue for all i.
IFAIL=14
On entry, YEVALMvalue=value was outside the bounding box.
Constraint: valueYEVALMjvalue for all j.
IFAIL=-999
Dynamic memory allocation failed.

7  Accuracy

E02JFF uses the de Casteljau algorithm and thus is numerically stable. See Farin and Hansford (2000) for details.

8  Further Comments

A real array of length O1 is dynamically allocated by each invocation of E02JFF.

9  Example

See Section 9 in E02JDF.

E02JFF (PDF version)
E02 Chapter Contents
E02 Chapter Introduction
NAG Library Manual

© The Numerical Algorithms Group Ltd, Oxford, UK. 2012