VbrMatrix: Variable block row matrix. More...

#include <Tpetra_VbrMatrix_decl.hpp>

Inheritance diagram for Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >:

Public Types
typedef Scalar	scalar_type
	The type of the entries of the input and output multivectors.
typedef LocalOrdinal	local_ordinal_type
	The local index type.
typedef GlobalOrdinal	global_ordinal_type
	The global index type.
typedef Node	node_type
	The Kokkos Node type.
Typedefs
typedef char	packet_type
	The type of each datum being sent or received in an Import or Export.
Public Member Functions
Constructor/Destructor Methods
	VbrMatrix (const Teuchos::RCP< const BlockMap< LocalOrdinal, GlobalOrdinal, Node > > &blkRowMap, size_t maxNumEntriesPerRow, ProfileType pftype=DynamicProfile)
	Constructor specifying the row-map and the max number of (block) non-zeros for all rows.
	VbrMatrix (const Teuchos::RCP< const BlockCrsGraph< LocalOrdinal, GlobalOrdinal, Node > > &blkGraph)
	Constructor specifying a pre-filled block-graph.
virtual	~VbrMatrix ()
	Destructor.
Advanced Mathematical operations
template<class DomainScalar , class RangeScalar >
void	multiply (const MultiVector< DomainScalar, LocalOrdinal, GlobalOrdinal, Node > &X, MultiVector< RangeScalar, LocalOrdinal, GlobalOrdinal, Node > &Y, Teuchos::ETransp trans, RangeScalar alpha, RangeScalar beta) const
	Multiply this matrix by a MultiVector.
template<class DomainScalar , class RangeScalar >
void	solve (const MultiVector< RangeScalar, LocalOrdinal, GlobalOrdinal, Node > &Y, MultiVector< DomainScalar, LocalOrdinal, GlobalOrdinal, Node > &X, Teuchos::ETransp trans) const
	Triangular Solve -- Matrix must be triangular.
Operator Methods
Teuchos::RCP< const Map < LocalOrdinal, GlobalOrdinal, Node > >	getDomainMap () const
	Returns the (point-entry) Map associated with the domain of this operator.
Teuchos::RCP< const Map < LocalOrdinal, GlobalOrdinal, Node > >	getRangeMap () const
	Returns the (point-entry) Map associated with the range of this operator.
void	apply (const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &X, MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &Y, Teuchos::ETransp trans=Teuchos::NO_TRANS, Scalar alpha=Teuchos::ScalarTraits< Scalar >::one(), Scalar beta=Teuchos::ScalarTraits< Scalar >::zero()) const
	Computes the operator-multivector application.
void	applyInverse (const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &Y, MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &X, Teuchos::ETransp trans) const
	Triangular Solve -- Matrix must be triangular.
bool	hasTransposeApply () const
	Indicates whether this operator supports applying the adjoint operator.
Attribute Query Methods
Teuchos::RCP< const BlockMap < LocalOrdinal, GlobalOrdinal, Node > >	getBlockRowMap () const
	Returns the block-row map.
Teuchos::RCP< const BlockMap < LocalOrdinal, GlobalOrdinal, Node > >	getBlockColMap () const
	Returns the block-column map.
Teuchos::RCP< const BlockMap < LocalOrdinal, GlobalOrdinal, Node > >	getBlockDomainMap () const
	Returns the block-domain map.
Teuchos::RCP< const BlockMap < LocalOrdinal, GlobalOrdinal, Node > >	getBlockRangeMap () const
	Returns the block-range map.
Teuchos::RCP< const Map < LocalOrdinal, GlobalOrdinal, Node > >	getPointRowMap () const
	Returns the point-row map.
Teuchos::RCP< const Map < LocalOrdinal, GlobalOrdinal, Node > >	getPointColMap () const
	Returns the point-column map.
bool	isFillComplete () const
	Return true if fillComplete has been called, false otherwise.
Insertion Methods
void	putScalar (Scalar s)
	Set the specified scalar throughout the matrix.
void	setGlobalBlockEntry (GlobalOrdinal globalBlockRow, GlobalOrdinal globalBlockCol, const Teuchos::SerialDenseMatrix< int, Scalar > &blockEntry)
	Copy the contents of the input block-entry into the matrix.
void	setLocalBlockEntry (LocalOrdinal localBlockRow, LocalOrdinal localBlockCol, const Teuchos::SerialDenseMatrix< int, Scalar > &blockEntry)
	Copy the contents of the input block-entry into the matrix.
void	sumIntoGlobalBlockEntry (GlobalOrdinal globalBlockRow, GlobalOrdinal globalBlockCol, const Teuchos::SerialDenseMatrix< int, Scalar > &blockEntry)
	Add the contents of the input block-entry into the matrix.
void	sumIntoLocalBlockEntry (LocalOrdinal localBlockRow, LocalOrdinal localBlockCol, const Teuchos::SerialDenseMatrix< int, Scalar > &blockEntry)
	Add the contents of the input block-entry into the matrix.
void	setGlobalBlockEntry (GlobalOrdinal globalBlockRow, GlobalOrdinal globalBlockCol, LocalOrdinal blkRowSize, LocalOrdinal blkColSize, LocalOrdinal LDA, const Teuchos::ArrayView< const Scalar > &blockEntry)
	Copy the contents of the input block-entry into the matrix.
void	setLocalBlockEntry (LocalOrdinal localBlockRow, LocalOrdinal localBlockCol, LocalOrdinal blkRowSize, LocalOrdinal blkColSize, LocalOrdinal LDA, const Teuchos::ArrayView< const Scalar > &blockEntry)
	Copy the contents of the input block-entry into the matrix.
void	sumIntoGlobalBlockEntry (GlobalOrdinal globalBlockRow, GlobalOrdinal globalBlockCol, LocalOrdinal blkRowSize, LocalOrdinal blkColSize, LocalOrdinal LDA, const Teuchos::ArrayView< const Scalar > &blockEntry)
	Add the contents of the input block-entry into the matrix.
void	sumIntoLocalBlockEntry (LocalOrdinal localBlockRow, LocalOrdinal localBlockCol, LocalOrdinal blkRowSize, LocalOrdinal blkColSize, LocalOrdinal LDA, const Teuchos::ArrayView< const Scalar > &blockEntry)
	Add the contents of the input block-entry into the matrix.
Transformational Methods
void	fillComplete (const Teuchos::RCP< const BlockMap< LocalOrdinal, GlobalOrdinal, Node > > &blockDomainMap, const Teuchos::RCP< const BlockMap< LocalOrdinal, GlobalOrdinal, Node > > &blockRangeMap)
	Transition the matrix to the packed, optimized-storage state.
