NOX::LAPACK::LinearSolver< T > Class Template Reference

A simple linear solver for use by NOX::LAPACK::Group. More...

`#include <NOX_LAPACK_LinearSolver.H>`

Inheritance diagram for NOX::LAPACK::LinearSolver< T >:

[legend]
Collaboration diagram for NOX::LAPACK::LinearSolver< T >:
[legend]
List of all members.

Public Member Functions

LinearSolver (int n)
Constructor.
LinearSolver (const LinearSolver< T > &s)
Copy constructor.
~LinearSolver ()
Destructor.
LinearSolveroperator= (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.

Detailed Description

template<typename T> class NOX::LAPACK::LinearSolver< T >

A simple linear solver for use by NOX::LAPACK::Group.

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

Member Function Documentation

 template void NOX::LAPACK::LinearSolver< T >::apply ( bool trans, int ncols, const T * input, T * output ) const
 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.

 template void NOX::LAPACK::LinearSolver< T >::reset ( )
 Reset. Resets the LU factorization indicating the matrix is updated.

 template bool NOX::LAPACK::LinearSolver< T >::solve ( bool trans, int ncols, T * output )
 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.

The documentation for this class was generated from the following file:
• NOX_LAPACK_LinearSolver.H

Generated on Thu Sep 18 12:38:39 2008 for NOX by  1.3.9.1