evaluate_5d_multi Subroutine

    pure subroutine evaluate_5d_multi(me,xval,yval,zval,qval,rval,idx,idy,idz,idq,idr,f,iflag)

    implicit none

    class(bspline_5d),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.
    real(wp),intent(in) :: zval  !! \(z\) coordinate of evaluation point.
    real(wp),intent(in) :: qval  !! \(q\) coordinate of evaluation point.
    real(wp),intent(in) :: rval  !! \(r\) 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.
    integer(ip),intent(in) :: idz  !! \(z\) derivative of piecewise polynomial to evaluate.
    integer(ip),intent(in) :: idq  !! \(q\) derivative of piecewise polynomial to evaluate.
    integer(ip),intent(in) :: idr  !! \(r\) 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 [[db5val]])

    if (me%initialized) then
        call db5val(xval,yval,zval,qval,rval,idx,idy,idz,idq,idr,me%tx,me%ty,me%tz,me%tq,me%tr,&
                    me%nx,me%ny,me%nz,me%nq,me%nr,me%kx,me%ky,me%kz,me%kq,me%kr,&
                    me%bcoef,f,iflag,me%inbvx,me%iloy,me%iloz,me%iloq,me%ilor,&
                    me%work_val_1,me%work_val_2,me%work_val_3,me%work_val_4,extrap=me%extrap)
    else
        f = 0.0_wp
        iflag = 1_ip
    end if
    me%iflag = iflag

    end subroutine evaluate_5d_multi