db6ink – bspline-fortran

Calls

Help

Called by

Help

private pure subroutine db6ink_default(x, nx, y, ny, z, nz, q, nq, r, nr, s, ns, fcn, kx, ky, kz, kq, kr, ks, iknot, tx, ty, tz, tq, tr, ts, bcoef, iflag)

Determines the parameters of a function that interpolates the six-dimensional gridded data:

Arguments

Type	Intent	Attributes		Name
real(kind=wp),	intent(in),	dimension(:)	::	x	`(nx)` array of $x$ abcissae. must be strictly increasing.
integer(kind=ip),	intent(in)		::	nx	number of $x$ abcissae ( $\ge 3$ )
real(kind=wp),	intent(in),	dimension(:)	::	y	`(ny)` array of $y$ abcissae. must be strictly increasing.
integer(kind=ip),	intent(in)		::	ny	number of $y$ abcissae ( $\ge 3$ )
real(kind=wp),	intent(in),	dimension(:)	::	z	`(nz)` array of $z$ abcissae. must be strictly increasing.
integer(kind=ip),	intent(in)		::	nz	number of $z$ abcissae ( $\ge 3$ )
real(kind=wp),	intent(in),	dimension(:)	::	q	`(nq)` array of $q$ abcissae. must be strictly increasing.
integer(kind=ip),	intent(in)		::	nq	number of $q$ abcissae ( $\ge 3$ )
real(kind=wp),	intent(in),	dimension(:)	::	r	`(nr)` array of $r$ abcissae. must be strictly increasing.
integer(kind=ip),	intent(in)		::	nr	number of $r$ abcissae ( $\ge 3$ )
real(kind=wp),	intent(in),	dimension(:)	::	s	`(ns)` array of $s$ abcissae. must be strictly increasing.
integer(kind=ip),	intent(in)		::	ns	number of $s$ abcissae ( $\ge 3$ )
real(kind=wp),	intent(in),	dimension(:,:,:,:,:,:)	::	fcn	`(nx,ny,nz,nq,nr,ns)` matrix of function values to interpolate. `fcn(i,j,k,q,r,s)` should contain the function value at the point (`x(i)`,`y(j)`,`z(k)`,`q(l)`,`r(m)`,`s(n)`)
integer(kind=ip),	intent(in)		::	kx	the order of spline pieces in $x$ ( $2 \le k_x < n_x$ ) (order = polynomial degree + 1)
integer(kind=ip),	intent(in)		::	ky	the order of spline pieces in $y$ ( $2 \le k_y < n_y$ ) (order = polynomial degree + 1)
integer(kind=ip),	intent(in)		::	kz	the order of spline pieces in $z$ ( $2 \le k_z < n_z$ ) (order = polynomial degree + 1)
integer(kind=ip),	intent(in)		::	kq	the order of spline pieces in $q$ ( $2 \le k_q < n_q$ ) (order = polynomial degree + 1)
integer(kind=ip),	intent(in)		::	kr	the order of spline pieces in $r$ ( $2 \le k_r < n_r$ ) (order = polynomial degree + 1)
integer(kind=ip),	intent(in)		::	ks	the order of spline pieces in $s$ ( $2 \le k_s < n_s$ ) (order = polynomial degree + 1)
integer(kind=ip),	intent(in)		::	iknot	knot sequence flag: Read more…
real(kind=wp),	intent(inout),	dimension(:)	::	tx	The `(nx+kx)` knots in the $x$ direction for the spline interpolant. Read more…
real(kind=wp),	intent(inout),	dimension(:)	::	ty	The `(ny+ky)` knots in the $y$ direction for the spline interpolant. Read more…
real(kind=wp),	intent(inout),	dimension(:)	::	tz	The `(nz+kz)` knots in the $z$ direction for the spline interpolant. Read more…
real(kind=wp),	intent(inout),	dimension(:)	::	tq	The `(nq+kq)` knots in the $q$ direction for the spline interpolant. Read more…
real(kind=wp),	intent(inout),	dimension(:)	::	tr	The `(nr+kr)` knots in the $r$ direction for the spline interpolant. Read more…
real(kind=wp),	intent(inout),	dimension(:)	::	ts	The `(ns+ks)` knots in the $s$ direction for the spline interpolant. Read more…
real(kind=wp),	intent(out),	dimension(:,:,:,:,:,:)	::	bcoef	`(nx,ny,nz,nq,nr,ns)` matrix of coefficients of the b-spline interpolant.
integer(kind=ip),	intent(out)		::	iflag	Read more…

