ergo
|
Normal matrix. More...
#include <MatrixBase.h>
Public Types | |
typedef VectorGeneral< Treal, typename Tmatrix::VectorType > | VectorType |
Public Member Functions | |
MatrixGeneral () | |
Default constructor More... | |
MatrixGeneral (const MatrixGeneral< Treal, Tmatrix > &matr) | |
Copy constructor More... | |
MatrixGeneral (const MatrixSymmetric< Treal, Tmatrix > &symm) | |
Copy from symmetric matrix constructor More... | |
MatrixGeneral (const MatrixTriangular< Treal, Tmatrix > &triang) | |
Copy from triangular matrix constructor More... | |
void | assignFromFull (std::vector< Treal > const &fullMat) |
void | fullMatrix (std::vector< Treal > &fullMat) const |
void | fullMatrix (std::vector< Treal > &fullMat, std::vector< int > const &rowInversePermutation, std::vector< int > const &colInversePermutation) const |
Save matrix as full matrix. More... | |
void | assign_from_sparse (std::vector< int > const &rowind, std::vector< int > const &colind, std::vector< Treal > const &values, SizesAndBlocks const &newRows, SizesAndBlocks const &newCols) |
Assign from sparse matrix given by three arrays. More... | |
void | assign_from_sparse (std::vector< int > const &rowind, std::vector< int > const &colind, std::vector< Treal > const &values, std::vector< int > const &rowPermutation, std::vector< int > const &colPermutation) |
Same as above, except taking two additional arguments specifying the permutation of rows and columns. More... | |
void | assign_from_sparse (std::vector< int > const &rowind, std::vector< int > const &colind, std::vector< Treal > const &values, SizesAndBlocks const &newRows, SizesAndBlocks const &newCols, std::vector< int > const &rowPermutation, std::vector< int > const &colPermutation) |
Same as above, except not assuming that sizes and blocks are set. More... | |
void | get_values (std::vector< int > const &rowind, std::vector< int > const &colind, std::vector< Treal > &values) const |
Get values given by row and column index lists. More... | |
void | get_values (std::vector< int > const &rowind, std::vector< int > const &colind, std::vector< Treal > &values, std::vector< int > const &rowPermutation, std::vector< int > const &colPermutation) const |
Same as above, except taking two additional arguments specifying the permutation of rows and columns. More... | |
void | get_all_values (std::vector< int > &rowind, std::vector< int > &colind, std::vector< Treal > &values) const |
Get all values and corresponding row and column index lists, in matrix. More... | |
void | get_all_values (std::vector< int > &rowind, std::vector< int > &colind, std::vector< Treal > &values, std::vector< int > const &rowInversePermutation, std::vector< int > const &colInversePermutation) const |
Same as above, except taking two additional arguments specifying the permutation of rows and columns. More... | |
MatrixGeneral< Treal, Tmatrix > & | operator= (const MatrixGeneral< Treal, Tmatrix > &mat) |
MatrixGeneral< Treal, Tmatrix > & | operator= (const Xtrans< MatrixGeneral< Treal, Tmatrix > > &mt) |
MatrixGeneral< Treal, Tmatrix > & | operator= (const MatrixSymmetric< Treal, Tmatrix > &symm) |
MatrixGeneral< Treal, Tmatrix > & | operator= (const MatrixTriangular< Treal, Tmatrix > &triang) |
MatrixGeneral< Treal, Tmatrix > & | operator= (int const k) |
Treal | frob () const |
Treal | eucl (Treal const requestedAccuracy, int maxIter=-1) const |
void | thresh (Treal const threshold, normType const norm) |
void | frob_thresh (Treal threshold) |
Treal | eucl_thresh (Treal const threshold) |
void | gershgorin (Treal &lmin, Treal &lmax) |
size_t | nnz () const |
size_t | nvalues () const |
void | write_to_buffer (void *buffer, const int n_bytes) const |
void | read_from_buffer (void *buffer, const int n_bytes) |
MatrixGeneral< Treal, Tmatrix > & | operator= (const XYZ< Treal, MatrixGeneral< Treal, Tmatrix >, MatrixGeneral< Treal, Tmatrix > > &smm) |
C = alpha * op(A) * op(B) More... | |
MatrixGeneral< Treal, Tmatrix > & | operator= (const XY< MatrixGeneral< Treal, Tmatrix >, MatrixGeneral< Treal, Tmatrix > > &mm) |
C = op(A) * op(B) More... | |
MatrixGeneral< Treal, Tmatrix > & | operator+= (const XYZ< Treal, MatrixGeneral< Treal, Tmatrix >, MatrixGeneral< Treal, Tmatrix > > &smm) |
C += alpha * op(A) * op(B) More... | |
MatrixGeneral< Treal, Tmatrix > & | operator= (const XYZpUV< Treal, MatrixGeneral< Treal, Tmatrix >, MatrixGeneral< Treal, Tmatrix >, Treal, MatrixGeneral< Treal, Tmatrix > > &smmpsm) |
C = alpha * op(A) * op(B) + beta * C More... | |
MatrixGeneral< Treal, Tmatrix > & | operator= (XpY< MatrixGeneral< Treal, Tmatrix >, MatrixGeneral< Treal, Tmatrix > > const &mpm) |
C = A + B. More... | |
MatrixGeneral< Treal, Tmatrix > & | operator+= (MatrixGeneral< Treal, Tmatrix > const &A) |
B += A. More... | |
MatrixGeneral< Treal, Tmatrix > & | operator-= (MatrixGeneral< Treal, Tmatrix > const &A) |
MatrixGeneral< Treal, Tmatrix > & | operator+= (XY< Treal, MatrixGeneral< Treal, Tmatrix > > const &sm) |
B += alpha * A. More... | |
MatrixGeneral< Treal, Tmatrix > & | operator= (const XYZ< Treal, MatrixSymmetric< Treal, Tmatrix >, MatrixGeneral< Treal, Tmatrix > > &smm) |
C = alpha * A * B : A is symmetric. More... | |
MatrixGeneral< Treal, Tmatrix > & | operator= (const XY< MatrixSymmetric< Treal, Tmatrix >, MatrixGeneral< Treal, Tmatrix > > &mm) |
C = A * B : A is symmetric. More... | |
MatrixGeneral< Treal, Tmatrix > & | operator+= (const XYZ< Treal, MatrixSymmetric< Treal, Tmatrix >, MatrixGeneral< Treal, Tmatrix > > &smm) |
C += alpha * A * B : A is symmetric. More... | |
MatrixGeneral< Treal, Tmatrix > & | operator= (const XYZpUV< Treal, MatrixSymmetric< Treal, Tmatrix >, MatrixGeneral< Treal, Tmatrix >, Treal, MatrixGeneral< Treal, Tmatrix > > &smmpsm) |
C = alpha * A * B + beta * C : A is symmetric. More... | |
MatrixGeneral< Treal, Tmatrix > & | operator= (const XYZ< Treal, MatrixGeneral< Treal, Tmatrix >, MatrixSymmetric< Treal, Tmatrix > > &smm) |
C = alpha * B * A : A is symmetric. More... | |
MatrixGeneral< Treal, Tmatrix > & | operator= (const XY< MatrixGeneral< Treal, Tmatrix >, MatrixSymmetric< Treal, Tmatrix > > &mm) |
C = B * A : A is symmetric. More... | |
MatrixGeneral< Treal, Tmatrix > & | operator+= (const XYZ< Treal, MatrixGeneral< Treal, Tmatrix >, MatrixSymmetric< Treal, Tmatrix > > &smm) |
C += alpha * B * A : A is symmetric. More... | |
MatrixGeneral< Treal, Tmatrix > & | operator= (const XYZpUV< Treal, MatrixGeneral< Treal, Tmatrix >, MatrixSymmetric< Treal, Tmatrix >, Treal, MatrixGeneral< Treal, Tmatrix > > &smmpsm) |
C = alpha * B * A + beta * C : A is symmetric. More... | |
MatrixGeneral< Treal, Tmatrix > & | operator= (const XYZ< Treal, MatrixSymmetric< Treal, Tmatrix >, MatrixSymmetric< Treal, Tmatrix > > &smm) |
C = alpha * A * B : A and B are symmetric. More... | |
MatrixGeneral< Treal, Tmatrix > & | operator= (const XY< MatrixSymmetric< Treal, Tmatrix >, MatrixSymmetric< Treal, Tmatrix > > &mm) |
C = A * B : A and B are symmetric. More... | |
MatrixGeneral< Treal, Tmatrix > & | operator+= (const XYZ< Treal, MatrixSymmetric< Treal, Tmatrix >, MatrixSymmetric< Treal, Tmatrix > > &smm) |
C += alpha * A * B : A and B are symmetric. More... | |
MatrixGeneral< Treal, Tmatrix > & | operator= (const XYZpUV< Treal, MatrixSymmetric< Treal, Tmatrix >, MatrixSymmetric< Treal, Tmatrix >, Treal, MatrixGeneral< Treal, Tmatrix > > &smmpsm) |
C = alpha * A * B + beta * C : A and B are symmetric. More... | |
MatrixGeneral< Treal, Tmatrix > & | operator= (const XYZ< Treal, MatrixTriangular< Treal, Tmatrix >, MatrixGeneral< Treal, Tmatrix > > &smm) |
B = alpha * op(A) * B : A is upper triangular. More... | |
MatrixGeneral< Treal, Tmatrix > & | operator= (const XYZ< Treal, MatrixGeneral< Treal, Tmatrix >, MatrixTriangular< Treal, Tmatrix > > &smm) |
B = alpha * B * op(A) : A is upper triangular. More... | |
void | random () |
void | randomZeroStructure (Treal probabilityBeingZero) |
template<typename TRule > | |
void | setElementsByRule (TRule &rule) |
std::string | obj_type_id () const |
Static Public Member Functions | |
static Treal | frob_diff (const MatrixGeneral< Treal, Tmatrix > &A, const MatrixGeneral< Treal, Tmatrix > &B) |
static Treal | trace_ab (const MatrixGeneral< Treal, Tmatrix > &A, const MatrixGeneral< Treal, Tmatrix > &B) |
static Treal | trace_aTb (const MatrixGeneral< Treal, Tmatrix > &A, const MatrixGeneral< Treal, Tmatrix > &B) |
Protected Member Functions | |
void | writeToFileProt (std::ofstream &file) const |
void | readFromFileProt (std::ifstream &file) |
Normal matrix.
This class belongs to the matrix API
Treal: Type for real numbers
Tmatrix: The matrix class
Tperm: Permutation used in the matrix class
typedef VectorGeneral<Treal, typename Tmatrix::VectorType> mat::MatrixGeneral< Treal, Tmatrix >::VectorType |
|
inline |
Default constructor
|
inlineexplicit |
Copy constructor
|
inlineexplicit |
Copy from symmetric matrix constructor
|
inlineexplicit |
Copy from triangular matrix constructor
|
inline |
Assign from sparse matrix given by three arrays.
The arrays contain row indices, column indices and values. The indices start at zero. nval is the length of the three arrays.
Referenced by mat::MatrixGeneral< Treal, Tmatrix >::assign_from_sparse(), and mainFun().
|
inline |
Same as above, except not assuming that sizes and blocks are set.
References mat::MatrixGeneral< Treal, Tmatrix >::assign_from_sparse().
|
inline |
Same as above, except taking two additional arguments specifying the permutation of rows and columns.
Also assuming that sizes and blocks are already known
References mat::MatrixBase< Treal, Tmatrix >::getPermutedIndexes().
|
inline |
Treal mat::MatrixGeneral< Treal, Tmatrix >::eucl | ( | Treal const | requestedAccuracy, |
int | maxIter = -1 |
||
) | const |
References cols, mat::arn::LanczosLargestMagnitudeEig< Treal, Tmatrix, Tvector >::getLargestMagnitudeEig(), mat::Interval< Treal >::length(), mat::Interval< Treal >::low(), mat::Interval< Treal >::midPoint(), mat::VectorGeneral< Treal, Tvector >::rand(), and mat::VectorGeneral< Treal, Tvector >::resetSizesAndBlocks().
Referenced by mainFun().
Treal mat::MatrixGeneral< Treal, Tmatrix >::eucl_thresh | ( | Treal const | threshold | ) |
References mat::EuclTruncationBase< Tmatrix, Treal >::run().
Referenced by mainFun(), and SCF_general::SCF_general().
|
inline |
Referenced by SCF_restricted::get_error_measure(), SCF_unrestricted::get_error_measure(), and mainFun().
|
inlinestatic |
|
inline |
Referenced by mainFun().
|
inline |
|
inline |
Save matrix as full matrix.
Whole matrix is written in columnwise order. Both lower and upper triangle. Permutation is used.
References mat::MatrixGeneral< Treal, Tmatrix >::get_all_values().
|
inline |
|
inline |
Get all values and corresponding row and column index lists, in matrix.
nval is the length of the three arrays and is preferably computed with nvalues() before hand. Returns the number of values.
Referenced by mat::MatrixGeneral< Treal, Tmatrix >::fullMatrix().
|
inline |
Same as above, except taking two additional arguments specifying the permutation of rows and columns.
Note, however, that this permutation is the inverse permutation compared to the permutations provided in the functions "assign_from_sparse", "add_values", and "get_values"
References mat::MatrixBase< Treal, Tmatrix >::getPermutedIndexes().
|
inline |
|
inline |
Same as above, except taking two additional arguments specifying the permutation of rows and columns.
References mat::MatrixBase< Treal, Tmatrix >::getPermutedIndexes().
|
inline |
|
inline |
|
inline |
MatrixGeneral< Treal, Tmatrix > & mat::MatrixGeneral< Treal, Tmatrix >::operator+= | ( | const XYZ< Treal, MatrixGeneral< Treal, Tmatrix >, MatrixGeneral< Treal, Tmatrix > > & | smm | ) |
C += alpha * op(A) * op(B)
References mat::gemm().
MatrixGeneral< Treal, Tmatrix > & mat::MatrixGeneral< Treal, Tmatrix >::operator+= | ( | const XYZ< Treal, MatrixGeneral< Treal, Tmatrix >, MatrixSymmetric< Treal, Tmatrix > > & | smm | ) |
C += alpha * B * A : A is symmetric.
References mat::symm().
MatrixGeneral< Treal, Tmatrix > & mat::MatrixGeneral< Treal, Tmatrix >::operator+= | ( | const XYZ< Treal, MatrixSymmetric< Treal, Tmatrix >, MatrixGeneral< Treal, Tmatrix > > & | smm | ) |
C += alpha * A * B : A is symmetric.
References mat::symm().
MatrixGeneral< Treal, Tmatrix > & mat::MatrixGeneral< Treal, Tmatrix >::operator+= | ( | const XYZ< Treal, MatrixSymmetric< Treal, Tmatrix >, MatrixSymmetric< Treal, Tmatrix > > & | smm | ) |
C += alpha * A * B : A and B are symmetric.
|
inline |
B += A.
References A.
|
inline |
B += alpha * A.
|
inline |
References A.
|
inline |
References mat::MatrixBase< Treal, Tmatrix >::operator=().
|
inline |
References mat::MatrixBase< Treal, Tmatrix >::operator=(), and mat::symm().
|
inline |
References mat::MatrixBase< Treal, Tmatrix >::operator=().
|
inline |
References mat::MatrixBase< Treal, Tmatrix >::operator=(), and mat::transpose().
MatrixGeneral< Treal, Tmatrix > & mat::MatrixGeneral< Treal, Tmatrix >::operator= | ( | const XY< MatrixGeneral< Treal, Tmatrix >, MatrixGeneral< Treal, Tmatrix > > & | mm | ) |
C = op(A) * op(B)
References mat::gemm().
MatrixGeneral< Treal, Tmatrix > & mat::MatrixGeneral< Treal, Tmatrix >::operator= | ( | const XY< MatrixGeneral< Treal, Tmatrix >, MatrixSymmetric< Treal, Tmatrix > > & | mm | ) |
C = B * A : A is symmetric.
References mat::symm().
MatrixGeneral< Treal, Tmatrix > & mat::MatrixGeneral< Treal, Tmatrix >::operator= | ( | const XY< MatrixSymmetric< Treal, Tmatrix >, MatrixGeneral< Treal, Tmatrix > > & | mm | ) |
C = A * B : A is symmetric.
References mat::symm().
MatrixGeneral< Treal, Tmatrix > & mat::MatrixGeneral< Treal, Tmatrix >::operator= | ( | const XY< MatrixSymmetric< Treal, Tmatrix >, MatrixSymmetric< Treal, Tmatrix > > & | mm | ) |
C = A * B : A and B are symmetric.
MatrixGeneral< Treal, Tmatrix > & mat::MatrixGeneral< Treal, Tmatrix >::operator= | ( | const XYZ< Treal, MatrixGeneral< Treal, Tmatrix >, MatrixGeneral< Treal, Tmatrix > > & | smm | ) |
C = alpha * op(A) * op(B)
References mat::gemm().
MatrixGeneral< Treal, Tmatrix > & mat::MatrixGeneral< Treal, Tmatrix >::operator= | ( | const XYZ< Treal, MatrixGeneral< Treal, Tmatrix >, MatrixSymmetric< Treal, Tmatrix > > & | smm | ) |
C = alpha * B * A : A is symmetric.
References mat::symm().
MatrixGeneral< Treal, Tmatrix > & mat::MatrixGeneral< Treal, Tmatrix >::operator= | ( | const XYZ< Treal, MatrixGeneral< Treal, Tmatrix >, MatrixTriangular< Treal, Tmatrix > > & | smm | ) |
B = alpha * B * op(A) : A is upper triangular.
References mat::trmm().
MatrixGeneral< Treal, Tmatrix > & mat::MatrixGeneral< Treal, Tmatrix >::operator= | ( | const XYZ< Treal, MatrixSymmetric< Treal, Tmatrix >, MatrixGeneral< Treal, Tmatrix > > & | smm | ) |
C = alpha * A * B : A is symmetric.
References mat::symm().
MatrixGeneral< Treal, Tmatrix > & mat::MatrixGeneral< Treal, Tmatrix >::operator= | ( | const XYZ< Treal, MatrixSymmetric< Treal, Tmatrix >, MatrixSymmetric< Treal, Tmatrix > > & | smm | ) |
C = alpha * A * B : A and B are symmetric.
MatrixGeneral< Treal, Tmatrix > & mat::MatrixGeneral< Treal, Tmatrix >::operator= | ( | const XYZ< Treal, MatrixTriangular< Treal, Tmatrix >, MatrixGeneral< Treal, Tmatrix > > & | smm | ) |
B = alpha * op(A) * B : A is upper triangular.
References mat::trmm().
MatrixGeneral< Treal, Tmatrix > & mat::MatrixGeneral< Treal, Tmatrix >::operator= | ( | const XYZpUV< Treal, MatrixGeneral< Treal, Tmatrix >, MatrixGeneral< Treal, Tmatrix >, Treal, MatrixGeneral< Treal, Tmatrix > > & | smmpsm | ) |
C = alpha * op(A) * op(B) + beta * C
References mat::gemm().
MatrixGeneral< Treal, Tmatrix > & mat::MatrixGeneral< Treal, Tmatrix >::operator= | ( | const XYZpUV< Treal, MatrixGeneral< Treal, Tmatrix >, MatrixSymmetric< Treal, Tmatrix >, Treal, MatrixGeneral< Treal, Tmatrix > > & | smmpsm | ) |
C = alpha * B * A + beta * C : A is symmetric.
References mat::symm().
MatrixGeneral< Treal, Tmatrix > & mat::MatrixGeneral< Treal, Tmatrix >::operator= | ( | const XYZpUV< Treal, MatrixSymmetric< Treal, Tmatrix >, MatrixGeneral< Treal, Tmatrix >, Treal, MatrixGeneral< Treal, Tmatrix > > & | smmpsm | ) |
C = alpha * A * B + beta * C : A is symmetric.
References mat::symm().
MatrixGeneral< Treal, Tmatrix > & mat::MatrixGeneral< Treal, Tmatrix >::operator= | ( | const XYZpUV< Treal, MatrixSymmetric< Treal, Tmatrix >, MatrixSymmetric< Treal, Tmatrix >, Treal, MatrixGeneral< Treal, Tmatrix > > & | smmpsm | ) |
C = alpha * A * B + beta * C : A and B are symmetric.
|
inline |
|
inline |
C = A + B.
|
inline |
|
inline |
Referenced by mainFun().
|
inline |
References mat::matrix_matr.
|
inlineprotected |
References mat::matrix_matr.
|
inline |
Referenced by mainFun().
void mat::MatrixGeneral< Treal, Tmatrix >::thresh | ( | Treal const | threshold, |
normType const | norm | ||
) |
References mat::frobNorm.
|
inlinestatic |
Referenced by mat::trace().
|
inlinestatic |
Referenced by mat::trace().
|
inline |
References mat::matrix_matr.
|
inlineprotected |
References mat::matrix_matr.