void e02jfc (Integer nxeval, Integer nyeval, const double xevalm[], const double yevalm[], const double coefs[], double fevalm[], const Integer iopts[], const double opts[], NagError *fail)

The function may be called by the names: e02jfc, nag_fit_dim2_spline_ts_evalm or nag_2d_spline_ts_eval_rect.

3 Description

e02jfc calculates values on a rectangular mesh of a bivariate spline computed by e02jdc. The points in the mesh are defined by

x

coordinates (

x_{i}

), for

i = 1, 2, \dots, n_{x}

, and

y

coordinates (

y_{j}

), for

j = 1, 2, \dots, n_{y}

. This function is derived from the TSFIT package of O. Davydov and F. Zeilfelder.

4 References

Davydov O, Morandi R and Sestini A (2006) Local hybrid approximation for scattered data fitting with bivariate splines Comput. Aided Geom. Design 23 703–721

Davydov O, Sestini A and Morandi R (2005) Local RBF approximation for scattered data fitting with bivariate splines Trends and Applications in Constructive Approximation M. G. de Bruin, D. H. Mache, and J. Szabados, Eds ISNM Vol. 151 Birkhauser 91–102

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 Arguments

1: $nxeval$ – Integer Input: On entry: $n_{x}$ , the number of values in the $x$ direction forming the mesh on which the spline is to be evaluated.

Constraint: $nxeval \geq 1$ .
2: $nyeval$ – Integer Input: On entry: $n_{y}$ , the number of values in the $y$ direction forming the mesh on which the spline is to be evaluated.

Constraint: $nyeval \geq 1$ .
3: $xevalm [nxeval]$ – const double Input: On entry: the $(x_{i})$ values forming the mesh on which the spline is to be evaluated.

Constraint: for all $i$ , $xevalm [i - 1]$ must lie inside, or on the boundary of, the spline's bounding box as determined by e02jdc.
4: $yevalm [nyeval]$ – const double Input: On entry: the $(y_{j})$ values forming the mesh on which the spline is to be evaluated.

Constraint: for all $j$ , $yevalm [j - 1]$ must lie inside, or on the boundary of, the spline's bounding box as determined by e02jdc.
5: $coefs [\dim]$ – const double Communication Array: Note: the dimension, $\dim$ , of this array is dictated by the requirements of associated functions that must have been previously called. This array MUST be the same array passed as argument coefs in the previous call to e02jdc.

On entry: the computed spline coefficients as output from e02jdc.
6: $fevalm [nxeval \times nyeval]$ – double Output: Note: the $(i, j)$ th element of the matrix is stored in $fevalm [(j - 1) \times nxeval + i - 1]$ .

On exit: if $fail . code =$ NE_NOERROR on exit $fevalm [(j - 1) \times nxeval + i - 1]$ contains the computed spline value at $(x_{i}, y_{j})$ .
7: $iopts [\dim]$ – const Integer Communication Array: Note: the dimension, $\dim$ , of this array is dictated by the requirements of associated functions that must have been previously called. This array MUST be the same array passed as argument iopts in the previous call to e02zkc.

On entry: the contents of the array MUST NOT have been modified either directly or indirectly, by a call to e02zkc, between calls to e02jdc and e02jfc.
8: $opts [\dim]$ – const double Communication Array: Note: the dimension, $\dim$ , of this array is dictated by the requirements of associated functions that must have been previously called. This array MUST be the same array passed as argument opts in the previous call to e02zkc.

On entry: the contents of the array MUST NOT have been modified either directly or indirectly, by a call to e02zkc, between calls to e02jdc and e02jfc.
9: $fail$ – NagError * Input/Output: The NAG error argument (see Section 7 in the Introduction to the NAG Library CL Interface).

6 Error Indicators and Warnings

NE_ALLOC_FAIL: Dynamic memory allocation failed.
See Section 3.1.2 in the Introduction to the NAG Library CL Interface for further information.
NE_BAD_PARAM: On entry, argument $⟨ value ⟩$ had an illegal value.
NE_INITIALIZATION: Option arrays are not initialized or are corrupted.
NE_INT: On entry, $nxeval = ⟨ value ⟩$ .
Constraint: $nxeval \geq 1$ .

On entry, $nyeval = ⟨ value ⟩$ .
Constraint: $nyeval \geq 1$ .
NE_INTERNAL_ERROR: An internal error has occurred in this function. Check the function call and any array sizes. If the call is correct then please contact NAG for assistance.
See Section 7.5 in the Introduction to the NAG Library CL Interface for further information.
NE_INVALID_SPLINE: The fitting routine has not been called, or the array of coefficients has been corrupted.
NE_NO_LICENCE: Your licence key may have expired or may not have been installed correctly.
See Section 8 in the Introduction to the NAG Library CL Interface for further information.
NE_POINT_OUTSIDE_RECT: On entry, $xevalm [⟨ value ⟩] = ⟨ value ⟩$ was outside the bounding box.
Constraint: $⟨ value ⟩ \leq xevalm [i - 1] \leq ⟨ value ⟩$ for all $i$ .

On entry, $yevalm [⟨ value ⟩] = ⟨ value ⟩$ was outside the bounding box.
Constraint: $⟨ value ⟩ \leq yevalm [j - 1] \leq ⟨ value ⟩$ for all $j$ .