void	fillComplete ()
	Transition the matrix to the packed, optimized-storage state.
void	globalAssemble ()
	Communicate non-local contributions to the processors that own those contributions.
Extraction Methods
void	getGlobalBlockRowView (GlobalOrdinal globalBlockRow, LocalOrdinal &numPtRows, Teuchos::ArrayView< const GlobalOrdinal > &blockCols, Teuchos::Array< LocalOrdinal > &ptColsPerBlockCol, Teuchos::Array< Teuchos::ArrayRCP< const Scalar > > &blockEntries) const
	Returns a const read-only view of the block-entries in the specified row.
void	getLocalBlockRowView (LocalOrdinal localBlockRow, LocalOrdinal &numPtRows, Teuchos::ArrayView< const LocalOrdinal > &blockCols, Teuchos::Array< LocalOrdinal > &ptColsPerBlockCol, Teuchos::ArrayRCP< const Scalar > &blockEntries) const
	Returns a const read-only view of the block-entries in the specified row.
void	getGlobalBlockEntryView (GlobalOrdinal globalBlockRow, GlobalOrdinal globalBlockCol, LocalOrdinal &numPtRows, LocalOrdinal &numPtCols, Teuchos::ArrayRCP< const Scalar > &blockEntry) const
	Returns a const read-only view of a block-entry.
void	getGlobalBlockEntryViewNonConst (GlobalOrdinal globalBlockRow, GlobalOrdinal globalBlockCol, LocalOrdinal &numPtRows, LocalOrdinal &numPtCols, Teuchos::ArrayRCP< Scalar > &blockEntry)
	Returns a non-const read-write view of a block-entry.
void	getLocalBlockEntryView (LocalOrdinal localBlockRow, LocalOrdinal localBlockCol, LocalOrdinal &numPtRows, LocalOrdinal &numPtCols, Teuchos::ArrayRCP< const Scalar > &blockEntry) const
	Returns a const read-only view of a block-entry.
void	getLocalBlockEntryViewNonConst (LocalOrdinal localBlockRow, LocalOrdinal localBlockCol, LocalOrdinal &numPtRows, LocalOrdinal &numPtCols, Teuchos::ArrayRCP< Scalar > &blockEntry)
	Returns a non-const read-write view of a block-entry.
void	getLocalDiagCopy (Tpetra::Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &diag) const
	Return a copy of the (point-entry) diagonal values.
Teuchos::RCP< const BlockCrsGraph< LocalOrdinal, GlobalOrdinal, Node > >	getBlockCrsGraph ()
Implementation of DistObject
virtual bool	checkSizes (const SrcDistObject &source)
	Compare the source and target (this) objects for compatibility.
virtual void	copyAndPermute (const SrcDistObject &source, size_t numSameIDs, const Teuchos::ArrayView< const LocalOrdinal > &permuteToLIDs, const Teuchos::ArrayView< const LocalOrdinal > &permuteFromLIDs)
	Perform copies and permutations that are local to this process.
virtual void	packAndPrepare (const SrcDistObject &source, const Teuchos::ArrayView< const LocalOrdinal > &exportLIDs, Teuchos::Array< char > &exports, const Teuchos::ArrayView< size_t > &numPacketsPerLID, size_t &constantNumPackets, Distributor &distor)
	Perform any packing or preparation required for communication.
virtual void	unpackAndCombine (const Teuchos::ArrayView< const LocalOrdinal > &importLIDs, const Teuchos::ArrayView< const char > &imports, const Teuchos::ArrayView< size_t > &numPacketsPerLID, size_t constantNumPackets, Distributor &distor, CombineMode CM)
	Perform any unpacking and combining after communication.
Implementation of Teuchos::Describable
std::string	description () const
	One-line descriptiion of this object.
void	describe (Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel=Teuchos::Describable::verbLevel_default) const
	Print the object with some verbosity level to a FancyOStream object.
