#include <deal.II/lac/trilinos_sparse_matrix.h>

Inherits Epetra_Operator.

Public Types
using	VectorType = Epetra_MultiVector

using	Range = VectorType

using	Domain = VectorType

Public Member Functions
Constructors / destructor
	TrilinosPayload ()

	TrilinosPayload (const TrilinosWrappers::SparseMatrix &matrix_exemplar, const TrilinosWrappers::SparseMatrix &matrix)

	TrilinosPayload (const TrilinosWrappers::SparseMatrix &matrix_exemplar, const TrilinosWrappers::PreconditionBase &preconditioner)

	TrilinosPayload (const TrilinosWrappers::PreconditionBase &preconditioner_exemplar, const TrilinosWrappers::PreconditionBase &preconditioner)

	TrilinosPayload (const TrilinosPayload &payload)

	TrilinosPayload (const TrilinosPayload &first_op, const TrilinosPayload &second_op)

virtual	~TrilinosPayload () override=default

TrilinosPayload	identity_payload () const

TrilinosPayload	null_payload () const

TrilinosPayload	transpose_payload () const

template<typename Solver , typename Preconditioner >
std::enable_if< std::is_base_of< TrilinosWrappers::SolverBase, Solver >::value &&std::is_base_of< TrilinosWrappers::PreconditionBase, Preconditioner >::value, TrilinosPayload >::type	inverse_payload (Solver &, const Preconditioner &) const

template<typename Solver , typename Preconditioner >
std::enable_if< !(std::is_base_of< TrilinosWrappers::SolverBase, Solver >::value &&std::is_base_of< TrilinosWrappers::PreconditionBase, Preconditioner >::value), TrilinosPayload >::type	inverse_payload (Solver &, const Preconditioner &) const

Core Epetra_Operator functionality
virtual bool	UseTranspose () const override

virtual int	SetUseTranspose (bool UseTranspose) override

virtual int	Apply (const VectorType &X, VectorType &Y) const override

virtual int	ApplyInverse (const VectorType &Y, VectorType &X) const override

Additional Epetra_Operator functionality
virtual const char *	Label () const override

virtual const Epetra_Comm &	Comm () const override

virtual const Epetra_Map &	OperatorDomainMap () const override

virtual const Epetra_Map &	OperatorRangeMap () const override

Private Member Functions
virtual bool	HasNormInf () const override

virtual double	NormInf () const override

Private Attributes
bool	use_transpose

Epetra_MpiComm	communicator

Epetra_Map	domain_map

Epetra_Map	range_map

LinearOperator functionality
std::function< void(VectorType &, const VectorType &)>	vmult

std::function< void(VectorType &, const VectorType &)>	Tvmult

std::function< void(VectorType &, const VectorType &)>	inv_vmult

std::function< void(VectorType &, const VectorType &)>	inv_Tvmult

IndexSet	locally_owned_domain_indices () const

IndexSet	locally_owned_range_indices () const

MPI_Comm	get_mpi_communicator () const

void	transpose ()

Detailed Description

This is an extension class to LinearOperators for Trilinos sparse matrix and preconditioner types. It provides the interface to performing basic operations (vmult and Tvmult) on Trilinos vector types. It fulfills the requirements necessary for wrapping a Trilinos solver, which calls Epetra_Operator functions, as a LinearOperator.

Note: The TrilinosWrappers::SparseMatrix or TrilinosWrappers::PreconditionBase that this payload wraps is passed by reference to the vmult and Tvmult functions. This object is not thread-safe when the transpose flag is set on it or the Trilinos object to which it refers. See the docuemtation for the TrilinosWrappers::internal::LinearOperatorImplementation::TrilinosPayload::SetUseTranspose() function for further details.

Author: Jean-Paul Pelteret, 2016

Definition at line 2241 of file trilinos_sparse_matrix.h.

Member Typedef Documentation

using TrilinosWrappers::internal::LinearOperatorImplementation::TrilinosPayload::VectorType = Epetra_MultiVector

Definition for the internally supported vector type.

Definition at line 2247 of file trilinos_sparse_matrix.h.

using TrilinosWrappers::internal::LinearOperatorImplementation::TrilinosPayload::Range = VectorType

Definition for the vector type for the domain space of the operator.

Definition at line 2252 of file trilinos_sparse_matrix.h.

using TrilinosWrappers::internal::LinearOperatorImplementation::TrilinosPayload::Domain = VectorType

Definition for the vector type for the range space of the operator.

Definition at line 2257 of file trilinos_sparse_matrix.h.

Constructor & Destructor Documentation

