43 #ifndef __Panzer_EpetraLinearObjContainer_hpp__ 44 #define __Panzer_EpetraLinearObjContainer_hpp__ 46 #include "PanzerDiscFE_config.hpp" 51 #include "Epetra_Vector.h" 52 #include "Epetra_CrsMatrix.h" 59 #include "Thyra_VectorBase.hpp" 60 #include "Thyra_LinearOpBase.hpp" 61 #include "Thyra_EpetraLinearOp.hpp" 62 #include "Thyra_EpetraThyraWrappers.hpp" 63 #include "Thyra_get_Epetra_Operator.hpp" 98 if(
get_x()!=Teuchos::null)
get_x()->PutScalar(0.0);
100 if(
get_f()!=Teuchos::null)
get_f()->PutScalar(0.0);
101 if(
get_A()!=Teuchos::null)
get_A()->PutScalar(0.0);
107 set_x(Teuchos::null);
109 set_f(Teuchos::null);
110 set_A(Teuchos::null);
126 {
A->PutScalar(value); }
129 {
x = (in==Teuchos::null) ? Teuchos::null : Thyra::get_Epetra_Vector(*
domainMap,in); }
131 {
return (
x==Teuchos::null) ? Teuchos::null : Thyra::create_Vector(
x,
domainSpace); }
134 {
dxdt = (in==Teuchos::null) ? Teuchos::null : Thyra::get_Epetra_Vector(*
domainMap,in); }
139 {
f = (in==Teuchos::null) ? Teuchos::null : Thyra::get_Epetra_Vector(*
rangeMap,in); }
141 {
return (
f==Teuchos::null) ? Teuchos::null : Thyra::create_Vector(
f,
rangeSpace); }
144 {
A = (in==Teuchos::null) ? Teuchos::null : Teuchos::rcp_dynamic_cast<Epetra_CrsMatrix>(Thyra::get_Epetra_Operator(*in),
true); }
146 {
return (
A==Teuchos::null) ? Teuchos::null : Thyra::nonconstEpetraLinearOp(
A); }
virtual void set_dxdt_th(const Teuchos::RCP< Thyra::VectorBase< double > > &in)
EpetraLinearObjContainer(const Teuchos::RCP< const Epetra_Map > &domain, const Teuchos::RCP< const Epetra_Map > &range)
void set_f(const Teuchos::RCP< Epetra_Vector > &in)
const Teuchos::RCP< Epetra_CrsMatrix > get_A() const
void clear()
Wipe out stored data.
virtual Teuchos::RCP< Thyra::VectorBase< double > > get_x_th() const
virtual Teuchos::RCP< Thyra::VectorBase< double > > get_f_th() const
virtual Teuchos::RCP< Thyra::LinearOpBase< double > > get_A_th() const
virtual void set_x_th(const Teuchos::RCP< Thyra::VectorBase< double > > &in)
void set_A(const Teuchos::RCP< Epetra_CrsMatrix > &in)
virtual void set_A_th(const Teuchos::RCP< Thyra::LinearOpBase< double > > &in)
Teuchos::RCP< Epetra_Vector > x
Teuchos::RCP< Epetra_Vector > dxdt
void set_dxdt(const Teuchos::RCP< Epetra_Vector > &in)
Epetra_CrsMatrix CrsMatrixType
Teuchos::RCP< const Thyra::VectorSpaceBase< double > > domainSpace
Teuchos::RCP< const Epetra_Map > domainMap
LinearObjContainer::Members Members
Teuchos::RCP< Epetra_Vector > f
const Teuchos::RCP< Epetra_Vector > get_x() const
virtual void set_f_th(const Teuchos::RCP< Thyra::VectorBase< double > > &in)
virtual void initialize()
virtual Teuchos::RCP< Thyra::VectorBase< double > > get_dxdt_th() const
Teuchos::RCP< const Thyra::VectorSpaceBase< double > > rangeSpace
Teuchos::RCP< Epetra_CrsMatrix > A
void set_x(const Teuchos::RCP< Epetra_Vector > &in)
EpetraLinearObjContainer(const Teuchos::RCP< const Epetra_Map > &domain, const Teuchos::RCP< const Thyra::VectorSpaceBase< double > > &domainS, const Teuchos::RCP< const Epetra_Map > &range, const Teuchos::RCP< const Thyra::VectorSpaceBase< double > > &rangeS)
Teuchos::RCP< const Epetra_Map > rangeMap
void initializeMatrix(double value)
Put a particular scalar in the matrix.
EpetraLinearObjContainer()
const Teuchos::RCP< Epetra_Vector > get_f() const
const Teuchos::RCP< Epetra_Vector > get_dxdt() const