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;
66class BlockedDOFManager;
77template <
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); }
110template <
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_;
206template <
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
Non-templated empty base class for template managers.
Wrapper to PHX::EvaluatorWithBaseImpl that implements Panzer-specific helpers.
ScatterDirichletResidual_BlockedTpetra(const Teuchos::RCP< const BlockedDOFManager > &indexer)
std::vector< PHX::MDField< const ScalarT, Cell, NODE > > scatterFields_
std::vector< int > productVectorBlockIndex_
panzer::Traits::Jacobian::ScalarT ScalarT
std::vector< PHX::MDField< const bool, Cell, NODE > > applyBC_
std::string globalDataKey_
Tpetra::CrsMatrix< RealType, LO, GO, NodeT > CrsMatrixType
std::size_t side_subcell_dim_
virtual Teuchos::RCP< CloneableEvaluator > clone(const Teuchos::ParameterList &pl) const
PHX::View< LO * > blockOffsets_
The offset values of the blocked DOFs per element. Size of number of blocks in the product vector + 1...
TRAITS::RealType RealType
std::vector< PHX::View< int * > > basisIndexForMDFieldOffsets_
BlockedTpetraLinearObjContainer< RealType, LO, GO, NodeT > ContainerType
Tpetra::Operator< RealType, LO, GO, NodeT > OperatorType
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.
Teuchos::RCP< Thyra::ProductVectorBase< double > > dirichletCounter_
std::size_t local_side_id_
Teuchos::RCP< const BlockedTpetraLinearObjContainer< RealType, LO, GO, NodeT > > blockedContainer_
Teuchos::RCP< const panzer::BlockedDOFManager > globalIndexer_
Thyra::TpetraLinearOp< RealType, LO, GO, NodeT > ThyraLinearOp
std::vector< int > fieldIds_
Teuchos::RCP< const std::map< std::string, std::string > > fieldMap_
Teuchos::RCP< PHX::FieldTag > scatterHolder_
Tpetra::Map< LO, GO, NodeT > MapType
PHX::View< LO ** > worksetLIDs_
Local indices for unknowns.
ScatterDirichletResidual_BlockedTpetra()
Tpetra::Export< LO, GO, NodeT > ExportType
std::size_t side_subcell_dim_
std::vector< Teuchos::RCP< const panzer::GlobalIndexer > > fieldGlobalIndexers_
Tpetra::Vector< RealType, LO, GO, NodeT > VectorType
std::vector< int > productVectorBlockIndex_
Teuchos::RCP< Thyra::ProductVectorBase< double > > dirichletCounter_
std::vector< PHX::View< int * > > fieldOffsets_
Offset into the cell lids for each field.
Tpetra::CrsMatrix< RealType, LO, GO, NodeT > CrsMatrixType
BlockedTpetraLinearObjContainer< RealType, LO, GO, NodeT > ContainerType
Teuchos::RCP< const std::map< std::string, std::string > > fieldMap_
Tpetra::Map< LO, GO, NodeT > MapType
panzer::Traits::Residual::ScalarT ScalarT
bool checkApplyBC_
If set to true, allows runtime disabling of dirichlet BCs on node-by-node basis.
TRAITS::RealType RealType
Teuchos::RCP< const BlockedTpetraLinearObjContainer< RealType, LO, GO, NodeT > > blockedContainer_
std::vector< PHX::MDField< const bool, Cell, NODE > > applyBC_
ScatterDirichletResidual_BlockedTpetra(const Teuchos::RCP< const BlockedDOFManager > &indexer)
Tpetra::CrsGraph< LO, GO, NodeT > CrsGraphType
std::size_t local_side_id_
Teuchos::RCP< PHX::FieldTag > scatterHolder_
std::vector< PHX::MDField< const ScalarT, Cell, NODE > > scatterFields_
Tpetra::Import< LO, GO, NodeT > ImportType
std::vector< int > fieldIds_
Field IDs in the local product vector block (not global field id)
std::vector< PHX::View< int * > > basisIndexForMDFieldOffsets_
PHX::View< LO ** > worksetLIDs_
Local indices for unknowns.
ScatterDirichletResidual_BlockedTpetra()
std::string globalDataKey_
Teuchos::RCP< const panzer::BlockedDOFManager > globalIndexer_
virtual Teuchos::RCP< CloneableEvaluator > clone(const Teuchos::ParameterList &pl) const
Pushes residual values into the residual vector for a Newton-based solve.
void evaluateFields(typename TRAITS::EvalData)
virtual Teuchos::RCP< CloneableEvaluator > clone(const Teuchos::ParameterList &pl) const
ScatterDirichletResidual_BlockedTpetra(const Teuchos::RCP< const BlockedDOFManager > &)
void postRegistrationSetup(typename TRAITS::SetupData, PHX::FieldManager< TRAITS > &)