43#ifndef PANZER_ASSEMBLY_ENGINE_IMPL_HPP
44#define PANZER_ASSEMBLY_ENGINE_IMPL_HPP
46#include "Phalanx_FieldManager.hpp"
54template <
typename EvalT>
56AssemblyEngine(
const Teuchos::RCP<panzer::FieldManagerBuilder>& fmb,
58 : m_field_manager_builder(fmb), m_lin_obj_factory(lof), countersInitialized_(false)
65template <
typename EvalT>
76 if ( flags.
getValue() & EvaluationFlags::Initialize ) {
77 PANZER_FUNC_TIME_MONITOR_DIFF(
"panzer::AssemblyEngine::evaluate_gather("+PHX::print<EvalT>()+
")", eval_gather);
91 if ( flags.
getValue() & EvaluationFlags::VolumetricFill) {
92 PANZER_FUNC_TIME_MONITOR_DIFF(
"panzer::AssemblyEngine::evaluate_volume("+PHX::print<EvalT>()+
")", eval_vol);
93 this->evaluateVolume(in);
103 if ( flags.
getValue() & EvaluationFlags::BoundaryFill) {
105 PANZER_FUNC_TIME_MONITOR_DIFF(
"panzer::AssemblyEngine::evaluate_neumannbcs("+PHX::print<EvalT>()+
")",eval_neumannbcs);
106 this->evaluateNeumannBCs(in);
110 PANZER_FUNC_TIME_MONITOR_DIFF(
"panzer::AssemblyEngine::evaluate_interfacebcs("+PHX::print<EvalT>()+
")",eval_interfacebcs);
111 this->evaluateInterfaceBCs(in);
116 PANZER_FUNC_TIME_MONITOR_DIFF(
"panzer::AssemblyEngine::evaluate_dirichletbcs("+PHX::print<EvalT>()+
")",eval_dirichletbcs);
117 this->evaluateDirichletBCs(in);
121 if ( flags.
getValue() & EvaluationFlags::Scatter) {
122 PANZER_FUNC_TIME_MONITOR_DIFF(
"panzer::AssemblyEngine::evaluate_scatter("+PHX::print<EvalT>()+
")",eval_scatter);
137template <
typename EvalT>
156 Teuchos::RCP<LOC> counter = this->evaluateDirichletBCs(in);
171template <
typename EvalT>
175 const std::vector< Teuchos::RCP< PHX::FieldManager<panzer::Traits> > > &
176 volume_field_managers = m_field_manager_builder->getVolumeFieldManagers();
177 const std::vector<WorksetDescriptor> & wkstDesc = m_field_manager_builder->getVolumeWorksetDescriptors();
179 Teuchos::RCP<panzer::WorksetContainer> wkstContainer = m_field_manager_builder->getWorksetContainer();
189 for (std::size_t block = 0; block < volume_field_managers.size(); ++block) {
191 Teuchos::RCP< PHX::FieldManager<panzer::Traits> > fm = volume_field_managers[block];
192 std::vector<panzer::Workset>& w = *wkstContainer->getWorksets(wd);
194 fm->template preEvaluate<EvalT>(ped);
197 for (std::size_t i = 0; i < w.size(); ++i) {
209 fm->template evaluateFields<EvalT>(workset);
217 fm->template postEvaluate<EvalT>(NULL);
223template <
typename EvalT>
232template <
typename EvalT>
241template <
typename EvalT>
247 if(!countersInitialized_) {
248 localCounter_ = m_lin_obj_factory->buildPrimitiveGhostedLinearObjContainer();
249 globalCounter_ = m_lin_obj_factory->buildPrimitiveLinearObjContainer();
250 summedGhostedCounter_ = m_lin_obj_factory->buildPrimitiveGhostedLinearObjContainer();
251 countersInitialized_ =
true;
260 summedGhostedCounter_->initialize();
261 globalCounter_->initialize();
269 m_lin_obj_factory->ghostToGlobalContainer(*localCounter_,*globalCounter_,LOC::F);
273 m_lin_obj_factory->globalToGhostContainer(*globalCounter_,*summedGhostedCounter_,LOC::F);
286 if(itr->second->requiresDirichletAdjustment()) {
287 Teuchos::RCP<LinearObjContainer> loc = Teuchos::rcp_dynamic_cast<LinearObjContainer>(itr->second);
288 if(loc!=Teuchos::null) {
289 m_lin_obj_factory->adjustForDirichletConditions(*localCounter_,*summedGhostedCounter_,*loc);
293 Teuchos::RCP<GlobalEvaluationData_BCAdjustment> bc_adjust = Teuchos::rcp_dynamic_cast<GlobalEvaluationData_BCAdjustment>(itr->second,
true);
294 bc_adjust->adjustForDirichletConditions(*localCounter_,*summedGhostedCounter_);
299 return globalCounter_;
304template <
typename EvalT>
308 const Teuchos::RCP<LinearObjContainer> preEval_loc)
310 Teuchos::RCP<panzer::WorksetContainer> wkstContainer = m_field_manager_builder->getWorksetContainer();
313 ped.
gedc->addDataObject(
"Dirichlet Counter",preEval_loc);
324 double betaValue = in.
beta;
331 std::map<unsigned,PHX::FieldManager<panzer::Traits> >,
333 m_field_manager_builder->getBCFieldManagers();
338 std::map<unsigned,PHX::FieldManager<panzer::Traits> >,
342 for (bcfm_it_type bcfm_it = bc_field_managers.begin();
343 bcfm_it != bc_field_managers.end(); ++bcfm_it) {
346 const std::map<unsigned,PHX::FieldManager<panzer::Traits> > bc_fm =
350 Teuchos::RCP<const std::map<unsigned,panzer::Workset> > bc_wkst_ptr = wkstContainer->getSideWorksets(desc);
351 TEUCHOS_TEST_FOR_EXCEPTION(bc_wkst_ptr == Teuchos::null, std::logic_error,
352 "Failed to find corresponding bc workset!");
353 const std::map<unsigned,panzer::Workset>& bc_wkst = *bc_wkst_ptr;
356 if (bc.
bcType() == bc_type) {
357 std::ostringstream timerName;
358 timerName <<
"panzer::AssemblyEngine::evaluateBCs: " << bc.
identifier();
359#ifdef PANZER_TEUCHOS_TIME_MONITOR
360 auto timer1 = Teuchos::TimeMonitor::getNewTimer(timerName.str());
361 Teuchos::TimeMonitor tm1(*timer1);
366 std::ostringstream timerSideName;
367 timerSideName <<
"panzer::AssemblyEngine::evaluateBCs: " << bc.
identifier() <<
", side=" << side->first;
368#ifdef PANZER_TEUCHOS_TIME_MONITOR
369 auto timer2 = Teuchos::TimeMonitor::getNewTimer(timerSideName.str());
370 Teuchos::TimeMonitor tm2(*timer2);
374 unsigned local_side_index = side->first;
379 std::map<unsigned,panzer::Workset>::const_iterator wkst_it =
380 bc_wkst.find(local_side_index);
382 TEUCHOS_TEST_FOR_EXCEPTION(wkst_it == bc_wkst.end(), std::logic_error,
383 "Failed to find corresponding bc workset side!");
389 local_side_fm.template preEvaluate<EvalT>(ped);
393 workset.
beta = betaValue;
398 local_side_fm.template evaluateFields<EvalT>(workset);
401 local_side_fm.template postEvaluate<EvalT>(NULL);
std::vector< double > gather_seeds
Teuchos::RCP< panzer::LinearObjContainer > ghostedContainer_
bool apply_dirichlet_beta
std::string second_sensitivities_name
bool evaluate_transient_terms
Teuchos::RCP< panzer::LinearObjContainer > container_
std::string first_sensitivities_name
void fillGlobalEvaluationDataContainer(GlobalEvaluationDataContainer &gedc) const
Using internal map fill the global evaluation data container object.
Teuchos::RCP< LinearObjContainer > evaluateOnlyDirichletBCs(const panzer::AssemblyEngineInArgs &input_arguments)
void evaluate(const panzer::AssemblyEngineInArgs &input_arguments, const EvaluationFlags flags=EvaluationFlags(EvaluationFlags::All))
void evaluateVolume(const panzer::AssemblyEngineInArgs &input_arguments)
void evaluateNeumannBCs(const panzer::AssemblyEngineInArgs &input_arguments)
void evaluateInterfaceBCs(const panzer::AssemblyEngineInArgs &input_arguments)
AssemblyEngine(const Teuchos::RCP< panzer::FieldManagerBuilder > &fmb, const Teuchos::RCP< const panzer::LinearObjFactory< panzer::Traits > > &lof)
Teuchos::RCP< LinearObjContainer > evaluateDirichletBCs(const panzer::AssemblyEngineInArgs &input_arguments)
This method returns the global counter used to indicate which rows are boundary conditions.
void evaluateBCs(const panzer::BCType bc_type, const panzer::AssemblyEngineInArgs &input_arguments, const Teuchos::RCP< LinearObjContainer > preEval_loc=Teuchos::null)
Stores input information for a boundary condition.
BCType bcType() const
Returns the boundary condition type (Dirichlet or Neumann or Interface).
std::string identifier() const
A unique string identifier for this boundary condition.
void addDataObject(const std::string &key, const Teuchos::RCP< GlobalEvaluationData > &ged)
void initialize()
Call initialize on all containers.
void globalToGhost(int p)
Call global to ghost on all the containers.
void ghostToGlobal(int p)
Call ghost to global on all the containers.
const_iterator begin() const
std::unordered_map< std::string, Teuchos::RCP< GlobalEvaluationData > >::iterator iterator
const_iterator end() const
virtual void initialize()=0
bool evaluate_transient_terms
std::vector< double > gather_seeds
WorksetDescriptor bcDescriptor(const panzer::BC &bc)
BCType
Type of boundary condition.
Teuchos::RCP< GlobalEvaluationDataContainer > gedc
std::string second_sensitivities_name
std::string first_sensitivities_name