Public methods for redistributing data
void	doImport (const SrcDistObject &source, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
	Import data into this object using an Import object ("forward mode").
void	doImport (const SrcDistObject &source, const Export< LocalOrdinal, GlobalOrdinal, Node > &exporter, CombineMode CM)
	Import data into this object using an Export object ("reverse mode").
void	doExport (const SrcDistObject &source, const Export< LocalOrdinal, GlobalOrdinal, Node > &exporter, CombineMode CM)
	Export data into this object using an Export object ("forward mode").
void	doExport (const SrcDistObject &source, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
	Export data into this object using an Import object ("reverse mode").
Attribute accessor methods
bool	isDistributed () const
	Whether this is a globally distributed object.
virtual Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > >	getMap () const
	The Map describing the parallel distribution of this object.
I/O methods
void	print (std::ostream &os) const
	Print this object to the given output stream.
Protected Member Functions
virtual size_t	constantNumberOfPackets () const
	Whether the implementation's instance promises always to have a constant number of packets per LID, and if so, how many packets per LID there are.
virtual void	doTransfer (const SrcDistObject &src, CombineMode CM, size_t numSameIDs, const Teuchos::ArrayView< const LocalOrdinal > &permuteToLIDs, const Teuchos::ArrayView< const LocalOrdinal > &permuteFromLIDs, const Teuchos::ArrayView< const LocalOrdinal > &remoteLIDs, const Teuchos::ArrayView< const LocalOrdinal > &exportLIDs, Distributor &distor, ReverseOption revOp)
	Redistribute data across memory images.
virtual void	createViews () const
	Hook for creating a const view.
virtual void	createViewsNonConst (KokkosClassic::ReadWriteOption rwo)
	Hook for creating a nonconst view.
virtual void	releaseViews () const
	Hook for releasing views.
Protected Attributes
Teuchos::RCP< const Map < LocalOrdinal, GlobalOrdinal, Node > >	map_
	The Map over which this object is distributed.
Teuchos::Array< char >	imports_
	Buffer into which packed data are imported (received from other processes).
Teuchos::Array< size_t >	numImportPacketsPerLID_
	Number of packets to receive for each receive operation.
Teuchos::Array< char >	exports_
	Buffer from which packed data are exported (sent to other processes).
Teuchos::Array< size_t >	numExportPacketsPerLID_
	Number of packets to send for each send operation.
Methods for use only by experts
virtual void	removeEmptyProcessesInPlace (const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &newMap)
	Remove processes which contain no elements in this object's Map.
void	removeEmptyProcessesInPlace (Teuchos::RCP< Tpetra::DistObject< PT, LO, GO, NT > > &input, const Teuchos::RCP< const Map< LO, GO, NT > > &newMap)
void	removeEmptyProcessesInPlace (Teuchos::RCP< Tpetra::DistObject< PT, LO, GO, NT > > &input)

Detailed Description

template<class Scalar = Operator<>::scalar_type, class LocalOrdinal = typename Operator<Scalar>::local_ordinal_type, class GlobalOrdinal = typename Operator<Scalar, LocalOrdinal>::global_ordinal_type, class Node = typename Operator<Scalar, LocalOrdinal, GlobalOrdinal>::node_type>
class Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >

VbrMatrix: Variable block row matrix.

The VbrMatrix class has two significant 'states', distinguished by whether or not storage has been optimized (packed) or not.

When the matrix is in the non-optimized-storage state, internal data storage is in a non-contiguous data-structure that allows for convenient insertion of data.

When the matrix is in the optimized-storage state, internal data is stored in contiguous (packed) arrays. When in this state, existing entries may be updated and replaced, but no new entries (indices and/or coefficients) may be inserted. In other words, the sparsity pattern or structure of the matrix may not be changed.

Use of the matrix as an Operator (performing matrix-vector multiplication) is only allowed when it is in the optimized-storage state.

VbrMatrix has two constructors, one which leaves the matrix in the optimized- storage state, and another which leaves the matrix in the non-optimized-storage state.

When the VbrMatrix is constructed in the non-optimized-storage state, (and then filled using methods such as setGlobalBlockEntry etc.), it can then be transformed to the optimized-storage state by calling the method fillComplete().

Once in the optimized-storage state, the VbrMatrix can not be returned to the non-optimized-storage state.

Definition at line 117 of file Tpetra_VbrMatrix_decl.hpp.

Member Typedef Documentation

template<class Scalar = Operator<>::scalar_type, class LocalOrdinal = typename Operator<Scalar>::local_ordinal_type, class GlobalOrdinal = typename Operator<Scalar, LocalOrdinal>::global_ordinal_type, class Node = typename Operator<Scalar, LocalOrdinal, GlobalOrdinal>::node_type>

typedef Scalar Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::scalar_type

The type of the entries of the input and output multivectors.

Reimplemented from Tpetra::Operator< Scalar, LocalOrdinal, GlobalOrdinal, Node >.

Definition at line 121 of file Tpetra_VbrMatrix_decl.hpp.

template<class Scalar = Operator<>::scalar_type, class LocalOrdinal = typename Operator<Scalar>::local_ordinal_type, class GlobalOrdinal = typename Operator<Scalar, LocalOrdinal>::global_ordinal_type, class Node = typename Operator<Scalar, LocalOrdinal, GlobalOrdinal>::node_type>

typedef LocalOrdinal Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::local_ordinal_type

The local index type.

Reimplemented from Tpetra::Operator< Scalar, LocalOrdinal, GlobalOrdinal, Node >.

Definition at line 122 of file Tpetra_VbrMatrix_decl.hpp.

template<class Scalar = Operator<>::scalar_type, class LocalOrdinal = typename Operator<Scalar>::local_ordinal_type, class GlobalOrdinal = typename Operator<Scalar, LocalOrdinal>::global_ordinal_type, class Node = typename Operator<Scalar, LocalOrdinal, GlobalOrdinal>::node_type>

typedef GlobalOrdinal Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::global_ordinal_type

The global index type.

Reimplemented from Tpetra::Operator< Scalar, LocalOrdinal, GlobalOrdinal, Node >.

Definition at line 123 of file Tpetra_VbrMatrix_decl.hpp.

template<class Scalar = Operator<>::scalar_type, class LocalOrdinal = typename Operator<Scalar>::local_ordinal_type, class GlobalOrdinal = typename Operator<Scalar, LocalOrdinal>::global_ordinal_type, class Node = typename Operator<Scalar, LocalOrdinal, GlobalOrdinal>::node_type>

typedef Node Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::node_type

The Kokkos Node type.

Reimplemented from Tpetra::Operator< Scalar, LocalOrdinal, GlobalOrdinal, Node >.

Definition at line 124 of file Tpetra_VbrMatrix_decl.hpp.

typedef char Tpetra::DistObject< char , LocalOrdinal, GlobalOrdinal, Node >::packet_type [inherited]

The type of each datum being sent or received in an Import or Export.

Note that this type does not always correspond to the Scalar template parameter of subclasses.

Definition at line 183 of file Tpetra_DistObject_decl.hpp.

Constructor & Destructor Documentation

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

Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::VbrMatrix	(	const Teuchos::RCP< const BlockMap< LocalOrdinal, GlobalOrdinal, Node > > &	blkRowMap,
		size_t	maxNumEntriesPerRow,
		ProfileType	pftype = `DynamicProfile`
	)

Constructor specifying the row-map and the max number of (block) non-zeros for all rows.

After this constructor completes, the VbrMatrix is in the non-packed, non-optimized-storage, isFillComplete()==false state. Block-entries (rectangular, dense submatrices) may be inserted using class methods such as setGlobalBlockEntry(...), declared below.

Definition at line 75 of file Tpetra_VbrMatrix_def.hpp.

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

Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::VbrMatrix ( const Teuchos::RCP< const BlockCrsGraph< LocalOrdinal, GlobalOrdinal, Node > > & blkGraph )

Constructor specifying a pre-filled block-graph.

Constructing a VbrMatrix with a pre-filled graph means that the matrix will start out in the optimized-storage state, i.e., isFillComplete()==true. The graph provided to this constructor must be already filled. (If blkGraph->isFillComplete() == false, an exception is thrown.)

Entries in the input BlockCrsGraph correspond to block-entries in the VbrMatrix. In other words, the VbrMatrix will have a block-row corresponding to each row in the graph, and a block-entry corresponding to each column- index in the graph.

Definition at line 101 of file Tpetra_VbrMatrix_def.hpp.

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

Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::~VbrMatrix ( ) [virtual]

Destructor.

Definition at line 139 of file Tpetra_VbrMatrix_def.hpp.

Member Function Documentation

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

template<class DomainScalar , class RangeScalar >

void Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::multiply	(	const MultiVector< DomainScalar, LocalOrdinal, GlobalOrdinal, Node > &	X,
		MultiVector< RangeScalar, LocalOrdinal, GlobalOrdinal, Node > &	Y,
		Teuchos::ETransp	trans,
		RangeScalar	alpha,
		RangeScalar	beta
	)		const

Multiply this matrix by a MultiVector.

X is required to be post-imported, i.e., described by the column map of the matrix. Y is required to be pre-exported, i.e., described by the row map of the matrix. See also the Operator::apply method which is implemented below.

Definition at line 172 of file Tpetra_VbrMatrix_def.hpp.

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

template<class DomainScalar , class RangeScalar >

void Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::solve	(	const MultiVector< RangeScalar, LocalOrdinal, GlobalOrdinal, Node > &	Y,
		MultiVector< DomainScalar, LocalOrdinal, GlobalOrdinal, Node > &	X,
		Teuchos::ETransp	trans
	)		const

Triangular Solve -- Matrix must be triangular.

Find X such that A*X = Y. X is required to be post-imported, i.e., described by the column map of the matrix. Y is required to be pre-exported, i.e., described by the row map of the matrix.

Both X and Y are required to have constant stride.

Note that if the diagonal block-entries are stored, they must be triangular. I.e., the matrix structure must be block-triangular, and any diagonal blocks must be "point"-triangular, meaning that coefficients on the "wrong" side of the point-diagonal must be zero.

Definition at line 194 of file Tpetra_VbrMatrix_def.hpp.

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

Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::getDomainMap ( ) const [virtual]

Returns the (point-entry) Map associated with the domain of this operator.

Note that this is a point-entry map, not a block-map.

Implements Tpetra::Operator< Scalar, LocalOrdinal, GlobalOrdinal, Node >.

Definition at line 146 of file Tpetra_VbrMatrix_def.hpp.

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

Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::getRangeMap ( ) const [virtual]

Returns the (point-entry) Map associated with the range of this operator.

Note that this is a point-entry map, not a block-map.

Implements Tpetra::Operator< Scalar, LocalOrdinal, GlobalOrdinal, Node >.

Definition at line 154 of file Tpetra_VbrMatrix_def.hpp.

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

void Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::apply	(	const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &	X,
		MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &	Y,
		Teuchos::ETransp	trans = `Teuchos::NO_TRANS`,
		Scalar	alpha = `Teuchos::ScalarTraits<Scalar>::one()`,
		Scalar	beta = `Teuchos::ScalarTraits<Scalar>::zero()`
	)		const `[virtual]`

Computes the operator-multivector application.

Loosely, performs $Y = \alpha \cdot A^{\textrm{trans}} \cdot X + \beta \cdot Y$ . However, the details of operation vary according to the values of alpha and beta. Specifically

if beta == 0, apply() must overwrite Y, so that any values in Y (including NaNs) are ignored.
if alpha == 0, apply() may short-circuit the operator, so that any values in X (including NaNs) are ignored.

Implements Tpetra::Operator< Scalar, LocalOrdinal, GlobalOrdinal, Node >.

Definition at line 253 of file Tpetra_VbrMatrix_def.hpp.

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

void Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::applyInverse	(	const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &	Y,
		MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &	X,
		Teuchos::ETransp	trans
	)		const

Triangular Solve -- Matrix must be triangular.

Find X such that A*X = Y. Both X and Y are required to have constant stride.

Definition at line 278 of file Tpetra_VbrMatrix_def.hpp.

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

bool Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::hasTransposeApply ( ) const [virtual]

Indicates whether this operator supports applying the adjoint operator.

VbrMatrix does support transpose-apply. (This method returns true.)

Reimplemented from Tpetra::Operator< Scalar, LocalOrdinal, GlobalOrdinal, Node >.

Definition at line 301 of file Tpetra_VbrMatrix_def.hpp.

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

Teuchos::RCP< const BlockMap< LocalOrdinal, GlobalOrdinal, Node > > Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::getBlockRowMap ( ) const

Returns the block-row map.

Definition at line 309 of file Tpetra_VbrMatrix_def.hpp.

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

Teuchos::RCP< const BlockMap< LocalOrdinal, GlobalOrdinal, Node > > Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::getBlockColMap ( ) const

Returns the block-column map.

Definition at line 325 of file Tpetra_VbrMatrix_def.hpp.

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

Teuchos::RCP< const BlockMap< LocalOrdinal, GlobalOrdinal, Node > > Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::getBlockDomainMap ( ) const

Returns the block-domain map.

Definition at line 333 of file Tpetra_VbrMatrix_def.hpp.

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

Teuchos::RCP< const BlockMap< LocalOrdinal, GlobalOrdinal, Node > > Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::getBlockRangeMap ( ) const

Returns the block-range map.

Definition at line 341 of file Tpetra_VbrMatrix_def.hpp.

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

Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::getPointRowMap ( ) const

Returns the point-row map.

Definition at line 317 of file Tpetra_VbrMatrix_def.hpp.

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

Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::getPointColMap ( ) const

Returns the point-column map.

Definition at line 349 of file Tpetra_VbrMatrix_def.hpp.

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

bool Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::isFillComplete ( ) const

Return true if fillComplete has been called, false otherwise.

Definition at line 162 of file Tpetra_VbrMatrix_def.hpp.

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

void Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::putScalar ( Scalar s )

Set the specified scalar throughout the matrix.

This method may be called any time (before or after fillComplete()).

Definition at line 981 of file Tpetra_VbrMatrix_def.hpp.

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

void Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::setGlobalBlockEntry	(	GlobalOrdinal	globalBlockRow,
		GlobalOrdinal	globalBlockCol,
		const Teuchos::SerialDenseMatrix< int, Scalar > &	blockEntry
	)

Copy the contents of the input block-entry into the matrix.

This method will create the specified block-entry if it doesn't already exist, but only if fillComplete() has not yet been called.

If the specified block-entry already exists in the matrix, it will be over-written (replaced) by the input block-entry.

Note that if globalBlockRow is not owned by the local processor (as indicated by getBlockRowMap()) then the block-entry is held in temporary storage until globalAssemble() is called (which is called internally by fillComplete()) and then globalAssemble() performans the communication needed to move the data to the owning processor.

This method may be called any time (before or after fillComplete()).

Definition at line 1032 of file Tpetra_VbrMatrix_def.hpp.

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

void Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::setLocalBlockEntry	(	LocalOrdinal	localBlockRow,
		LocalOrdinal	localBlockCol,
		const Teuchos::SerialDenseMatrix< int, Scalar > &	blockEntry
	)

Copy the contents of the input block-entry into the matrix.

This method will throw an exception if fillComplete() has not yet been called, or if the specified block-entry doesn't already exist in the matrix.

The coefficients of the specified block-entry will be over-written (replaced) by the input block-entry.

Definition at line 1049 of file Tpetra_VbrMatrix_def.hpp.

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

void Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::sumIntoGlobalBlockEntry	(	GlobalOrdinal	globalBlockRow,
		GlobalOrdinal	globalBlockCol,
		const Teuchos::SerialDenseMatrix< int, Scalar > &	blockEntry
	)

Add the contents of the input block-entry into the matrix.

This method will create the specified block-entry if it doesn't already exist, but only if fillComplete() has not yet been called.

If the specified block-entry already exists in the matrix, the contents of the input block-entry will be added to the values that are already present.

Note that if globalBlockRow is not owned by the local processor (as indicated by getBlockRowMap()) then the block-entry is held in temporary storage until globalAssemble() is called (which is called internally by fillComplete()) and then globalAssemble() performans the communication needed to move the data to the owning processor.

This method may be called any time (before or after fillComplete()).

Definition at line 1066 of file Tpetra_VbrMatrix_def.hpp.

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

void Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::sumIntoLocalBlockEntry	(	LocalOrdinal	localBlockRow,
		LocalOrdinal	localBlockCol,
		const Teuchos::SerialDenseMatrix< int, Scalar > &	blockEntry
	)

Add the contents of the input block-entry into the matrix.

This method will throw an exception if fillComplete() has not yet been called, or if the specified block-entry doesn't already exist in the matrix.

The contents of the input block-entry will be added to the values that are already present in the matrix.

Definition at line 1083 of file Tpetra_VbrMatrix_def.hpp.

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

void Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::setGlobalBlockEntry	(	GlobalOrdinal	globalBlockRow,
		GlobalOrdinal	globalBlockCol,
		LocalOrdinal	blkRowSize,
		LocalOrdinal	blkColSize,
		LocalOrdinal	LDA,
		const Teuchos::ArrayView< const Scalar > &	blockEntry
	)

Copy the contents of the input block-entry into the matrix.

This method will create the specified block-entry if it doesn't already exist, but only if fillComplete() has not yet been called.

If the specified block-entry already exists in the matrix, it will be over-written (replaced) by the input block-entry.

Note that if globalBlockRow is not owned by the local processor (as indicated by getBlockRowMap()) then the block-entry is held in temporary storage until globalAssemble() is called (which is called internally by fillComplete()) and then globalAssemble() performans the communication needed to move the data to the owning processor.

This method may be called any time (before or after fillComplete()).

Definition at line 1100 of file Tpetra_VbrMatrix_def.hpp.

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

void Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::setLocalBlockEntry	(	LocalOrdinal	localBlockRow,
		LocalOrdinal	localBlockCol,
		LocalOrdinal	blkRowSize,
		LocalOrdinal	blkColSize,
		LocalOrdinal	LDA,
		const Teuchos::ArrayView< const Scalar > &	blockEntry
	)

Copy the contents of the input block-entry into the matrix.

This method will throw an exception if fillComplete() has not yet been called, or if the specified block-entry doesn't already exist in the matrix.

The coefficients for the specified block-entry will be over-written (replaced) by the input block-entry.

Definition at line 1117 of file Tpetra_VbrMatrix_def.hpp.

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

void Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::sumIntoGlobalBlockEntry	(	GlobalOrdinal	globalBlockRow,
		GlobalOrdinal	globalBlockCol,
		LocalOrdinal	blkRowSize,
		LocalOrdinal	blkColSize,
		LocalOrdinal	LDA,
		const Teuchos::ArrayView< const Scalar > &	blockEntry
	)

Add the contents of the input block-entry into the matrix.

This method will create the specified block-entry if it doesn't already exist, but only if fillComplete() has not yet been called.

If the specified block-entry already exists in the matrix, the contents of the input block-entry will be added to the values that are already present.

Note that if globalBlockRow is not owned by the local processor (as indicated by getBlockRowMap()) then the block-entry is held in temporary storage until globalAssemble() is called (which is called internally by fillComplete()) and then globalAssemble() performans the communication needed to move the data to the owning processor.

This method may be called any time (before or after fillComplete()).

Definition at line 1134 of file Tpetra_VbrMatrix_def.hpp.

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

void Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::sumIntoLocalBlockEntry	(	LocalOrdinal	localBlockRow,
		LocalOrdinal	localBlockCol,
		LocalOrdinal	blkRowSize,
		LocalOrdinal	blkColSize,
		LocalOrdinal	LDA,
		const Teuchos::ArrayView< const Scalar > &	blockEntry
	)

Add the contents of the input block-entry into the matrix.

This method will throw an exception if fillComplete() has not yet been called, or if the specified block-entry doesn't already exist in the matrix.

The contents of the input block-entry will be added to the values that are already present in the matrix.

Definition at line 1151 of file Tpetra_VbrMatrix_def.hpp.

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

void Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::fillComplete	(	const Teuchos::RCP< const BlockMap< LocalOrdinal, GlobalOrdinal, Node > > &	blockDomainMap,
		const Teuchos::RCP< const BlockMap< LocalOrdinal, GlobalOrdinal, Node > > &	blockRangeMap
	)

Transition the matrix to the packed, optimized-storage state.

This method also sets the domain and range maps. This method internally calls globalAssemble().

Definition at line 1323 of file Tpetra_VbrMatrix_def.hpp.

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

void Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::fillComplete ( )

Transition the matrix to the packed, optimized-storage state.

This method internally calls fillComplete(getBlockRowMap(),getBlockRowMap()).

Definition at line 1343 of file Tpetra_VbrMatrix_def.hpp.

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

void Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::globalAssemble ( )

Communicate non-local contributions to the processors that own those contributions.

Definition at line 1353 of file Tpetra_VbrMatrix_def.hpp.

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

void Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::getGlobalBlockRowView	(	GlobalOrdinal	globalBlockRow,
		LocalOrdinal &	numPtRows,
		Teuchos::ArrayView< const GlobalOrdinal > &	blockCols,
		Teuchos::Array< LocalOrdinal > &	ptColsPerBlockCol,
		Teuchos::Array< Teuchos::ArrayRCP< const Scalar > > &	blockEntries
	)		const

Returns a const read-only view of the block-entries in the specified row.

Can only be called if isFillComplete()==false

Definition at line 358 of file Tpetra_VbrMatrix_def.hpp.

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

void Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::getLocalBlockRowView	(	LocalOrdinal	localBlockRow,
		LocalOrdinal &	numPtRows,
		Teuchos::ArrayView< const LocalOrdinal > &	blockCols,
		Teuchos::Array< LocalOrdinal > &	ptColsPerBlockCol,
		Teuchos::ArrayRCP< const Scalar > &	blockEntries
	)		const

Returns a const read-only view of the block-entries in the specified row.

Can only be called if isFillComplete()==true

Definition at line 385 of file Tpetra_VbrMatrix_def.hpp.

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

void Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::getGlobalBlockEntryView	(	GlobalOrdinal	globalBlockRow,
		GlobalOrdinal	globalBlockCol,
		LocalOrdinal &	numPtRows,
		LocalOrdinal &	numPtCols,
		Teuchos::ArrayRCP< const Scalar > &	blockEntry
	)		const

Returns a const read-only view of a block-entry.

The arguments numPtRows and numPtCols are set to the dimensions of the block- entry on output. The stride (LDA in Blas terminology) is equal to numPtRows.

This method may be called any time (before or after fillComplete()), but will throw an exception if the specified block-entry doesn't already exist.

Definition at line 485 of file Tpetra_VbrMatrix_def.hpp.

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

void Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::getGlobalBlockEntryViewNonConst	(	GlobalOrdinal	globalBlockRow,
		GlobalOrdinal	globalBlockCol,
		LocalOrdinal &	numPtRows,
		LocalOrdinal &	numPtCols,
		Teuchos::ArrayRCP< Scalar > &	blockEntry
	)

Returns a non-const read-write view of a block-entry.

Creates the block-entry if it doesn't already exist, and if:

the arguments numPtRows and numPtCols are set on entry (and nonzero),
and if fillComplete() has not yet been called.

Important Note: Be very careful managing the lifetime of this view. If fillComplete() has been called, and if you are running on a GPU, this view may be a copy of memory from the GPU, and your changes to the view won't be copied back to the GPU until your ArrayRCP is destroyed or set to Teuchos::null.

Definition at line 419 of file Tpetra_VbrMatrix_def.hpp.

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

void Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::getLocalBlockEntryView	(	LocalOrdinal	localBlockRow,
		LocalOrdinal	localBlockCol,
		LocalOrdinal &	numPtRows,
		LocalOrdinal &	numPtCols,
		Teuchos::ArrayRCP< const Scalar > &	blockEntry
	)		const

Returns a const read-only view of a block-entry.

The arguments numPtRows and numPtCols are set to the dimensions of the block- entry on output. The stride (LDA in Blas terminology) is equal to numPtRows. Throws an exception if fillComplete() has not yet been called, or if the specified block-entry doesn't exist.

This method may only be called after fillComplete() has been called, and will throw an exception if the specified block-entry doesn't already exist.

Definition at line 929 of file Tpetra_VbrMatrix_def.hpp.

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

void Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::getLocalBlockEntryViewNonConst	(	LocalOrdinal	localBlockRow,
		LocalOrdinal	localBlockCol,
		LocalOrdinal &	numPtRows,
		LocalOrdinal &	numPtCols,
		Teuchos::ArrayRCP< Scalar > &	blockEntry
	)

Returns a non-const read-write view of a block-entry.

The arguments numPtRows and numPtCols are set to the dimensions of the block- entry on output. The stride (LDA in Blas terminology) is equal to numPtRows. Throws an exception if fillComplete() has not yet been called, or if the specified block-entry doesn't exist.

Important Note: Be very careful managing the lifetime of this view. If fillComplete() has been called, and if you are running on a GPU, this view may be a copy of memory from the GPU, and your changes to the view won't be copied back to the GPU until your ArrayRCP is destroyed or set to Teuchos::null.

This method may only be called after fillComplete() has been called, and will throw an exception if the specified block-entry doesn't already exist.

Definition at line 533 of file Tpetra_VbrMatrix_def.hpp.

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

void Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::getLocalDiagCopy ( Tpetra::Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > & diag ) const

Return a copy of the (point-entry) diagonal values.

Throws an exception if the input Vector's Map is not the same as getBlockRowMap()->getPointMap().

Definition at line 580 of file Tpetra_VbrMatrix_def.hpp.

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

bool Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::checkSizes ( const SrcDistObject & source ) [virtual]

Compare the source and target (this) objects for compatibility.

Returns:: True if they are compatible, else false.

Implements Tpetra::DistObject< char, LocalOrdinal, GlobalOrdinal, Node >.

Definition at line 607 of file Tpetra_VbrMatrix_def.hpp.

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

void Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::copyAndPermute	(	const SrcDistObject &	source,
		size_t	numSameIDs,
		const Teuchos::ArrayView< const LocalOrdinal > &	permuteToLIDs,
		const Teuchos::ArrayView< const LocalOrdinal > &	permuteFromLIDs
	)		`[virtual]`

Perform copies and permutations that are local to this process.

Parameters:

source	[in] On entry, the source object, from which we are distributing. We distribute to the destination object, which is `*this` object.
numSameIDs	[in] The umber of elements that are the same on the source and destination (this) objects. These elements are owned by the same process in both the source and destination objects. No permutation occurs.
numPermuteIDs	[in] The number of elements that are locally permuted between the source and destination objects.
permuteToLIDs	[in] List of the elements that are permuted. They are listed by their LID in the destination object.
permuteFromLIDs	[in] List of the elements that are permuted. They are listed by their LID in the source object.

Implements Tpetra::DistObject< char, LocalOrdinal, GlobalOrdinal, Node >.

Definition at line 638 of file Tpetra_VbrMatrix_def.hpp.

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

void Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::packAndPrepare	(	const SrcDistObject &	source,
		const Teuchos::ArrayView< const LocalOrdinal > &	exportLIDs,
		Teuchos::Array< char > &	exports,
		const Teuchos::ArrayView< size_t > &	numPacketsPerLID,
		size_t &	constantNumPackets,
		Distributor &	distor
	)		`[virtual]`

Perform any packing or preparation required for communication.

Parameters:

source	[in] Source object for the redistribution.
exportLIDs	[in] List of the entries (as local IDs in the source object) we will be sending to other images.
exports	[out] On exit, the buffer for data to send.
numPacketsPerLID	[out] On exit, the implementation of this method must do one of two things: set numPacketsPerLID[i] to contain the number of packets to be exported for exportLIDs[i] and set constantNumPackets to zero, or set constantNumPackets to a nonzero value. If the latter, the implementation need not fill numPacketsPerLID.
constantNumPackets	[out] On exit, 0 if numPacketsPerLID has variable contents (different size for each LID). If nonzero, then it is expected that the number of packets per LID is constant, and that constantNumPackets is that value.
distor	[in] The Distributor object we are using.

Implements Tpetra::DistObject< char, LocalOrdinal, GlobalOrdinal, Node >.

Definition at line 724 of file Tpetra_VbrMatrix_def.hpp.

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

void Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::unpackAndCombine	(	const Teuchos::ArrayView< const LocalOrdinal > &	importLIDs,
		const Teuchos::ArrayView< const char > &	imports,
		const Teuchos::ArrayView< size_t > &	numPacketsPerLID,
		size_t	constantNumPackets,
		Distributor &	distor,
		CombineMode	CM
	)		`[virtual]`

Perform any unpacking and combining after communication.

Parameters:

importLIDs	[in] List of the entries (as LIDs in the destination object) we received from other images.
imports	[in] Buffer containing data we received.
numPacketsPerLID	[in] If constantNumPackets is zero, then numPacketsPerLID[i] contains the number of packets imported for importLIDs[i].
constantNumPackets	[in] If nonzero, then numPacketsPerLID is constant (same value in all entries) and constantNumPackets is that value. If zero, then numPacketsPerLID[i] is the number of packets imported for importLIDs[i].
distor	[in] The Distributor object we are using.
CM	[in] The combine mode to use when combining the imported entries with existing entries.

Implements Tpetra::DistObject< char, LocalOrdinal, GlobalOrdinal, Node >.

Definition at line 878 of file Tpetra_VbrMatrix_def.hpp.

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

std::string Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::description ( ) const [virtual]

One-line descriptiion of this object.

We declare this method virtual so that subclasses of DistObject may override it.

Reimplemented from Tpetra::DistObject< char, LocalOrdinal, GlobalOrdinal, Node >.

Definition at line 1392 of file Tpetra_VbrMatrix_def.hpp.

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

void Tpetra::VbrMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node >::describe	(	Teuchos::FancyOStream &	out,
		const Teuchos::EVerbosityLevel	verbLevel = `Teuchos::Describable::verbLevel_default`
	)		const `[virtual]`

Print the object with some verbosity level to a FancyOStream object.

Reimplemented from Tpetra::DistObject< char, LocalOrdinal, GlobalOrdinal, Node >.

Definition at line 1410 of file Tpetra_VbrMatrix_def.hpp.

void Tpetra::DistObject< char , LocalOrdinal, GlobalOrdinal, Node >::doImport	(	const SrcDistObject< char, LocalOrdinal, GlobalOrdinal, Node > &	source,
		const Import< LocalOrdinal, GlobalOrdinal, Node > &	importer,
		CombineMode	CM
	)		`[inherited]`

Import data into this object using an Import object ("forward mode").

The input DistObject is always the source of the data redistribution operation, and the *this object is always the target.

If you don't know the difference between forward and reverse mode, then you probably want forward mode. Use this method with your precomputed Import object if you want to do an Import, else use doExport() with a precomputed Export object.

Parameters:

source	[in] The "source" object for redistribution.
importer	[in] Precomputed data redistribution plan. Its source Map must be the same as the input DistObject's Map, and its target Map must be the same as `this->getMap()`.
CM	[in] How to combine incoming data with the same global index.

void Tpetra::DistObject< char , LocalOrdinal, GlobalOrdinal, Node >::doImport	(	const SrcDistObject< char, LocalOrdinal, GlobalOrdinal, Node > &	source,
		const Export< LocalOrdinal, GlobalOrdinal, Node > &	exporter,
		CombineMode	CM
	)		`[inherited]`

Import data into this object using an Export object ("reverse mode").

The input DistObject is always the source of the data redistribution operation, and the *this object is always the target.

If you don't know the difference between forward and reverse mode, then you probably want forward mode. Use the version of doImport() that takes a precomputed Import object in that case.

Parameters:

source	[in] The "source" object for redistribution.
exporter	[in] Precomputed data redistribution plan. Its target Map must be the same as the input DistObject's Map, and its source Map must be the same as `this->getMap()`. (Note the difference from forward mode.)
CM	[in] How to combine incoming data with the same global index.

void Tpetra::DistObject< char , LocalOrdinal, GlobalOrdinal, Node >::doExport	(	const SrcDistObject< char, LocalOrdinal, GlobalOrdinal, Node > &	source,
		const Export< LocalOrdinal, GlobalOrdinal, Node > &	exporter,
		CombineMode	CM
	)		`[inherited]`

Export data into this object using an Export object ("forward mode").

The input DistObject is always the source of the data redistribution operation, and the *this object is always the target.

If you don't know the difference between forward and reverse mode, then you probably want forward mode. Use this method with your precomputed Export object if you want to do an Export, else use doImport() with a precomputed Import object.

Parameters:

source	[in] The "source" object for redistribution.
exporter	[in] Precomputed data redistribution plan. Its source Map must be the same as the input DistObject's Map, and its target Map must be the same as `this->getMap()`.
CM	[in] How to combine incoming data with the same global index.

void Tpetra::DistObject< char , LocalOrdinal, GlobalOrdinal, Node >::doExport	(	const SrcDistObject< char, LocalOrdinal, GlobalOrdinal, Node > &	source,
		const Import< LocalOrdinal, GlobalOrdinal, Node > &	importer,
		CombineMode	CM
	)		`[inherited]`

Export data into this object using an Import object ("reverse mode").

The input DistObject is always the source of the data redistribution operation, and the *this object is always the target.

If you don't know the difference between forward and reverse mode, then you probably want forward mode. Use the version of doExport() that takes a precomputed Export object in that case.

Parameters:

source	[in] The "source" object for redistribution.
importer	[in] Precomputed data redistribution plan. Its target Map must be the same as the input DistObject's Map, and its source Map must be the same as `this->getMap()`. (Note the difference from forward mode.)
CM	[in] How to combine incoming data with the same global index.

bool Tpetra::DistObject< char , LocalOrdinal, GlobalOrdinal, Node >::isDistributed ( ) const [inherited]

Whether this is a globally distributed object.

For a definition of "globally distributed" (and its opposite, "locally replicated"), see the documentation of Map's isDistributed() method.

virtual Teuchos::RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> > Tpetra::DistObject< char , LocalOrdinal, GlobalOrdinal, Node >::getMap ( ) const [inline, virtual, inherited]

The Map describing the parallel distribution of this object.

Note that some Tpetra objects might be distributed using multiple Map objects. For example, CrsMatrix has both a row Map and a column Map. It is up to the subclass to decide which Map to use when invoking the DistObject constructor.

Definition at line 316 of file Tpetra_DistObject_decl.hpp.

void Tpetra::DistObject< char , LocalOrdinal, GlobalOrdinal, Node >::print ( std::ostream & os ) const [inherited]

Print this object to the given output stream.

We generally assume that all MPI processes can print to the given stream.

virtual void Tpetra::DistObject< char , LocalOrdinal, GlobalOrdinal, Node >::removeEmptyProcessesInPlace ( const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > & newMap ) [virtual, inherited]

Remove processes which contain no elements in this object's Map.

Warning:: This method is ONLY for use by experts. We highly recommend using the nonmember function of the same name defined in this file.; We make NO promises of backwards compatibility. This method may change or disappear at any time.

On input, this object is distributed over the Map returned by getMap() (the "original Map," with its communicator, the "original communicator"). The input newMap of this method must be the same as the result of calling getMap()->removeEmptyProcesses(). On processes in the original communicator which contain zero elements ("excluded processes," as opposed to "included processes"), the input newMap must be Teuchos::null (which is what getMap()->removeEmptyProcesses() returns anyway).

On included processes, reassign this object's Map (that would be returned by getMap()) to the input newMap, and do any work that needs to be done to restore correct semantics. On excluded processes, free any data that needs freeing, and do any other work that needs to be done to restore correct semantics.

This method has collective semantics over the original communicator. On exit, the only method of this object which is safe to call on excluded processes is the destructor. This implies that subclasses' destructors must not contain communication operations.

Returns:: The object's new Map. Its communicator is a new communicator, distinct from the old Map's communicator, which contains a subset of the processes in the old communicator.

Note:: The name differs from Map's method removeEmptyProcesses(), in order to emphasize that the operation on DistObject happens in place, modifying the input, whereas the operation removeEmptyProcess() on Map does not modify the input.; To implementers of DistObject subclasses: The default implementation of this class throws std::logic_error.

virtual size_t Tpetra::DistObject< char , LocalOrdinal, GlobalOrdinal, Node >::constantNumberOfPackets ( ) const [protected, virtual, inherited]

Whether the implementation's instance promises always to have a constant number of packets per LID, and if so, how many packets per LID there are.

If this method returns zero, the instance says that it might possibly have a different number of packets for each LID to send or receive. If it returns nonzero, the instance promises that the number of packets is the same for all LIDs, and that the return value is this number of packets per LID.

The default implementation of this method returns zero. This does not affect the behavior of doTransfer() in any way. If a nondefault implementation returns nonzero, doTransfer() will use this information to avoid unnecessary allocation and / or resizing of arrays.

virtual void Tpetra::DistObject< char , LocalOrdinal, GlobalOrdinal, Node >::doTransfer	(	const SrcDistObject< char, LocalOrdinal, GlobalOrdinal, Node > &	src,
		CombineMode	CM,
		size_t	numSameIDs,
		const Teuchos::ArrayView< const LocalOrdinal > &	permuteToLIDs,
		const Teuchos::ArrayView< const LocalOrdinal > &	permuteFromLIDs,
		const Teuchos::ArrayView< const LocalOrdinal > &	remoteLIDs,
		const Teuchos::ArrayView< const LocalOrdinal > &	exportLIDs,
		Distributor &	distor,
		ReverseOption	revOp
	)		`[protected, virtual, inherited]`

Redistribute data across memory images.

Parameters:

src	[in] The source object, to redistribute into the target object, which is `*this` object.
CM	[in] The combine mode that describes how to combine values that map to the same global ID on the same process.
permuteToLIDs	[in] See copyAndPermute().
permuteFromLIDs	[in] See copyAndPermute().
remoteLIDs	[in] List of entries (as local IDs) in the destination object to receive from other processes.
exportLIDs	[in] See packAndPrepare().
distor	[in/out] The Distributor object that knows how to redistribute data.
revOp	[in] Whether to do a forward or reverse mode redistribution.

virtual void Tpetra::DistObject< char , LocalOrdinal, GlobalOrdinal, Node >::createViews ( ) const [protected, virtual, inherited]

Hook for creating a const view.

doTransfer() calls this on the source object. By default, it does nothing, but the source object can use this as a hint to fetch data from a compute buffer on an off-CPU device (such as a GPU) into host memory.

virtual void Tpetra::DistObject< char , LocalOrdinal, GlobalOrdinal, Node >::createViewsNonConst ( KokkosClassic::ReadWriteOption rwo ) [protected, virtual, inherited]

Hook for creating a nonconst view.

doTransfer() calls this on the destination (*this) object. By default, it does nothing, but the destination object can use this as a hint to fetch data from a compute buffer on an off-CPU device (such as a GPU) into host memory.

Parameters:

rwo [in] Whether to create a write-only or a read-and-write view. For Kokkos Node types where compute buffers live in a separate memory space (e.g., in the device memory of a discrete accelerator like a GPU), a write-only view only requires copying from host memory to the compute buffer, whereas a read-and-write view requires copying both ways (once to read, from the compute buffer to host memory, and once to write, back to the compute buffer).

virtual void Tpetra::DistObject< char , LocalOrdinal, GlobalOrdinal, Node >::releaseViews ( ) const [protected, virtual, inherited]

Hook for releasing views.

doTransfer() calls this on both the source and destination objects, once it no longer needs to access that object's data. By default, this method does nothing. Implementations may use this as a hint to free host memory which is a view of a compute buffer, once the host memory view is no longer needed. Some implementations may prefer to mirror compute buffers in host memory; for these implementations, releaseViews() may do nothing.

Member Data Documentation

Teuchos::RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> > Tpetra::DistObject< char , LocalOrdinal, GlobalOrdinal, Node >::map_ [protected, inherited]

The Map over which this object is distributed.

Definition at line 611 of file Tpetra_DistObject_decl.hpp.

Teuchos::Array<char > Tpetra::DistObject< char , LocalOrdinal, GlobalOrdinal, Node >::imports_ [protected, inherited]

Buffer into which packed data are imported (received from other processes).

Definition at line 615 of file Tpetra_DistObject_decl.hpp.

Teuchos::Array<size_t> Tpetra::DistObject< char , LocalOrdinal, GlobalOrdinal, Node >::numImportPacketsPerLID_ [protected, inherited]

Number of packets to receive for each receive operation.

This array is used in Distributor::doPosts() (and doReversePosts()) when starting the ireceive operation.

This may be ignored in doTransfer() if constantNumPackets is nonzero, indicating a constant number of packets per LID. (For example, MultiVector sets the constantNumPackets output argument of packAndPrepare() to the number of columns in the multivector.)

Definition at line 627 of file Tpetra_DistObject_decl.hpp.

Teuchos::Array<char > Tpetra::DistObject< char , LocalOrdinal, GlobalOrdinal, Node >::exports_ [protected, inherited]

Buffer from which packed data are exported (sent to other processes).

Definition at line 630 of file Tpetra_DistObject_decl.hpp.

Teuchos::Array<size_t> Tpetra::DistObject< char , LocalOrdinal, GlobalOrdinal, Node >::numExportPacketsPerLID_ [protected, inherited]

Number of packets to send for each send operation.

This array is used in Distributor::doPosts() (and doReversePosts()) for preparing for the send operation.

This may be ignored in doTransfer() if constantNumPackets is nonzero, indicating a constant number of packets per LID. (For example, MultiVector sets the constantNumPackets output argument of packAndPrepare() to the number of columns in the multivector.)

Definition at line 642 of file Tpetra_DistObject_decl.hpp.

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

Public Types

Public Member Functions

Protected Member Functions

Protected Attributes

Methods for use only by experts

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation