QUISPL

Use

The function can be called by any user-defined subroutine.

Format

Fortran Calling Syntax

CALL QUISPL (XVAL, ZVAL, ID, IORD, ARRAY, ERRFLG)

C/C++ Calling Syntax

c_quispl (xval, zval, id, iord, array, errflg)

Python Calling Syntax

[array, errflg] = py_quispl (xval, zval, id, iord)

Compose/MATLAB Calling Syntax

[array, errflg] = m_quispl (xval, zval, id, iord)

Description

QUISPL is a quintic, spline-based interpolation (5th order) that requires a minimum of six points on a curve. It provides a continuous and smooth curve for the zeroth, first, and second derivative. The smoother a derivative is, the easier it is for the solver to converge. However, a noisy dataset with sudden changes in value can led to oscillatory results.

Attributes

XVAL

[double]

The first independent variable of a curve or surface at which the spline function is supposed to interpolate the output value.

ZVAL

[double]

The second independent variable of a surface at which the spline function is supposed to interpolate the output value. Set this value to 0 if only a curve interpolation is needed.

ID

[integer]

The ID of the Reference_Spline element.

IORD

The order of the derivative to be returned by the spline function. Only 0, 1, or 2 are valid entities.

Output

ARRAY

[double]

The vector output value of dimension 3.

• If IORD is 0, then ARRAY(1) contains the computed interpolated value y.
• If IORD is 1, then ARRAY(1) and ARRAY(2) contain the derivatives of y with respect to x and z.
• If IORD is 2, then ARRAY(1), ARRAY(2), and ARRAY(3) contain the second derivatives of y with respect to xx, xz, and zz, respectively.

ERRFLG

[logical]

Logical variable that is returned to the calling subroutine and that indicates the success of the spline function call.

Example

def REQSUB(id, time, par, npar, iflag):
  results = [0.0]*8
  results[0]  = QUISPL(time,0,100,0)
  results[1]  = QUISPL(time,0,100,1)
  results[2]  = QUISPL(time,0,100,2)
  return results