43 #ifndef PANZER_PHYSICS_BLOCK_HPP 44 #define PANZER_PHYSICS_BLOCK_HPP 50 #include "Teuchos_RCP.hpp" 51 #include "Phalanx_FieldManager.hpp" 71 class IntegrationRule;
72 struct EquationSetFactory;
83 void buildPhysicsBlocks(
const std::map<std::string,std::string>& block_ids_to_physics_ids,
84 const std::map<std::string,Teuchos::RCP<const shards::CellTopology> >& block_ids_to_cell_topo,
85 const Teuchos::RCP<Teuchos::ParameterList>& physics_blocks_plist,
86 const int default_integration_order,
87 const std::size_t workset_size,
88 const Teuchos::RCP<const panzer::EquationSetFactory>& eqset_factory,
89 const Teuchos::RCP<panzer::GlobalData>& global_data,
90 const bool build_transient_support,
91 std::vector<Teuchos::RCP<panzer::PhysicsBlock> > & physicsBlocks,
92 const std::vector<std::string>& tangent_param_names = std::vector<std::string>());
100 void readPhysicsBlocks(
const std::map<std::string,std::string>& block_ids_to_physics_ids,
101 const Teuchos::RCP<Teuchos::ParameterList>& physics_blocks_plist,
102 std::vector<Teuchos::RCP<panzer::PhysicsBlock> > & physicsBlocks);
111 Teuchos::RCP<panzer::PhysicsBlock> findPhysicsBlock(
const std::string element_block_id,
112 const std::vector<Teuchos::RCP<panzer::PhysicsBlock> > & physics_blocks,
113 bool throw_on_failure =
true);
123 { std::cout <<
"WARNING: Default constructor for panzer::PhysicsBlock is for testing purposes only!" << std::endl; }
127 PhysicsBlock(
const Teuchos::RCP<Teuchos::ParameterList>& physics_block_plist,
128 const std::string & element_block_id,
129 const int default_integration_order,
131 const Teuchos::RCP<const panzer::EquationSetFactory>& factory,
132 const Teuchos::RCP<panzer::GlobalData>& global_data,
133 const bool build_transient_support,
134 const std::vector<std::string>& tangent_param_names = std::vector<std::string>());
141 PhysicsBlock(
const Teuchos::RCP<Teuchos::ParameterList>& physics_block_plist,
142 const std::string & element_block_id);
152 const std::string & physics_block_id,
153 const int integration_order,
155 const Teuchos::RCP<panzer::GlobalData>& global_data,
156 const Teuchos::RCP<panzer::PureBasis> & fields);
161 void initialize(
const int default_integration_order,
162 const bool build_transient_support,
164 const Teuchos::RCP<const panzer::EquationSetFactory>& factory,
165 const Teuchos::RCP<panzer::GlobalData>& global_data,
166 const std::vector<std::string>& tangent_param_names = std::vector<std::string>());
175 const Teuchos::ParameterList& user_data)
const;
179 const Teuchos::ParameterList& user_data)
const;
183 const Teuchos::ParameterList& user_data)
const;
187 const Teuchos::ParameterList& user_data)
const;
191 const Teuchos::ParameterList& models,
192 const Teuchos::ParameterList& user_data)
const;
196 const std::string& model_name,
197 const Teuchos::ParameterList& models,
199 const Teuchos::ParameterList& user_data)
const;
203 const std::string& model_name,
204 const Teuchos::ParameterList& models,
205 const Teuchos::ParameterList& user_data)
const;
207 template<
typename EvalT>
209 const Teuchos::ParameterList& user_data)
const;
211 template<
typename EvalT>
214 const Teuchos::ParameterList& user_data)
const;
216 template<
typename EvalT>
219 const Teuchos::ParameterList& user_data)
const;
221 template<
typename EvalT>
224 const Teuchos::ParameterList& user_data)
const;
226 template<
typename EvalT>
229 const Teuchos::ParameterList& models,
230 const Teuchos::ParameterList& user_data)
const;
232 template<
typename EvalT>
235 const std::string& model_name,
236 const Teuchos::ParameterList& models,
238 const Teuchos::ParameterList& user_data)
const;
240 const std::vector<std::string>&
getDOFNames()
const;
253 const std::map<std::string,Teuchos::RCP<panzer::PureBasis> >&
getBases()
const;
271 Teuchos::RCP<panzer::GlobalData>
globalData()
const;
284 void initialize(
const Teuchos::RCP<Teuchos::ParameterList>& input_parameters,
285 const int& default_integration_order,
286 const std::string & element_block_id,
288 const bool build_transient_support,
289 const std::vector<std::string>& tangent_param_names = std::vector<std::string>());
306 std::map<std::string,Teuchos::RCP<panzer::PureBasis> >
m_bases;
310 std::vector< Teuchos::RCP<panzer::EquationSet_TemplateManager<panzer::Traits> > >
m_equation_sets;
324 template<
typename EvalT>
326 const Teuchos::ParameterList& user_data)
const 333 vector< RCP<EquationSet_TemplateManager<panzer::Traits> > >::const_iterator
347 template<
typename EvalT>
350 const Teuchos::ParameterList& user_data)
const 357 vector< RCP<EquationSet_TemplateManager<panzer::Traits> > >::const_iterator
358 eq_set = m_equation_sets.begin();
360 for (;eq_set != m_equation_sets.end(); ++eq_set,++idx) {
361 if (m_active_evaluation_types[idx]) {
365 const int di = eqstm.getAsObject<EvalT>()->setDetailsIndex(this->getDetailsIndex());
366 eqstm.getAsObject<EvalT>()->buildAndRegisterGatherAndOrientationEvaluators(fm,*m_field_lib,lof,user_data);
367 eqstm.getAsObject<EvalT>()->setDetailsIndex(di);
372 template<
typename EvalT>
375 const Teuchos::ParameterList& user_data)
const 382 vector< RCP<EquationSet_TemplateManager<panzer::Traits> > >::const_iterator
383 eq_set = m_equation_sets.begin();
384 for (;eq_set != m_equation_sets.end(); ++eq_set) {
388 for (std::map<
int,Teuchos::RCP<panzer::IntegrationRule> >::const_iterator ir_iter = m_integration_rules.begin();
389 ir_iter != m_integration_rules.end(); ++ ir_iter) {
391 Teuchos::RCP<panzer::IntegrationRule> ir = ir_iter->second;
393 const int di = eqstm.getAsObject<EvalT>()->setDetailsIndex(this->getDetailsIndex());
394 eqstm.getAsObject<EvalT>()->buildAndRegisterDOFProjectionsToIPEvaluators(fm,*m_field_lib->buildFieldLayoutLibrary(*ir),ir,lof,user_data);
395 eqstm.getAsObject<EvalT>()->setDetailsIndex(di);
401 template<
typename EvalT>
404 const Teuchos::ParameterList& user_data)
const 411 vector< RCP<EquationSet_TemplateManager<panzer::Traits> > >::const_iterator
412 eq_set = m_equation_sets.begin();
413 for (;eq_set != m_equation_sets.end(); ++eq_set) {
417 const int di = eqstm.getAsObject<EvalT>()->setDetailsIndex(this->getDetailsIndex());
418 eqstm.getAsObject<EvalT>()->buildAndRegisterScatterEvaluators(fm,*m_field_lib,lof,user_data);
419 eqstm.getAsObject<EvalT>()->setDetailsIndex(di);
423 template<
typename EvalT>
426 const Teuchos::ParameterList& models,
427 const Teuchos::ParameterList& user_data)
const 434 vector< RCP<EquationSet_TemplateManager<panzer::Traits> > >::const_iterator
435 eq_set = m_equation_sets.begin();
436 for (;eq_set != m_equation_sets.end(); ++eq_set) {
441 for (std::map<
int,Teuchos::RCP<panzer::IntegrationRule> >::const_iterator ir_iter = m_integration_rules.begin();
442 ir_iter != m_integration_rules.end(); ++ ir_iter) {
444 Teuchos::RCP<panzer::IntegrationRule> ir = ir_iter->second;
446 const int di = eqstm.getAsObject<EvalT>()->setDetailsIndex(this->getDetailsIndex());
447 eqstm.getAsObject<EvalT>()->buildAndRegisterClosureModelEvaluators(fm,*m_field_lib->buildFieldLayoutLibrary(*ir),ir,factory,models,user_data);
448 eqstm.getAsObject<EvalT>()->setDetailsIndex(di);
454 template<
typename EvalT>
457 const std::string& model_name,
458 const Teuchos::ParameterList& models,
460 const Teuchos::ParameterList& user_data)
const 467 vector< RCP<EquationSet_TemplateManager<panzer::Traits> > >::const_iterator
468 eq_set = m_equation_sets.begin();
469 for (;eq_set != m_equation_sets.end(); ++eq_set) {
470 std::vector<StrBasisPair> providedDOFs;
474 const int di = eqstm.getAsObject<EvalT>()->setDetailsIndex(this->getDetailsIndex());
475 eqstm.getAsObject<EvalT>()->buildAndRegisterInitialConditionEvaluators(fm, *m_field_lib, factory, model_name, models, lof, user_data);
476 eqstm.getAsObject<EvalT>()->setDetailsIndex(di);
std::vector< StrPureBasisPair > m_provided_dofs
Teuchos::RCP< panzer::GlobalData > m_global_data
void buildAndRegisterDOFProjectionsToIPEvaluatorsForType(PHX::FieldManager< panzer::Traits > &fm, const Teuchos::Ptr< const panzer::LinearObjFactory< panzer::Traits > > &lof, const Teuchos::ParameterList &user_data) const
Object that contains information on the physics and discretization of a block of elements with the SA...
Teuchos::RCP< panzer::GlobalData > globalData() const
int getDetailsIndex() const
Get the WorksetDetails index.
Teuchos::RCP< PhysicsBlock > copyWithCellData(const panzer::CellData &cell_data) const
int setDetailsIndex(const int details_index)
Teuchos::RCP< const panzer::EquationSetFactory > m_eqset_factory
void buildAndRegisterClosureModelEvaluatorsForType(PHX::FieldManager< panzer::Traits > &fm, const panzer::ClosureModelFactory_TemplateManager< panzer::Traits > &factory, const Teuchos::ParameterList &models, const Teuchos::ParameterList &user_data) const
std::string elementBlockID() const
void buildAndRegisterClosureModelEvaluators(PHX::FieldManager< panzer::Traits > &fm, const panzer::ClosureModelFactory_TemplateManager< panzer::Traits > &factory, const Teuchos::ParameterList &models, const Teuchos::ParameterList &user_data) const
void buildAndRegisterGatherAndOrientationEvaluatorsForType(PHX::FieldManager< panzer::Traits > &fm, const LinearObjFactory< panzer::Traits > &lof, const Teuchos::ParameterList &user_data) const
void buildAndRegisterInitialConditionEvaluatorsForType(PHX::FieldManager< panzer::Traits > &fm, const panzer::ClosureModelFactory_TemplateManager< panzer::Traits > &factory, const std::string &model_name, const Teuchos::ParameterList &models, const panzer::LinearObjFactory< panzer::Traits > &lof, const Teuchos::ParameterList &user_data) const
const std::vector< std::string > & getDOFNames() const
int m_default_integration_order
void buildAndRegisterScatterEvaluatorsForType(PHX::FieldManager< panzer::Traits > &fm, const LinearObjFactory< panzer::Traits > &lof, const Teuchos::ParameterList &user_data) const
void setActiveEvaluationTypes(const std::vector< bool > &aet)
Used to save memory by disabling unneeded evaluation types.
std::vector< Teuchos::RCP< panzer::EquationSet_TemplateManager< panzer::Traits > > > m_equation_sets
void buildAndRegisterInitialConditionEvaluators(PHX::FieldManager< panzer::Traits > &fm, const panzer::ClosureModelFactory_TemplateManager< panzer::Traits > &factory, const std::string &model_name, const Teuchos::ParameterList &models, const panzer::LinearObjFactory< panzer::Traits > &lof, const Teuchos::ParameterList &user_data) const
std::vector< std::vector< std::string > > m_coordinate_dofs
Data for determining cell topology and dimensionality.
panzer::CellData m_cell_data
std::string physicsBlockID() const
void buildAndRegisterGatherAndOrientationEvaluators(PHX::FieldManager< panzer::Traits > &fm, const panzer::LinearObjFactory< panzer::Traits > &lof, const Teuchos::ParameterList &user_data) const
std::map< std::string, Teuchos::RCP< panzer::PureBasis > > m_bases
map of unique bases, key is the panzer::PureBasis::name() corresponding to its value ...
Teuchos::RCP< const FieldLibraryBase > getFieldLibraryBase() const
Teuchos::RCP< Teuchos::ParameterList > m_input_parameters
store the input parameter list for copy ctors
std::map< int, Teuchos::RCP< panzer::IntegrationRule > > m_integration_rules
map of unique integration rules, key is panzer::IntegrationRule::order() corresponding to its value ...
void buildAndRegisterDOFProjectionsToIPEvaluators(PHX::FieldManager< panzer::Traits > &fm, const Teuchos::Ptr< const panzer::LinearObjFactory< panzer::Traits > > &lof, const Teuchos::ParameterList &user_data) const
const std::vector< std::vector< std::string > > & getCoordinateDOFs() const
const panzer::CellData & cellData() const
Teuchos::RCP< const FieldLibrary > getFieldLibrary() const
void buildAndRegisterEquationSetEvaluators(PHX::FieldManager< panzer::Traits > &fm, const Teuchos::ParameterList &user_data) const
void buildAndRegisterEquationSetEvaluatorsForType(PHX::FieldManager< panzer::Traits > &fm, const Teuchos::ParameterList &user_data) const
const std::vector< StrPureBasisPair > & getProvidedDOFs() const
void activateAllEvaluationTypes()
Used to reactivate all evaluation types if some were temporarily disabled with a call to setActiveEva...
const std::map< std::string, Teuchos::RCP< panzer::PureBasis > > & getBases() const
Returns the unique set of bases, key is the unique panzer::PureBasis::name() of the basis...
WorksetNeeds getWorksetNeeds() const
void buildAndRegisterScatterEvaluators(PHX::FieldManager< panzer::Traits > &fm, const panzer::LinearObjFactory< panzer::Traits > &lof, const Teuchos::ParameterList &user_data) const
void initialize(const int default_integration_order, const bool build_transient_support, const panzer::CellData &cell_data, const Teuchos::RCP< const panzer::EquationSetFactory > &factory, const Teuchos::RCP< panzer::GlobalData > &global_data, const std::vector< std::string > &tangent_param_names=std::vector< std::string >())
bool m_build_transient_support
Teuchos::RCP< FieldLibrary > m_field_lib
std::vector< bool > m_active_evaluation_types
Returns true for evaluation types that are active.
std::string m_element_block_id
const std::vector< StrPureBasisPair > & getTangentFields() const
Returns list of tangent fields from DOFs and tangent param names.
std::vector< std::string > m_dof_names
const std::map< int, Teuchos::RCP< panzer::IntegrationRule > > & getIntegrationRules() const
Returns the unique set of point rules, key is the unique panzer::PointRule::name() ...
const shards::CellTopology getBaseCellTopology() const
std::vector< StrPureBasisPair > m_tangent_fields
Teuchos::RCP< const Teuchos::ParameterList > getParameterList() const