db5val – bspline-fortran

public interface db5val

5D evaluation routines.

Calls

Help

Called by

Help

Module Procedures

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

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

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.
real(kind=wp),	intent(in)			::	rval	$r$ 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.
integer(kind=ip),	intent(in)			::	idr	$r$ 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 db5ink)
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 db5ink)
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 db5ink)
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 db5ink)
real(kind=wp),	intent(in),		dimension(nr+kr)	::	tr	sequence of knots defining the piecewise polynomial in the $r$ direction. (same as in last call to db5ink)
integer(kind=ip),	intent(in)			::	nx	the number of interpolation points in $x$ . (same as in last call to db5ink)
integer(kind=ip),	intent(in)			::	ny	the number of interpolation points in $y$ . (same as in last call to db5ink)
integer(kind=ip),	intent(in)			::	nz	the number of interpolation points in $z$ . (same as in last call to db5ink)
integer(kind=ip),	intent(in)			::	nq	the number of interpolation points in $q$ . (same as in last call to db5ink)
integer(kind=ip),	intent(in)			::	nr	the number of interpolation points in $r$ . (same as in last call to db5ink)
integer(kind=ip),	intent(in)			::	kx	order of polynomial pieces in $x$ . (same as in last call to db5ink)
integer(kind=ip),	intent(in)			::	ky	order of polynomial pieces in $y$ . (same as in last call to db5ink)
integer(kind=ip),	intent(in)			::	kz	order of polynomial pieces in $z$ . (same as in last call to db5ink)
integer(kind=ip),	intent(in)			::	kq	order of polynomial pieces in $q$ . (same as in last call to db5ink)
integer(kind=ip),	intent(in)			::	kr	order of polynomial pieces in $r$ . (same as in last call to db5ink)
real(kind=wp),	intent(in),		dimension(nx,ny,nz,nq,nr)	::	bcoef	the b-spline coefficients computed by db5ink.
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)			::	inbvr	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.
integer(kind=ip),	intent(inout)			::	ilor	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,kr)	::	w4	work array
real(kind=wp),	intent(inout),		dimension(kz,kq,kr)	::	w3	work array
real(kind=wp),	intent(inout),		dimension(kq,kr)	::	w2	work array
real(kind=wp),	intent(inout),		dimension(kr)	::	w1	work array
real(kind=wp),	intent(inout),		dimension(3_ip*max(kx,ky,kz,kq,kr))	::	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 db5val_multi(xval, yval, zval, qval, rval, idx, idy, idz, idq, idr, tx, ty, tz, tq, tr, nx, ny, nz, nq, nr, kx, ky, kz, kq, kr, bcoef, f, iflag, inbvx, iloy, iloz, iloq, ilor, w4, w3, w2, w1, extrap)

Multiple-output version of db5val: evaluates all nf outputs of a spline from db5ink (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.
real(kind=wp),	intent(in)			::	rval	$r$ 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.
integer(kind=ip),	intent(in)			::	idr	$r$ 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.
real(kind=wp),	intent(in),		dimension(nr+kr)	::	tr	knots in the $r$ 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)			::	nr	the number of interpolation points in $r$ .
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$ .
integer(kind=ip),	intent(in)			::	kr	order of polynomial pieces in $r$ .
real(kind=wp),	intent(in),		dimension(:,:,:,:,:,:)	::	bcoef	`(nx,ny,nz,nq,nr,nf)` b-spline coefficients from db5ink (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.
integer(kind=ip),	intent(inout)			::	ilor	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,kr)	::	w4	work array of dimension `(ky,kz,kq,kr)`
real(kind=wp),	intent(inout),		dimension(kz,kq,kr)	::	w3	work array of dimension `(kz,kq,kr)`
real(kind=wp),	intent(inout),		dimension(kq,kr)	::	w2	work array of dimension `(kq,kr)`
real(kind=wp),	intent(inout),		dimension(kr)	::	w1	work array of dimension `(kr)`
logical,	intent(in),	optional		::	extrap	if extrapolation is allowed (if not present, default is False)

db5val Interface

Contents

Module Procedures

public interface db5val

Calls

Called by

Module Procedures

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

Arguments

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

Arguments