43 #ifndef PANZER_EVALUATOR_SCATTER_DIRICHLET_RESIDUAL_BLOCKEDTPETRA_HPP 44 #define PANZER_EVALUATOR_SCATTER_DIRICHLET_RESIDUAL_BLOCKEDTPETRA_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;
66 class BlockedDOFManager;
77 template <
typename EvalT,
typename TRAITS,
typename LO,
typename GO,
typename NodeT=panzer::TpetraNodeType>
80 public PHX::EvaluatorDerived<panzer::Traits::Residual, TRAITS>,
88 const Teuchos::ParameterList& p);
90 virtual Teuchos::RCP<CloneableEvaluator>
clone(
const Teuchos::ParameterList & pl)
const 96 { std::cout <<
"unspecialized version of \"ScatterDirichletResidual_BlockedTpetra::evaluateFields\" on \""+PHX::print<EvalT>()+
"\" should not be used!" << std::endl;
97 TEUCHOS_ASSERT(
false); }
110 template <
typename TRAITS,
typename LO,
typename GO,
typename NodeT>
113 public PHX::EvaluatorDerived<panzer::Traits::Residual, TRAITS>,
118 : globalIndexer_(indexer) {}
121 const Teuchos::ParameterList& p);
126 void preEvaluate(
typename TRAITS::PreEvalData d);
130 virtual Teuchos::RCP<CloneableEvaluator>
clone(
const Teuchos::ParameterList & pl)
const 171 Teuchos::RCP<const std::map<std::string,std::string> >
fieldMap_;
198 std::vector< PHX::MDField<const bool,Cell,NODE> >
applyBC_;
206 template <
typename TRAITS,
typename LO,
typename GO,
typename NodeT>
209 public PHX::EvaluatorDerived<panzer::Traits::Jacobian, TRAITS>,
214 : globalIndexer_(indexer) {}
217 const Teuchos::ParameterList& p);
219 void preEvaluate(
typename TRAITS::PreEvalData d);
226 virtual Teuchos::RCP<CloneableEvaluator>
clone(
const Teuchos::ParameterList & pl)
const 260 Teuchos::RCP<const std::map<std::string,std::string> >
fieldMap_;
287 std::vector< PHX::MDField<const bool,Cell,NODE> >
applyBC_;
290 static constexpr
int maxDerivativeArraySize_ = 256;
297 #ifdef Panzer_BUILD_HESSIAN_SUPPORT
PHX::View< LO ** > worksetLIDs_
Local indices for unknowns.
virtual Teuchos::RCP< CloneableEvaluator > clone(const Teuchos::ParameterList &pl) const
TRAITS::RealType RealType
TRAITS::RealType RealType
std::size_t side_subcell_dim_
Tpetra::Export< LO, GO, NodeT > ExportType
Teuchos::RCP< const std::map< std::string, std::string > > fieldMap_
std::vector< PHX::MDField< const ScalarT, Cell, NODE > > scatterFields_
Teuchos::RCP< const panzer::BlockedDOFManager > globalIndexer_
BlockedTpetraLinearObjContainer< RealType, LO, GO, NodeT > ContainerType
virtual Teuchos::RCP< CloneableEvaluator > clone(const Teuchos::ParameterList &pl) const
ScatterDirichletResidual_BlockedTpetra(const Teuchos::RCP< const BlockedDOFManager > &)
std::vector< int > productVectorBlockIndex_
Teuchos::RCP< Thyra::ProductVectorBase< double > > dirichletCounter_
std::size_t side_subcell_dim_
BlockedTpetraLinearObjContainer< RealType, LO, GO, NodeT > ContainerType
Teuchos::RCP< Thyra::ProductVectorBase< double > > dirichletCounter_
std::string globalDataKey_
std::size_t local_side_id_
std::vector< PHX::MDField< const ScalarT, Cell, NODE > > scatterFields_
bool checkApplyBC_
If set to true, allows runtime disabling of dirichlet BCs on node-by-node basis.
Tpetra::Operator< RealType, LO, GO, NodeT > OperatorType
Tpetra::CrsMatrix< RealType, LO, GO, NodeT > CrsMatrixType
Wrapper to PHX::EvaluatorWithBaseImpl that implements Panzer-specific helpers.
PHX::View< LO ** > worksetLIDs_
Local indices for unknowns.
Thyra::TpetraLinearOp< RealType, LO, GO, NodeT > ThyraLinearOp
virtual Teuchos::RCP< CloneableEvaluator > clone(const Teuchos::ParameterList &pl) const
Teuchos::RCP< const BlockedTpetraLinearObjContainer< RealType, LO, GO, NodeT > > blockedContainer_
Tpetra::Import< LO, GO, NodeT > ImportType
std::vector< PHX::View< int * > > basisIndexForMDFieldOffsets_
ScatterDirichletResidual_BlockedTpetra(const Teuchos::RCP< const BlockedDOFManager > &indexer)
Teuchos::RCP< const std::map< std::string, std::string > > fieldMap_
std::vector< PHX::MDField< const bool, Cell, NODE > > applyBC_
Tpetra::Vector< RealType, LO, GO, NodeT > VectorType
Tpetra::CrsMatrix< RealType, LO, GO, NodeT > CrsMatrixType
Tpetra::Map< LO, GO, NodeT > MapType
void evaluateFields(typename TRAITS::EvalData)
Tpetra::Map< LO, GO, NodeT > MapType
std::vector< int > fieldIds_
Field IDs in the local product vector block (not global field id)
Teuchos::RCP< const panzer::BlockedDOFManager > globalIndexer_
PHX::View< LO * > blockOffsets_
The offset values of the blocked DOFs per element. Size of number of blocks in the product vector + 1...
Teuchos::RCP< PHX::FieldTag > scatterHolder_
bool checkApplyBC_
If set to true, allows runtime disabling of dirichlet BCs on node-by-node basis.
std::vector< PHX::View< int * > > fieldOffsets_
Offset into the cell lids for each field. Size of number of fields to scatter.
std::vector< PHX::View< int * > > basisIndexForMDFieldOffsets_
std::vector< int > productVectorBlockIndex_
std::string globalDataKey_
Tpetra::CrsGraph< LO, GO, NodeT > CrsGraphType
std::vector< PHX::MDField< const bool, Cell, NODE > > applyBC_
Pushes residual values into the residual vector for a Newton-based solve.
panzer::Traits::Residual::ScalarT ScalarT
std::size_t local_side_id_
Teuchos::RCP< PHX::FieldTag > scatterHolder_
ScatterDirichletResidual_BlockedTpetra()
std::vector< int > fieldIds_
Teuchos::RCP< const BlockedTpetraLinearObjContainer< RealType, LO, GO, NodeT > > blockedContainer_
ScatterDirichletResidual_BlockedTpetra(const Teuchos::RCP< const BlockedDOFManager > &indexer)
panzer::Traits::Jacobian::ScalarT ScalarT
void postRegistrationSetup(typename TRAITS::SetupData, PHX::FieldManager< TRAITS > &)
std::vector< Teuchos::RCP< const panzer::GlobalIndexer > > fieldGlobalIndexers_
Non-templated empty base class for template managers.
std::vector< PHX::View< int * > > fieldOffsets_
Offset into the cell lids for each field.