evaluate_2d Subroutine

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

Evaluate a bspline_2d interpolate. This is a wrapper for db2val.

Type Bound

bspline_2d

Arguments

Type IntentOptional Attributes Name
class(bspline_2d), intent(inout) :: me
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) :: f

interpolated value
integer(kind=ip), intent(out) :: iflag

status flag (see db2val)

Calls

proc~~evaluate_2d~~CallsGraph proc~evaluate_2d bspline_oo_module::bspline_2d%evaluate_2d proc~db2val bspline_sub_module::db2val proc~evaluate_2d->proc~db2val proc~check_value bspline_sub_module::check_value proc~db2val->proc~check_value proc~dbvalu bspline_sub_module::dbvalu proc~db2val->proc~dbvalu proc~dintrv bspline_sub_module::dintrv proc~db2val->proc~dintrv proc~dbvalu->proc~dintrv proc~get_temp_x_for_extrap bspline_sub_module::get_temp_x_for_extrap proc~dintrv->proc~get_temp_x_for_extrap

Source Code

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

    implicit none

    class(bspline_2d),intent(inout) :: me
    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),intent(out)            :: f     !! interpolated value
    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%inbvy,me%iloy,&
                    me%work_val_1,me%work_val_2,&
                    extrap=me%extrap)
    else
        iflag = 1_ip
    end if

    me%iflag = iflag

    end subroutine evaluate_2d