43#ifndef PANZER_EVALUATOR_SCATTER_RESIDUAL_BLOCKEDEPETRA_DECL_HPP
44#define PANZER_EVALUATOR_SCATTER_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>
class ProductVectorBase;
62 template <
typename>
class BlockedLinearOpBase;
75 public PHX::EvaluatorDerived<panzer::Traits::Residual, TRAITS>,
82 virtual Teuchos::RCP<CloneableEvaluator>
clone(
const Teuchos::ParameterList & pl)
const
88 { std::cout <<
"unspecialized version of \"ScatterResidual_BlockedEpetra::evaluateFields\" on "+PHX::print<EvalT>()+
" \" should not be used!" << std::endl;
89 TEUCHOS_ASSERT(
false); }
102template<
typename TRAITS,
typename LO,
typename GO>
105 public PHX::EvaluatorDerived<panzer::Traits::Residual, TRAITS>,
111 const std::vector<Teuchos::RCP<const GlobalIndexer> > & ,
113 : rowIndexers_(rIndexers) {}
116 const std::vector<Teuchos::RCP<const GlobalIndexer> > & cIndexers,
117 const Teuchos::ParameterList& p,
118 bool useDiscreteAdjoint=
false);
123 void preEvaluate(
typename TRAITS::PreEvalData d);
127 virtual Teuchos::RCP<CloneableEvaluator>
clone(
const Teuchos::ParameterList & pl)
const
149 Teuchos::RCP<const std::map<std::string,std::string> >
fieldMap_;
153 Teuchos::RCP<Thyra::ProductVectorBase<double> >
r_;
161template<
typename TRAITS,
typename LO,
typename GO>
164 public PHX::EvaluatorDerived<panzer::Traits::Tangent, TRAITS>,
170 const std::vector<Teuchos::RCP<const GlobalIndexer> > & ,
172 : rowIndexers_(rIndexers) {}
175 const std::vector<Teuchos::RCP<const GlobalIndexer> > & cIndexers,
176 const Teuchos::ParameterList& p,
177 bool useDiscreteAdjoint=
false);
182 void preEvaluate(
typename TRAITS::PreEvalData d);
186 virtual Teuchos::RCP<CloneableEvaluator>
clone(
const Teuchos::ParameterList & pl)
const
208 Teuchos::RCP<const std::map<std::string,std::string> >
fieldMap_;
212 Teuchos::RCP<Thyra::ProductVectorBase<double> >
r_;
221template<
typename TRAITS,
typename LO,
typename GO>
224 public PHX::EvaluatorDerived<panzer::Traits::Jacobian, TRAITS>,
250 const std::vector<Teuchos::RCP<const GlobalIndexer> > & cIndexers,
251 bool useDiscreteAdjoint=
false)
252 : rowIndexers_(rIndexers), colIndexers_(cIndexers), useDiscreteAdjoint_(useDiscreteAdjoint) {}
255 const std::vector<Teuchos::RCP<const GlobalIndexer> > & cIndexers,
256 const Teuchos::ParameterList& p,
257 bool useDiscreteAdjoint=
false);
262 void preEvaluate(
typename TRAITS::PreEvalData d);
266 virtual Teuchos::RCP<CloneableEvaluator>
clone(
const Teuchos::ParameterList & pl)
const
289 Teuchos::RCP<const std::map<std::string,std::string> >
fieldMap_;
294 Teuchos::RCP<Thyra::ProductVectorBase<double> >
r_;
295 Teuchos::RCP<Thyra::BlockedLinearOpBase<double> >
Jac_;
303#ifdef Panzer_BUILD_HESSIAN_SUPPORT
Non-templated empty base class for template managers.
Wrapper to PHX::EvaluatorWithBaseImpl that implements Panzer-specific helpers.
Teuchos::RCP< PHX::FieldTag > scatterHolder_
std::vector< Teuchos::RCP< const GlobalIndexer > > colIndexers_
Teuchos::RCP< Thyra::ProductVectorBase< double > > r_
Teuchos::RCP< const std::map< std::string, std::string > > fieldMap_
std::string globalDataKey_
ScatterResidual_BlockedEpetra()
std::vector< int > indexerIds_
ScatterResidual_BlockedEpetra(const std::vector< Teuchos::RCP< const GlobalIndexer > > &rIndexers, const std::vector< Teuchos::RCP< const GlobalIndexer > > &cIndexers, bool useDiscreteAdjoint=false)
Teuchos::RCP< Thyra::BlockedLinearOpBase< double > > Jac_
virtual Teuchos::RCP< CloneableEvaluator > clone(const Teuchos::ParameterList &pl) const
std::vector< Teuchos::RCP< const GlobalIndexer > > rowIndexers_
panzer::Traits::Jacobian::ScalarT ScalarT
std::vector< PHX::MDField< const ScalarT, Cell, NODE > > scatterFields_
std::vector< int > subFieldIds_
Teuchos::RCP< const std::map< std::string, std::string > > fieldMap_
std::vector< int > indexerIds_
Teuchos::RCP< Thyra::ProductVectorBase< double > > r_
std::vector< Teuchos::RCP< const GlobalIndexer > > colIndexers_
ScatterResidual_BlockedEpetra(const std::vector< Teuchos::RCP< const GlobalIndexer > > &rIndexers, const std::vector< Teuchos::RCP< const GlobalIndexer > > &, bool)
std::vector< PHX::MDField< const ScalarT, Cell, NODE > > scatterFields_
Teuchos::RCP< PHX::FieldTag > scatterHolder_
std::vector< Teuchos::RCP< const GlobalIndexer > > rowIndexers_
std::vector< int > subFieldIds_
panzer::Traits::Residual::ScalarT ScalarT
std::string globalDataKey_
virtual Teuchos::RCP< CloneableEvaluator > clone(const Teuchos::ParameterList &pl) const
ScatterResidual_BlockedEpetra()
Teuchos::RCP< PHX::FieldTag > scatterHolder_
Teuchos::RCP< const std::map< std::string, std::string > > fieldMap_
std::vector< Teuchos::RCP< const GlobalIndexer > > rowIndexers_
std::vector< PHX::MDField< const ScalarT, Cell, NODE > > scatterFields_
Teuchos::RCP< Thyra::ProductVectorBase< double > > r_
ScatterResidual_BlockedEpetra()
std::vector< Teuchos::RCP< const GlobalIndexer > > colIndexers_
ScatterResidual_BlockedEpetra(const std::vector< Teuchos::RCP< const GlobalIndexer > > &rIndexers, const std::vector< Teuchos::RCP< const GlobalIndexer > > &, bool)
std::vector< int > subFieldIds_
std::string globalDataKey_
std::vector< int > indexerIds_
panzer::Traits::Tangent::ScalarT ScalarT
virtual Teuchos::RCP< CloneableEvaluator > clone(const Teuchos::ParameterList &pl) const
Pushes residual values into the residual vector for a Newton-based solve.
virtual Teuchos::RCP< CloneableEvaluator > clone(const Teuchos::ParameterList &pl) const
ScatterResidual_BlockedEpetra(const Teuchos::ParameterList &p)
void postRegistrationSetup(typename TRAITS::SetupData d, PHX::FieldManager< TRAITS > &vm)
void evaluateFields(typename TRAITS::EvalData d)