evaluate_2d_multi Subroutine

private pure subroutine evaluate_2d_multi(me, xval, yval, idx, idy, f, iflag)

Evaluate all the outputs of a bspline_2d with multiple outputs. This is a wrapper for db2val (multiple-output version).

Type Bound

bspline_2d

Arguments

Type IntentOptional Attributes Name
class(bspline_2d), intent(inout) :: me

the spline
real(kind=wp), intent(in) :: xval

coordinate of evaluation point.
real(kind=wp), intent(in) :: yval

coordinate of evaluation point.
integer(kind=ip), intent(in) :: idx

derivative of piecewise polynomial to evaluate.
integer(kind=ip), intent(in) :: idy

derivative of piecewise polynomial to evaluate.
real(kind=wp), intent(out), dimension(:) :: f

(nf) interpolated values of the outputs (see number_of_outputs)
integer(kind=ip), intent(out) :: iflag

status flag (see db2val)

Calls

proc~~evaluate_2d_multi~~CallsGraph proc~evaluate_2d_multi bspline_2d%evaluate_2d_multi interface~db2val db2val proc~evaluate_2d_multi->interface~db2val proc~db2val_default db2val_default interface~db2val->proc~db2val_default proc~db2val_multi db2val_multi interface~db2val->proc~db2val_multi proc~check_value check_value proc~db2val_default->proc~check_value proc~dbvalu_weights dbvalu_weights proc~db2val_default->proc~dbvalu_weights proc~db2val_multi->proc~check_value proc~db2val_multi->proc~dbvalu_weights proc~dintrv dintrv proc~dbvalu_weights->proc~dintrv proc~get_temp_x_for_extrap get_temp_x_for_extrap proc~dintrv->proc~get_temp_x_for_extrap

Source Code

    pure subroutine evaluate_2d_multi(me,xval,yval,idx,idy,f,iflag)

    implicit none

    class(bspline_2d),intent(inout) :: me  !! the spline
    real(wp),intent(in) :: xval  !! \(x\) coordinate of evaluation point.
    real(wp),intent(in) :: yval  !! \(y\) coordinate of evaluation point.
    integer(ip),intent(in) :: idx  !! \(x\) derivative of piecewise polynomial to evaluate.
    integer(ip),intent(in) :: idy  !! \(y\) derivative of piecewise polynomial to evaluate.
    real(wp),dimension(:),intent(out) :: f  !! `(nf)` interpolated values of the outputs
                                            !! (see [[number_of_outputs]])
    integer(ip),intent(out) :: iflag  !! status flag (see [[db2val]])

    if (me%initialized) then
        call db2val(xval,yval,idx,idy,me%tx,me%ty,&
                    me%nx,me%ny,me%kx,me%ky,&
                    me%bcoef,f,iflag,me%inbvx,me%iloy,&
                    me%work_val_1,extrap=me%extrap)
    else
        f = 0.0_wp
        iflag = 1_ip
    end if
    me%iflag = iflag

    end subroutine evaluate_2d_multi