#include <NOX_LAPACK_LinearSolver.H>
Public Member Functions | |
LinearSolver (int n) | |
Constructor. | |
LinearSolver (const LinearSolver< T > &s) | |
Copy constructor. | |
~LinearSolver () | |
Destructor. | |
LinearSolver & | operator= (const LinearSolver< T > &s) |
Assignment. | |
Matrix< T > & | getMatrix () |
Return matrix. | |
const Matrix< T > & | getMatrix () const |
Return matrix. | |
void | reset () |
Reset. | |
void | apply (bool trans, int ncols, const T *input, T *output) const |
Apply matrix. | |
bool | solve (bool trans, int ncols, T *output) |
Solve inverse. | |
Protected Attributes | |
Matrix< T > | mat |
Matrix. | |
Matrix< T > | lu |
LU factorization of matrix. | |
std::vector< int > | pivots |
Pivots. | |
bool | isValidLU |
Is LU factorization valid. | |
Teuchos::BLAS< int, T > | blas |
BLAS wrappers. | |
Teuchos::LAPACK< int, T > | lapack |
LAPACK wrappers. |
This class provides a simple linear solver class that stores a NOX::LAPACK::Matrix and provides routines to apply the matrix and solve it using BLAS and LAPACK routines. It is templated so that it can be used to solve both real and complex matrices. It also stores an LU factorization of the matrix so repeated solves are more efficient. The group should signal that the matrix has changed by calling reset().
Definition at line 66 of file NOX_LAPACK_LinearSolver.H.
NOX::LAPACK::LinearSolver< T >::LinearSolver | ( | int | n | ) | [inline] |
NOX::LAPACK::LinearSolver< T >::LinearSolver | ( | const LinearSolver< T > & | s | ) | [inline] |
NOX::LAPACK::LinearSolver< T >::~LinearSolver | ( | ) | [inline] |
NOX::LAPACK::LinearSolver< T > & NOX::LAPACK::LinearSolver< T >::operator= | ( | const LinearSolver< T > & | s | ) | [inline] |
Assignment.
Definition at line 166 of file NOX_LAPACK_LinearSolver.H.
References NOX::LAPACK::LinearSolver< T >::isValidLU, NOX::LAPACK::LinearSolver< T >::lu, NOX::LAPACK::LinearSolver< T >::mat, and NOX::LAPACK::LinearSolver< T >::pivots.
NOX::LAPACK::Matrix< T > & NOX::LAPACK::LinearSolver< T >::getMatrix | ( | ) | [inline] |
Return matrix.
Definition at line 180 of file NOX_LAPACK_LinearSolver.H.
References NOX::LAPACK::LinearSolver< T >::mat.
Referenced by LOCA::LAPACK::Group::applyComplex(), LOCA::LAPACK::Group::applyComplexInverseMultiVector(), LOCA::LAPACK::Group::applyComplexMultiVector(), LOCA::LAPACK::Group::applyComplexTranspose(), LOCA::LAPACK::Group::applyComplexTransposeInverseMultiVector(), LOCA::LAPACK::Group::applyComplexTransposeMultiVector(), NOX::LAPACK::Group::applyJacobianInverseMultiVector(), LOCA::LAPACK::Group::applyJacobianTransposeInverseMultiVector(), LOCA::LAPACK::Group::applyShiftedMatrixInverseMultiVector(), LOCA::LAPACK::Group::applyShiftedMatrixMultiVector(), LOCA::LAPACK::Group::augmentJacobianForHomotopy(), LOCA::LAPACK::Group::computeComplex(), NOX::LAPACK::Group::computeJacobian(), LOCA::LAPACK::Group::computeShiftedMatrix(), LOCA::LAPACK::Group::getComplexMatrix(), LOCA::LAPACK::Group::getJacobianMatrix(), and LOCA::LAPACK::Group::getShiftedMatrix().
const NOX::LAPACK::Matrix< T > & NOX::LAPACK::LinearSolver< T >::getMatrix | ( | ) | const [inline] |
Return matrix.
Definition at line 187 of file NOX_LAPACK_LinearSolver.H.
References NOX::LAPACK::LinearSolver< T >::mat.
void NOX::LAPACK::LinearSolver< T >::reset | ( | ) | [inline] |
Reset.
Resets the LU factorization indicating the matrix is updated.
Definition at line 194 of file NOX_LAPACK_LinearSolver.H.
References NOX::LAPACK::LinearSolver< T >::isValidLU.
Referenced by NOX::LAPACK::Group::resetIsValid(), and LOCA::LAPACK::Group::resetIsValid().
void NOX::LAPACK::LinearSolver< T >::apply | ( | bool | trans, | |
int | ncols, | |||
const T * | input, | |||
T * | output | |||
) | const [inline] |
Apply matrix.
Set trans
to true
to apply the transpose. ncols
is the number of columns in input
and output
, which should be stored column-wise.
Definition at line 201 of file NOX_LAPACK_LinearSolver.H.
References NOX::LAPACK::LinearSolver< T >::blas, and NOX::LAPACK::LinearSolver< T >::mat.
Referenced by LOCA::LAPACK::Group::applyComplex(), LOCA::LAPACK::Group::applyComplexMultiVector(), LOCA::LAPACK::Group::applyComplexTranspose(), LOCA::LAPACK::Group::applyComplexTransposeMultiVector(), LOCA::LAPACK::Group::applyShiftedMatrix(), LOCA::LAPACK::Group::applyShiftedMatrixMultiVector(), and NOX::LAPACK::Group::computeGradient().
bool NOX::LAPACK::LinearSolver< T >::solve | ( | bool | trans, | |
int | ncols, | |||
T * | output | |||
) | [inline] |
Solve inverse.
Set trans
to true
to solve the transpose. ncols
is the number of columns in output
, which should be stored column- wise. The right-hand-side is passed through output
, which is then overwritten with the result.
Definition at line 219 of file NOX_LAPACK_LinearSolver.H.
References NOX::LAPACK::LinearSolver< T >::isValidLU, NOX::LAPACK::LinearSolver< T >::lapack, NOX::LAPACK::LinearSolver< T >::lu, NOX::LAPACK::LinearSolver< T >::mat, and NOX::LAPACK::LinearSolver< T >::pivots.
Referenced by LOCA::LAPACK::Group::applyComplexInverseMultiVector(), LOCA::LAPACK::Group::applyComplexTransposeInverseMultiVector(), NOX::LAPACK::Group::applyJacobianInverseMultiVector(), LOCA::LAPACK::Group::applyJacobianTransposeInverse(), LOCA::LAPACK::Group::applyJacobianTransposeInverseMultiVector(), and LOCA::LAPACK::Group::applyShiftedMatrixInverseMultiVector().
Matrix<T> NOX::LAPACK::LinearSolver< T >::mat [protected] |
Definition at line 114 of file NOX_LAPACK_LinearSolver.H.
Referenced by NOX::LAPACK::LinearSolver< T >::apply(), NOX::LAPACK::LinearSolver< T >::getMatrix(), NOX::LAPACK::LinearSolver< T >::operator=(), and NOX::LAPACK::LinearSolver< T >::solve().
Matrix<T> NOX::LAPACK::LinearSolver< T >::lu [protected] |
LU factorization of matrix.
Definition at line 117 of file NOX_LAPACK_LinearSolver.H.
Referenced by NOX::LAPACK::LinearSolver< T >::operator=(), and NOX::LAPACK::LinearSolver< T >::solve().
std::vector<int> NOX::LAPACK::LinearSolver< T >::pivots [protected] |
Pivots.
Definition at line 120 of file NOX_LAPACK_LinearSolver.H.
Referenced by NOX::LAPACK::LinearSolver< T >::operator=(), and NOX::LAPACK::LinearSolver< T >::solve().
bool NOX::LAPACK::LinearSolver< T >::isValidLU [protected] |
Is LU factorization valid.
Definition at line 123 of file NOX_LAPACK_LinearSolver.H.
Referenced by NOX::LAPACK::LinearSolver< T >::operator=(), NOX::LAPACK::LinearSolver< T >::reset(), and NOX::LAPACK::LinearSolver< T >::solve().
Teuchos::BLAS<int,T> NOX::LAPACK::LinearSolver< T >::blas [protected] |
BLAS wrappers.
Definition at line 126 of file NOX_LAPACK_LinearSolver.H.
Referenced by NOX::LAPACK::LinearSolver< T >::apply().
Teuchos::LAPACK<int,T> NOX::LAPACK::LinearSolver< T >::lapack [protected] |
LAPACK wrappers.
Definition at line 129 of file NOX_LAPACK_LinearSolver.H.
Referenced by NOX::LAPACK::LinearSolver< T >::solve().