naginterfaces.library.interp.dimn_​scat_​shep_​eval

naginterfaces.library.interp.dimn_scat_shep_eval(x, f, iq, rq, xe)

dimn_scat_shep_eval evaluates the multidimensional interpolating function generated by dimn_scat_shep() and its first partial derivatives.

For full information please refer to the NAG Library document for e01zn

https://support.nag.com/numeric/nl/nagdoc_30.3/flhtml/e01/e01znf.html

Parameters

xfloat, array-like, shape $(d,m)$

Note: the $i$th ordinate of the point $x_j$ is stored in $x[i-1,j-1]$.

Must be the same array supplied as argument $\text{x}$ in the preceding call to dimn_scat_shep(). It must remain unchanged between calls.

ffloat, array-like, shape $\left(m\right)$

Must be the same array supplied as argument $\text{f}$ in the preceding call to dimn_scat_shep(). It must remain unchanged between calls.

iqint, array-like, shape $(2\times m+1)$

Must be the same array returned as argument $\text{iq}$ in the preceding call to dimn_scat_shep(). It must remain unchanged between calls.

rqfloat, array-like, shape $\left(\text{lrq}\right)$

Must be the same array returned as argument $\text{rq}$ in the preceding call to dimn_scat_shep(). It must remain unchanged between calls.

xefloat, array-like, shape $(d,n)$

Note: the $i$th ordinate of the point $x_j$ is stored in $xe[i-1,j-1]$.

$xe[0,\text{j}-1],\dots ,xe[d-1,\text{j}-1]$ must be set to the evaluation point $x_{\text{j}}$, for $\text{j}=1,2,\dots ,n$.

Returns

qfloat, ndarray, shape $\left(n\right)$

$q[\text{i}-1]$ contains the value of the interpolant, at $x_{\text{i}}$, for $\text{i}=1,2,\dots ,n$. If any of these evaluation points lie outside the region of definition of the interpolant the corresponding entries in $q$ are set to an extrapolated approximation, and dimn_scat_shep_eval returns with $errno$ = 3.

qxfloat, ndarray, shape $(d,n)$

$qx[i-1,j-1]$ contains the value of the partial derivatives with respect to the $i$th independent variable (dimension) of the interpolant $Q\left(x\right)$ at $x_{\text{j}}$, for $\text{j}=1,2,\dots ,n$, and for each of the partial derivatives $i=1,2,\dots ,d$. If any of these evaluation points lie outside the region of definition of the interpolant, the corresponding entries in $qx$ are set to extrapolated approximations to the partial derivatives, and dimn_scat_shep_eval returns with $errno$ = 3.

Raises

NagValueError

(errno $1$)

On entry, $n=⟨value⟩$.

Constraint: $n\ge 1$.

(errno $1$)

On entry, $m=⟨value⟩$ and $d=⟨value⟩$.

Constraint: $m\ge (d+1)\times (d+2)/2+2$.

(errno $1$)

On entry, $((d+1)\times (d+2)/2)\times m+2\times d+1$ exceeds the largest machine integer.

$d=⟨value⟩$ and $m=⟨value⟩$.

(errno $1$)

On entry, $d=⟨value⟩$.

Constraint: $d\ge 2$.

(errno $2$)

On entry, values in $rq$ appear to be invalid. Check that $rq$ has not been corrupted between calls to dimn_scat_shep() and dimn_scat_shep_eval.

(errno $2$)

On entry, values in $iq$ appear to be invalid. Check that $iq$ has not been corrupted between calls to dimn_scat_shep() and dimn_scat_shep_eval.

Warns

NagAlgorithmicWarning

(errno $3$)

On entry, at least one evaluation point lies outside the region of definition of the interpolant. At such points the corresponding values in $q$ and $qx$ contain extrapolated approximations. Points should be evaluated one by one to identify extrapolated values.

Notes

dimn_scat_shep_eval takes as input the interpolant $Q\left(x\right)$, $x\in R^d$ of a set of scattered data points $(x_{\text{r}},f_{\text{r}})$, for $\text{r}=1,2,\dots ,m$, as computed by dimn_scat_shep(), and evaluates the interpolant and its first partial derivatives at the set of points $x_{\text{i}}$, for $\text{i}=1,2,\dots ,n$.

dimn_scat_shep_eval must only be called after a call to dimn_scat_shep().

dimn_scat_shep_eval is derived from the new implementation of QS3GRD described by Renka (1988). It uses the modification for high-dimensional interpolation described by Berry and Minser (1999).

References

Berry, M W, Minser, K S, 1999, Algorithm 798: high-dimensional interpolation using the modified Shepard method, ACM Trans. Math. Software (25), 353–366

Renka, R J, 1988, Algorithm 661: QSHEP3D: Quadratic Shepard method for trivariate interpolation of scattered data, ACM Trans. Math. Software (14), 151–152