Thyra::DefaultScaledAdjointLinearOp< Scalar > Class Template Reference
[Development of linear Abstract Numerical Algorithms (ANAs)]

Concrete decorator LinearOpBase subclass that wraps a LinearOpBase object and adds on an extra scaling factor and/or a transpose enum. More...

#include <Thyra_DefaultScaledAdjointLinearOpDecl.hpp>

Inheritance diagram for Thyra::DefaultScaledAdjointLinearOp< Scalar >:

[legend]List of all members.


Constructors/initializers/accessors
	DefaultScaledAdjointLinearOp ()
	Constructs to uninitialized.
	DefaultScaledAdjointLinearOp (const Scalar &scalar, const ETransp &transp, const RCP< LinearOpBase< Scalar > > &Op)
	Calls `initialize()`.
	DefaultScaledAdjointLinearOp (const Scalar &scalar, const ETransp &transp, const RCP< const LinearOpBase< Scalar > > &Op)
	Calls `initialize()`.
void	initialize (const Scalar &scalar, const ETransp &transp, const RCP< LinearOpBase< Scalar > > &Op)
	Initialize with an operator with by defining adjoint (transpose) and scaling arguments.
void	initialize (const Scalar &scalar, const ETransp &transp, const RCP< const LinearOpBase< Scalar > > &Op)
	Initialize with an operator with by defining adjoint (transpose) and scaling arguments.
RCP< LinearOpBase< Scalar > >	getNonconstOp ()
	Return the non-const linear operator passed into `initialize()`.
RCP< const LinearOpBase< Scalar > >	getOp () const
	Return the const linear operator passed into `initialize()`.
void	uninitialize ()
	Set to uninitialized and (optionally) extract the objects passed into `initialize()`.
Overridden from Teuchos::Describable
std::string	description () const
	Outputs `DefaultScaledAdjointLinearOp<Scalar>{this->getOrigOp().description())` along with the dimensions.
void	describe (Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel) const
	Prints out the original operator as well as all of the scalings and transpositions in the order that they occurred.
Overridden from OpBase
RCP< const VectorSpaceBase< Scalar > >	range () const
	Return the range space of the logical linear operator.
RCP< const VectorSpaceBase< Scalar > >	domain () const
	Return the domain space of the logical linear operator.
RCP< const LinearOpBase< Scalar > >	clone () const

bool	opSupported (ETransp M_trans) const
Overridden from LinearOpBase
void	apply (const ETransp M_trans, const MultiVectorBase< Scalar > &X, MultiVectorBase< Scalar > *Y, const Scalar alpha, const Scalar beta) const
	Apply the linear operator (or its transpose) to a multi-vector : `Y = alphaop(M)X + beta*Y`.
Overridden from ScaledAdointLinearOpBase
Scalar	overallScalar () const

ETransp	overallTransp () const

RCP< LinearOpBase< Scalar > >	getNonconstOrigOp ()

RCP< const LinearOpBase< Scalar > >	getOrigOp () const

Public Member Functions
	TEUCHOS_GET_RCP (LinearOpBase< Scalar >)

Related Functions
(Note that these are not member functions.)
RCP< LinearOpBase< Scalar > >	nonconstScale (const Scalar &scalar, const RCP< LinearOpBase< Scalar > > &Op)
	Build an implicit non-`const` scaled linear operator.
RCP< const LinearOpBase< Scalar > >	scale (const Scalar &scalar, const RCP< const LinearOpBase< Scalar > > &Op)
	Build an implicit `const` scaled linear operator.
RCP< LinearOpBase< Scalar > >	nonconstAdjoint (const RCP< LinearOpBase< Scalar > > &Op)
	Build an implicit non-`const` adjoined linear operator.
RCP< const LinearOpBase< Scalar > >	adjoint (const RCP< const LinearOpBase< Scalar > > &Op)
	Build an implicit `const` adjoined linear operator.
RCP< LinearOpBase< Scalar > >	nonconstTranspose (const RCP< LinearOpBase< Scalar > > &Op)
	Build an implicit non-`const` transposed linear operator.
RCP< const LinearOpBase< Scalar > >	transpose (const RCP< const LinearOpBase< Scalar > > &Op)
	Build an implicit `const` transposed linear operator.
RCP< LinearOpBase< Scalar > >	nonconstScaleAndAdjoint (const Scalar &scalar, const ETransp &transp, const RCP< LinearOpBase< Scalar > > &Op)
	Build an implicit non-`const` scaled and/or adjoined (transposed) linear operator.
RCP< const LinearOpBase< Scalar > >	scaleAndAdjoint (const Scalar &scalar, const ETransp &transp, const RCP< const LinearOpBase< Scalar > > &Op)
	Build an implicit `const` scaled and/or adjoined (transposed) linear operator.

Detailed Description

template<class Scalar>
class Thyra::DefaultScaledAdjointLinearOp< Scalar >

Concrete decorator LinearOpBase subclass that wraps a LinearOpBase object and adds on an extra scaling factor and/or a transpose enum.

This class represents a scaled, adjointed (transposed) linear operator M of the form:

  M = scalar * op(Op)

where Op is another LinearOpBase object, scalar is a Scalar, and the operation op(Op) is specified by a ETransp and is given as op(Op) = Op (NOTRANS), or op(Op) = Op^T (TRANS), or op(Op) = Op^H (CONJTRANS). Of course the operator M is not constructed explicitly but instead just applies the decorated operator Op by modifying the apply() function that calls Op.apply().

This subclass is designed to allow the efficient handling of multiple implicit scalings and/or adjoints (transposes) and allow these implicit transformations to be reversed. A sequence of scalings/adjoints from some original LinearOpBase object origOp is shown as:

 M = scalar_n * op_n( ... scalar_2 * op_2( scalar_1 * op_1( origOp ) ) ... )

   =>

 M  = overallScalar * overall_op(origOp)

where overallScalar and overall_op(...) are the cumulative transformations given as:

 overallScalar = scalar_n * ... * scalar_2 * ... * scalar_1

 overall_op(origOp) = op_n( ... op_2( op_1( ... op_n( origOp ) ... ) ) )

Each individual transformation pair (scalar_i,op_i(...)) is specified with arguments Scalar scalar and ETransp transp. The overall scaling is returned using this->overallScalar(), the overall adjoint enum is returned using this->overallTransp(), and the original linear operator is returned from <tt>this->getOrigOp().

The operator most recently wrapped is returned from <tt>this->getOp(). The individual scalings and transformations are not exposed from this interface but can be viewed by calling this->description() with a verbosity level of ???. The arguments passed into the constructor DefaultScaledAdjointLinearOp() or initialize() can always be extracted using uninitialize().

This subclass keeps track of all of the individual scalings scalar_i and adjoining operations op_i(...) so that the description() function can print this out and also so that the operations can be reversed.

The copy constructor and assignment operators are declared private since some thought needs to be put into what they should mean.

Note: This class does not maintain any specific cached information about the original operator so it is safe, in some cases, to manipulate the original operator and still keep this intact and automatically updated for any changes that are made.

Examples:: test_composite_linear_ops.cpp.

Definition at line 115 of file Thyra_DefaultScaledAdjointLinearOpDecl.hpp.

Constructor & Destructor Documentation

template<class Scalar>

Thyra::DefaultScaledAdjointLinearOp< Scalar >::DefaultScaledAdjointLinearOp ( ) [inline]

Constructs to uninitialized.

Postconditions:

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

Definition at line 548 of file Thyra_DefaultScaledAdjointLinearOpDecl.hpp.

template<class Scalar>

Thyra::DefaultScaledAdjointLinearOp< Scalar >::DefaultScaledAdjointLinearOp	(	const Scalar &	scalar,
		const ETransp &	transp,
		const RCP< LinearOpBase< Scalar > > &	Op
	)			`[inline]`

Calls initialize().

Note, instead of calling this constructor directly consider using the non-member functions described belos which create dynamically allocated RCP-wrapped objects.

Definition at line 556 of file Thyra_DefaultScaledAdjointLinearOpDecl.hpp.

template<class Scalar>

Thyra::DefaultScaledAdjointLinearOp< Scalar >::DefaultScaledAdjointLinearOp	(	const Scalar &	scalar,
		const ETransp &	transp,
		const RCP< const LinearOpBase< Scalar > > &	Op
	)			`[inline]`

Calls initialize().

Note, instead of calling this constructor directly consider using the non-member functions described belos which create dynamically allocated RCP-wrapped objects.

Definition at line 570 of file Thyra_DefaultScaledAdjointLinearOpDecl.hpp.

Member Function Documentation

template<class Scalar>

Thyra::DefaultScaledAdjointLinearOp< Scalar >::TEUCHOS_GET_RCP ( LinearOpBase< Scalar > )

template<class Scalar>

void Thyra::DefaultScaledAdjointLinearOp< Scalar >::initialize	(	const Scalar &	scalar,
		const ETransp &	transp,
		const RCP< LinearOpBase< Scalar > > &	Op
	)

Initialize with an operator with by defining adjoint (transpose) and scaling arguments.

Parameters:

	scalar	[in] Scalar argument defining `scalar_0` (see introduction).
	transp	[in] Transpose argument defining `op_0(...)` (see introduction).
	Op	[in] Smart pointer to linear operator (persisting relationship).

Preconditions:

Op.get() != NULL

Postconditions:

ToDo: Fill these in!!!!

Definition at line 40 of file Thyra_DefaultScaledAdjointLinearOp.hpp.

template<class Scalar>

void Thyra::DefaultScaledAdjointLinearOp< Scalar >::initialize	(	const Scalar &	scalar,
		const ETransp &	transp,
		const RCP< const LinearOpBase< Scalar > > &	Op
	)

Initialize with an operator with by defining adjoint (transpose) and scaling arguments.

Parameters:

	scalar	[in] Scalar argument defining `scalar_0` (see introduction).
	transp	[in] Transpose argument defining `op_0(...)` (see introduction).
	Op	[in] Smart pointer to linear operator (persisting relationship).

Preconditions:

Op.get() != NULL

Postconditions:

ToDo: Fill these in!!!!

Definition at line 50 of file Thyra_DefaultScaledAdjointLinearOp.hpp.

template<class Scalar>

Teuchos::RCP< LinearOpBase< Scalar > > Thyra::DefaultScaledAdjointLinearOp< Scalar >::getNonconstOp ( )

Return the non-const linear operator passed into initialize().

Definition at line 61 of file Thyra_DefaultScaledAdjointLinearOp.hpp.

template<class Scalar>

Teuchos::RCP< const LinearOpBase< Scalar > > Thyra::DefaultScaledAdjointLinearOp< Scalar >::getOp ( ) const

Return the const linear operator passed into initialize().

Definition at line 68 of file Thyra_DefaultScaledAdjointLinearOp.hpp.

template<class Scalar>

void Thyra::DefaultScaledAdjointLinearOp< Scalar >::uninitialize ( )

Set to uninitialized and (optionally) extract the objects passed into initialize().

Postconditions:

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

Definition at line 74 of file Thyra_DefaultScaledAdjointLinearOp.hpp.

template<class Scalar>

std::string Thyra::DefaultScaledAdjointLinearOp< Scalar >::description ( ) const [virtual]

Outputs DefaultScaledAdjointLinearOp<Scalar>{this->getOrigOp().description()) along with the dimensions.

Reimplemented from Thyra::LinearOpDefaultBase< Scalar >.

Definition at line 87 of file Thyra_DefaultScaledAdjointLinearOp.hpp.

template<class Scalar>

void Thyra::DefaultScaledAdjointLinearOp< Scalar >::describe	(	Teuchos::FancyOStream &	out,
		const Teuchos::EVerbosityLevel	verbLevel
	)			const `[virtual]`

Prints out the original operator as well as all of the scalings and transpositions in the order that they occurred.

This function outputs different levels of detail based on the value passed in for verbLevel:

ToDo: Finish documentation!

Reimplemented from Thyra::LinearOpDefaultBase< Scalar >.

Definition at line 99 of file Thyra_DefaultScaledAdjointLinearOp.hpp.

template<class Scalar>

Teuchos::RCP< const VectorSpaceBase< Scalar > > Thyra::DefaultScaledAdjointLinearOp< Scalar >::range ( ) const [virtual]

Return the range space of the logical linear operator.

Simply returns:

   return ( this->overallTransp()==NOTRANS ? this->getOrigOp()->range() : this->getOrigOp()->domain() );

Implements Thyra::LinearOpBase< Scalar >.

Definition at line 154 of file Thyra_DefaultScaledAdjointLinearOp.hpp.

template<class Scalar>

Teuchos::RCP< const VectorSpaceBase< Scalar > > Thyra::DefaultScaledAdjointLinearOp< Scalar >::domain ( ) const [virtual]

Return the domain space of the logical linear operator.

Simply returns:

   return ( this->overallTransp()==NOTRANS ? this->getOrigOp()->domain() : this->getOrigOp()->range() );

Implements Thyra::LinearOpBase< Scalar >.

Definition at line 162 of file Thyra_DefaultScaledAdjointLinearOp.hpp.

template<class Scalar>

Teuchos::RCP< const LinearOpBase< Scalar > > Thyra::DefaultScaledAdjointLinearOp< Scalar >::clone ( ) const [virtual]

Reimplemented from Thyra::LinearOpBase< Scalar >.

Definition at line 170 of file Thyra_DefaultScaledAdjointLinearOp.hpp.

template<class Scalar>

bool Thyra::DefaultScaledAdjointLinearOp< Scalar >::opSupported ( ETransp M_trans ) const [virtual]

Return if the operation is supported on the logical linear operator.

Simply returns:

   return this->getOrigOp()->opSupported(trans_trans(this->overallTransp(),M_trans));

Implements Thyra::SingleScalarLinearOpBase< Scalar >.

Definition at line 176 of file Thyra_DefaultScaledAdjointLinearOp.hpp.

template<class Scalar>

void Thyra::DefaultScaledAdjointLinearOp< Scalar >::apply	(	const ETransp	M_trans,
		const MultiVectorBase< Scalar > &	X,
		MultiVectorBase< Scalar > *	Y,
		const Scalar	alpha,
		const Scalar	beta
	)			const `[virtual]`

Apply the linear operator (or its transpose) to a multi-vector : Y = alpha*op(M)*X + beta*Y.

Simply calls:

   this->getOrigOp()->apply(trans_trans(M_trans,this->overallTransp()),X,Y,(this->overallScalar()*alpha),beta)

Implements Thyra::SingleScalarLinearOpBase< Scalar >.

Definition at line 185 of file Thyra_DefaultScaledAdjointLinearOp.hpp.

template<class Scalar>

Scalar Thyra::DefaultScaledAdjointLinearOp< Scalar >::overallScalar ( ) const [virtual]

Implements Thyra::ScaledAdjointLinearOpBase< Scalar >.

Definition at line 203 of file Thyra_DefaultScaledAdjointLinearOp.hpp.

template<class Scalar>

ETransp Thyra::DefaultScaledAdjointLinearOp< Scalar >::overallTransp ( ) const [virtual]

Implements Thyra::ScaledAdjointLinearOpBase< Scalar >.

Definition at line 209 of file Thyra_DefaultScaledAdjointLinearOp.hpp.

template<class Scalar>

Teuchos::RCP< LinearOpBase< Scalar > > Thyra::DefaultScaledAdjointLinearOp< Scalar >::getNonconstOrigOp ( ) [virtual]

Implements Thyra::ScaledAdjointLinearOpBase< Scalar >.

Definition at line 216 of file Thyra_DefaultScaledAdjointLinearOp.hpp.

template<class Scalar>

Teuchos::RCP< const LinearOpBase< Scalar > > Thyra::DefaultScaledAdjointLinearOp< Scalar >::getOrigOp ( ) const [virtual]

Implements Thyra::ScaledAdjointLinearOpBase< Scalar >.

Definition at line 223 of file Thyra_DefaultScaledAdjointLinearOp.hpp.

Friends And Related Function Documentation

template<class Scalar>

RCP< LinearOpBase< Scalar > > nonconstScale	(	const Scalar &	scalar,
		const RCP< LinearOpBase< Scalar > > &	Op
	)			`[related]`

Build an implicit non-const scaled linear operator.

Returns Teuchos::rcp(new DefaultScaledAdjointLinearOp<Scalar>(scalar,NOTRANS,Op).

Preconditions:

Op.get()!=NULL

Postconditions:

return.get()!=NULL

Definition at line 601 of file Thyra_DefaultScaledAdjointLinearOpDecl.hpp.

template<class Scalar>

RCP< const LinearOpBase< Scalar > > scale	(	const Scalar &	scalar,
		const RCP< const LinearOpBase< Scalar > > &	Op
	)			`[related]`

Build an implicit const scaled linear operator.

Returns Teuchos::rcp(new DefaultScaledAdjointLinearOp<Scalar>(scalar,NOTRANS,Op).

Preconditions:

Op.get()!=NULL

Postconditions:

return.get()!=NULL

Definition at line 611 of file Thyra_DefaultScaledAdjointLinearOpDecl.hpp.

template<class Scalar>

RCP< LinearOpBase< Scalar > > nonconstAdjoint ( const RCP< LinearOpBase< Scalar > > & Op ) [related]

Build an implicit non-const adjoined linear operator.

Returns Teuchos::rcp(new DefaultScaledAdjointLinearOp<Scalar>(Teuchos::ScalarTraits<Scalar>::one(),CONJTRANS,Op).

Preconditions:

Op.get()!=NULL

Postconditions:

return.get()!=NULL

Definition at line 623 of file Thyra_DefaultScaledAdjointLinearOpDecl.hpp.

template<class Scalar>

RCP< const LinearOpBase< Scalar > > adjoint ( const RCP< const LinearOpBase< Scalar > > & Op ) [related]

Build an implicit const adjoined linear operator.

Returns Teuchos::rcp(new DefaultScaledAdjointLinearOp<Scalar>(Teuchos::ScalarTraits<Scalar>::one(),CONJTRANS,Op).

Preconditions:

Op.get()!=NULL

Postconditions:

return.get()!=NULL

Definition at line 635 of file Thyra_DefaultScaledAdjointLinearOpDecl.hpp.

template<class Scalar>

RCP< LinearOpBase< Scalar > > nonconstTranspose ( const RCP< LinearOpBase< Scalar > > & Op ) [related]

Build an implicit non-const transposed linear operator.

Returns Teuchos::rcp(new DefaultScaledAdjointLinearOp<Scalar>(Teuchos::ScalarTraits<Scalar>::one(),TRANS,Op).

Preconditions:

Op.get()!=NULL

Postconditions:

return.get()!=NULL

Definition at line 647 of file Thyra_DefaultScaledAdjointLinearOpDecl.hpp.

template<class Scalar>

RCP< const LinearOpBase< Scalar > > transpose ( const RCP< const LinearOpBase< Scalar > > & Op ) [related]

Build an implicit const transposed linear operator.

Returns Teuchos::rcp(new DefaultScaledAdjointLinearOp<Scalar>(Teuchos::ScalarTraits<Scalar>::one(),TRANS,Op).

Preconditions:

Op.get()!=NULL

Postconditions:

return.get()!=NULL

Definition at line 659 of file Thyra_DefaultScaledAdjointLinearOpDecl.hpp.

template<class Scalar>

RCP< LinearOpBase< Scalar > > nonconstScaleAndAdjoint	(	const Scalar &	scalar,
		const ETransp &	transp,
		const RCP< LinearOpBase< Scalar > > &	Op
	)			`[related]`

Build an implicit non-const scaled and/or adjoined (transposed) linear operator.

Returns Teuchos::rcp(new DefaultScaledAdjointLinearOp<Scalar>(scale,transp,Op).

Preconditions:

Op.get()!=NULL

Postconditions:

return.get()!=NULL

Definition at line 671 of file Thyra_DefaultScaledAdjointLinearOpDecl.hpp.

template<class Scalar>

RCP< const LinearOpBase< Scalar > > scaleAndAdjoint	(	const Scalar &	scalar,
		const ETransp &	transp,
		const RCP< const LinearOpBase< Scalar > > &	Op
	)			`[related]`

Build an implicit const scaled and/or adjoined (transposed) linear operator.

Returns Teuchos::rcp(new DefaultScaledAdjointLinearOp<Scalar>(scale,transp,Op).

Preconditions:

Op.get()!=NULL

Postconditions:

return.get()!=NULL

Definition at line 682 of file Thyra_DefaultScaledAdjointLinearOpDecl.hpp.

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

Generated on Tue Oct 20 12:47:00 2009 for Thyra Operator/Vector Support by

1.4.7

Thyra::DefaultScaledAdjointLinearOp< Scalar > Class Template Reference [Development of linear Abstract Numerical Algorithms (ANAs)]

Constructors/initializers/accessors

Overridden from Teuchos::Describable

Overridden from OpBase

Overridden from LinearOpBase

Overridden from ScaledAdointLinearOpBase

Public Member Functions

Related Functions

Detailed Description

template<class Scalar> class Thyra::DefaultScaledAdjointLinearOp< Scalar >

Constructor & Destructor Documentation

Member Function Documentation

Friends And Related Function Documentation

Thyra::DefaultScaledAdjointLinearOp< Scalar > Class Template Reference
[Development of linear Abstract Numerical Algorithms (ANAs)]

template<class Scalar>
class Thyra::DefaultScaledAdjointLinearOp< Scalar >