@@ -3120,29 +3120,29 @@ pure subroutine dbspvn(t,jhigh,k,index,x,ileft,vnikx,work,iwork,iflag)
31203120
31213121 implicit none
31223122
3123- real (wp),dimension (* ),intent (in ) :: t ! ! knot vector of length `n+k`, where
3124- ! ! `n` = number of b-spline basis functions
3125- ! ! `n` = sum of knot multiplicities-`k`
3126- ! ! dimension `t(ileft+jhigh)`
3127- integer (ip),intent (in ) :: jhigh ! ! order of b-spline, `1 <= jhigh <= k`
3128- integer (ip),intent (in ) :: k ! ! highest possible order
3129- integer (ip),intent (in ) :: index ! ! index = 1 gives basis functions of order `jhigh`
3130- ! ! = 2 denotes previous entry with `work`, `iwork`
3131- ! ! values saved for subsequent calls to
3132- ! ! dbspvn.
3133- real (wp),intent (in ) :: x ! ! argument of basis functions, `t(k) <= x <= t(n+1)`
3134- integer (ip),intent (in ) :: ileft ! ! largest integer such that `t(ileft) <= x < t(ileft+1)`
3135- real (wp),dimension (k),intent (out ) :: vnikx ! ! vector of length `k` for spline values.
3136- real (wp),dimension (* ),intent (out ) :: work ! ! a work vector of length `2*k`
3137- integer (ip),intent (out ) :: iwork ! ! a work parameter. both `work` and `iwork` contain
3138- ! ! information necessary to continue for `index = 2`.
3139- ! ! when `index = 1` exclusively, these are scratch
3140- ! ! variables and can be used for other purposes.
3141- integer (ip),intent (out ) :: iflag ! ! * 0: no errors
3142- ! ! * 201: `k` does not satisfy `k>=1`
3143- ! ! * 202: `jhigh` does not satisfy `1<=jhigh<=k`
3144- ! ! * 203: `index` is not 1 or 2
3145- ! ! * 204: `x` does not satisfy `t(ileft)<=x<=t(ileft+1)`
3123+ real (wp),dimension (* ),intent (in ) :: t ! ! knot vector of length `n+k`, where
3124+ ! ! `n` = number of b-spline basis functions
3125+ ! ! `n` = sum of knot multiplicities-`k`
3126+ ! ! dimension `t(ileft+jhigh)`
3127+ integer (ip),intent (in ) :: jhigh ! ! order of b-spline, `1 <= jhigh <= k`
3128+ integer (ip),intent (in ) :: k ! ! highest possible order
3129+ integer (ip),intent (in ) :: index ! ! index = 1 gives basis functions of order `jhigh`
3130+ ! ! = 2 denotes previous entry with `work`, `iwork`
3131+ ! ! values saved for subsequent calls to
3132+ ! ! dbspvn.
3133+ real (wp),intent (in ) :: x ! ! argument of basis functions, `t(k) <= x <= t(n+1)`
3134+ integer (ip),intent (in ) :: ileft ! ! largest integer such that `t(ileft) <= x < t(ileft+1)`
3135+ real (wp),dimension (k),intent (out ) :: vnikx ! ! vector of length `k` for spline values.
3136+ real (wp),dimension (* ),intent (inout ) :: work ! ! a work vector of length `2*k`
3137+ integer (ip),intent (inout ) :: iwork ! ! a work parameter. both `work` and `iwork` contain
3138+ ! ! information necessary to continue for `index = 2`.
3139+ ! ! when `index = 1` exclusively, these are scratch
3140+ ! ! variables and can be used for other purposes.
3141+ integer (ip),intent (out ) :: iflag ! ! * 0: no errors
3142+ ! ! * 201: `k` does not satisfy `k>=1`
3143+ ! ! * 202: `jhigh` does not satisfy `1<=jhigh<=k`
3144+ ! ! * 203: `index` is not 1 or 2
3145+ ! ! * 204: `x` does not satisfy `t(ileft)<=x<=t(ileft+1)`
31463146
31473147 integer (ip) :: imjp1, ipj, jp1, jp1ml, l
31483148 real (wp) :: vm, vmprev
0 commit comments