Concrete LinearOpBase adapter subclass for Epetra_Operator object. More...

#include <Thyra_EpetraLinearOp.hpp>

Inheritance diagram for Thyra::EpetraLinearOp:

Private Member Functions
const Epetra_Map &	getRangeMap () const

const Epetra_Map &	getDomainMap () const

void	computeAbsRowSum (Epetra_Vector &rowStatVec_in) const
	Compute the absolute row sum for this matrix. More...

Private Attributes
bool	isFullyInitialized_

RCP< Epetra_Operator >	op_

RCP< Epetra_RowMatrix >	rowMatrix_

EOpTransp	opTrans_

EApplyEpetraOpAs	applyAs_

EAdjointEpetraOp	adjointSupport_

RCP< const SpmdVectorSpaceBase < double > >	range_

RCP< const SpmdVectorSpaceBase < double > >	domain_

Related Functions
(Note that these are not member functions.)
RCP< EpetraLinearOp >	nonconstEpetraLinearOp ()
	Default nonmember constructor. More...

RCP< EpetraLinearOp >	partialNonconstEpetraLinearOp (const RCP< const VectorSpaceBase< double > > &range, const RCP< const VectorSpaceBase< double > > &domain, const RCP< Epetra_Operator > &op, EOpTransp opTrans=NOTRANS, EApplyEpetraOpAs applyAs=EPETRA_OP_APPLY_APPLY, EAdjointEpetraOp adjointSupport=EPETRA_OP_ADJOINT_SUPPORTED)
	Partially initialized EpetraLinearOp. More...

RCP< EpetraLinearOp >	nonconstEpetraLinearOp (const RCP< Epetra_Operator > &op, EOpTransp opTrans=NOTRANS, EApplyEpetraOpAs applyAs=EPETRA_OP_APPLY_APPLY, EAdjointEpetraOp adjointSupport=EPETRA_OP_ADJOINT_SUPPORTED, const RCP< const VectorSpaceBase< double > > &range=Teuchos::null, const RCP< const VectorSpaceBase< double > > &domain=Teuchos::null)
	Dynamically allocate an const EpetraLinearOp to wrap a const Epetra_Operator object. More...

RCP< const EpetraLinearOp >	epetraLinearOp (const RCP< const Epetra_Operator > &op, EOpTransp opTrans=NOTRANS, EApplyEpetraOpAs applyAs=EPETRA_OP_APPLY_APPLY, EAdjointEpetraOp adjointSupport=EPETRA_OP_ADJOINT_SUPPORTED, const RCP< const VectorSpaceBase< double > > &range=Teuchos::null, const RCP< const VectorSpaceBase< double > > &domain=Teuchos::null)
	Dynamically allocate a nonconst EpetraLinearOp to wrap a const Epetra_Operator object. More...

RCP< EpetraLinearOp >	nonconstEpetraLinearOp (const RCP< Epetra_Operator > &op, const std::string &label, EOpTransp opTrans=NOTRANS, EApplyEpetraOpAs applyAs=EPETRA_OP_APPLY_APPLY, EAdjointEpetraOp adjointSupport=EPETRA_OP_ADJOINT_SUPPORTED, const RCP< const VectorSpaceBase< double > > &range=Teuchos::null, const RCP< const VectorSpaceBase< double > > &domain=Teuchos::null)
	Dynamically allocate an const EpetraLinearOp to wrap a const Epetra_Operator object and give it a string label. More...

RCP< const EpetraLinearOp >	epetraLinearOp (const RCP< const Epetra_Operator > &op, const std::string &label, EOpTransp opTrans=NOTRANS, EApplyEpetraOpAs applyAs=EPETRA_OP_APPLY_APPLY, EAdjointEpetraOp adjointSupport=EPETRA_OP_ADJOINT_SUPPORTED, const RCP< const SpmdVectorSpaceBase< double > > &range=Teuchos::null, const RCP< const SpmdVectorSpaceBase< double > > &domain=Teuchos::null)
	Dynamically allocate a nonconst EpetraLinearOp to wrap a const Epetra_Operator object. More...

Constructors / initializers / accessors
	EpetraLinearOp ()
	Construct to uninitialized. More...

void	initialize (const RCP< Epetra_Operator > &op, EOpTransp opTrans=NOTRANS, EApplyEpetraOpAs applyAs=EPETRA_OP_APPLY_APPLY, EAdjointEpetraOp adjointSupport=EPETRA_OP_ADJOINT_SUPPORTED, const RCP< const VectorSpaceBase< double > > &range=Teuchos::null, const RCP< const VectorSpaceBase< double > > &domain=Teuchos::null)
	Fully initialize. More...

void	partiallyInitialize (const RCP< const VectorSpaceBase< double > > &range, const RCP< const VectorSpaceBase< double > > &domain, const RCP< Epetra_Operator > &op, EOpTransp opTrans=NOTRANS, EApplyEpetraOpAs applyAs=EPETRA_OP_APPLY_APPLY, EAdjointEpetraOp adjointSupport=EPETRA_OP_ADJOINT_SUPPORTED)
	Partially initialize. More...

void	setFullyInitialized (bool isFullyInitialized=true)
	Set to fully initialized. More...

void	uninitialize (RCP< Epetra_Operator > op=NULL, EOpTransp opTrans=NULL, EApplyEpetraOpAs applyAs=NULL, EAdjointEpetraOp adjointSupport=NULL, RCP< const VectorSpaceBase< double > > range=NULL, RCP< const VectorSpaceBase< double > > domain=NULL)
	Set to uninitialized and optionally return the current state. More...

RCP< const SpmdVectorSpaceBase < double > >	spmdRange () const
	Return a smart pointer to the SpmdVectorSpaceBase object for the range. More...

RCP< const SpmdVectorSpaceBase < double > >	spmdDomain () const
	Return a smart pointer to the SpmdVectorSpaceBase object for the domain. More...

RCP< Epetra_Operator >	epetra_op ()

RCP< const Epetra_Operator >	epetra_op () const

Overridden from EpetraLinearOpBase
void	getNonconstEpetraOpView (const Ptr< RCP< Epetra_Operator > > &epetraOp, const Ptr< EOpTransp > &epetraOpTransp, const Ptr< EApplyEpetraOpAs > &epetraOpApplyAs, const Ptr< EAdjointEpetraOp > &epetraOpAdjointSupport)

void	getEpetraOpView (const Ptr< RCP< const Epetra_Operator > > &epetraOp, const Ptr< EOpTransp > &epetraOpTransp, const Ptr< EApplyEpetraOpAs > &epetraOpApplyAs, const Ptr< EAdjointEpetraOp > &epetraOpAdjointSupport) const

Overridden from LinearOpBase
RCP< const VectorSpaceBase < double > >	range () const

RCP< const VectorSpaceBase < double > >	domain () const

RCP< const LinearOpBase< double > >	clone () const

Overridden from Teuchos::Describable
std::string	description () const

void	describe (FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel) const

Protected member functions overridden from LinearOpBase.
bool	opSupportedImpl (EOpTransp M_trans) const

void	applyImpl (const EOpTransp M_trans, const MultiVectorBase< double > &X, const Ptr< MultiVectorBase< double > > &Y, const double alpha, const double beta) const

Protected member functions overridden from ScaledLinearOpBase.
virtual bool	supportsScaleLeftImpl () const

virtual bool	supportsScaleRightImpl () const

virtual void	scaleLeftImpl (const VectorBase< double > &row_scaling)

virtual void	scaleRightImpl (const VectorBase< double > &col_scaling)

Protected member functions overridden from RowStatLinearOpBase.
virtual bool	rowStatIsSupportedImpl (const RowStatLinearOpBaseUtils::ERowStat rowStat) const

virtual void	getRowStatImpl (const RowStatLinearOpBaseUtils::ERowStat rowStat, const Ptr< VectorBase< double > > &rowStatVec) const

Allocators for domain and range spaces
virtual RCP< const SpmdVectorSpaceBase< double > >	allocateDomain (const RCP< Epetra_Operator > &op, EOpTransp op_trans) const
	Allocate the domain space of the operator. More...

virtual RCP< const SpmdVectorSpaceBase< double > >	allocateRange (const RCP< Epetra_Operator > &op, EOpTransp op_trans) const
	Allocate the range space of the operator. More...

Additional Inherited Members
Public Member Functions inherited from Thyra::EpetraLinearOpBase
Public Member Functions inherited from Teuchos::Describable
DescribableStreamManipulatorState	describe (const Describable &describable, const EVerbosityLevel verbLevel=Describable::verbLevel_default)

std::ostream &	operator<< (std::ostream &os, const DescribableStreamManipulatorState &d)

virtual void	describe (FancyOStream &out, const EVerbosityLevel verbLevel=verbLevel_default) const

void	describe (std::ostream &out, const EVerbosityLevel verbLevel=verbLevel_default) const

virtual	~Describable ()

	LabeledObject ()

virtual	~LabeledObject ()

virtual void	setObjectLabel (const std::string &objectLabel)

virtual std::string	getObjectLabel () const

Static Public Attributes inherited from Teuchos::Describable
static const EVerbosityLevel	verbLevel_default

Detailed Description

Concrete LinearOpBase adapter subclass for Epetra_Operator object.

This subclass can be used to represent the non-transposed operator or transposed operator defined by an Epetra_Operator object. This class can implement apply() using either Epetra_Operator::Apply() or Epetra_Operator::ApplyInverse(). In addition, the user can specify whether adjoints are supported or not.

Partial Automatic Change Propagation: This class shall maintain no state with respect to the values of the internally stored Epetra_Operator object. Therefore, as long as the domain and range maps do not change, the the Epetra_Operator can be changed and this will automatically update *this object. This simplifies some types of update operations. Since this is a simple concrete class, this is harmless. However, if the range and domain maps change, then one must call the this->initialize() function.

Definition at line 78 of file Thyra_EpetraLinearOp.hpp.

Constructor & Destructor Documentation

Thyra::EpetraLinearOp::EpetraLinearOp ( )

Construct to uninitialized.

See the postconditions for uninitialize()

Definition at line 65 of file Thyra_EpetraLinearOp.cpp.

References NOTRANS.

Member Function Documentation

void Thyra::EpetraLinearOp::initialize	(	const RCP< Epetra_Operator > &	op,
		EOpTransp	opTrans = `NOTRANS`,
		EApplyEpetraOpAs	applyAs = `EPETRA_OP_APPLY_APPLY`,
		EAdjointEpetraOp	adjointSupport = `EPETRA_OP_ADJOINT_SUPPORTED`,
		const RCP< const VectorSpaceBase< double > > &	range = `Teuchos::null`,
		const RCP< const VectorSpaceBase< double > > &	domain = `Teuchos::null`
	)

Fully initialize.

Parameters

op	[in] The `Epetra_Operator` this `*this` will wrap.
opTrans	[in] If `opTrans==NOTRANS` then `op` will be viewed as `op` and if `opTrans==TRANS` then `op` will be viewed as its transpose `op'` for the behavior of `apply()`.
applyAs	[in] If `applyAs==APPLY_APPLY` then `op->Apply()` will be used and if `applyAs==APPLY_APPLY_INVERSE` then `op->ApplyInverse()` is used instead.
adjointSupport	[in] Determines if it is to be assumed that adjoints are supported on the underlying `Epetra_Operator` object `op`. If `adjointSupport==EPETRA_OP_ADJOINT_SUPPORTED` then `this->opSupported(TRANS)` will return `true`. If `adjointSupport==EPETRA_OP_ADJOINT_UNSUPPORTED` then `this->opSupported(TRANS)` will return `false`.
range	[in] Smart pointer to the range space for the `Epetra_Operator`. The default value is `Teuchos::null` in which case `*this` will allocate a new `SpmdVectorSpace` given range map from `op`. A client may only bother to specify this space if one wants to override the defintion of the scalar product.
domain	[in] Smart pointer to the domain space for the `Epetra_Operator`. The default value is `Teuchos::null` in which case `*this` will allocate a new `DefaultSpmdVectorSpace` given map from `op`. A client may only bother to specify this space if one wants to override the defintion of the scalar product.

Preconditions:

!is_null(op)

Postconditions:

this->epetra_op().get() == op.get()
[range.get() != NULL] this->range().get() == range.get()
[domain.get() != NULL] this->domain().get() == domain.get()
[range.get() == NULL] this->range().get() != NULL
[domain.get() == NULL] this->domain().get() != NULL
this->opSupported(NOTRANS) == true
this->opSupported(TRNAS) == adjointSupport==EPETRA_OP_ADJOINT_SUPPORTED

After this function is called, this will be fully initialized and ready to go.

Definition at line 73 of file Thyra_EpetraLinearOp.cpp.

References adjointSupport_, allocateDomain(), allocateRange(), applyAs_, domain_, Thyra::EPETRA_OP_APPLY_APPLY, is_null(), isFullyInitialized_, op_, opTrans_, range_, rowMatrix_, and TEUCHOS_TEST_FOR_EXCEPTION.

Referenced by Thyra::EpetraExtAddTransformer::transform(), Thyra::EpetraExtDiagScalingTransformer::transform(), and Thyra::EpetraExtDiagScaledMatProdTransformer::transform().

void Thyra::EpetraLinearOp::partiallyInitialize	(	const RCP< const VectorSpaceBase< double > > &	range,
		const RCP< const VectorSpaceBase< double > > &	domain,
		const RCP< Epetra_Operator > &	op,
		EOpTransp	opTrans = `NOTRANS`,
		EApplyEpetraOpAs	applyAs = `EPETRA_OP_APPLY_APPLY`,
		EAdjointEpetraOp	adjointSupport = `EPETRA_OP_ADJOINT_SUPPORTED`
	)

Partially initialize.

Parameters

range	[in] Smart pointer to the range space for the `Epetra_Operator`.
domain	[in] Smart pointer to the domain space for the `Epetra_Operator`.
op	[in] The `Epetra_Operator` this `*this` will wrap. This object is assumed to not be fully unitialized.
opTrans	[in] If `opTrans==NOTRANS` then `op` will be viewed as `op` and if `opTrans==TRANS` then `op` will be viewed as its transpose `op'` for the behavior of `apply()`.
applyAs	[in] If `applyAs==APPLY_APPLY` then `op->Apply()` will be used and if `applyAs==APPLY_APPLY_INVERSE` then `op->ApplyInverse()` is used instead.
adjointSupport	[in] Determines if it is to be assumed that adjoints are supported on the underlying `Epetra_Operator` object `op`. If `adjointSupport==EPETRA_OP_ADJOINT_SUPPORTED` then `this->opSupported(TRANS)` will return `true`. If `adjointSupport==EPETRA_OP_ADJOINT_UNSUPPORTED` then `this->opSupported(TRANS)` will return `false`.

Preconditions:

!is_null(range)
!is_null(domain)
!is_null(op)

Postconditions:

this->epetra_op().get() == op.get()
[range.get() != NULL] this->range().get() == range.get()
[domain.get() != NULL] this->domain().get() == domain.get()
[range.get() == NULL] this->range().get() != NULL
[domain.get() == NULL] this->domain().get() != NULL
this->opSupported(NOTRANS) == true
this->opSupported(TRNAS) == adjointSupport==EPETRA_OP_ADJOINT_SUPPORTED

After this function is called, only the range and domain spaces will be supported and this must be followed up by a call to setFullyInitialized().

Definition at line 120 of file Thyra_EpetraLinearOp.cpp.

References adjointSupport_, applyAs_, domain_, is_null(), isFullyInitialized_, op_, opTrans_, range_, rowMatrix_, and TEUCHOS_TEST_FOR_EXCEPTION.

void Thyra::EpetraLinearOp::setFullyInitialized ( bool isFullyInitialized = true )

Set to fully initialized.

In debug mode, asserts will be performed to ensure that everything matches up as it should.

The functions initialize() ore partiallyInitialize() must have been called prior to calling this function.

Definition at line 161 of file Thyra_EpetraLinearOp.cpp.

References isFullyInitialized_.

void Thyra::EpetraLinearOp::uninitialize	(	RCP< Epetra_Operator > *	op = `NULL`,
		EOpTransp *	opTrans = `NULL`,
		EApplyEpetraOpAs *	applyAs = `NULL`,
		EAdjointEpetraOp *	adjointSupport = `NULL`,
		RCP< const VectorSpaceBase< double > > *	range = `NULL`,
		RCP< const VectorSpaceBase< double > > *	domain = `NULL`
	)

Set to uninitialized and optionally return the current state.

Postconditions:

this->domain().get() == NULL
this->range().get() == NULL

Definition at line 168 of file Thyra_EpetraLinearOp.cpp.

References adjointSupport_, applyAs_, domain_, Thyra::EPETRA_OP_ADJOINT_SUPPORTED, Thyra::EPETRA_OP_APPLY_APPLY, isFullyInitialized_, NOTRANS, Teuchos::null, op_, opTrans_, range_, and rowMatrix_.

RCP< const SpmdVectorSpaceBase< double > > Thyra::EpetraLinearOp::spmdRange ( ) const

Return a smart pointer to the SpmdVectorSpaceBase object for the range.

Postconditions:

[this->range().get() != NULL] return.get() != NULL
[this->range().get() == NULL] return.get() == NULL

Definition at line 198 of file Thyra_EpetraLinearOp.cpp.

References range_.

RCP< const SpmdVectorSpaceBase< double > > Thyra::EpetraLinearOp::spmdDomain ( ) const

Return a smart pointer to the SpmdVectorSpaceBase object for the domain.

Postconditions:

[this->domain().get() != NULL] return.get() != NULL
[this->domain().get() == NULL] return.get() == NULL

Definition at line 205 of file Thyra_EpetraLinearOp.cpp.

References domain_.

RCP< Epetra_Operator > Thyra::EpetraLinearOp::epetra_op ( )

Definition at line 212 of file Thyra_EpetraLinearOp.cpp.

References op_.

Referenced by Thyra::get_Epetra_Operator(), Thyra::EpetraExtAddTransformer::transform(), Thyra::EpetraExtDiagScalingTransformer::transform(), and Thyra::EpetraExtDiagScaledMatProdTransformer::transform().

RCP< const Epetra_Operator > Thyra::EpetraLinearOp::epetra_op ( ) const

Definition at line 219 of file Thyra_EpetraLinearOp.cpp.

References op_.

void Thyra::EpetraLinearOp::getNonconstEpetraOpView	(	const Ptr< RCP< Epetra_Operator > > &	epetraOp,
		const Ptr< EOpTransp > &	epetraOpTransp,
		const Ptr< EApplyEpetraOpAs > &	epetraOpApplyAs,
		const Ptr< EAdjointEpetraOp > &	epetraOpAdjointSupport
	)

virtual

Implements Thyra::EpetraLinearOpBase.

Definition at line 228 of file Thyra_EpetraLinearOp.cpp.

References adjointSupport_, applyAs_, op_, and opTrans_.

void Thyra::EpetraLinearOp::getEpetraOpView	(	const Ptr< RCP< const Epetra_Operator > > &	epetraOp,
		const Ptr< EOpTransp > &	epetraOpTransp,
		const Ptr< EApplyEpetraOpAs > &	epetraOpApplyAs,
		const Ptr< EAdjointEpetraOp > &	epetraOpAdjointSupport
	)		const

virtual

Implements Thyra::EpetraLinearOpBase.

Definition at line 242 of file Thyra_EpetraLinearOp.cpp.

References adjointSupport_, applyAs_, op_, and opTrans_.

RCP< const VectorSpaceBase< double > > Thyra::EpetraLinearOp::range ( ) const

Definition at line 260 of file Thyra_EpetraLinearOp.cpp.

References range_.

Referenced by describe(), and description().

RCP< const VectorSpaceBase< double > > Thyra::EpetraLinearOp::domain ( ) const

Definition at line 267 of file Thyra_EpetraLinearOp.cpp.

References domain_.

Referenced by describe(), and description().

RCP< const LinearOpBase< double > > Thyra::EpetraLinearOp::clone ( ) const

Definition at line 274 of file Thyra_EpetraLinearOp.cpp.

References Teuchos::null.

std::string Thyra::EpetraLinearOp::description ( ) const

virtual

Reimplemented from Teuchos::Describable.

Definition at line 284 of file Thyra_EpetraLinearOp.cpp.

References Teuchos::Describable::description(), domain(), Teuchos::RCP< T >::get(), op_, range(), and typeName().

Referenced by describe().

void Thyra::EpetraLinearOp::describe	(	FancyOStream &	out,
		const Teuchos::EVerbosityLevel	verbLevel
	)		const

Definition at line 301 of file Thyra_EpetraLinearOp.cpp.

References adjointSupport_, applyAs_, Teuchos::as(), Teuchos::Describable::description(), description(), domain(), Teuchos::RCP< T >::get(), Teuchos::includesVerbLevel(), includesVerbLevel(), is_null(), op_, opTrans_, range(), Thyra::toString(), typeName(), Teuchos::VERB_EXTREME, Teuchos::VERB_HIGH, Teuchos::VERB_LOW, and Teuchos::VERB_MEDIUM.

bool Thyra::EpetraLinearOp::opSupportedImpl ( EOpTransp M_trans ) const

protected

Definition at line 352 of file Thyra_EpetraLinearOp.cpp.

References adjointSupport_, Thyra::EPETRA_OP_ADJOINT_SUPPORTED, isFullyInitialized_, and NOTRANS.

void Thyra::EpetraLinearOp::applyImpl	(	const EOpTransp	M_trans,
		const MultiVectorBase< double > &	X,
		const Ptr< MultiVectorBase< double > > &	Y,
		const double	alpha,
		const double	beta
	)		const

protected

bool Thyra::EpetraLinearOp::supportsScaleLeftImpl ( ) const

protectedvirtual

Definition at line 524 of file Thyra_EpetraLinearOp.cpp.

References nonnull(), and rowMatrix_.

bool Thyra::EpetraLinearOp::supportsScaleRightImpl ( ) const

protectedvirtual

Definition at line 530 of file Thyra_EpetraLinearOp.cpp.

References nonnull(), and rowMatrix_.

void Thyra::EpetraLinearOp::scaleLeftImpl ( const VectorBase< double > & row_scaling )

protectedvirtual

Definition at line 536 of file Thyra_EpetraLinearOp.cpp.

References Thyra::get_Epetra_Vector(), getRangeMap(), and rowMatrix_.

void Thyra::EpetraLinearOp::scaleRightImpl ( const VectorBase< double > & col_scaling )

protectedvirtual

Definition at line 545 of file Thyra_EpetraLinearOp.cpp.

References Thyra::get_Epetra_Vector(), getDomainMap(), and rowMatrix_.

bool Thyra::EpetraLinearOp::rowStatIsSupportedImpl ( const RowStatLinearOpBaseUtils::ERowStat rowStat ) const

protectedvirtual

Definition at line 557 of file Thyra_EpetraLinearOp.cpp.

References is_null(), rowMatrix_, TEUCHOS_TEST_FOR_EXCEPT, and TEUCHOS_UNREACHABLE_RETURN.

void Thyra::EpetraLinearOp::getRowStatImpl	(	const RowStatLinearOpBaseUtils::ERowStat	rowStat,
		const Ptr< VectorBase< double > > &	rowStatVec
	)		const

protectedvirtual

Definition at line 574 of file Thyra_EpetraLinearOp.cpp.

References computeAbsRowSum(), Thyra::get_Epetra_Vector(), getRangeMap(), rowMatrix_, and TEUCHOS_TEST_FOR_EXCEPT.

RCP< const SpmdVectorSpaceBase< double > > Thyra::EpetraLinearOp::allocateDomain	(	const RCP< Epetra_Operator > &	op,
		EOpTransp	op_trans
	)		const

protectedvirtual

Allocate the domain space of the operator.

Purpose: In TSFExtended, both EpetraLinearOp and EpetraVectorSpace are extended from the Thyra versions by inheritance, and the TSFExtended operator subclasses expect to work with an extended vector space subclass. Thus, it is necessary for the base operator class to never directly allocate vector space objects, and allocation is delegated to a virtual allocator function.

Definition at line 600 of file Thyra_EpetraLinearOp.cpp.

References Thyra::create_VectorSpace(), and Teuchos::rcp().

Referenced by initialize().

RCP< const SpmdVectorSpaceBase< double > > Thyra::EpetraLinearOp::allocateRange	(	const RCP< Epetra_Operator > &	op,
		EOpTransp	op_trans
	)		const

protectedvirtual

Allocate the range space of the operator.

Purpose: In TSFExtended, both EpetraLinearOp and EpetraVectorSpace are extended from the Thyra versions by inheritance, and the TSFExtended operator subclasses expect to work with an extended vector space subclass. Thus, it is necessary for the base operator class to never directly allocate vector space objects, and allocation is delegated to a virtual allocator function.

Definition at line 613 of file Thyra_EpetraLinearOp.cpp.

References Thyra::create_VectorSpace(), and Teuchos::rcp().

Referenced by initialize().

const Epetra_Map & Thyra::EpetraLinearOp::getRangeMap ( ) const

private

Definition at line 627 of file Thyra_EpetraLinearOp.cpp.

References applyAs_, Thyra::EPETRA_OP_APPLY_APPLY, and op_.

Referenced by applyImpl(), getRowStatImpl(), and scaleLeftImpl().

const Epetra_Map & Thyra::EpetraLinearOp::getDomainMap ( ) const

private

Definition at line 635 of file Thyra_EpetraLinearOp.cpp.

References applyAs_, Thyra::EPETRA_OP_APPLY_APPLY, and op_.

Referenced by applyImpl(), and scaleRightImpl().

void Thyra::EpetraLinearOp::computeAbsRowSum ( Epetra_Vector & rowStatVec_in ) const

private

Compute the absolute row sum for this matrix.

A concrete implementation is required because Epetra does not support absolute row sums.

Note: This only works for Epetra_CrsMatrix objects.

Definition at line 642 of file Thyra_EpetraLinearOp.cpp.

References is_null(), rowMatrix_, TEUCHOS_ASSERT, TEUCHOS_TEST_FOR_EXCEPT, and TEUCHOS_TEST_FOR_EXCEPTION.

Referenced by getRowStatImpl().

Friends And Related Function Documentation

RCP< EpetraLinearOp > nonconstEpetraLinearOp ( )

RCP< EpetraLinearOp > partialNonconstEpetraLinearOp	(	const RCP< const VectorSpaceBase< double > > &	range,
		const RCP< const VectorSpaceBase< double > > &	domain,
		const RCP< Epetra_Operator > &	op,
		EOpTransp	opTrans = `NOTRANS`,
		EApplyEpetraOpAs	applyAs = `EPETRA_OP_APPLY_APPLY`,
		EAdjointEpetraOp	adjointSupport = `EPETRA_OP_ADJOINT_SUPPORTED`
	)

RCP< EpetraLinearOp > nonconstEpetraLinearOp	(	const RCP< Epetra_Operator > &	op,
		EOpTransp	opTrans = `NOTRANS`,
		EApplyEpetraOpAs	applyAs = `EPETRA_OP_APPLY_APPLY`,
		EAdjointEpetraOp	adjointSupport = `EPETRA_OP_ADJOINT_SUPPORTED`,
		const RCP< const VectorSpaceBase< double > > &	range = `Teuchos::null`,
		const RCP< const VectorSpaceBase< double > > &	domain = `Teuchos::null`
	)

RCP< const EpetraLinearOp > epetraLinearOp	(	const RCP< const Epetra_Operator > &	op,
		EOpTransp	opTrans = `NOTRANS`,
		EApplyEpetraOpAs	applyAs = `EPETRA_OP_APPLY_APPLY`,
		EAdjointEpetraOp	adjointSupport = `EPETRA_OP_ADJOINT_SUPPORTED`,
		const RCP< const VectorSpaceBase< double > > &	range = `Teuchos::null`,
		const RCP< const VectorSpaceBase< double > > &	domain = `Teuchos::null`
	)

RCP< EpetraLinearOp > nonconstEpetraLinearOp	(	const RCP< Epetra_Operator > &	op,
		const std::string &	label,
		EOpTransp	opTrans = `NOTRANS`,
		EApplyEpetraOpAs	applyAs = `EPETRA_OP_APPLY_APPLY`,
		EAdjointEpetraOp	adjointSupport = `EPETRA_OP_ADJOINT_SUPPORTED`,
		const RCP< const VectorSpaceBase< double > > &	range = `Teuchos::null`,
		const RCP< const VectorSpaceBase< double > > &	domain = `Teuchos::null`
	)

RCP< const EpetraLinearOp > epetraLinearOp	(	const RCP< const Epetra_Operator > &	op,
		const std::string &	label,
		EOpTransp	opTrans = `NOTRANS`,
		EApplyEpetraOpAs	applyAs = `EPETRA_OP_APPLY_APPLY`,
		EAdjointEpetraOp	adjointSupport = `EPETRA_OP_ADJOINT_SUPPORTED`,
		const RCP< const SpmdVectorSpaceBase< double > > &	range = `Teuchos::null`,
		const RCP< const SpmdVectorSpaceBase< double > > &	domain = `Teuchos::null`
	)