db4val – bspline-fortran

public interface db4val

4D evaluation routines.

Calls

Help

Called by

Help

Module Procedures

private pure subroutine db4val_default(xval, yval, zval, qval, idx, idy, idz, idq, tx, ty, tz, tq, nx, ny, nz, nq, kx, ky, kz, kq, bcoef, f, iflag, inbvx, inbvy, inbvz, inbvq, iloy, iloz, iloq, w3, w2, w1, w0, extrap)

Evaluates the tensor product piecewise polynomial interpolant constructed by the routine db4ink or one of its derivatives at the point (xval,yval,zval,qval).

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=wp),	intent(in)			::	xval	$x$ coordinate of evaluation point.
real(kind=wp),	intent(in)			::	yval	$y$ coordinate of evaluation point.
real(kind=wp),	intent(in)			::	zval	$z$ coordinate of evaluation point.
real(kind=wp),	intent(in)			::	qval	$q$ coordinate of evaluation point.
integer(kind=ip),	intent(in)			::	idx	$x$ derivative of piecewise polynomial to evaluate.
integer(kind=ip),	intent(in)			::	idy	$y$ derivative of piecewise polynomial to evaluate.
integer(kind=ip),	intent(in)			::	idz	$z$ derivative of piecewise polynomial to evaluate.
integer(kind=ip),	intent(in)			::	idq	$q$ derivative of piecewise polynomial to evaluate.
real(kind=wp),	intent(in),		dimension(nx+kx)	::	tx	sequence of knots defining the piecewise polynomial in the $x$ direction. (same as in last call to db4ink)
real(kind=wp),	intent(in),		dimension(ny+ky)	::	ty	sequence of knots defining the piecewise polynomial in the $y$ direction. (same as in last call to db4ink)
real(kind=wp),	intent(in),		dimension(nz+kz)	::	tz	sequence of knots defining the piecewise polynomial in the $z$ direction. (same as in last call to db4ink)
real(kind=wp),	intent(in),		dimension(nq+kq)	::	tq	sequence of knots defining the piecewise polynomial in the $q$ direction. (same as in last call to db4ink)
integer(kind=ip),	intent(in)			::	nx	the number of interpolation points in $x$ . (same as in last call to db4ink)
integer(kind=ip),	intent(in)			::	ny	the number of interpolation points in $y$ . (same as in last call to db4ink)
integer(kind=ip),	intent(in)			::	nz	the number of interpolation points in $z$ . (same as in last call to db4ink)
integer(kind=ip),	intent(in)			::	nq	the number of interpolation points in $q$ . (same as in last call to db4ink)
integer(kind=ip),	intent(in)			::	kx	order of polynomial pieces in $x$ . (same as in last call to db4ink)
integer(kind=ip),	intent(in)			::	ky	order of polynomial pieces in $y$ . (same as in last call to db4ink)
integer(kind=ip),	intent(in)			::	kz	order of polynomial pieces in $z$ . (same as in last call to db4ink)
integer(kind=ip),	intent(in)			::	kq	order of polynomial pieces in $q$ . (same as in last call to db4ink)
real(kind=wp),	intent(in),		dimension(nx,ny,nz,nq)	::	bcoef	the b-spline coefficients computed by db4ink.
real(kind=wp),	intent(out)			::	f	interpolated value
integer(kind=ip),	intent(out)			::	iflag	status flag: Read more…
integer(kind=ip),	intent(inout)			::	inbvx	initialization parameter which must be set to 1 the first time this routine is called, and must not be changed by the user.
integer(kind=ip),	intent(inout)			::	inbvy	not used (kept for backward compatibility)
integer(kind=ip),	intent(inout)			::	inbvz	not used (kept for backward compatibility)
integer(kind=ip),	intent(inout)			::	inbvq	not used (kept for backward compatibility)
integer(kind=ip),	intent(inout)			::	iloy	initialization parameter which must be set to 1 the first time this routine is called, and must not be changed by the user.
integer(kind=ip),	intent(inout)			::	iloz	initialization parameter which must be set to 1 the first time this routine is called, and must not be changed by the user.
integer(kind=ip),	intent(inout)			::	iloq	initialization parameter which must be set to 1 the first time this routine is called, and must not be changed by the user.
real(kind=wp),	intent(inout),		dimension(ky,kz,kq)	::	w3	work array
real(kind=wp),	intent(inout),		dimension(kz,kq)	::	w2	work array
real(kind=wp),	intent(inout),		dimension(kq)	::	w1	work array
real(kind=wp),	intent(inout),		dimension(3_ip*max(kx,ky,kz,kq))	::	w0	not used (kept for backward compatibility)
logical,	intent(in),	optional		::	extrap	if extrapolation is allowed (if not present, default is False)

private pure subroutine db4val_multi(xval, yval, zval, qval, idx, idy, idz, idq, tx, ty, tz, tq, nx, ny, nz, nq, kx, ky, kz, kq, bcoef, f, iflag, inbvx, iloy, iloz, iloq, w3, w2, w1, extrap)

Multiple-output version of db4val: evaluates all nf outputs of a spline from db4ink (multiple-output version) at one point. The interval search and the basis functions in each direction are only computed once, for all of the outputs.

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=wp),	intent(in)			::	xval	$x$ coordinate of evaluation point.
real(kind=wp),	intent(in)			::	yval	$y$ coordinate of evaluation point.
real(kind=wp),	intent(in)			::	zval	$z$ coordinate of evaluation point.
real(kind=wp),	intent(in)			::	qval	$q$ coordinate of evaluation point.
integer(kind=ip),	intent(in)			::	idx	$x$ derivative of piecewise polynomial to evaluate.
integer(kind=ip),	intent(in)			::	idy	$y$ derivative of piecewise polynomial to evaluate.
integer(kind=ip),	intent(in)			::	idz	$z$ derivative of piecewise polynomial to evaluate.
integer(kind=ip),	intent(in)			::	idq	$q$ derivative of piecewise polynomial to evaluate.
real(kind=wp),	intent(in),		dimension(nx+kx)	::	tx	knots in the $x$ direction.
real(kind=wp),	intent(in),		dimension(ny+ky)	::	ty	knots in the $y$ direction.
real(kind=wp),	intent(in),		dimension(nz+kz)	::	tz	knots in the $z$ direction.
real(kind=wp),	intent(in),		dimension(nq+kq)	::	tq	knots in the $q$ direction.
integer(kind=ip),	intent(in)			::	nx	the number of interpolation points in $x$ .
integer(kind=ip),	intent(in)			::	ny	the number of interpolation points in $y$ .
integer(kind=ip),	intent(in)			::	nz	the number of interpolation points in $z$ .
integer(kind=ip),	intent(in)			::	nq	the number of interpolation points in $q$ .
integer(kind=ip),	intent(in)			::	kx	order of polynomial pieces in $x$ .
integer(kind=ip),	intent(in)			::	ky	order of polynomial pieces in $y$ .
integer(kind=ip),	intent(in)			::	kz	order of polynomial pieces in $z$ .
integer(kind=ip),	intent(in)			::	kq	order of polynomial pieces in $q$ .
real(kind=wp),	intent(in),		dimension(:,:,:,:,:)	::	bcoef	`(nx,ny,nz,nq,nf)` b-spline coefficients from db4ink (multiple-output version).
real(kind=wp),	intent(out),		dimension(:)	::	f	`(nf)` interpolated values of the outputs
integer(kind=ip),	intent(out)			::	iflag	status flag: Read more…
integer(kind=ip),	intent(inout)			::	inbvx	initialization parameter which must be set to 1 the first time this routine is called, and must not be changed by the user.
integer(kind=ip),	intent(inout)			::	iloy	initialization parameter which must be set to 1 the first time this routine is called, and must not be changed by the user.
integer(kind=ip),	intent(inout)			::	iloz	initialization parameter which must be set to 1 the first time this routine is called, and must not be changed by the user.
integer(kind=ip),	intent(inout)			::	iloq	initialization parameter which must be set to 1 the first time this routine is called, and must not be changed by the user.
real(kind=wp),	intent(inout),		dimension(ky,kz,kq)	::	w3	work array of dimension `(ky,kz,kq)`
real(kind=wp),	intent(inout),		dimension(kz,kq)	::	w2	work array of dimension `(kz,kq)`
real(kind=wp),	intent(inout),		dimension(kq)	::	w1	work array of dimension `(kq)`
logical,	intent(in),	optional		::	extrap	if extrapolation is allowed (if not present, default is False)

db4val Interface

Contents

Module Procedures

public interface db4val

Calls

Called by

Module Procedures

private pure subroutine db4val_default(xval, yval, zval, qval, idx, idy, idz, idq, tx, ty, tz, tq, nx, ny, nz, nq, kx, ky, kz, kq, bcoef, f, iflag, inbvx, inbvy, inbvz, inbvq, iloy, iloz, iloq, w3, w2, w1, w0, extrap)

Arguments

private pure subroutine db4val_multi(xval, yval, zval, qval, idx, idy, idz, idq, tx, ty, tz, tq, nx, ny, nz, nq, kx, ky, kz, kq, bcoef, f, iflag, inbvx, iloy, iloz, iloq, w3, w2, w1, extrap)

Arguments