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 Teuchos::RCP< Thyra::VectorBase< double > > get_x_th() const
void clear()
Wipe out stored data.
virtual void initialize()
void set_f(const Teuchos::RCP< Epetra_Vector > &in)
virtual void set_f_th(const Teuchos::RCP< Thyra::VectorBase< double > > &in)
Teuchos::RCP< Epetra_CrsMatrix > A
Teuchos::RCP< Epetra_Vector > x
void set_x(const Teuchos::RCP< Epetra_Vector > &in)
EpetraLinearObjContainer(const Teuchos::RCP< const Epetra_Map > &domain, const Teuchos::RCP< const Epetra_Map > &range)
const Teuchos::RCP< Epetra_CrsMatrix > get_A() const
virtual void set_x_th(const Teuchos::RCP< Thyra::VectorBase< double > > &in)
virtual void set_dxdt_th(const Teuchos::RCP< Thyra::VectorBase< double > > &in)
EpetraLinearObjContainer()
Epetra_CrsMatrix CrsMatrixType
Teuchos::RCP< const Epetra_Map > domainMap
const Teuchos::RCP< Epetra_Vector > get_dxdt() const
Teuchos::RCP< const Thyra::VectorSpaceBase< double > > domainSpace
void initializeMatrix(double value)
Put a particular scalar in the matrix.
Teuchos::RCP< const Epetra_Map > rangeMap
virtual Teuchos::RCP< Thyra::LinearOpBase< double > > get_A_th() const
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< Epetra_Vector > dxdt
void set_A(const Teuchos::RCP< Epetra_CrsMatrix > &in)
const Teuchos::RCP< Epetra_Vector > get_f() const
virtual void set_A_th(const Teuchos::RCP< Thyra::LinearOpBase< double > > &in)
Teuchos::RCP< const Thyra::VectorSpaceBase< double > > rangeSpace
virtual Teuchos::RCP< Thyra::VectorBase< double > > get_dxdt_th() const
const Teuchos::RCP< Epetra_Vector > get_x() const
virtual Teuchos::RCP< Thyra::VectorBase< double > > get_f_th() const
void set_dxdt(const Teuchos::RCP< Epetra_Vector > &in)
LinearObjContainer::Members Members
Teuchos::RCP< Epetra_Vector > f