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;
79 public PHX::EvaluatorDerived<panzer::Traits::Residual, TRAITS>,
86 virtual Teuchos::RCP<CloneableEvaluator>
clone(
const Teuchos::ParameterList & pl)
const
92 { std::cout <<
"unspecialized version of \"ScatterDirichletResidual_BlockedEpetra::evaluateFields\" on \""+PHX::print<EvalT>()+
"\" should not be used!" << std::endl;
93 TEUCHOS_ASSERT(
false); }
106template<
typename TRAITS,
typename LO,
typename GO>
109 public PHX::EvaluatorDerived<panzer::Traits::Residual, TRAITS>,
115 const std::vector<Teuchos::RCP<const GlobalIndexer> > & )
116 : rowIndexers_(rIndexers) {}
119 const std::vector<Teuchos::RCP<const GlobalIndexer> > & cIndexers,
120 const Teuchos::ParameterList& p,
121 bool useDiscreteAdjoint=
false);
126 void preEvaluate(
typename TRAITS::PreEvalData d);
130 virtual Teuchos::RCP<CloneableEvaluator>
clone(
const Teuchos::ParameterList & pl)
const
153 Teuchos::RCP<const std::map<std::string,std::string> >
fieldMap_;
161 Teuchos::RCP<Thyra::ProductVectorBase<double> >
r_;
179template<
typename TRAITS,
typename LO,
typename GO>
182 public PHX::EvaluatorDerived<panzer::Traits::Tangent, TRAITS>,
188 const std::vector<Teuchos::RCP<const GlobalIndexer> > & )
189 : rowIndexers_(rIndexers) {}
192 const std::vector<Teuchos::RCP<const GlobalIndexer> > & cIndexers,
193 const Teuchos::ParameterList& p,
194 bool useDiscreteAdjoint=
false);
199 void preEvaluate(
typename TRAITS::PreEvalData d);
203 virtual Teuchos::RCP<CloneableEvaluator>
clone(
const Teuchos::ParameterList & pl)
const
226 Teuchos::RCP<const std::map<std::string,std::string> >
fieldMap_;
234 Teuchos::RCP<Thyra::ProductVectorBase<double> >
r_;
252template<
typename TRAITS,
typename LO,
typename GO>
255 public PHX::EvaluatorDerived<panzer::Traits::Jacobian, TRAITS>,
261 const std::vector<Teuchos::RCP<const GlobalIndexer> > & cIndexers)
262 : rowIndexers_(rIndexers), colIndexers_(cIndexers) {}
265 const std::vector<Teuchos::RCP<const GlobalIndexer> > & cIndexers,
266 const Teuchos::ParameterList& p,
267 bool useDiscreteAdjoint=
false);
269 void preEvaluate(
typename TRAITS::PreEvalData d);
276 virtual Teuchos::RCP<CloneableEvaluator>
clone(
const Teuchos::ParameterList & pl)
const
300 Teuchos::RCP<const std::map<std::string,std::string> >
fieldMap_;
311 Teuchos::RCP<Thyra::ProductVectorBase<double> >
r_;
312 Teuchos::RCP<Thyra::BlockedLinearOpBase<double> >
Jac_;
326#ifdef Panzer_BUILD_HESSIAN_SUPPORT
Non-templated empty base class for template managers.
Wrapper to PHX::EvaluatorWithBaseImpl that implements Panzer-specific helpers.
Teuchos::RCP< Thyra::BlockedLinearOpBase< double > > Jac_
bool checkApplyBC_
If set to true, allows runtime disabling of dirichlet BCs on node-by-node basis.
std::vector< Teuchos::RCP< const GlobalIndexer > > rowIndexers_
virtual Teuchos::RCP< CloneableEvaluator > clone(const Teuchos::ParameterList &pl) const
std::vector< BCFieldType > applyBC_
std::vector< int > subFieldIds_
std::string globalDataKey_
Teuchos::RCP< const std::map< std::string, std::string > > fieldMap_
Teuchos::RCP< PHX::FieldTag > scatterHolder_
std::size_t local_side_id_
Teuchos::RCP< Thyra::ProductVectorBase< double > > r_
ScatterDirichletResidual_BlockedEpetra()
std::size_t side_subcell_dim_
PHX::MDField< const bool, Cell, NODE > BCFieldType
Teuchos::RCP< Thyra::ProductVectorBase< double > > dirichletCounter_
panzer::Traits::Jacobian::ScalarT ScalarT
std::vector< PHX::MDField< const ScalarT, Cell, NODE > > scatterFields_
ScatterDirichletResidual_BlockedEpetra(const std::vector< Teuchos::RCP< const GlobalIndexer > > &rIndexers, const std::vector< Teuchos::RCP< const GlobalIndexer > > &cIndexers)
std::vector< int > indexerIds_
std::vector< Teuchos::RCP< const GlobalIndexer > > colIndexers_
std::vector< BCFieldType > applyBC_
std::vector< PHX::MDField< const ScalarT, Cell, NODE > > scatterFields_
panzer::Traits::Residual::ScalarT ScalarT
std::vector< int > indexerIds_
std::size_t side_subcell_dim_
Teuchos::RCP< Thyra::ProductVectorBase< double > > dirichletCounter_
std::vector< int > subFieldIds_
std::size_t local_side_id_
std::string globalDataKey_
Teuchos::RCP< const std::map< std::string, std::string > > fieldMap_
bool checkApplyBC_
If set to true, allows runtime disabling of dirichlet BCs on node-by-node basis.
ScatterDirichletResidual_BlockedEpetra(const std::vector< Teuchos::RCP< const GlobalIndexer > > &rIndexers, const std::vector< Teuchos::RCP< const GlobalIndexer > > &)
std::vector< Teuchos::RCP< const GlobalIndexer > > colIndexers_
virtual Teuchos::RCP< CloneableEvaluator > clone(const Teuchos::ParameterList &pl) const
Teuchos::RCP< Thyra::ProductVectorBase< double > > r_
PHX::MDField< const bool, Cell, NODE > BCFieldType
ScatterDirichletResidual_BlockedEpetra()
std::vector< Teuchos::RCP< const GlobalIndexer > > rowIndexers_
Teuchos::RCP< PHX::FieldTag > scatterHolder_
std::vector< int > subFieldIds_
PHX::MDField< const bool, Cell, NODE > BCFieldType
bool checkApplyBC_
If set to true, allows runtime disabling of dirichlet BCs on node-by-node basis.
std::vector< Teuchos::RCP< const GlobalIndexer > > colIndexers_
std::vector< Teuchos::RCP< const GlobalIndexer > > rowIndexers_
std::vector< PHX::MDField< const ScalarT, Cell, NODE > > scatterFields_
Teuchos::RCP< PHX::FieldTag > scatterHolder_
std::vector< BCFieldType > applyBC_
ScatterDirichletResidual_BlockedEpetra(const std::vector< Teuchos::RCP< const GlobalIndexer > > &rIndexers, const std::vector< Teuchos::RCP< const GlobalIndexer > > &)
Teuchos::RCP< Thyra::ProductVectorBase< double > > dirichletCounter_
Teuchos::RCP< Thyra::ProductVectorBase< double > > r_
std::string globalDataKey_
std::vector< int > indexerIds_
virtual Teuchos::RCP< CloneableEvaluator > clone(const Teuchos::ParameterList &pl) const
std::size_t side_subcell_dim_
Teuchos::RCP< const std::map< std::string, std::string > > fieldMap_
panzer::Traits::Tangent::ScalarT ScalarT
ScatterDirichletResidual_BlockedEpetra()
std::size_t local_side_id_
Pushes residual values into the residual vector for a Newton-based solve.
ScatterDirichletResidual_BlockedEpetra(const Teuchos::ParameterList &p)
virtual Teuchos::RCP< CloneableEvaluator > clone(const Teuchos::ParameterList &pl) const
void evaluateFields(typename TRAITS::EvalData d)
void postRegistrationSetup(typename TRAITS::SetupData d, PHX::FieldManager< TRAITS > &vm)