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" 57 #include "Teuchos_RCP.hpp" 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" 79 const Teuchos::RCP<const Epetra_Map> & range)
88 const Teuchos::RCP<const Epetra_Map> & range,
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);
113 inline void set_x(
const Teuchos::RCP<Epetra_Vector> & in) {
x = in; }
114 inline const Teuchos::RCP<Epetra_Vector>
get_x()
const {
return x; }
116 inline void set_dxdt(
const Teuchos::RCP<Epetra_Vector> & in) {
dxdt = in; }
117 inline const Teuchos::RCP<Epetra_Vector>
get_dxdt()
const {
return dxdt; }
119 inline void set_f(
const Teuchos::RCP<Epetra_Vector> & in) {
f = in; }
120 inline const Teuchos::RCP<Epetra_Vector>
get_f()
const {
return f; }
122 inline void set_A(
const Teuchos::RCP<Epetra_CrsMatrix> & in) {
A = in; }
123 inline const Teuchos::RCP<Epetra_CrsMatrix>
get_A()
const {
return A; }
126 {
A->PutScalar(value); }
129 {
x = (in==Teuchos::null) ? Teuchos::null : Thyra::get_Epetra_Vector(*
domainMap,in); }
130 virtual Teuchos::RCP<Thyra::VectorBase<double> >
get_x_th()
const 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); }
135 virtual Teuchos::RCP<Thyra::VectorBase<double> >
get_dxdt_th()
const 139 {
f = (in==Teuchos::null) ? Teuchos::null : Thyra::get_Epetra_Vector(*
rangeMap,in); }
140 virtual Teuchos::RCP<Thyra::VectorBase<double> >
get_f_th()
const 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); }
145 virtual Teuchos::RCP<Thyra::LinearOpBase<double> >
get_A_th()
const 146 {
return (
A==Teuchos::null) ? Teuchos::null : Thyra::nonconstEpetraLinearOp(
A); }
152 Teuchos::RCP<const Thyra::VectorSpaceBase<double> >
rangeSpace;
154 Teuchos::RCP<Epetra_CrsMatrix>
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