db6val – bspline-fortran

public interface db6val

6D evaluation routines.

Calls

Help

Called by

Help

Module Procedures

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

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

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

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

db6val Interface

Contents

Module Procedures

public interface db6val

Calls

Called by

Module Procedures

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

Arguments

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

Arguments