Go to the documentation of this file.
7 #ifndef __IPEXPANSIONMATRIX_HPP__
8 #define __IPEXPANSIONMATRIX_HPP__
17 class ExpansionMatrixSpace;
50 const Index* ExpandedPosIndices()
const;
61 const Index* CompressedPosIndices()
const;
110 const std::string& name,
112 const std::string& prefix
115 PrintImplOffset(jnlst, level, category, name, indent, prefix, 1, 1);
123 const std::string& name,
125 const std::string& prefix,
130 friend class ParExpansionMatrix;
185 delete[] compressed_pos_;
186 delete[] expanded_pos_;
198 return MakeNewExpansionMatrix();
211 return expanded_pos_;
225 return compressed_pos_;
MatrixSpace base class, corresponding to the Matrix base class.
virtual void AddMSinvZImpl(Number alpha, const Vector &S, const Vector &Z, Vector &X) const
X = X + alpha*(Matrix S^{-1} Z).
virtual void PrintImpl(const Journalist &jnlst, EJournalLevel level, EJournalCategory category, const std::string &name, Index indent, const std::string &prefix) const
Print detailed information about the matrix.
This file contains a base class for all exceptions and a set of macros to help with exceptions.
double Number
Type of all numbers.
virtual void ComputeColAMaxImpl(Vector &cols_norms, bool init) const
Compute the max-norm of the columns in the matrix.
const Index * CompressedPosIndices() const
Return the vector of indices marking the compressed position.
EJournalLevel
Print Level Enum.
virtual void MultVectorImpl(Number alpha, const Vector &x, Number beta, Vector &y) const
Matrix-vector multiply.
void PrintImplOffset(const Journalist &jnlst, EJournalLevel level, EJournalCategory category, const std::string &name, Index indent, const std::string &prefix, Index row_offset, Index col_offset) const
~ExpansionMatrix()
Destructor.
int Index
Type of all indices of vectors, matrices etc.
ExpansionMatrixSpace(Index NLargeVec, Index NSmallVec, const Index *ExpPos, const int offset=0)
Constructor, given the list of elements of the large vector (of size NLargeVec) to be filtered into t...
This is the matrix space for ExpansionMatrix.
virtual Matrix * MakeNew() const
Pure virtual method for creating a new Matrix of the corresponding type.
virtual void SinvBlrmZMTdBrImpl(Number alpha, const Vector &S, const Vector &R, const Vector &Z, const Vector &D, Vector &X) const
X = S^{-1} (r + alpha*Z*M^Td).
EJournalCategory
Category Selection Enum.
virtual void TransMultVectorImpl(Number alpha, const Vector &x, Number beta, Vector &y) const
Matrix(transpose) vector multiply.
void operator=(const ExpansionMatrix &)
Default Assignment Operator.
virtual void ComputeRowAMaxImpl(Vector &rows_norms, bool init) const
Compute the max-norm of the rows in the matrix.
const Index * CompressedPosIndices() const
Accessor Method to obtain the Index array (of length NLargeVec=NRows()) that stores the mapping from ...
~ExpansionMatrixSpace()
Destructor.
Class responsible for all message output.
ExpansionMatrix(const ExpansionMatrix &)
Copy Constructor.
const Index * ExpandedPosIndices() const
Accessor Method to obtain the Index array (of length NSmallVec=NCols()) that stores the mapping from ...
const ExpansionMatrixSpace * owner_space_
ExpansionMatrix(const ExpansionMatrixSpace *owner_space)
Constructor, taking the owner_space.
ExpansionMatrix()
Default Constructor.
Class for expansion/projection matrices.
const Index * ExpandedPosIndices() const
Return the vector of indices marking the expanded position.
ExpansionMatrix * MakeNewExpansionMatrix() const
Method for creating a new matrix of this specific type.