from * a practical guide to splines * by c. de boor to be called in s m o o t h put delx = x(.+1) - x(.) into v(.,4), put the three bands of q-transp*d into v(.,1-3), and put the three bands of (d*q)-transp*(d*q) at and above the diagonal
into v(.,5-7) . here, q is the tridiagonal matrix of order (npoint-2,npoint) with general row 1/delx(i) , -1/delx(i) - 1/delx(i+1) , 1/delx(i+1) and d is the diagonal matrix with general row dx(i) .
Nemo Release 3.1 | setupq (3) | June 29, 2025 |