The LUSOL function is used in conjunction with the LUDC procedure to solve a set of n linear equations in n unknowns Ax = b. The parameter A is input not as the original array, but as its LU decomposition, created by the routine LUDC.
LUSOL is based on the routine
lubksb described in section 2.3 of Numerical Recipes in C: The Art of Scientific Computing (Second Edition), published by Cambridge University Press, and is used by permission.
Result = LUSOL(A, Index, B [, /COLUMN] [, /DOUBLE])
The result is an n-element vector whose type is identical to A.
The n by n LU decomposition of an array created by the LUDC procedure.
An input vector, created by the LUDC procedure, containing a record of the row permutations which occurred as a result of partial pivoting.
An n-element vector containing the right-hand side of the linear system Ax = b.
Set this keyword if the input array A is in column-major format (composed of column vectors) rather than in row-major format (composed of row vectors).
Set this keyword to force the computation to be done in double-precision arithmetic.
This example solves the linear system Ax = b using LU decomposition and back substitution:
; Define array A: A = [[ 2.0, 1.0, 1.0], $ [ 4.0, -6.0, 0.0], $ [-2.0, 7.0, 2.0]] ; Define right-hand side vector B: B = [3.0, -8.0, 10.0] ; Decompose A: LUDC, A, INDEX ; Compute the solution using back substitution: result = LUSOL(A, INDEX, B) ; Print the result: PRINT, result
1.00000 2.00000 -1.00000
This is the exact solution vector.
Introduced: Pre 4.0
CHOLSOL, CRAMER, GS_ITER, LA_LUSOL, LU_COMPLEX, LUDC, SVSOL, TRISOL