************80
! SBBLOK solves a linear system that was factored by FCBLOK.
Discussion:
The routine supervises the solution, by forward and backward
substitution, of the linear system
A * x = b
for X, with the PLU factorization of A already generated in FCBLOK.
Individual blocks of equations are solved via SUBFOR and SUBBAK.
Modified:
14 February 2007
Author:
Carl DeBoor
Reference:
Carl DeBoor,
A Practical Guide to Splines,
Springer, 2001,
ISBN: 0387953663,
LC: QA1.A647.v27.
Parameters:
Input, real ( kind = 8 ) BLOKS(*), integer INTEGS(3,NBLOKS), integer
NBLOKS, integer IPIVOT(*), are as on return from FCBLOK.
Input, real ( kind = 8 ) B(*), the right hand side, stored corresponding
to the storage of the equations. See comments in SLVBLK for details.
Output, real ( kind = 8 ) X(*), the solution vector.
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=8) | :: | bloks(*) | ||||
integer(kind=4) | :: | integs(3,nbloks) | ||||
integer(kind=4) | :: | nbloks | ||||
integer(kind=4) | :: | ipivot(*) | ||||
real(kind=8) | :: | b(*) | ||||
real(kind=8) | :: | x(*) |