TrilinosWrappers::internal::LinearOperatorImplementation::TrilinosPayload::TrilinosPayload ( )

Default constructor

Note: By design, the resulting object is inoperable since there is insufficient information with which to construct the domain and range maps.

Definition at line 2627 of file trilinos_sparse_matrix.cc.

TrilinosWrappers::internal::LinearOperatorImplementation::TrilinosPayload::TrilinosPayload	(	const TrilinosWrappers::SparseMatrix &	matrix_exemplar,
		const TrilinosWrappers::SparseMatrix &	matrix
	)

Constructor for a sparse matrix based on an exemplary matrix

Definition at line 2666 of file trilinos_sparse_matrix.cc.

TrilinosWrappers::internal::LinearOperatorImplementation::TrilinosPayload::TrilinosPayload	(	const TrilinosWrappers::SparseMatrix &	matrix_exemplar,
		const TrilinosWrappers::PreconditionBase &	preconditioner
	)

Constructor for a preconditioner based on an exemplary matrix

Definition at line 2751 of file trilinos_sparse_matrix.cc.

TrilinosWrappers::internal::LinearOperatorImplementation::TrilinosPayload::TrilinosPayload	(	const TrilinosWrappers::PreconditionBase &	preconditioner_exemplar,
		const TrilinosWrappers::PreconditionBase &	preconditioner
	)

Constructor for a preconditioner based on an exemplary preconditioner

Definition at line 2871 of file trilinos_sparse_matrix.cc.

TrilinosWrappers::internal::LinearOperatorImplementation::TrilinosPayload::TrilinosPayload ( const TrilinosPayload & payload )

Default copy constructor

Definition at line 2990 of file trilinos_sparse_matrix.cc.

TrilinosWrappers::internal::LinearOperatorImplementation::TrilinosPayload::TrilinosPayload	(	const TrilinosPayload &	first_op,
		const TrilinosPayload &	second_op
	)

Composite copy constructor

This is required for PackagedOperations as it sets up the domain and range maps, and composite vmult and Tvmult operations based on the combined operation of both operations

Definition at line 3005 of file trilinos_sparse_matrix.cc.

virtual TrilinosWrappers::internal::LinearOperatorImplementation::TrilinosPayload::~TrilinosPayload ( )

overridevirtualdefault

Destructor

Member Function Documentation

TrilinosPayload TrilinosWrappers::internal::LinearOperatorImplementation::TrilinosPayload::identity_payload ( ) const

Return a payload configured for identity operations

Definition at line 3018 of file trilinos_sparse_matrix.cc.

TrilinosPayload TrilinosWrappers::internal::LinearOperatorImplementation::TrilinosPayload::null_payload ( ) const

Return a payload configured for null operations

Definition at line 3044 of file trilinos_sparse_matrix.cc.

TrilinosPayload TrilinosWrappers::internal::LinearOperatorImplementation::TrilinosPayload::transpose_payload ( ) const

Return a payload configured for transpose operations

Definition at line 3084 of file trilinos_sparse_matrix.cc.

template<typename Solver , typename Preconditioner >

std::enable_if< std::is_base_of<TrilinosWrappers::SolverBase, Solver>::value && std::is_base_of<TrilinosWrappers::PreconditionBase, Preconditioner>::value, TrilinosPayload>::type TrilinosWrappers::internal::LinearOperatorImplementation::TrilinosPayload::inverse_payload	(	Solver &	,
		const Preconditioner &
	)		const

Return a payload configured for inverse operations

Invoking this factory function will configure two additional functions, namely inv_vmult and inv_Tvmult, both of which wrap inverse operations. The vmult and Tvmult operations retain the standard definitions inherited from op.

Note

This function is enabled only if the solver and preconditioner derive from the respective TrilinosWrappers base classes. The C++ compiler will therefore only consider this function if the following criterion are satisfied:

the Solver derives from TrilinosWrappers::SolverBase, and
the Preconditioner derives from TrilinosWrappers::PreconditionBase.

template<typename Solver , typename Preconditioner >

std::enable_if< !(std::is_base_of<TrilinosWrappers::SolverBase, Solver>::value && std::is_base_of<TrilinosWrappers::PreconditionBase, Preconditioner>::value), TrilinosPayload>::type TrilinosWrappers::internal::LinearOperatorImplementation::TrilinosPayload::inverse_payload	(	Solver &	,
		const Preconditioner &
	)		const

Return a payload configured for inverse operations

Invoking this factory function will configure two additional functions, namely inv_vmult and inv_Tvmult, both of which are disabled because the Solver or Preconditioner are not compatible with Epetra_MultiVector. The vmult and Tvmult operations retain the standard definitions inherited from op.

