linalg
1.4.3
A linear algebra library that provides a user-friendly interface to several BLAS and LAPACK routines.
|
Solves a triangular system of equations. More...
Private Member Functions | |
subroutine | solve_tri_mtx (lside, upper, trans, nounit, alpha, a, b, err) |
Solves one of the matrix equations: op(A) * X = alpha * B, or X * op(A) = alpha * B, where A is a triangular matrix. More... | |
subroutine | solve_tri_vec (upper, trans, nounit, a, x, err) |
Solves the system of equations: op(A) * X = B, where A is a triangular matrix. More... | |
Solves a triangular system of equations.
Definition at line 92 of file linalg_solve.f90.
|
private |
Solves one of the matrix equations: op(A) * X = alpha * B, or X * op(A) = alpha * B, where A is a triangular matrix.
[in] | lside | Set to true to solve op(A) * X = alpha * B; else, set to false to solve X * op(A) = alpha * B. |
[in] | upper | Set to true if A is an upper triangular matrix; else, set to false if A is a lower triangular matrix. |
[in] | trans | Set to true if op(A) = A**T; else, set to false if op(A) = A. |
[in] | nounit | Set to true if A is not a unit-diagonal matrix (ones on every diagonal element); else, set to false if A is a unit-diagonal matrix. |
[in] | alpha | The scalar multiplier to B. |
[in] | a | If lside is true, the M-by-M triangular matrix on which to operate; else, if lside is false, the N-by-N triangular matrix on which to operate. |
[in,out] | b | On input, the M-by-N right-hand-side. On output, the M-by-N solution. |
[out] | err | An optional errors-based object that if provided can be used to retrieve information relating to any errors encountered during execution. If not provided, a default implementation of the errors class is used internally to provide error handling. Possible errors and warning messages that may be encountered are as follows.
|
Definition at line 496 of file linalg_solve.f90.
|
private |
Solves the system of equations: op(A) * X = B, where A is a triangular matrix.
[in] | upper | Set to true if A is an upper triangular matrix; else, set to false if A is a lower triangular matrix. |
[in] | trans | Set to true if op(A) = A**T; else, set to false if op(A) = A. |
[in] | nounit | Set to true if A is not a unit-diagonal matrix (ones on every diagonal element); else, set to false if A is a unit-diagonal matrix. |
[in] | a | The N-by-N triangular matrix. |
[in,out] | x | On input, the N-element right-hand-side array. On output, the N-element solution array. |
[out] | err | An optional errors-based object that if provided can be used to retrieve information relating to any errors encountered during execution. If not provided, a default implementation of the errors class is used internally to provide error handling. Possible errors and warning messages that may be encountered are as follows.
|
Definition at line 602 of file linalg_solve.f90.