|
|
|
****** |
i n p u t ******
|
|
|
t |
knot sequence, of size left+kpm (at least)
|
|
kpm |
order of spline
|
|
left |
integer indicating interval of interest, viz the interval
(t(left), t(left+1))
|
|
nrow |
number of rows in block to be put together
|
|
|
|
****** |
w o r k a r e a ******
|
|
|
scrtch |
used in bsplvd, of size (kpm,kpm)
|
|
dbiatx |
used to contain derivatives of b-splines, of size (kpm,m+1)
with dbiatx(j,i+1) containing the i-th derivative of the
j-th b-spline of interest
|
|
|
|
****** |
o u t p u t ******
|
|
|
q |
the block, of size (nrow,kpm)
|
|
b |
the corresponding piece of the right side, of size (nrow)
|
|
|
|
****** |
m e t h o d ******
|
|
|
the |
k collocation equations for the interval (t(left),t(left+1))
|
|
are constructed with the aid of the subroutine | | |
d i f e q u ( 2, .,
[char46] ) and interspersed (in order) with the side conditions (if any) in
|
|
this interval, using | | |
d i f e q u ( 3, ., . ) for the information.
|
|
|
the block | | |
q has kpm columns, corresponding to the kpm b-
|
|
|
|
splines of order | | |
kpm which have the interval (t(left),t(left+1))
in their support. the block’s diagonal is part of the diagonal of the
total system. the first equation in this block not overlapped by the
|
|
preceding block is therefore equation | | |
l o w r o w , with lowrow = 1+
number of side conditions in preceding intervals (or blocks).
|
|
|
|