private pure subroutine db6ink_multi(x, nx, y, ny, z, nz, q, nq, r, nr, s, ns, fcn, kx, ky, kz, kq, kr, ks, iknot, tx, ty, tz, tq, tr, ts, bcoef, iflag)

Multiple-output version of db6ink: interpolates nf functions on the same grid, with the same knots. The last dimension of fcn and bcoef is the output index. Evaluate them with db6val (multiple-output version).

Arguments

Type	Intent	Attributes		Name
real(kind=wp),	intent(in),	dimension(:)	::	x	`(nx)` array of $x$ abcissae. Must be strictly increasing.
integer(kind=ip),	intent(in)		::	nx	Number of $x$ abcissae
real(kind=wp),	intent(in),	dimension(:)	::	y	`(ny)` array of $y$ abcissae. Must be strictly increasing.
integer(kind=ip),	intent(in)		::	ny	Number of $y$ abcissae
real(kind=wp),	intent(in),	dimension(:)	::	z	`(nz)` array of $z$ abcissae. Must be strictly increasing.
integer(kind=ip),	intent(in)		::	nz	Number of $z$ abcissae
real(kind=wp),	intent(in),	dimension(:)	::	q	`(nq)` array of $q$ abcissae. Must be strictly increasing.
integer(kind=ip),	intent(in)		::	nq	Number of $q$ abcissae
real(kind=wp),	intent(in),	dimension(:)	::	r	`(nr)` array of $r$ abcissae. Must be strictly increasing.
integer(kind=ip),	intent(in)		::	nr	Number of $r$ abcissae
real(kind=wp),	intent(in),	dimension(:)	::	s	`(ns)` array of $s$ abcissae. Must be strictly increasing.
integer(kind=ip),	intent(in)		::	ns	Number of $s$ abcissae
real(kind=wp),	intent(in),	dimension(:,:,:,:,:,:,:)	::	fcn	`(nx,ny,nz,nq,nr,ns,nf)` array of function values to interpolate: `fcn(i1,i2,i3,i4,i5,i6,j)` is output `j` at grid point (`x(i1)`,`y(i2)`,`z(i3)`,`q(i4)`,`r(i5)`,`s(i6)`)
integer(kind=ip),	intent(in)		::	kx	The order of spline pieces in $x$
integer(kind=ip),	intent(in)		::	ky	The order of spline pieces in $y$
integer(kind=ip),	intent(in)		::	kz	The order of spline pieces in $z$
integer(kind=ip),	intent(in)		::	kq	The order of spline pieces in $q$
integer(kind=ip),	intent(in)		::	kr	The order of spline pieces in $r$
integer(kind=ip),	intent(in)		::	ks	The order of spline pieces in $s$
integer(kind=ip),	intent(in)		::	iknot	knot sequence flag (see db6ink_default)
real(kind=wp),	intent(inout),	dimension(:)	::	tx	The `(nx+kx)` knots in the $x$ direction
real(kind=wp),	intent(inout),	dimension(:)	::	ty	The `(ny+ky)` knots in the $y$ direction
real(kind=wp),	intent(inout),	dimension(:)	::	tz	The `(nz+kz)` knots in the $z$ direction
real(kind=wp),	intent(inout),	dimension(:)	::	tq	The `(nq+kq)` knots in the $q$ direction
real(kind=wp),	intent(inout),	dimension(:)	::	tr	The `(nr+kr)` knots in the $r$ direction
real(kind=wp),	intent(inout),	dimension(:)	::	ts	The `(ns+ks)` knots in the $s$ direction
real(kind=wp),	intent(out),	dimension(:,:,:,:,:,:,:)	::	bcoef	`(nx,ny,nz,nq,nr,ns,nf)` array of b-spline coefficients
integer(kind=ip),	intent(out)		::	iflag	status flag (see db6ink_default), or 6001 if `fcn` and `bcoef` do not have the same number of outputs (or it is zero)

db6ink Interface

Contents

Module Procedures

public interface db6ink

Calls

Called by

Module Procedures

private pure subroutine db6ink_default(x, nx, y, ny, z, nz, q, nq, r, nr, s, ns, fcn, kx, ky, kz, kq, kr, ks, iknot, tx, ty, tz, tq, tr, ts, bcoef, iflag)

Arguments

private pure subroutine db6ink_multi(x, nx, y, ny, z, nz, q, nq, r, nr, s, ns, fcn, kx, ky, kz, kq, kr, ks, iknot, tx, ty, tz, tq, tr, ts, bcoef, iflag)

Arguments