Teuchos Package Browser (Single Doxygen Collection)
Version of the Day
|
This class creates and provides basic support for symmetric, positive-definite dense matrices of templated type. More...
#include <Teuchos_SerialSymDenseMatrix.hpp>
Public Types | |
typedef OrdinalType | ordinalType |
Typedef for ordinal type. More... | |
typedef ScalarType | scalarType |
Typedef for scalar type. More... | |
Public Types inherited from Teuchos::DefaultBLASImpl< OrdinalType, ScalarType > | |
typedef details::GivensRotator < ScalarType >::c_type | rotg_c_type |
The type used for c in ROTG. More... | |
Protected Member Functions | |
void | scale (const ScalarType alpha) |
void | copyMat (bool inputUpper, ScalarType *inputMatrix, OrdinalType inputStride, OrdinalType numRowCols, bool outputUpper, ScalarType *outputMatrix, OrdinalType outputStride, OrdinalType startRowCol, ScalarType alpha=ScalarTraits< ScalarType >::zero()) |
void | copyUPLOMat (bool inputUpper, ScalarType *inputMatrix, OrdinalType inputStride, OrdinalType inputRows) |
void | deleteArrays () |
void | checkIndex (OrdinalType rowIndex, OrdinalType colIndex=0) const |
Static Protected Member Functions | |
static ScalarType * | allocateValues (const OrdinalType numRows, const OrdinalType numCols) |
Protected Attributes | |
OrdinalType | numRowCols_ = 0 |
OrdinalType | stride_ = 0 |
bool | valuesCopied_ = false |
ScalarType * | values_ = nullptr |
bool | upper_ = false |
char UPLO_ | L |
Protected Attributes inherited from Teuchos::CompObject | |
Flops * | flopCounter_ |
Related Functions | |
(Note that these are not member functions.) | |
template<typename OrdinalType , typename ScalarType > | |
void | symMatTripleProduct (ETransp transw, const ScalarType alpha, const SerialSymDenseMatrix< OrdinalType, ScalarType > &A, const SerialDenseMatrix< OrdinalType, ScalarType > &W, SerialSymDenseMatrix< OrdinalType, ScalarType > &B) |
A templated, non-member, helper function for computing the matrix triple-product: B = alpha*W^T*A*W or B = alpha*W*A*W^T. More... | |
Constructor/Destructor Methods | |
SerialSymDenseMatrix ()=default | |
Default constructor; defines a zero size object. More... | |
SerialSymDenseMatrix (OrdinalType numRowsCols, bool zeroOut=true) | |
Basic constructor; defines a matrix of numRowsCols size and (optionally) initializes it. More... | |
SerialSymDenseMatrix (DataAccess CV, bool upper, ScalarType *values, OrdinalType stride, OrdinalType numRowsCols) | |
Set object values from two-dimensional array. More... | |
SerialSymDenseMatrix (const SerialSymDenseMatrix< OrdinalType, ScalarType > &Source) | |
Teuchos::SerialSymDenseMatrix copy constructor. More... | |
SerialSymDenseMatrix (DataAccess CV, const SerialSymDenseMatrix< OrdinalType, ScalarType > &Source, OrdinalType numRowCols, OrdinalType startRowCol=0) | |
Submatrix Copy Constructor. More... | |
virtual | ~SerialSymDenseMatrix () |
Teuchos::SerialSymDenseMatrix destructor. More... | |
Shaping Methods | |
int | shape (OrdinalType numRowsCols) |
Set dimensions of a Teuchos::SerialSymDenseMatrix object; init values to zero. More... | |
int | shapeUninitialized (OrdinalType numRowsCols) |
Set dimensions of a Teuchos::SerialSymDenseMatrix object; don't initialize values. More... | |
int | reshape (OrdinalType numRowsCols) |
Reshape a Teuchos::SerialSymDenseMatrix object. More... | |
void | setLower () |
Specify that the lower triangle of the this matrix should be used. More... | |
void | setUpper () |
Specify that the upper triangle of the this matrix should be used. More... | |
Set methods. | |
SerialSymDenseMatrix < OrdinalType, ScalarType > & | operator= (const SerialSymDenseMatrix< OrdinalType, ScalarType > &Source) |
Copies values from one matrix to another. More... | |
SerialSymDenseMatrix < OrdinalType, ScalarType > & | assign (const SerialSymDenseMatrix< OrdinalType, ScalarType > &Source) |
Copies values from one matrix to another. More... | |
SerialSymDenseMatrix < OrdinalType, ScalarType > & | operator= (const ScalarType value) |
Set all values in the matrix to a constant value. More... | |
int | putScalar (const ScalarType value=Teuchos::ScalarTraits< ScalarType >::zero(), bool fullMatrix=false) |
Set all values in the matrix to a constant value. More... | |
void | swap (SerialSymDenseMatrix< OrdinalType, ScalarType > &B) |
Swap values between this matrix and incoming matrix. More... | |
int | random (const ScalarType bias=0.1 *Teuchos::ScalarTraits< ScalarType >::one()) |
Set all values in the active area (upper/lower triangle) of this matrix to be random numbers. More... | |
Accessor methods. | |
ScalarType & | operator() (OrdinalType rowIndex, OrdinalType colIndex) |
Element access method (non-const). More... | |
const ScalarType & | operator() (OrdinalType rowIndex, OrdinalType colIndex) const |
Element access method (const). More... | |
ScalarType * | values () const |
Returns the pointer to the ScalarType data array contained in the object. More... | |
Query methods | |
bool | upper () const |
Returns true if upper triangular part of this matrix has and will be used. More... | |
char | UPLO () const |
Returns character value of UPLO used by LAPACK routines. More... | |
Mathematical Methods | |
SerialSymDenseMatrix < OrdinalType, ScalarType > & | operator*= (const ScalarType alpha) |
Inplace scalar-matrix product A = alpha*A . More... | |
SerialSymDenseMatrix < OrdinalType, ScalarType > & | operator+= (const SerialSymDenseMatrix< OrdinalType, ScalarType > &Source) |
Add another matrix to this matrix. More... | |
SerialSymDenseMatrix < OrdinalType, ScalarType > & | operator-= (const SerialSymDenseMatrix< OrdinalType, ScalarType > &Source) |
Subtract another matrix from this matrix. More... | |
Comparison methods. | |
bool | operator== (const SerialSymDenseMatrix< OrdinalType, ScalarType > &Operand) const |
Equality of two matrices. More... | |
bool | operator!= (const SerialSymDenseMatrix< OrdinalType, ScalarType > &Operand) const |
Inequality of two matrices. More... | |
Attribute methods. | |
OrdinalType | numRows () const |
Returns the row dimension of this matrix. More... | |
OrdinalType | numCols () const |
Returns the column dimension of this matrix. More... | |
OrdinalType | stride () const |
Returns the stride between the columns of this matrix in memory. More... | |
bool | empty () const |
Returns whether this matrix is empty. More... | |
Norm methods. | |
ScalarTraits< ScalarType > ::magnitudeType | normOne () const |
Returns the 1-norm of the matrix. More... | |
ScalarTraits< ScalarType > ::magnitudeType | normInf () const |
Returns the Infinity-norm of the matrix. More... | |
ScalarTraits< ScalarType > ::magnitudeType | normFrobenius () const |
Returns the Frobenius-norm of the matrix. More... | |
I/O methods. | |
virtual std::ostream & | print (std::ostream &os) const |
Print method. Defines the behavior of the std::ostream << operator. More... | |
This class creates and provides basic support for symmetric, positive-definite dense matrices of templated type.
The Teuchos_SerialSymDenseMatrix class enables the construction and use of symmetric, positive-definite, dense matrices of templated type. The Teuchos::SerialSymDenseMatrix class is intended to provide full-featured support for solving linear and eigen system problems for symmetric positive-definite (SPD) matrices. It is written on top of BLAS and LAPACK and thus has excellent performance and numerical capabilities. Using this class, one can either perform simple factorizations and solves or apply all the tricks available in LAPACK to get the best possible solution for very ill-conditioned problems. <b>Teuchos::SerialSymDenseMatrix vs. Teuchos::LAPACK</b> The Teuchos::LAPACK class provides access to most of the same functionality as Teuchos::SerialSymDenseMatrix. The primary difference is that Teuchos::LAPACK is a "thin" layer on top of LAPACK and Teuchos::SerialSymDenseMatrix attempts to provide easy access to the more sophisticated aspects of solving dense linear and eigensystems.
Constructing Teuchos::SerialSymDenseMatrix Objects
There are three Teuchos::SerialSymDenseMatrix constructors. The first constructs a zero-sized object which should be made to appropriate length using the Shape() or Reshape() functions and then filled with the [] or () operators. The second is a constructor that accepts user data as a 2D array, the third is a copy constructor. The second constructor has two data access modes (specified by the Teuchos::DataAccess argument):
Extracting Data from Teuchos::SerialSymDenseMatrix Objects
Once a Teuchos::SerialSymDenseMatrix is constructed, it is possible to view the data via access functions.
Vector and Utility Functions
Once a Teuchos::SerialSymDenseMatrix is constructed, several mathematical functions can be applied to the object. Specifically:
Definition at line 91 of file Teuchos_SerialSymDenseMatrix.hpp.
typedef OrdinalType Teuchos::SerialSymDenseMatrix< OrdinalType, ScalarType >::ordinalType |
Typedef for ordinal type.
Definition at line 96 of file Teuchos_SerialSymDenseMatrix.hpp.
typedef ScalarType Teuchos::SerialSymDenseMatrix< OrdinalType, ScalarType >::scalarType |
Typedef for scalar type.
Definition at line 98 of file Teuchos_SerialSymDenseMatrix.hpp.
|
default |
Default constructor; defines a zero size object.
Teuchos::SerialSymDenseMatrix objects defined by the default constructor should be sized with the Shape() or Reshape() functions. Values should be defined by using the [] or ()operators.
Note: By default the active part of the matrix is assumed to be the lower triangular part. To set the upper part as active, call SetUpper(). See Detailed Description section for further discussion.
Teuchos::SerialSymDenseMatrix< OrdinalType, ScalarType >::SerialSymDenseMatrix | ( | OrdinalType | numRowsCols, |
bool | zeroOut = true |
||
) |
Basic constructor; defines a matrix of numRowsCols
size and (optionally) initializes it.
numRowsCols | - Number of rows and columns in the matrix. |
zeroOut | - Initializes values to 0 if true (default) |
Creates a shaped matrix with numRowsCols
rows and cols. All values are initialized to 0 when zeroOut
is true. Values of this matrix should be set using the [] or the () operators.
Definition at line 419 of file Teuchos_SerialSymDenseMatrix.hpp.
Teuchos::SerialSymDenseMatrix< OrdinalType, ScalarType >::SerialSymDenseMatrix | ( | DataAccess | CV, |
bool | upper, | ||
ScalarType * | values, | ||
OrdinalType | stride, | ||
OrdinalType | numRowsCols | ||
) |
Set object values from two-dimensional array.
CV | - Enumerated type set to Teuchos::Copy or Teuchos::View. |
values | - Pointer to an array of ScalarType. The first column starts at values , the second at values+stride , etc. |
stride | - The stride between the columns of the matrix in memory. |
numRowsCols | - Number of rows and columns in the matrix. |
Definition at line 431 of file Teuchos_SerialSymDenseMatrix.hpp.
Teuchos::SerialSymDenseMatrix< OrdinalType, ScalarType >::SerialSymDenseMatrix | ( | const SerialSymDenseMatrix< OrdinalType, ScalarType > & | Source | ) |
Teuchos::SerialSymDenseMatrix copy constructor.
Definition at line 452 of file Teuchos_SerialSymDenseMatrix.hpp.
Teuchos::SerialSymDenseMatrix< OrdinalType, ScalarType >::SerialSymDenseMatrix | ( | DataAccess | CV, |
const SerialSymDenseMatrix< OrdinalType, ScalarType > & | Source, | ||
OrdinalType | numRowCols, | ||
OrdinalType | startRowCol = 0 |
||
) |
Submatrix Copy Constructor.
CV | - Enumerated type set to Teuchos::Copy or Teuchos::View. |
Source | - Reference to another dense matrix from which values are to be copied. |
numRowCols | - The number of rows and columns in this matrix. |
startRowCol | - The row and column of Source from which the submatrix copy should start. |
Creates a shaped matrix with numRowCols
rows and columns, which is a submatrix of Source
. If startRowCol
are not given, then the submatrix is the leading submatrix of Source
.
Definition at line 478 of file Teuchos_SerialSymDenseMatrix.hpp.
|
virtual |
Teuchos::SerialSymDenseMatrix destructor.
Definition at line 498 of file Teuchos_SerialSymDenseMatrix.hpp.
int Teuchos::SerialSymDenseMatrix< OrdinalType, ScalarType >::shape | ( | OrdinalType | numRowsCols | ) |
Set dimensions of a Teuchos::SerialSymDenseMatrix object; init values to zero.
numRowsCols | - Number of rows and columns in object. |
Allows user to define the dimensions of a Teuchos::SerialSymDenseMatrix at any point. This function can be called at any point after construction. Any values that were previously in this object are destroyed and the resized matrix starts off with all zero values.
Definition at line 508 of file Teuchos_SerialSymDenseMatrix.hpp.
int Teuchos::SerialSymDenseMatrix< OrdinalType, ScalarType >::shapeUninitialized | ( | OrdinalType | numRowsCols | ) |
Set dimensions of a Teuchos::SerialSymDenseMatrix object; don't initialize values.
numRowsCols | - Number of rows and columns in object. |
Allows user to define the dimensions of a Teuchos::SerialSymDenseMatrix at any point. This function can be called at any point after construction. Any values that were previously in this object are destroyed. The resized matrix has uninitialized values.
Definition at line 520 of file Teuchos_SerialSymDenseMatrix.hpp.
int Teuchos::SerialSymDenseMatrix< OrdinalType, ScalarType >::reshape | ( | OrdinalType | numRowsCols | ) |
Reshape a Teuchos::SerialSymDenseMatrix object.
numRowsCols | - Number of rows and columns in object. |
Allows user to define the dimensions of a Teuchos::SerialSymDenseMatrix at any point. This function can be called at any point after construction. Any values that were previously in this object are copied into the new shape. If the new shape is smaller than the original, the upper left portion of the original matrix (the principal submatrix) is copied to the new matrix.
Definition at line 531 of file Teuchos_SerialSymDenseMatrix.hpp.
void Teuchos::SerialSymDenseMatrix< OrdinalType, ScalarType >::setLower | ( | ) |
Specify that the lower triangle of the this matrix should be used.
Definition at line 558 of file Teuchos_SerialSymDenseMatrix.hpp.
void Teuchos::SerialSymDenseMatrix< OrdinalType, ScalarType >::setUpper | ( | ) |
Specify that the upper triangle of the this matrix should be used.
Definition at line 569 of file Teuchos_SerialSymDenseMatrix.hpp.
SerialSymDenseMatrix< OrdinalType, ScalarType > & Teuchos::SerialSymDenseMatrix< OrdinalType, ScalarType >::operator= | ( | const SerialSymDenseMatrix< OrdinalType, ScalarType > & | Source | ) |
Copies values from one matrix to another.
The operator= copies the values from one existing SerialSymDenseMatrix to another. If Source
is a view (i.e. CV = Teuchos::View), then this method will return a view. Otherwise, it will return a copy of Source
. this object will be resized if it is not large enough to copy Source
into.
Definition at line 659 of file Teuchos_SerialSymDenseMatrix.hpp.
SerialSymDenseMatrix< OrdinalType, ScalarType > & Teuchos::SerialSymDenseMatrix< OrdinalType, ScalarType >::assign | ( | const SerialSymDenseMatrix< OrdinalType, ScalarType > & | Source | ) |
Copies values from one matrix to another.
Copies the values from one existing SerialSymDenseMatrix to another if the dimension of both matrices are the same. If not, this matrix will be returned unchanged.
Definition at line 751 of file Teuchos_SerialSymDenseMatrix.hpp.
|
inline |
Set all values in the matrix to a constant value.
value | - Value to use; |
Definition at line 234 of file Teuchos_SerialSymDenseMatrix.hpp.
int Teuchos::SerialSymDenseMatrix< OrdinalType, ScalarType >::putScalar | ( | const ScalarType | value = Teuchos::ScalarTraits<ScalarType>::zero() , |
bool | fullMatrix = false |
||
) |
Set all values in the matrix to a constant value.
value | - Value to use; zero if none specified. |
fullMatrix | - set full matrix entries to value , not just active portion of symmetric matrix. |
Definition at line 580 of file Teuchos_SerialSymDenseMatrix.hpp.
void Teuchos::SerialSymDenseMatrix< OrdinalType, ScalarType >::swap | ( | SerialSymDenseMatrix< OrdinalType, ScalarType > & | B | ) |
Swap values between this matrix and incoming matrix.
Swaps pointers and associated state without copying the matrix data.
Definition at line 613 of file Teuchos_SerialSymDenseMatrix.hpp.
int Teuchos::SerialSymDenseMatrix< OrdinalType, ScalarType >::random | ( | const ScalarType | bias = 0.1*Teuchos::ScalarTraits<ScalarType>::one() | ) |
Set all values in the active area (upper/lower triangle) of this matrix to be random numbers.
Definition at line 625 of file Teuchos_SerialSymDenseMatrix.hpp.
|
inline |
Element access method (non-const).
Returns the element in the ith row and jth column if A(i,j) is specified.
rowIndex
row and colIndex
column. rowIndex
and colIndex
will only be checked if Teuchos is configured with –enable-teuchos-abc. Definition at line 773 of file Teuchos_SerialSymDenseMatrix.hpp.
|
inline |
Element access method (const).
Returns the element in the ith row and jth column if A(i,j) is specified.
rowIndex
row and colIndex
column. rowIndex
and colIndex
will only be checked if Teuchos is configured with –enable-teuchos-abc. Definition at line 795 of file Teuchos_SerialSymDenseMatrix.hpp.
|
inline |
Returns the pointer to the ScalarType data array contained in the object.
Definition at line 283 of file Teuchos_SerialSymDenseMatrix.hpp.
|
inline |
Returns true if upper triangular part of this matrix has and will be used.
Definition at line 291 of file Teuchos_SerialSymDenseMatrix.hpp.
|
inline |
Returns character value of UPLO used by LAPACK routines.
Definition at line 294 of file Teuchos_SerialSymDenseMatrix.hpp.
SerialSymDenseMatrix< OrdinalType, ScalarType > & Teuchos::SerialSymDenseMatrix< OrdinalType, ScalarType >::operator*= | ( | const ScalarType | alpha | ) |
SerialSymDenseMatrix< OrdinalType, ScalarType > & Teuchos::SerialSymDenseMatrix< OrdinalType, ScalarType >::operator+= | ( | const SerialSymDenseMatrix< OrdinalType, ScalarType > & | Source | ) |
Add another matrix to this matrix.
Add Source
to this if the dimension of both matrices are the same. If not, this matrix will be returned unchanged.
Definition at line 727 of file Teuchos_SerialSymDenseMatrix.hpp.
SerialSymDenseMatrix< OrdinalType, ScalarType > & Teuchos::SerialSymDenseMatrix< OrdinalType, ScalarType >::operator-= | ( | const SerialSymDenseMatrix< OrdinalType, ScalarType > & | Source | ) |
Subtract another matrix from this matrix.
Subtract Source
from this if the dimension of both matrices are the same. If not, this matrix will be returned unchanged.
Definition at line 739 of file Teuchos_SerialSymDenseMatrix.hpp.
bool Teuchos::SerialSymDenseMatrix< OrdinalType, ScalarType >::operator== | ( | const SerialSymDenseMatrix< OrdinalType, ScalarType > & | Operand | ) | const |
Equality of two matrices.
Operand
are of the same shape (rows / columns) and have the same entries in the active (upper / lower triangular) area of the matrix, else False will be returned. Definition at line 903 of file Teuchos_SerialSymDenseMatrix.hpp.
bool Teuchos::SerialSymDenseMatrix< OrdinalType, ScalarType >::operator!= | ( | const SerialSymDenseMatrix< OrdinalType, ScalarType > & | Operand | ) | const |
Inequality of two matrices.
Operand
of not of the same shape (rows / columns) or don't have the same entries in the active (upper / lower triangular), else False will be returned. Definition at line 923 of file Teuchos_SerialSymDenseMatrix.hpp.
|
inline |
Returns the row dimension of this matrix.
Definition at line 344 of file Teuchos_SerialSymDenseMatrix.hpp.
|
inline |
Returns the column dimension of this matrix.
Definition at line 347 of file Teuchos_SerialSymDenseMatrix.hpp.
|
inline |
Returns the stride between the columns of this matrix in memory.
Definition at line 350 of file Teuchos_SerialSymDenseMatrix.hpp.
|
inline |
Returns whether this matrix is empty.
Definition at line 353 of file Teuchos_SerialSymDenseMatrix.hpp.
ScalarTraits< ScalarType >::magnitudeType Teuchos::SerialSymDenseMatrix< OrdinalType, ScalarType >::normOne | ( | ) | const |
Returns the 1-norm of the matrix.
Definition at line 821 of file Teuchos_SerialSymDenseMatrix.hpp.
ScalarTraits< ScalarType >::magnitudeType Teuchos::SerialSymDenseMatrix< OrdinalType, ScalarType >::normInf | ( | ) | const |
Returns the Infinity-norm of the matrix.
Definition at line 827 of file Teuchos_SerialSymDenseMatrix.hpp.
ScalarTraits< ScalarType >::magnitudeType Teuchos::SerialSymDenseMatrix< OrdinalType, ScalarType >::normFrobenius | ( | ) | const |
Returns the Frobenius-norm of the matrix.
Definition at line 870 of file Teuchos_SerialSymDenseMatrix.hpp.
|
virtual |
Print method. Defines the behavior of the std::ostream << operator.
Definition at line 982 of file Teuchos_SerialSymDenseMatrix.hpp.
|
protected |
Definition at line 933 of file Teuchos_SerialSymDenseMatrix.hpp.
|
protected |
Definition at line 1034 of file Teuchos_SerialSymDenseMatrix.hpp.
|
protected |
Definition at line 1117 of file Teuchos_SerialSymDenseMatrix.hpp.
|
protected |
Definition at line 1023 of file Teuchos_SerialSymDenseMatrix.hpp.
|
inlineprotected |
Definition at line 1013 of file Teuchos_SerialSymDenseMatrix.hpp.
|
inlinestaticprotected |
Definition at line 396 of file Teuchos_SerialSymDenseMatrix.hpp.
|
related |
A templated, non-member, helper function for computing the matrix triple-product: B = alpha*W^T*A*W or B = alpha*W*A*W^T.
transw | - [in] Compute B = alpha*W^T*A*W if transw = Teuchos::TRANS, else compute B = alpha*W*A*W^T if transw = Teuchos::NO_TRANS. |
alpha | - [in] The scaling factor. |
A | - [in] SerialSymDenseMatrix |
W | - [in] SerialDenseMatrix |
B | - [out] SerialSymDenseMatrix |
Teuchos::symMatTripleProduct<int,double>( Teuchos::TRANS, alpha, A, W, B )
Definition at line 44 of file Teuchos_SerialDenseHelpers.hpp.
|
protected |
Definition at line 406 of file Teuchos_SerialSymDenseMatrix.hpp.
|
protected |
Definition at line 407 of file Teuchos_SerialSymDenseMatrix.hpp.
|
protected |
Definition at line 408 of file Teuchos_SerialSymDenseMatrix.hpp.
|
protected |
Definition at line 409 of file Teuchos_SerialSymDenseMatrix.hpp.
|
protected |
Definition at line 410 of file Teuchos_SerialSymDenseMatrix.hpp.
|
protected |
Definition at line 411 of file Teuchos_SerialSymDenseMatrix.hpp.