Welcome to the Exelis VIS product documentation center! Here you will find reference guides, help documents, and product libraries. Discover the products ENVI, IDL, ENVI LiDAR, and ESE, developed by Exelis VIS.
﻿

LUSOL

LUSOL

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.

Note: If you are working with complex inputs, use the LA_LUSOL function instead.

Examples

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`

IDL prints:

`  1.00000   2.00000   -1.00000`

This is the exact solution vector.

Syntax

Result = LUSOL(A, Index, B [, /COLUMN] [, /DOUBLE])

Return Value

The result is an n-element vector whose type is identical to A.

Arguments

A

The n by n LU decomposition of an array created by the LUDC procedure.

Note: If LUSOL is complex then only the real part is used for the computation.

Index

An input vector, created by the LUDC procedure, containing a record of the row permutations which occurred as a result of partial pivoting.

B

An n-element vector containing the right-hand side of the linear system Ax = b.

Keywords

COLUMN

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).

DOUBLE

Set this keyword to force the computation to be done in double-precision arithmetic.

Version History

 Pre 4.0 Introduced

Resources and References

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.

﻿

Notes

This page has no user notes yet. Be the first one!

This information is not subject to the controls of the International Traffic in Arms Regulations (ITAR) or the Export Administration Regulations (EAR). However, it may be restricted from transfer to various embargoed countries under U.S. laws and regulations.
 PRODUCTS & SERVICES ENVI IDL Professional Services

 INDUSTRIES Defense & Intelligence Environmental Monitoring Academic
 LEARN Code Videos Events & Training Case Studies Whitepapers Resources
 SUPPORT Forums Help Articles Reference Guides Updates & Maintenance
 COMPANY Mission Careers Press Room Legal Exelis, Inc
© 2014 Exelis Visual Information Solutions