This section describes functions from file numerics.ct.
[y] = intpol(A...)
intpol(A,index1,index2...) is a general interpolation
function. A must be an array from which values are interpolated.
The rank of A must equal the number of index arguments.
Each index argument may be a real scalar or real array.
All index arguments must mutually agree in type and rank.
The array A may also be complex. The result y is of same
rank and size as each of the index arguments.
intpol(A,i,j,...) is a generalization of mapped indexing
A<[i,j,...]> for non-integral indices. The function benefits
from vectorization even more than most other Tela functions.
Currently intpol uses linear interpolation.
Error codes:
-1: First arg not a numerical array
-2: Rank of first arg does not match number of index args
-3: Non-real index arg
-4: Dissimilar index args
-6: Range overflow
[Lu] = stencil2d_4(u,ap0,am0,a0p,a0m)
stencil2d_4(u,ap0,am0,a0p,a0m) computes the two-dimensional
five-point "molecule" where the coefficient of the central
term is unity:
Lu = u[i,j]
+ ap0*u[i+1,j] + am0*u[i-1,j]
+ a0p*u[i,j+1] + a0m*u[i,j-1];
where the indices i and j run from 2..nx and 2..ny where
[nx,ny] = size(u). The size of ap0,am0,a0p,a0m must be two
less than the size of u in both directions.
Error codes:
-1: One of the arguments is not a real matrix
-2: One of the coefficient args has bad size
Next Chapter, Previous Chapter
Table of contents of this chapter, General table of contents
Top of the document, Beginning of this Chapter