Note

The C++ compiler will only consider this function if the following criterion are satisfied:

the Solver does not derive from TrilinosWrappers::SolverBase, and
the Preconditioner does not derive from TrilinosWrappers::PreconditionBase.

IndexSet TrilinosWrappers::internal::LinearOperatorImplementation::TrilinosPayload::locally_owned_domain_indices ( ) const

Return an IndexSet that defines the partitioning of the domain space of this matrix, i.e., the partitioning of the vectors this matrix has to be multiplied with / operate on.

Definition at line 3094 of file trilinos_sparse_matrix.cc.

IndexSet TrilinosWrappers::internal::LinearOperatorImplementation::TrilinosPayload::locally_owned_range_indices ( ) const

Return an IndexSet that defines the partitioning of the range space of this matrix, i.e., the partitioning of the vectors that result from matrix-vector products.

Definition at line 3102 of file trilinos_sparse_matrix.cc.

MPI_Comm TrilinosWrappers::internal::LinearOperatorImplementation::TrilinosPayload::get_mpi_communicator ( ) const

Return the MPI communicator object in use with this Payload.

Definition at line 3110 of file trilinos_sparse_matrix.cc.

void TrilinosWrappers::internal::LinearOperatorImplementation::TrilinosPayload::transpose ( )

Sets an internal flag so that all operations performed by the matrix, i.e., multiplications, are done in transposed order.

Note: This does not reshape the matrix to transposed form directly, so care should be taken when using this flag.

Definition at line 3122 of file trilinos_sparse_matrix.cc.

bool TrilinosWrappers::internal::LinearOperatorImplementation::TrilinosPayload::UseTranspose ( ) const

overridevirtual

Return the status of the transpose flag for this operator

This overloads the same function from the Trilinos class Epetra_Operator.

Definition at line 3130 of file trilinos_sparse_matrix.cc.

int TrilinosWrappers::internal::LinearOperatorImplementation::TrilinosPayload::SetUseTranspose ( bool UseTranspose )

overridevirtual

Sets an internal flag so that all operations performed by the matrix, i.e., multiplications, are done in transposed order.

This overloads the same function from the Trilinos class Epetra_Operator.

Note: This does not reshape the matrix to transposed form directly, so care should be taken when using this flag. When the flag is set to true (either here or directly on the underlying Trilinos object itself), this object is no longer thread-safe. In essence, it is not possible ensure that the transposed state of the LinearOperator and the underlying Trilinos object remain synchronized throughout all operations that may occur on different threads simultaneously.

Definition at line 3138 of file trilinos_sparse_matrix.cc.

int TrilinosWrappers::internal::LinearOperatorImplementation::TrilinosPayload::Apply	(	const VectorType &	X,
		VectorType &	Y
	)		const

overridevirtual

Apply the vmult operation on a vector X (of internally defined type VectorType) and store the result in the vector Y.

This overloads the same function from the Trilinos class Epetra_Operator.

Note: The intended operation depends on the status of the internal transpose flag. If this flag is set to true, the result will be the equivalent of performing a Tvmult operation.

Definition at line 3153 of file trilinos_sparse_matrix.cc.

int TrilinosWrappers::internal::LinearOperatorImplementation::TrilinosPayload::ApplyInverse	(	const VectorType &	Y,
		VectorType &	X
	)		const

overridevirtual

Apply the vmult inverse operation on a vector X (of internally defined type VectorType) and store the result in the vector Y.

In practise, this function is only called from a Trilinos solver if the wrapped object is to act as a preconditioner.

This overloads the same function from the Trilinos class Epetra_Operator.

Note: This function will only be operable if the payload has been initialized with an InverseOperator, or is a wrapper to a preconditioner. If not, then using this function will lead to an error being thrown.; The intended operation depends on the status of the internal transpose flag. If this flag is set to true, the result will be the equivalent of performing a Tvmult operation.

Definition at line 3164 of file trilinos_sparse_matrix.cc.

const char * TrilinosWrappers::internal::LinearOperatorImplementation::TrilinosPayload::Label ( ) const

overridevirtual

Return a label to describe this class.

This overloads the same function from the Trilinos class Epetra_Operator.

Definition at line 3175 of file trilinos_sparse_matrix.cc.

const Epetra_Comm & TrilinosWrappers::internal::LinearOperatorImplementation::TrilinosPayload::Comm ( ) const

overridevirtual

Return a reference to the underlying MPI communicator for this object.

This overloads the same function from the Trilinos class Epetra_Operator.

Definition at line 3183 of file trilinos_sparse_matrix.cc.

