43#ifndef PANZER_STK_MODEL_EVALUATOR_FACTORY_HPP
44#define PANZER_STK_MODEL_EVALUATOR_FACTORY_HPP
50#include "Teuchos_RCP.hpp"
51#include "Teuchos_Ptr.hpp"
52#include "Teuchos_Comm.hpp"
53#include "Teuchos_DefaultMpiComm.hpp"
54#include "Teuchos_ParameterList.hpp"
55#include "Teuchos_ParameterListAcceptorDefaultBase.hpp"
57#include "PanzerAdaptersSTK_config.hpp"
71#include "Thyra_EpetraModelEvaluator.hpp"
73#ifdef PANZER_HAVE_TEKO
74#include "Teko_RequestHandler.hpp"
79#ifdef PANZER_HAVE_TEMPUS
80 template <
typename ScalarT>
class TempusSolverForwardOnly;
85 template<
typename ScalarT>
class ModelEvaluator;
86 template<
typename ScalarT>
class LinearOpWithSolveFactoryBase;
92 template <
typename>
class LinearObjFactory;
94 class BlockedDOFManager;
101 class STKConnManager;
102 class NOXObserverFactory;
103 class RythmosObserverFactory;
104#ifdef PANZER_HAVE_TEMPUS
105 class TempusObserverFactory;
107 class WorksetFactory;
109 template<
typename ScalarT>
116 void setParameterList(Teuchos::RCP<Teuchos::ParameterList>
const& paramList);
128 void buildObjects(
const Teuchos::RCP<
const Teuchos::Comm<int> >& comm,
129 const Teuchos::RCP<panzer::GlobalData>& global_data,
130 const Teuchos::RCP<const panzer::EquationSetFactory>& eqset_factory,
133 bool meConstructionOn=
true);
140 void setNOXObserverFactory(
const Teuchos::RCP<const panzer_stk::NOXObserverFactory>& nox_observer_factory);
144#ifdef PANZER_HAVE_TEMPUS
145 void setTempusObserverFactory(
const Teuchos::RCP<const panzer_stk::TempusObserverFactory>& tempus_observer_factory);
148 template <
typename BuilderT>
149 int addResponse(
const std::string & responseName,
const std::vector<panzer::WorksetDescriptor> & wkstDesc,
const BuilderT & builder);
152 const bool write_graphviz_file=
false,
153 const std::string& graphviz_file_prefix=
"");
157 Teuchos::RCP<Thyra::ModelEvaluator<ScalarT> >
159 const Teuchos::RCP<panzer::GlobalData>& global_data,
161#ifdef PANZER_HAVE_TEMPUS
162 const Teuchos::RCP<Piro::TempusSolverForwardOnly<ScalarT> > tempusSolver = Teuchos::null,
164 const Teuchos::Ptr<const panzer_stk::NOXObserverFactory> & in_nox_observer_factory=Teuchos::null,
165 const Teuchos::Ptr<const panzer_stk::RythmosObserverFactory> & in_rythmos_observer_factory=Teuchos::null
166#ifdef PANZER_HAVE_TEMPUS
167 ,
const Teuchos::Ptr<const panzer_stk::TempusObserverFactory> & in_tempus_observer_factory=Teuchos::null
178 const std::vector<Teuchos::RCP<panzer::PhysicsBlock> > &
getPhysicsBlocks()
const;
181 Teuchos::RCP<panzer_stk::STK_Interface>
getMesh()
const
206 Teuchos::RCP<Thyra::ModelEvaluator<double> >
208 const Teuchos::RCP<Teuchos::ParameterList> & physics_block_plist,
209 const Teuchos::RCP<const panzer::EquationSetFactory>& eqset_factory,
212 bool is_transient,
bool is_explicit,
213 const Teuchos::Ptr<const Teuchos::ParameterList> & bc_list=Teuchos::null,
221 const std::vector<Teuchos::RCP<panzer::PhysicsBlock> >& physicsBlocks,
224 const Teuchos::ParameterList & closure_pl,
225 const Teuchos::ParameterList & initial_cond_pl,
226 const Teuchos::ParameterList & user_data_pl,
227 bool write_dot_files,
const std::string & dot_file_prefix)
const;
233 const std::vector<Teuchos::RCP<panzer::PhysicsBlock> >& physicsBlocks,
234 const Teuchos::RCP<panzer::WorksetContainer> & wc,
235 const Teuchos::RCP<const panzer::GlobalIndexer> & ugi,
237 const Teuchos::RCP<panzer_stk::STK_Interface> & mesh,
239 const Teuchos::ParameterList & closure_model_pl,
240 const Teuchos::ParameterList & user_data_pl,
241 int workset_size)
const;
245 Teuchos::RCP<Thyra::ModelEvaluatorDefaultBase<double> >
247 const Teuchos::RCP<panzer::FieldManagerBuilder> & fmb,
250 const std::vector<Teuchos::RCP<Teuchos::Array<std::string> > > & p_names,
251 const std::vector<Teuchos::RCP<Teuchos::Array<double> > > & p_values,
253 const Teuchos::RCP<panzer::GlobalData> & global_data,
254 bool is_transient,
double t_init)
const;
265 double getInitialTime(Teuchos::ParameterList& transient_ic_params,
268 Teuchos::RCP<Thyra::LinearOpWithSolveFactoryBase<double> >
270 const Teuchos::RCP<const panzer::GlobalIndexer> & globalIndexer,
271 const Teuchos::RCP<panzer::ConnManager> & conn_manager,
272 const Teuchos::RCP<panzer_stk::STK_Interface> & mesh,
274 #ifdef PANZER_HAVE_TEKO
275 ,
const Teuchos::RCP<Teko::RequestHandler> & req_handler=Teuchos::null
287 Teuchos::RCP<STK_MeshFactory>
buildSTKMeshFactory(
const Teuchos::ParameterList & mesh_params)
const;
290 const std::vector<Teuchos::RCP<panzer::PhysicsBlock> > & physicsBlocks,
296 Teuchos::RCP<panzer::FieldManagerBuilder>
298 const std::vector<Teuchos::RCP<panzer::PhysicsBlock> >& physicsBlocks,
299 const std::vector<panzer::BC> & bcs,
304 const Teuchos::ParameterList& closure_models,
306 const Teuchos::ParameterList& user_data,
307 bool writeGraph,
const std::string & graphPrefix,
308 bool write_field_managers,
const std::string & field_manager_prefix)
const;
313 const Teuchos::RCP<panzer::WorksetContainer> & wc,
314 const Teuchos::RCP<const panzer::GlobalIndexer> & ugi,
316 const Teuchos::RCP<panzer_stk::STK_Interface> & mesh)
const;
321 const std::vector<Teuchos::RCP<panzer::PhysicsBlock> > & physicsBlocks,
323 const Teuchos::ParameterList & closure_models,
324 int workset_size, Teuchos::ParameterList & user_data)
const;
330 Teuchos::RCP<Thyra::ModelEvaluator<ScalarT> >
m_rome_me;
335 Teuchos::RCP<panzer_stk::STK_Interface>
m_mesh;
347#ifdef PANZER_HAVE_TEMPUS
348 Teuchos::RCP<const panzer_stk::TempusObserverFactory> m_tempus_observer_factory;
356template<
typename ScalarT>
357template <
typename BuilderT>
359addResponse(
const std::string & responseName,
const std::vector<panzer::WorksetDescriptor> & wkstDesc,
const BuilderT & builder)
363 Teuchos::RCP<Thyra::EpetraModelEvaluator> thyra_ep_me = Teuchos::rcp_dynamic_cast<Thyra::EpetraModelEvaluator>(m_physics_me);
364 Teuchos::RCP<PanzerME> panzer_me = Teuchos::rcp_dynamic_cast<PanzerME>(m_physics_me);
366 if(thyra_ep_me!=Teuchos::null && panzer_me==Teuchos::null) {
368 Teuchos::RCP<EpetraExt::ModelEvaluator> ep_me = Teuchos::rcp_const_cast<EpetraExt::ModelEvaluator>(thyra_ep_me->getEpetraModel());
369 Teuchos::RCP<panzer::ModelEvaluator_Epetra> ep_panzer_me = Teuchos::rcp_dynamic_cast<panzer::ModelEvaluator_Epetra>(ep_me);
371 return ep_panzer_me->addResponse(responseName,wkstDesc,builder);
373 else if(panzer_me!=Teuchos::null && thyra_ep_me==Teuchos::null) {
374 return panzer_me->addResponse(responseName,wkstDesc,builder);
377 TEUCHOS_ASSERT(
false);
Class that provides access to worksets on each element block and side set.
Teuchos::RCP< Thyra::LinearOpWithSolveFactoryBase< double > > buildLOWSFactory(bool blockedAssembly, const Teuchos::RCP< const panzer::GlobalIndexer > &globalIndexer, const Teuchos::RCP< panzer::ConnManager > &conn_manager, const Teuchos::RCP< panzer_stk::STK_Interface > &mesh, const Teuchos::RCP< const Teuchos::MpiComm< int > > &mpi_comm) const
Teuchos::RCP< panzer::ResponseLibrary< panzer::Traits > > getResponseLibrary()
Teuchos::RCP< panzer::ResponseLibrary< panzer::Traits > > initializeSolnWriterResponseLibrary(const Teuchos::RCP< panzer::WorksetContainer > &wc, const Teuchos::RCP< const panzer::GlobalIndexer > &ugi, const Teuchos::RCP< const panzer::LinearObjFactory< panzer::Traits > > &lof, const Teuchos::RCP< panzer_stk::STK_Interface > &mesh) const
void setRythmosObserverFactory(const Teuchos::RCP< const panzer_stk::RythmosObserverFactory > &rythmos_observer_factory)
Teuchos::RCP< panzer_stk::WorksetFactory > m_user_wkst_factory
Teuchos::RCP< panzer::FieldManagerBuilder > buildFieldManagerBuilder(const Teuchos::RCP< panzer::WorksetContainer > &wc, const std::vector< Teuchos::RCP< panzer::PhysicsBlock > > &physicsBlocks, const std::vector< panzer::BC > &bcs, const panzer::EquationSetFactory &eqset_factory, const panzer::BCStrategyFactory &bc_factory, const panzer::ClosureModelFactory_TemplateManager< panzer::Traits > &volume_cm_factory, const panzer::ClosureModelFactory_TemplateManager< panzer::Traits > &bc_cm_factory, const Teuchos::ParameterList &closure_models, const panzer::LinearObjFactory< panzer::Traits > &lo_factory, const Teuchos::ParameterList &user_data, bool writeGraph, const std::string &graphPrefix, bool write_field_managers, const std::string &field_manager_prefix) const
Teuchos::RCP< const panzer_stk::NOXObserverFactory > m_nox_observer_factory
Teuchos::RCP< panzer_stk::STK_Interface > m_mesh
Teuchos::RCP< Thyra::ModelEvaluator< ScalarT > > buildResponseOnlyModelEvaluator(const Teuchos::RCP< Thyra::ModelEvaluator< ScalarT > > &thyra_me, const Teuchos::RCP< panzer::GlobalData > &global_data, const Teuchos::RCP< Piro::RythmosSolver< ScalarT > > rythmosSolver=Teuchos::null, const Teuchos::Ptr< const panzer_stk::NOXObserverFactory > &in_nox_observer_factory=Teuchos::null, const Teuchos::Ptr< const panzer_stk::RythmosObserverFactory > &in_rythmos_observer_factory=Teuchos::null)
double getInitialTime(Teuchos::ParameterList &transient_ic_params, const panzer_stk::STK_Interface &mesh) const
Gets the initial time from either the input parameter list or an exodus file.
void finalizeMeshConstruction(const STK_MeshFactory &mesh_factory, const std::vector< Teuchos::RCP< panzer::PhysicsBlock > > &physicsBlocks, const Teuchos::MpiComm< int > mpi_comm, STK_Interface &mesh) const
Teuchos::RCP< Thyra::ModelEvaluator< ScalarT > > getResponseOnlyModelEvaluator()
void buildObjects(const Teuchos::RCP< const Teuchos::Comm< int > > &comm, const Teuchos::RCP< panzer::GlobalData > &global_data, const Teuchos::RCP< const panzer::EquationSetFactory > &eqset_factory, const panzer::BCStrategyFactory &bc_factory, const panzer::ClosureModelFactory_TemplateManager< panzer::Traits > &cm_factory, bool meConstructionOn=true)
Builds the model evaluators for a panzer assembly.
void addUserFieldsToMesh(panzer_stk::STK_Interface &mesh, const Teuchos::ParameterList &output_list) const
Add the user fields specified by output_list to the mesh.
Teuchos::RCP< panzer::GlobalIndexer > getGlobalIndexer() const
Get global indexer used to build model evaluator.
Teuchos::RCP< const panzer::EquationSetFactory > m_eqset_factory
Teuchos::RCP< Thyra::ModelEvaluator< double > > cloneWithNewPhysicsBlocks(const Teuchos::RCP< Thyra::LinearOpWithSolveFactoryBase< ScalarT > > &solverFactory, const Teuchos::RCP< Teuchos::ParameterList > &physics_block_plist, const Teuchos::RCP< const panzer::EquationSetFactory > &eqset_factory, const panzer::BCStrategyFactory &bc_factory, const panzer::ClosureModelFactory_TemplateManager< panzer::Traits > &user_cm_factory, bool is_transient, bool is_explicit, const Teuchos::Ptr< const Teuchos::ParameterList > &bc_list=Teuchos::null, const Teuchos::RCP< Thyra::ModelEvaluator< ScalarT > > &physics_me=Teuchos::null) const
Teuchos::RCP< panzer::ResponseLibrary< panzer::Traits > > m_response_library
Teuchos::RCP< panzer::WorksetContainer > getWorksetContainer() const
Get the workset container associated with the mesh database.
Teuchos::RCP< panzer::ConnManager > m_conn_manager
bool useDynamicCoordinates_
void buildResponses(const panzer::ClosureModelFactory_TemplateManager< panzer::Traits > &cm_factory, const bool write_graphviz_file=false, const std::string &graphviz_file_prefix="")
void setParameterList(Teuchos::RCP< Teuchos::ParameterList > const ¶mList)
bool isBlockedAssembly() const
Is blocked assembly?
void setupInitialConditions(Thyra::ModelEvaluator< ScalarT > &model, panzer::WorksetContainer &wkstContainer, const std::vector< Teuchos::RCP< panzer::PhysicsBlock > > &physicsBlocks, const panzer::ClosureModelFactory_TemplateManager< panzer::Traits > &cm_factory, const panzer::LinearObjFactory< panzer::Traits > &lof, const Teuchos::ParameterList &closure_pl, const Teuchos::ParameterList &initial_cond_pl, const Teuchos::ParameterList &user_data_pl, bool write_dot_files, const std::string &dot_file_prefix) const
Setup the initial conditions in a model evaluator. Note that this is entirely self contained.
Teuchos::RCP< panzer::WorksetContainer > m_wkstContainer
bool useDynamicCoordinates() const
void setUserWorksetFactory(Teuchos::RCP< panzer_stk::WorksetFactory > &user_wkst_factory)
Set user defined workset factory.
const std::vector< Teuchos::RCP< panzer::PhysicsBlock > > & getPhysicsBlocks() const
std::vector< Teuchos::RCP< panzer::PhysicsBlock > > m_physics_blocks
Teuchos::RCP< Thyra::ModelEvaluator< ScalarT > > m_rome_me
Teuchos::RCP< Thyra::ModelEvaluator< ScalarT > > getPhysicsModelEvaluator()
Teuchos::RCP< panzer_stk::STK_Interface > getMesh() const
Get mesh object used to build model evaluator.
Teuchos::RCP< panzer::LinearObjFactory< panzer::Traits > > m_lin_obj_factory
Teuchos::RCP< panzer::GlobalIndexer > m_global_indexer
Teuchos::RCP< const panzer_stk::RythmosObserverFactory > m_rythmos_observer_factory
Teuchos::RCP< panzer::GlobalData > m_global_data
Teuchos::RCP< panzer::LinearObjFactory< panzer::Traits > > getLinearObjFactory() const
Get linear object factory used to build model evaluator.
Teuchos::RCP< STK_MeshFactory > buildSTKMeshFactory(const Teuchos::ParameterList &mesh_params) const
build STK mesh factory from a mesh parameter list
void setNOXObserverFactory(const Teuchos::RCP< const panzer_stk::NOXObserverFactory > &nox_observer_factory)
void finalizeSolnWriterResponseLibrary(panzer::ResponseLibrary< panzer::Traits > &rl, const std::vector< Teuchos::RCP< panzer::PhysicsBlock > > &physicsBlocks, const panzer::ClosureModelFactory_TemplateManager< panzer::Traits > &cm_factory, const Teuchos::ParameterList &closure_models, int workset_size, Teuchos::ParameterList &user_data) const
Teuchos::RCP< Thyra::ModelEvaluatorDefaultBase< double > > buildPhysicsModelEvaluator(bool buildThyraME, const Teuchos::RCP< panzer::FieldManagerBuilder > &fmb, const Teuchos::RCP< panzer::ResponseLibrary< panzer::Traits > > &rLibrary, const Teuchos::RCP< panzer::LinearObjFactory< panzer::Traits > > &lof, const std::vector< Teuchos::RCP< Teuchos::Array< std::string > > > &p_names, const std::vector< Teuchos::RCP< Teuchos::Array< double > > > &p_values, const Teuchos::RCP< Thyra::LinearOpWithSolveFactoryBase< ScalarT > > &solverFactory, const Teuchos::RCP< panzer::GlobalData > &global_data, bool is_transient, double t_init) const
Teuchos::RCP< panzer::ConnManager > getConnManager() const
Get connection manager.
int addResponse(const std::string &responseName, const std::vector< panzer::WorksetDescriptor > &wkstDesc, const BuilderT &builder)
Teuchos::RCP< const Teuchos::ParameterList > getValidParameters() const
Teuchos::RCP< Thyra::ModelEvaluator< ScalarT > > m_physics_me
void writeInitialConditions(const Thyra::ModelEvaluator< ScalarT > &model, const std::vector< Teuchos::RCP< panzer::PhysicsBlock > > &physicsBlocks, const Teuchos::RCP< panzer::WorksetContainer > &wc, const Teuchos::RCP< const panzer::GlobalIndexer > &ugi, const Teuchos::RCP< const panzer::LinearObjFactory< panzer::Traits > > &lof, const Teuchos::RCP< panzer_stk::STK_Interface > &mesh, const panzer::ClosureModelFactory_TemplateManager< panzer::Traits > &cm_factory, const Teuchos::ParameterList &closure_model_pl, const Teuchos::ParameterList &user_data_pl, int workset_size) const
Write the initial conditions to exodus. Note that this is entirely self contained.
Interface for constructing a BCStrategy_TemplateManager.
Allocates and initializes an equation set template manager.