NAG Library Routine Document

d03pzf (dim1_parab_fd_interp)

1

Purpose

2

Specification

3

Description

4

References

5

Arguments

1:     $\mathbf{npde}$ – IntegerInput

On entry: the number of PDEs.

Constraint: $\mathbf{npde}\ge 1$.

2:     $\mathbf{m}$ – IntegerInput

On entry: the coordinate system used. If the call to d03pzf follows one of the finite difference routines then m must be the same argument m as used in that call. For the Keller box scheme only Cartesian coordinate systems are valid and so m must be set to zero. No check will be made by d03pzf in this case.

$\mathbf{m}=0$

Indicates Cartesian coordinates.

$\mathbf{m}=1$

Indicates cylindrical polar coordinates.

$\mathbf{m}=2$

Indicates spherical polar coordinates.

Constraints:

• $0\le \mathbf{m}\le 2$ following a finite difference routine;
• $\mathbf{m}=0$ following a Keller box scheme routine.

3:     $\mathbf{u}\left(\mathbf{npde},\mathbf{npts}\right)$ – Real (Kind=nag_wp) arrayInput

On entry: the PDE part of the original solution returned in the argument u by the PDE routine.

Constraint: $\mathbf{npde}\ge 1$.

4:     $\mathbf{npts}$ – IntegerInput

On entry: the number of mesh points.

Constraint: $\mathbf{npts}\ge 3$.

5:     $\mathbf{x}\left(\mathbf{npts}\right)$ – Real (Kind=nag_wp) arrayInput

On entry: $\mathbf{x}\left(\mathit{i}\right)$, for $\mathit{i}=1,2,\dots ,\mathbf{npts}$, must contain the mesh points as used by the PDE routine.

6:     $\mathbf{xp}\left(\mathbf{intpts}\right)$ – Real (Kind=nag_wp) arrayInput

On entry: $\mathbf{xp}\left(\mathit{i}\right)$, for $\mathit{i}=1,2,\dots ,\mathbf{intpts}$, must contain the spatial interpolation points.

Constraint: $\mathbf{x}\left(1\right)\le \mathbf{xp}\left(1\right)<\mathbf{xp}\left(2\right)<\cdots <\mathbf{xp}\left(\mathbf{intpts}\right)\le \mathbf{x}\left(\mathbf{npts}\right)$.

7:     $\mathbf{intpts}$ – IntegerInput

On entry: the number of interpolation points.

Constraint: $\mathbf{intpts}\ge 1$.

8:     $\mathbf{itype}$ – IntegerInput

On entry: specifies the interpolation to be performed.

$\mathbf{itype}=1$

The solutions at the interpolation points are computed.

$\mathbf{itype}=2$

Both the solutions and their first derivatives at the interpolation points are computed.

Constraint: $\mathbf{itype}=1$ or $2$.

9:     $\mathbf{up}\left(\mathbf{npde},\mathbf{intpts},\mathbf{itype}\right)$ – Real (Kind=nag_wp) arrayOutput

On exit: if $\mathbf{itype}=1$, $\mathbf{up}\left(\mathit{i},\mathit{j},1\right)$, contains the value of the solution $U_{\mathit{i}}\left(x_{\mathit{j}},t_{\mathrm{out}}\right)$, at the interpolation points $x_{\mathit{j}}=\mathbf{xp}\left(\mathit{j}\right)$, for $\mathit{j}=1,2,\dots ,\mathbf{intpts}$ and $\mathit{i}=1,2,\dots ,\mathbf{npde}$.

If $\mathbf{itype}=2$, $\mathbf{up}\left(\mathit{i},\mathit{j},1\right)$ contains $U_{\mathit{i}}\left(x_{\mathit{j}},t_{\mathrm{out}}\right)$ and $\mathbf{up}\left(\mathit{i},\mathit{j},2\right)$ contains $\frac{\partial U_{\mathit{i}}}{\partial x}$ at these points.

10:   $\mathbf{ifail}$ – IntegerInput/Output

On entry: ifail must be set to $0$, $-1\text{ or }1$. If you are unfamiliar with this argument you should refer to Section 3.4 in How to Use the NAG Library and its Documentation for details.

For environments where it might be inappropriate to halt program execution when an error is detected, the value $-1\text{ 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 argument, the recommended value is $0$. When the value $-\mathbf{1}\text{ or }\mathbf{1}$ is used it is essential to test the value of ifail on exit.

On exit: $\mathbf{ifail}=\mathbf{0}$ unless the routine detects an error or a warning has been flagged (see Section 6).

6

Error Indicators and Warnings

$\mathbf{ifail}=1$

On entry, $\mathit{i}=〈\mathit{\text{value}}〉$, $\mathbf{x}\left(\mathit{i}\right)=〈\mathit{\text{value}}〉$, $\mathit{j}=〈\mathit{\text{value}}〉$ and $\mathbf{x}\left(\mathit{j}\right)=〈\mathit{\text{value}}〉$.
Constraint: $\mathbf{x}\left(1\right)<\mathbf{x}\left(2\right)<\cdots <\mathbf{x}\left(\mathbf{npts}\right)$.

On entry, $\mathbf{intpts}=〈\mathit{\text{value}}〉$.
Constraint: $\mathbf{intpts}\ge 1$.

On entry, $\mathbf{itype}=〈\mathit{\text{value}}〉$.
Constraint: $\mathbf{itype}=1$ or $2$.

On entry, $\mathbf{m}=〈\mathit{\text{value}}〉$.
Constraint: $\mathbf{m}=0$, $1$ or $2$.

On entry, $\mathbf{npde}=〈\mathit{\text{value}}〉$.
Constraint: $\mathbf{npde}>0$.

On entry, $\mathbf{npts}=〈\mathit{\text{value}}〉$.
Constraint: $\mathbf{npts}>2$.

$\mathbf{ifail}=2$

On entry, $\mathit{i}=〈\mathit{\text{value}}〉$, $\mathbf{xp}\left(\mathit{i}\right)=〈\mathit{\text{value}}〉$, $\mathit{j}=〈\mathit{\text{value}}〉$ and $\mathbf{xp}\left(\mathit{j}\right)=〈\mathit{\text{value}}〉$.
Constraint: $\mathbf{x}\left(1\right)\le \mathbf{xp}\left(1\right)<\mathbf{xp}\left(2\right)<\cdots <\mathbf{xp}\left(\mathbf{intpts}\right)\le \mathbf{x}\left(\mathbf{npts}\right)$.

$\mathbf{ifail}=3$

On entry, interpolating point $〈\mathit{\text{value}}〉$ with the value $〈\mathit{\text{value}}〉$ is outside the x range.

$\mathbf{ifail}=-99$

An unexpected error has been triggered by this routine. Please contact NAG.

See Section 3.9 in How to Use the NAG Library and its Documentation for further information.

$\mathbf{ifail}=-399$

Your licence key may have expired or may not have been installed correctly.

See Section 3.8 in How to Use the NAG Library and its Documentation for further information.

$\mathbf{ifail}=-999$

Dynamic memory allocation failed.

See Section 3.7 in How to Use the NAG Library and its Documentation for further information.

7

Accuracy

8

Parallelism and Performance

9

Further Comments

10

Example