43 #ifndef PANZER_EVALUATOR_SCATTER_DIRICHLET_RESIDUAL_BLOCKEDEPETRA_DECL_HPP 44 #define PANZER_EVALUATOR_SCATTER_DIRICHLET_RESIDUAL_BLOCKEDEPETRA_DECL_HPP 46 #include "Phalanx_config.hpp" 47 #include "Phalanx_Evaluator_Macros.hpp" 48 #include "Phalanx_MDField.hpp" 50 #include "Teuchos_ParameterList.hpp" 52 #include "PanzerDiscFE_config.hpp" 61 template <
typename ScalarT>
class ProductVectorBase;
67 template <
typename LocalOrdinalT,
typename GlobalOrdinalT>
68 class UniqueGlobalIndexer;
70 template <
typename LocalOrdinalT,
typename GlobalOrdinalT>
71 class BlockedDOFManager;
83 public PHX::EvaluatorDerived<panzer::Traits::Residual, TRAITS>,
90 virtual Teuchos::RCP<CloneableEvaluator>
clone(
const Teuchos::ParameterList & pl)
const 96 { std::cout <<
"unspecialized version of \"ScatterDirichletResidual_BlockedEpetra::evaluateFields\" on \""+PHX::typeAsString<EvalT>()+
"\" should not be used!" << std::endl;
97 TEUCHOS_ASSERT(
false); }
110 template<
typename TRAITS,
typename LO,
typename GO>
113 public PHX::EvaluatorDerived<panzer::Traits::Residual, TRAITS>,
120 : rowIndexers_(rIndexers) {}
124 const Teuchos::ParameterList& p,
125 bool useDiscreteAdjoint=
false);
130 void preEvaluate(
typename TRAITS::PreEvalData d);
134 virtual Teuchos::RCP<CloneableEvaluator>
clone(
const Teuchos::ParameterList & pl)
const 146 std::vector<Teuchos::RCP<const UniqueGlobalIndexer<LO,int> > >
rowIndexers_;
147 std::vector<Teuchos::RCP<const UniqueGlobalIndexer<LO,int> > >
colIndexers_;
156 Teuchos::RCP<const std::map<std::string,std::string> >
fieldMap_;
164 Teuchos::RCP<Thyra::ProductVectorBase<double> >
r_;
174 std::vector< PHX::MDField<const bool,Cell,NODE> >
applyBC_;
182 template<
typename TRAITS,
typename LO,
typename GO>
185 public PHX::EvaluatorDerived<panzer::Traits::Tangent, TRAITS>,
192 : rowIndexers_(rIndexers) {}
196 const Teuchos::ParameterList& p,
197 bool useDiscreteAdjoint=
false);
202 void preEvaluate(
typename TRAITS::PreEvalData d);
206 virtual Teuchos::RCP<CloneableEvaluator>
clone(
const Teuchos::ParameterList & pl)
const 218 std::vector<Teuchos::RCP<const UniqueGlobalIndexer<LO,int> > >
rowIndexers_;
219 std::vector<Teuchos::RCP<const UniqueGlobalIndexer<LO,int> > >
colIndexers_;
228 Teuchos::RCP<const std::map<std::string,std::string> >
fieldMap_;
236 Teuchos::RCP<Thyra::ProductVectorBase<double> >
r_;
246 std::vector< PHX::MDField<const bool,Cell,NODE> >
applyBC_;
254 template<
typename TRAITS,
typename LO,
typename GO>
257 public PHX::EvaluatorDerived<panzer::Traits::Jacobian, TRAITS>,
264 : rowIndexers_(rIndexers), colIndexers_(cIndexers) {}
268 const Teuchos::ParameterList& p,
269 bool useDiscreteAdjoint=
false);
271 void preEvaluate(
typename TRAITS::PreEvalData d);
278 virtual Teuchos::RCP<CloneableEvaluator>
clone(
const Teuchos::ParameterList & pl)
const 291 std::vector<Teuchos::RCP<const UniqueGlobalIndexer<LO,int> > >
rowIndexers_;
292 std::vector<Teuchos::RCP<const UniqueGlobalIndexer<LO,int> > >
colIndexers_;
301 Teuchos::RCP<const std::map<std::string,std::string> >
fieldMap_;
312 Teuchos::RCP<Thyra::ProductVectorBase<double> >
r_;
313 Teuchos::RCP<Thyra::BlockedLinearOpBase<double> >
Jac_;
319 std::vector< PHX::MDField<const bool,Cell,NODE> >
applyBC_;
327 #ifdef Panzer_BUILD_HESSIAN_SUPPORT
Teuchos::RCP< Thyra::BlockedLinearOpBase< double > > Jac_
std::size_t side_subcell_dim_
std::vector< Teuchos::RCP< const UniqueGlobalIndexer< LO, int > > > colIndexers_
ScatterDirichletResidual_BlockedEpetra(const std::vector< Teuchos::RCP< const UniqueGlobalIndexer< LO, int > > > &rIndexers, const std::vector< Teuchos::RCP< const UniqueGlobalIndexer< LO, int > > > &cIndexers)
std::vector< int > subFieldIds_
std::size_t local_side_id_
std::vector< Teuchos::RCP< const UniqueGlobalIndexer< LO, int > > > rowIndexers_
ScatterDirichletResidual_BlockedEpetra(const std::vector< Teuchos::RCP< const UniqueGlobalIndexer< LO, int > > > &rIndexers, const std::vector< Teuchos::RCP< const UniqueGlobalIndexer< LO, int > > > &cIndexers)
bool checkApplyBC_
If set to true, allows runtime disabling of dirichlet BCs on node-by-node basis.
std::vector< Teuchos::RCP< const UniqueGlobalIndexer< LO, int > > > rowIndexers_
Teuchos::RCP< const std::map< std::string, std::string > > fieldMap_
Teuchos::RCP< PHX::FieldTag > scatterHolder_
std::size_t side_subcell_dim_
std::vector< int > subFieldIds_
void postRegistrationSetup(typename TRAITS::SetupData d, PHX::FieldManager< TRAITS > &vm)
Teuchos::RCP< Thyra::ProductVectorBase< double > > dirichletCounter_
std::vector< PHX::MDField< const bool, Cell, NODE > > applyBC_
Teuchos::RCP< Thyra::ProductVectorBase< double > > r_
std::vector< PHX::MDField< const ScalarT, Cell, NODE > > scatterFields_
panzer::Traits::Tangent::ScalarT ScalarT
std::vector< Teuchos::RCP< const UniqueGlobalIndexer< LO, int > > > colIndexers_
PHX::MDField< ScalarT > vector
std::vector< PHX::MDField< const bool, Cell, NODE > > applyBC_
Pushes residual values into the residual vector for a Newton-based solve.
std::vector< int > subFieldIds_
std::string globalDataKey_
bool checkApplyBC_
If set to true, allows runtime disabling of dirichlet BCs on node-by-node basis.
std::string globalDataKey_
Wrapper to PHX::EvaluatorWithBaseImpl that implements Panzer-specific helpers.
ScatterDirichletResidual_BlockedEpetra()
std::vector< int > indexerIds_
std::size_t side_subcell_dim_
void evaluateFields(typename TRAITS::EvalData d)
std::string globalDataKey_
Teuchos::RCP< PHX::FieldTag > scatterHolder_
bool checkApplyBC_
If set to true, allows runtime disabling of dirichlet BCs on node-by-node basis.
ScatterDirichletResidual_BlockedEpetra(const Teuchos::ParameterList &p)
Teuchos::RCP< const std::map< std::string, std::string > > fieldMap_
std::vector< int > indexerIds_
ScatterDirichletResidual_BlockedEpetra(const std::vector< Teuchos::RCP< const UniqueGlobalIndexer< LO, int > > > &rIndexers, const std::vector< Teuchos::RCP< const UniqueGlobalIndexer< LO, int > > > &cIndexers)
Teuchos::RCP< Thyra::ProductVectorBase< double > > r_
virtual Teuchos::RCP< CloneableEvaluator > clone(const Teuchos::ParameterList &pl) const
panzer::Traits::Jacobian::ScalarT ScalarT
Teuchos::RCP< Thyra::ProductVectorBase< double > > dirichletCounter_
std::vector< int > indexerIds_
std::vector< PHX::MDField< const ScalarT, Cell, NODE > > scatterFields_
std::vector< Teuchos::RCP< const UniqueGlobalIndexer< LO, int > > > rowIndexers_
std::vector< PHX::MDField< const bool, Cell, NODE > > applyBC_
virtual Teuchos::RCP< CloneableEvaluator > clone(const Teuchos::ParameterList &pl) const
Teuchos::RCP< Thyra::ProductVectorBase< double > > r_
std::vector< PHX::MDField< const ScalarT, Cell, NODE > > scatterFields_
Teuchos::RCP< const std::map< std::string, std::string > > fieldMap_
Teuchos::RCP< Thyra::ProductVectorBase< double > > dirichletCounter_
std::size_t local_side_id_
virtual Teuchos::RCP< CloneableEvaluator > clone(const Teuchos::ParameterList &pl) const
ScatterDirichletResidual_BlockedEpetra()
Non-templated empty base class for template managers.
std::size_t local_side_id_
std::vector< Teuchos::RCP< const UniqueGlobalIndexer< LO, int > > > colIndexers_
virtual Teuchos::RCP< CloneableEvaluator > clone(const Teuchos::ParameterList &pl) const
Teuchos::RCP< PHX::FieldTag > scatterHolder_
panzer::Traits::Residual::ScalarT ScalarT