|
static void | Multiply (const Matrix &A, bool transposeA, const Matrix &B, bool transposeB, Matrix &C, bool call_FillComplete_on_result=true, bool doOptimizeStorage=true, const std::string &label=std::string(), const RCP< ParameterList > ¶ms=null) |
|
static RCP< Matrix > | Multiply (const Matrix &A, bool transposeA, const Matrix &B, bool transposeB, RCP< Matrix > C_in, Teuchos::FancyOStream &fos, bool doFillComplete=true, bool doOptimizeStorage=true, const std::string &label=std::string(), const RCP< ParameterList > ¶ms=null) |
| Helper function to do matrix-matrix multiply. More...
|
|
static RCP< Matrix > | Multiply (const Matrix &A, bool transposeA, const Matrix &B, bool transposeB, Teuchos::FancyOStream &fos, bool callFillCompleteOnResult=true, bool doOptimizeStorage=true, const std::string &label=std::string(), const RCP< ParameterList > ¶ms=null) |
| Helper function to do matrix-matrix multiply. More...
|
|
static RCP< Epetra_CrsMatrix > | MLTwoMatrixMultiply (const Epetra_CrsMatrix &epA, const Epetra_CrsMatrix &epB, Teuchos::FancyOStream &fos) |
|
static RCP< BlockedCrsMatrix > | TwoMatrixMultiplyBlock (const BlockedCrsMatrix &A, bool transposeA, const BlockedCrsMatrix &B, bool transposeB, Teuchos::FancyOStream &fos, bool doFillComplete=true, bool doOptimizeStorage=true) |
| Helper function to do matrix-matrix multiply "in-place". More...
|
|
static void | TwoMatrixAdd (const Matrix &A, bool transposeA, SC alpha, Matrix &B, SC beta) |
| Helper function to calculate B = alpha*A + beta*B. More...
|
|
static void | TwoMatrixAdd (const Matrix &A, bool transposeA, const SC &alpha, const Matrix &B, bool transposeB, const SC &beta, RCP< Matrix > &C, Teuchos::FancyOStream &fos, bool AHasFixedNnzPerRow=false) |
| Helper function to calculate C = alpha*A + beta*B. More...
|
|
template<>
class Xpetra::MatrixMatrix< double, int, int, EpetraNode >
Definition at line 767 of file Xpetra_MatrixMatrix.hpp.
static void Xpetra::MatrixMatrix< double, int, int, EpetraNode >::Multiply |
( |
const Matrix & |
A, |
|
|
bool |
transposeA, |
|
|
const Matrix & |
B, |
|
|
bool |
transposeB, |
|
|
Matrix & |
C, |
|
|
bool |
call_FillComplete_on_result = true , |
|
|
bool |
doOptimizeStorage = true , |
|
|
const std::string & |
label = std::string() , |
|
|
const RCP< ParameterList > & |
params = null |
|
) |
| |
|
inlinestatic |
Given CrsMatrix objects A, B and C, form the product C = A*B. In a parallel setting, A and B need not have matching distributions, but C needs to have the same row-map as A (if transposeA is false). At this time C=AT*B and C=A*BT are known to not work. However, C=A*B and C=AT*BT are known to work, Kurtis Nusbaum 03/24/2011
- Parameters
-
A | Input, must already have had 'FillComplete()' called. |
transposeA | Input, whether to use transpose of matrix A. |
B | Input, must already have had 'FillComplete()' called. |
transposeB | Input, whether to use transpose of matrix B. |
C | Result. On entry to this method, it doesn't matter whether FillComplete() has already been called on C or not. If it has, then C's graph must already contain all nonzero locations that will be produced when forming the product A*B. On exit, C.FillComplete() will have been called, unless the last argument to this function is specified to be false. |
call_FillComplete_on_result | Optional argument, defaults to true. Power users may specify this argument to be false if they DON'T want this function to call C.FillComplete. (It is often useful to allow this function to call C.FillComplete, in cases where one or both of the input matrices are rectangular and it is not trivial to know which maps to use for the domain- and range-maps.) |
Definition at line 800 of file Xpetra_MatrixMatrix.hpp.
static RCP<Matrix> Xpetra::MatrixMatrix< double, int, int, EpetraNode >::Multiply |
( |
const Matrix & |
A, |
|
|
bool |
transposeA, |
|
|
const Matrix & |
B, |
|
|
bool |
transposeB, |
|
|
RCP< Matrix > |
C_in, |
|
|
Teuchos::FancyOStream & |
fos, |
|
|
bool |
doFillComplete = true , |
|
|
bool |
doOptimizeStorage = true , |
|
|
const std::string & |
label = std::string() , |
|
|
const RCP< ParameterList > & |
params = null |
|
) |
| |
|
inlinestatic |
Helper function to do matrix-matrix multiply.
Given CrsMatrix objects A, B and C, form the product C = A*B. In a parallel setting, A and B need not have matching distributions, but C needs to have the same row-map as A (if transposeA is false). At this time C=AT*B and C=A*BT are known to not work. However, C=A*B and C=AT*BT are known to work, Kurtis Nusbaum 03/24/2011
- Parameters
-
A | Input, must already have had 'FillComplete()' called. |
transposeA | Input, whether to use transpose of matrix A. |
B | Input, must already have had 'FillComplete()' called. |
transposeB | Input, whether to use transpose of matrix B. |
C | Result. If Teuchos::null, a new CrsMatrix is created with optimal number of nnz per row. |
call_FillComplete_on_result | Optional argument, defaults to true. Power users may specify this argument to be false if they DON'T want this function to call C.FillComplete. (It is often useful to allow this function to call C.FillComplete, in cases where one or both of the input matrices are rectangular and it is not trivial to know which maps to use for the domain- and range-maps.) |
Definition at line 899 of file Xpetra_MatrixMatrix.hpp.