const Epetra_Map & TrilinosWrappers::internal::LinearOperatorImplementation::TrilinosPayload::OperatorDomainMap ( ) const

overridevirtual

Return the partitioning of the domain space of this matrix, i.e., the partitioning of the vectors this matrix has to be multiplied with.

This overloads the same function from the Trilinos class Epetra_Operator.

Definition at line 3191 of file trilinos_sparse_matrix.cc.

const Epetra_Map & TrilinosWrappers::internal::LinearOperatorImplementation::TrilinosPayload::OperatorRangeMap ( ) const

overridevirtual

Return the partitioning of the range space of this matrix, i.e., the partitioning of the vectors that are result from matrix-vector products.

This overloads the same function from the Trilinos class Epetra_Operator.

Definition at line 3199 of file trilinos_sparse_matrix.cc.

bool TrilinosWrappers::internal::LinearOperatorImplementation::TrilinosPayload::HasNormInf ( ) const

overrideprivatevirtual

Return a flag that describes whether this operator can return the computation of the infinity norm. Since in general this is not the case, this always returns a negetive result.

This overloads the same function from the Trilinos class Epetra_Operator.

Definition at line 3207 of file trilinos_sparse_matrix.cc.

double TrilinosWrappers::internal::LinearOperatorImplementation::TrilinosPayload::NormInf ( ) const

overrideprivatevirtual

Return the infinity norm of this operator. Throws an error since, in general, we cannot compute this value.

This overloads the same function from the Trilinos class Epetra_Operator.

Definition at line 3215 of file trilinos_sparse_matrix.cc.

Member Data Documentation

std::function<void(VectorType &, const VectorType &)> TrilinosWrappers::internal::LinearOperatorImplementation::TrilinosPayload::vmult

The standard matrix-vector operation to be performed by the payload when Apply is called.

Note: This is not called by a LinearOperator, but rather by Trilinos functions that expect this to mimic the action of the LinearOperator.

Definition at line 2425 of file trilinos_sparse_matrix.h.

std::function<void(VectorType &, const VectorType &)> TrilinosWrappers::internal::LinearOperatorImplementation::TrilinosPayload::Tvmult

The standard transpose matrix-vector operation to be performed by the payload when Apply is called.

Note: This is not called by a LinearOperator, but rather by Trilinos functions that expect this to mimic the action of the LinearOperator.

Definition at line 2434 of file trilinos_sparse_matrix.h.

std::function<void(VectorType &, const VectorType &)> TrilinosWrappers::internal::LinearOperatorImplementation::TrilinosPayload::inv_vmult

The inverse matrix-vector operation to be performed by the payload when ApplyInverse is called.

Note: This is not called by a LinearOperator, but rather by Trilinos functions that expect this to mimic the action of the InverseOperator.

Definition at line 2444 of file trilinos_sparse_matrix.h.

std::function<void(VectorType &, const VectorType &)> TrilinosWrappers::internal::LinearOperatorImplementation::TrilinosPayload::inv_Tvmult

The inverse transpose matrix-vector operation to be performed by the payload when ApplyInverse is called.

Note: This is not called by a LinearOperator, but rather by Trilinos functions that expect this to mimic the action of the InverseOperator.

Definition at line 2454 of file trilinos_sparse_matrix.h.

bool TrilinosWrappers::internal::LinearOperatorImplementation::TrilinosPayload::use_transpose

private

A flag recording whether the operator is to perform standard matrix-vector multiplication, or the transpose operation.

Definition at line 2577 of file trilinos_sparse_matrix.h.

Epetra_MpiComm TrilinosWrappers::internal::LinearOperatorImplementation::TrilinosPayload::communicator

private

Internal communication pattern in case the matrix needs to be copied from deal.II format.

Definition at line 2584 of file trilinos_sparse_matrix.h.

Epetra_Map TrilinosWrappers::internal::LinearOperatorImplementation::TrilinosPayload::domain_map

private

Epetra_Map that sets the partitioning of the domain space of this operator.

Definition at line 2593 of file trilinos_sparse_matrix.h.

Epetra_Map TrilinosWrappers::internal::LinearOperatorImplementation::TrilinosPayload::range_map

private

Epetra_Map that sets the partitioning of the range space of this operator.

Definition at line 2599 of file trilinos_sparse_matrix.h.

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

deal.II/lac/trilinos_sparse_matrix.h
/mnt/data/darndt/Sources/dealii-clang-6.0.0/source/lac/trilinos_sparse_matrix.cc

Public Types

Public Member Functions

Private Member Functions

Private Attributes

LinearOperator functionality

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation