42 #ifndef THYRA_EPETRA_LINEAR_OP_HPP 43 #define THYRA_EPETRA_LINEAR_OP_HPP 45 #include "Thyra_LinearOpBase.hpp" 47 #include "Thyra_ScaledLinearOpBase.hpp" 48 #include "Thyra_RowStatLinearOpBase.hpp" 49 #include "Thyra_SpmdVectorSpaceBase.hpp" 51 #include "Epetra_RowMatrix.h" 79 :
virtual public LinearOpBase<double>,
80 virtual public ScaledLinearOpBase<double>,
81 virtual public RowStatLinearOpBase<double>,
148 const RCP<Epetra_Operator> &op,
152 const RCP<
const VectorSpaceBase<double> > &
range = Teuchos::null,
153 const RCP<
const VectorSpaceBase<double> > &
domain = Teuchos::null
205 const RCP<
const VectorSpaceBase<double> > &
range,
206 const RCP<
const VectorSpaceBase<double> > &
domain,
207 const RCP<Epetra_Operator> &op,
231 RCP<Epetra_Operator> *op= NULL,
232 EOpTransp *opTrans = NULL,
235 RCP<
const VectorSpaceBase<double> > *
range = NULL,
236 RCP<
const VectorSpaceBase<double> > *
domain = NULL
247 RCP<const SpmdVectorSpaceBase<double> >
spmdRange()
const;
257 RCP<const SpmdVectorSpaceBase<double> >
spmdDomain()
const;
263 RCP<const Epetra_Operator>
epetra_op()
const;
272 const Ptr<RCP<Epetra_Operator> > &epetraOp,
273 const Ptr<EOpTransp> &epetraOpTransp,
274 const Ptr<EApplyEpetraOpAs> &epetraOpApplyAs,
275 const Ptr<EAdjointEpetraOp> &epetraOpAdjointSupport
279 const Ptr<RCP<const Epetra_Operator> > &epetraOp,
280 const Ptr<EOpTransp> &epetraOpTransp,
281 const Ptr<EApplyEpetraOpAs> &epetraOpApplyAs,
282 const Ptr<EAdjointEpetraOp> &epetraOpAdjointSupport
291 RCP<const VectorSpaceBase<double> >
range()
const;
294 RCP<const VectorSpaceBase<double> >
domain()
const;
297 RCP<const LinearOpBase<double> >
clone()
const;
309 const Teuchos::EVerbosityLevel verbLevel
324 const EOpTransp M_trans,
325 const MultiVectorBase<double> &X,
326 const Ptr<MultiVectorBase<double> > &Y,
343 virtual void scaleLeftImpl(
const VectorBase<double> &row_scaling);
346 virtual void scaleRightImpl(
const VectorBase<double> &col_scaling);
355 const RowStatLinearOpBaseUtils::ERowStat rowStat)
const;
359 const RowStatLinearOpBaseUtils::ERowStat rowStat,
360 const Ptr<VectorBase<double> > &rowStatVec)
const;
377 virtual RCP< const SpmdVectorSpaceBase<double> >
379 const RCP<Epetra_Operator> &op,
393 virtual RCP< const SpmdVectorSpaceBase<double> >
395 const RCP<Epetra_Operator> &op,
412 RCP<const SpmdVectorSpaceBase<double> >
range_;
413 RCP<const SpmdVectorSpaceBase<double> >
domain_;
437 RCP<EpetraLinearOp> nonconstEpetraLinearOp();
445 partialNonconstEpetraLinearOp(
446 const RCP<
const VectorSpaceBase<double> > &range,
447 const RCP<
const VectorSpaceBase<double> > &domain,
448 const RCP<Epetra_Operator> &op,
449 EOpTransp opTrans = NOTRANS,
461 nonconstEpetraLinearOp(
462 const RCP<Epetra_Operator> &op,
463 EOpTransp opTrans = NOTRANS,
466 const RCP<
const VectorSpaceBase<double> > &range = Teuchos::null,
467 const RCP<
const VectorSpaceBase<double> > &domain = Teuchos::null
476 RCP<const EpetraLinearOp>
478 const RCP<const Epetra_Operator> &op,
479 EOpTransp opTrans = NOTRANS,
482 const RCP<
const VectorSpaceBase<double> > &range = Teuchos::null,
483 const RCP<
const VectorSpaceBase<double> > &domain = Teuchos::null
493 nonconstEpetraLinearOp(
494 const RCP<Epetra_Operator> &op,
495 const std::string &label,
496 EOpTransp opTrans = NOTRANS,
499 const RCP<
const VectorSpaceBase<double> > &range = Teuchos::null,
500 const RCP<
const VectorSpaceBase<double> > &domain = Teuchos::null
509 RCP<const EpetraLinearOp>
511 const RCP<const Epetra_Operator> &op,
512 const std::string &label,
513 EOpTransp opTrans = NOTRANS,
516 const RCP<
const SpmdVectorSpaceBase<double> > &range = Teuchos::null,
517 const RCP<
const SpmdVectorSpaceBase<double> > &domain = Teuchos::null
524 #endif // THYRA_EPETRA_LINEAR_OP_HPP const Epetra_Map & getRangeMap() const
Concrete LinearOpBase adapter subclass for Epetra_Operator object.
virtual RCP< const SpmdVectorSpaceBase< double > > allocateRange(const RCP< Epetra_Operator > &op, EOpTransp op_trans) const
Allocate the range space of the operator.
virtual bool supportsScaleLeftImpl() const
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.
virtual void getRowStatImpl(const RowStatLinearOpBaseUtils::ERowStat rowStat, const Ptr< VectorBase< double > > &rowStatVec) const
Abstract base class for all LinearOpBase objects that can return an Epetra_Operator view of themselve...
RCP< const SpmdVectorSpaceBase< double > > spmdDomain() const
Return a smart pointer to the SpmdVectorSpaceBase object for the domain.
void computeAbsRowSum(Epetra_Vector &rowStatVec_in) const
Compute the absolute row sum for this matrix.
void getEpetraOpView(const Ptr< RCP< const Epetra_Operator > > &epetraOp, const Ptr< EOpTransp > &epetraOpTransp, const Ptr< EApplyEpetraOpAs > &epetraOpApplyAs, const Ptr< EAdjointEpetraOp > &epetraOpAdjointSupport) const
RCP< const SpmdVectorSpaceBase< double > > domain_
RCP< const SpmdVectorSpaceBase< double > > range_
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.
void describe(FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel) const
virtual void scaleLeftImpl(const VectorBase< double > &row_scaling)
EAdjointEpetraOp adjointSupport_
RCP< const LinearOpBase< double > > clone() const
EAdjointEpetraOp
Determine if adjoints are supported on Epetra_Opeator or not.
EpetraLinearOp()
Construct to uninitialized.
void setFullyInitialized(bool isFullyInitialized=true)
Set to fully initialized.
void getNonconstEpetraOpView(const Ptr< RCP< Epetra_Operator > > &epetraOp, const Ptr< EOpTransp > &epetraOpTransp, const Ptr< EApplyEpetraOpAs > &epetraOpApplyAs, const Ptr< EAdjointEpetraOp > &epetraOpAdjointSupport)
virtual bool rowStatIsSupportedImpl(const RowStatLinearOpBaseUtils::ERowStat rowStat) const
virtual RCP< const SpmdVectorSpaceBase< double > > allocateDomain(const RCP< Epetra_Operator > &op, EOpTransp op_trans) const
Allocate the domain space of the operator.
RCP< const VectorSpaceBase< double > > range() const
RCP< const VectorSpaceBase< double > > domain() const
RCP< const SpmdVectorSpaceBase< double > > spmdRange() const
Return a smart pointer to the SpmdVectorSpaceBase object for the range.
EApplyEpetraOpAs
Determine how the apply an Epetra_Operator as a linear operator.
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.
bool opSupportedImpl(EOpTransp M_trans) const
Apply using Epetra_Operator::Apply(...)
virtual void scaleRightImpl(const VectorBase< double > &col_scaling)
virtual bool supportsScaleRightImpl() const
EApplyEpetraOpAs applyAs_
const Epetra_Map & getDomainMap() const
RCP< Epetra_Operator > op_
RCP< Epetra_RowMatrix > rowMatrix_
std::string description() const
void applyImpl(const EOpTransp M_trans, const MultiVectorBase< double > &X, const Ptr< MultiVectorBase< double > > &Y, const double alpha, const double beta) const
RCP< Epetra_Operator > epetra_op()