#include <Xpetra_MatrixMatrix.hpp>

Static Public Member Functions
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 > &params=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 > &params=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 > &params=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 = alphaA + betaB. 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 = alphaA + betaB. More...

Detailed Description

template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
class Xpetra::MatrixMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >

Definition at line 220 of file Xpetra_MatrixMatrix.hpp.

Member Function Documentation

template<class Scalar , class LocalOrdinal , class GlobalOrdinal , class Node >

static void Xpetra::MatrixMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::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 250 of file Xpetra_MatrixMatrix.hpp.

template<class Scalar , class LocalOrdinal , class GlobalOrdinal , class Node >

static RCP<Matrix> Xpetra::MatrixMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::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 325 of file Xpetra_MatrixMatrix.hpp.

template<class Scalar , class LocalOrdinal , class GlobalOrdinal , class Node >

static RCP<Matrix> Xpetra::MatrixMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::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 > &	params = `null`
	)

inlinestatic

Helper function to do matrix-matrix multiply.

Returns C = AB.

Parameters

A	left matrix
transposeA	if true, use the transpose of A
B	right matrix
transposeB	if true, use the transpose of B
callFillCompleteOnResult	if true, the resulting matrix should be fillComplete'd

Definition at line 381 of file Xpetra_MatrixMatrix.hpp.

template<class Scalar , class LocalOrdinal , class GlobalOrdinal , class Node >

static RCP<Epetra_CrsMatrix> Xpetra::MatrixMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::MLTwoMatrixMultiply	(	const Epetra_CrsMatrix &	epA,
		const Epetra_CrsMatrix &	epB,
		Teuchos::FancyOStream &	fos
	)

inlinestatic

Definition at line 389 of file Xpetra_MatrixMatrix.hpp.

template<class Scalar , class LocalOrdinal , class GlobalOrdinal , class Node >

static RCP<BlockedCrsMatrix> Xpetra::MatrixMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::TwoMatrixMultiplyBlock	(	const BlockedCrsMatrix &	A,
		bool	transposeA,
		const BlockedCrsMatrix &	B,
		bool	transposeB,
		Teuchos::FancyOStream &	fos,
		bool	doFillComplete = `true`,
		bool	doOptimizeStorage = `true`
	)

inlinestatic

Helper function to do matrix-matrix multiply "in-place".

Returns RCP to non-constant Xpetra::BlockedCrsMatrix.

Parameters

A	left matrix
transposeA	if true, use the transpose of A
B	right matrix
transposeB	if true, use the transpose of B
doOptimizeStorage	if true, the resulting matrix should be fillComplete'd

Definition at line 407 of file Xpetra_MatrixMatrix.hpp.

template<class Scalar , class LocalOrdinal , class GlobalOrdinal , class Node >

static void Xpetra::MatrixMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::TwoMatrixAdd	(	const Matrix &	A,
		bool	transposeA,
		SC	alpha,
		Matrix &	B,
		SC	beta
	)

inlinestatic

Helper function to calculate B = alpha*A + beta*B.

Parameters

A	left matrix operand
transposeA	indicate whether to use transpose of A
alpha	scalar multiplier for A
B	right matrix operand
beta	scalar multiplier for B

Returns: sum in B.

Note that B does not have to be fill-completed.

Definition at line 519 of file Xpetra_MatrixMatrix.hpp.

template<class Scalar , class LocalOrdinal , class GlobalOrdinal , class Node >

static void Xpetra::MatrixMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::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`
	)

inlinestatic

Helper function to calculate C = alpha*A + beta*B.

Parameters

A	left matrix operand
transposeA	indicate whether to use transpose of A
alpha	scalar multiplier for A, defaults to 1.0
B	right matrix operand
transposeB	indicate whether to use transpose of B
beta	scalar multiplier for B, defaults to 1.0
C	resulting sum

It is up to the caller to ensure that the resulting matrix sum is fillComplete'd.

Definition at line 550 of file Xpetra_MatrixMatrix.hpp.

The documentation for this class was generated from the following file:

Xpetra_MatrixMatrix.hpp

Static Public Member Functions

Detailed Description

template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node> class Xpetra::MatrixMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >

Member Function Documentation

template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
class Xpetra::MatrixMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >