evaluate_1d_multi Subroutine

private pure subroutine evaluate_1d_multi(me, xval, idx, f, iflag)

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

Type Bound

bspline_1d

Arguments

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

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

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

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 db1val)

Calls

proc~~evaluate_1d_multi~~CallsGraph proc~evaluate_1d_multi bspline_1d%evaluate_1d_multi interface~db1val db1val proc~evaluate_1d_multi->interface~db1val proc~db1val_alt db1val_alt interface~db1val->proc~db1val_alt proc~db1val_default db1val_default interface~db1val->proc~db1val_default proc~db1val_multi db1val_multi interface~db1val->proc~db1val_multi proc~check_value check_value proc~db1val_alt->proc~check_value proc~dbvalu dbvalu proc~db1val_alt->proc~dbvalu proc~db1val_default->proc~check_value proc~db1val_default->proc~dbvalu proc~db1val_multi->proc~check_value proc~dbvalu_weights dbvalu_weights proc~db1val_multi->proc~dbvalu_weights proc~dintrv dintrv proc~dbvalu->proc~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_1d_multi(me,xval,idx,f,iflag)

    implicit none

    class(bspline_1d),intent(inout) :: me  !! the spline
    real(wp),intent(in) :: xval  !! \(x\) coordinate of evaluation point.
    integer(ip),intent(in) :: idx  !! \(x\) 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 [[db1val]])

    if (me%initialized) then
        call db1val(xval,idx,me%tx,&
                    me%nx,me%kx,&
                    me%bcoef,f,iflag,me%inbvx,extrap=me%extrap)
    else
        f = 0.0_wp
        iflag = 1_ip
    end if
    me%iflag = iflag

    end subroutine evaluate_1d_multi