Panzer  Version of the Day
Panzer_ModelEvaluator.hpp
Go to the documentation of this file.
1 // @HEADER
2 // ***********************************************************************
3 //
4 // Panzer: A partial differential equation assembly
5 // engine for strongly coupled complex multiphysics systems
6 // Copyright (2011) Sandia Corporation
7 //
8 // Under the terms of Contract DE-AC04-94AL85000 with Sandia Corporation,
9 // the U.S. Government retains certain rights in this software.
10 //
11 // Redistribution and use in source and binary forms, with or without
12 // modification, are permitted provided that the following conditions are
13 // met:
14 //
15 // 1. Redistributions of source code must retain the above copyright
16 // notice, this list of conditions and the following disclaimer.
17 //
18 // 2. Redistributions in binary form must reproduce the above copyright
19 // notice, this list of conditions and the following disclaimer in the
20 // documentation and/or other materials provided with the distribution.
21 //
22 // 3. Neither the name of the Corporation nor the names of the
23 // contributors may be used to endorse or promote products derived from
24 // this software without specific prior written permission.
25 //
26 // THIS SOFTWARE IS PROVIDED BY SANDIA CORPORATION "AS IS" AND ANY
27 // EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
28 // IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
29 // PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL SANDIA CORPORATION OR THE
30 // CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
31 // EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
32 // PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
33 // PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
34 // LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
35 // NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
36 // SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
37 //
38 // Questions? Contact Roger P. Pawlowski (rppawlo@sandia.gov) and
39 // Eric C. Cyr (eccyr@sandia.gov)
40 // ***********************************************************************
41 // @HEADER
42 
43 #ifndef PANZER_MODEL_EVALUATOR_DECL_HPP
44 #define PANZER_MODEL_EVALUATOR_DECL_HPP
45 
46 #include "PanzerDiscFE_config.hpp"
47 
48 #include "Panzer_Traits.hpp"
55 
56 #include "Teuchos_RCP.hpp"
57 #include "Teuchos_AbstractFactory.hpp"
58 
59 #include "Thyra_VectorBase.hpp"
60 #include "Thyra_VectorSpaceBase.hpp"
61 #include "Thyra_StateFuncModelEvaluatorBase.hpp"
62 #include "Thyra_LinearOpWithSolveFactoryBase.hpp"
63 
64 #include <Panzer_NodeType.hpp>
65 
66 namespace panzer {
67 
68 class FieldManagerBuilder;
69 template<typename> class LinearObjFactory;
70 struct GlobalData;
71 
72 template<typename Scalar>
74  : public Thyra::StateFuncModelEvaluatorBase<Scalar>
75 {
76 public:
77 
78 public:
79 
82 
86  const std::vector<Teuchos::RCP<Teuchos::Array<std::string> > >& p_names,
87  const std::vector<Teuchos::RCP<Teuchos::Array<double> > >& p_values,
88  const Teuchos::RCP<const Thyra::LinearOpWithSolveFactoryBase<Scalar> > & solverFactory,
89  const Teuchos::RCP<panzer::GlobalData>& global_data,
90  bool build_transient_support,double t_init);
91 
93  const Teuchos::RCP<const Thyra::LinearOpWithSolveFactoryBase<Scalar> > & solverFactory,
94  const Teuchos::RCP<panzer::GlobalData>& global_data,
95  bool build_transient_support,double t_init);
96 
99 
101 
104 
107 
110 
113 
116 
119 
121  const std::string & get_g_name(int i) const;
122 
125 
128 
131 
134 
136  Thyra::ModelEvaluatorBase::InArgs<Scalar> createInArgs() const;
137 
138  Thyra::ModelEvaluatorBase::InArgs<Scalar> getNominalValues() const;
139 
141 
143  const std::vector<Teuchos::RCP<panzer::PhysicsBlock> >& physicsBlocks,
144  const std::vector<panzer::BC> & bcs,
145  const panzer::EquationSetFactory & eqset_factory,
146  const panzer::BCStrategyFactory& bc_factory,
149  const Teuchos::ParameterList& closure_models,
150  const Teuchos::ParameterList& user_data,
151  bool writeGraph=false,const std::string & graphPrefix="",
152  const Teuchos::ParameterList& me_params = Teuchos::ParameterList());
153 
165  int addParameter(const std::string & name,const Scalar & initial);
166 
177  int addParameter(const Teuchos::Array<std::string> & names,
178  const Teuchos::Array<Scalar> & initialValues);
179 
196  int addDistributedParameter(const std::string & name,
197  const Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> > & vs,
199  const Teuchos::RCP<const Thyra::VectorBase<Scalar> > & initial,
200  const Teuchos::RCP<const UniqueGlobalIndexerBase> & ugi=Teuchos::null);
201 
210  void addNonParameterGlobalEvaluationData(const std::string & name,
212 
229  int addFlexibleResponse(const std::string & responseName,
230  const std::vector<WorksetDescriptor> & wkst_desc,
232 
248  template <typename ResponseEvaluatorFactory_BuilderT>
249  int addResponse(const std::string & responseName,
250  const std::vector<WorksetDescriptor> & wkst_desc,
251  const ResponseEvaluatorFactory_BuilderT & builder);
252 
257  const std::vector<Teuchos::RCP<panzer::PhysicsBlock> >& physicsBlocks,
258  const panzer::EquationSetFactory & eqset_factory,
260  const Teuchos::ParameterList& closure_models,
261  const Teuchos::ParameterList& user_data,
262  const bool write_graphviz_file=false,
263  const std::string& graphviz_file_prefix="")
264  { responseLibrary_->buildResponseEvaluators(physicsBlocks,eqset_factory,cm_factory,closure_models,user_data,write_graphviz_file,graphviz_file_prefix);
267 
268  typedef Thyra::ModelEvaluatorBase MEB;
269  MEB::OutArgsSetup<Scalar> outArgs;
270  outArgs.setModelEvalDescription(this->description());
271  outArgs.set_Np_Ng(num_me_parameters_, responses_.size());
272  outArgs.setSupports(MEB::OUT_ARG_f);
273  outArgs.setSupports(MEB::OUT_ARG_W_op);
274  prototypeOutArgs_ = outArgs; }
275 
280  const std::vector<Teuchos::RCP<panzer::PhysicsBlock> >& physicsBlocks,
282  const Teuchos::ParameterList& closure_models,
283  const Teuchos::ParameterList& user_data,
284  const bool write_graphviz_file=false,
285  const std::string& graphviz_file_prefix="")
286  { responseLibrary_->buildResponseEvaluators(physicsBlocks,cm_factory,closure_models,user_data,write_graphviz_file,graphviz_file_prefix);
289 
290  typedef Thyra::ModelEvaluatorBase MEB;
291  MEB::OutArgsSetup<Scalar> outArgs;
292  outArgs.setModelEvalDescription(this->description());
293  outArgs.set_Np_Ng(num_me_parameters_, responses_.size());
294  outArgs.setSupports(MEB::OUT_ARG_f);
295  outArgs.setSupports(MEB::OUT_ARG_W_op);
296  prototypeOutArgs_ = outArgs; }
297 
305  const std::vector<Teuchos::RCP<panzer::PhysicsBlock> >& physicsBlocks,
306  const std::vector<panzer::BC> & bcs,
307  const panzer::EquationSetFactory & eqset_factory,
308  const panzer::BCStrategyFactory& bc_factory,
310  const Teuchos::ParameterList& closure_models,
311  const Teuchos::ParameterList& user_data,
312  const bool write_graphviz_file=false,
313  const std::string& graphviz_file_prefix="");
314 
323  const std::vector<Teuchos::RCP<panzer::PhysicsBlock> >& physicsBlocks,
324  const std::vector<panzer::BC> & bcs,
325  const panzer::EquationSetFactory & eqset_factory,
326  const panzer::BCStrategyFactory& bc_factory,
328  const Teuchos::ParameterList& closure_models,
329  const Teuchos::ParameterList& user_data,
330  const bool write_graphviz_file=false,
331  const std::string& graphviz_file_prefix="");
332 
340  void setOneTimeDirichletBeta(const Scalar & beta) const;
341 
346  const Teuchos::RCP<Thyra::VectorBase<Scalar> > & f) const;
347 
353  void setupAssemblyInArgs(const Thyra::ModelEvaluatorBase::InArgs<Scalar> & inArgs,
354  panzer::AssemblyEngineInArgs & ae_inargs) const;
355 
356 
362  return ae_tm_;
363  }
364 
366  { return responseLibrary_; }
367 
370  int getXTangentVectorIndex(const int index) const {
371  int v_index = 0;
372  for (int i=0; i<index; i++) {
373  if (!parameters_[i]->is_distributed)
374  ++v_index;
375  }
376  return v_index + parameters_.size();
377  }
378 
381  int getXDotTangentVectorIndex(const int index) const {
382  int v_index = 0;
383  for (int i=0; i<index; i++) {
384  if (!parameters_[i]->is_distributed)
385  ++v_index;
386  }
387  return v_index + parameters_.size() + tangent_space_.size();
388  }
389 
393  return parameters_[index]->initial_value;
394  }
395 
403  void evalModel_D2gDx2(int rIndex,
404  const Thyra::ModelEvaluatorBase::InArgs<Scalar> & inArgs,
405  const Teuchos::RCP<const Thyra::VectorBase<Scalar> > & delta_x,
406  const Teuchos::RCP<Thyra::VectorBase<Scalar> > & D2gDx2) const;
407 
416  void evalModel_D2gDp2(int rIndex,
417  int pIndex,
418  const Thyra::ModelEvaluatorBase::InArgs<Scalar> & inArgs,
419  const Teuchos::RCP<const Thyra::VectorBase<Scalar> > & delta_x,
420  const Teuchos::RCP<Thyra::VectorBase<Scalar> > & D2gDp2) const;
421 
430  void evalModel_D2gDpDx(int rIndex,
431  int pIndex,
432  const Thyra::ModelEvaluatorBase::InArgs<Scalar> & inArgs,
433  const Teuchos::RCP<const Thyra::VectorBase<Scalar> > & delta_x,
434  const Teuchos::RCP<Thyra::VectorBase<Scalar> > & D2gDpDx) const;
435 
444  void evalModel_D2gDxDp(int rIndex,
445  int pIndex,
446  const Thyra::ModelEvaluatorBase::InArgs<Scalar> & inArgs,
447  const Teuchos::RCP<const Thyra::VectorBase<Scalar> > & delta_p,
448  const Teuchos::RCP<Thyra::VectorBase<Scalar> > & D2gDxDp) const;
449 
457  void evalModel_D2fDx2(const Thyra::ModelEvaluatorBase::InArgs<Scalar> & inArgs,
458  const Teuchos::RCP<const Thyra::VectorBase<Scalar> > & delta_x,
459  const Teuchos::RCP<Thyra::LinearOpBase<Scalar> > & D2fDx2) const;
460 
469  void evalModel_D2fDp2(int pIndex,
470  const Thyra::ModelEvaluatorBase::InArgs<Scalar> & inArgs,
471  const Teuchos::RCP<const Thyra::VectorBase<Scalar> > & delta_x,
472  const Teuchos::RCP<Thyra::LinearOpBase<Scalar> > & D2fDp2) const;
473 
482  void evalModel_D2fDpDx(int pIndex,
483  const Thyra::ModelEvaluatorBase::InArgs<Scalar> & inArgs,
484  const Teuchos::RCP<const Thyra::VectorBase<Scalar> > & delta_x,
485  const Teuchos::RCP<Thyra::LinearOpBase<Scalar> > & D2fDpDx) const;
486 
495  void evalModel_D2fDxDp(int pIndex,
496  const Thyra::ModelEvaluatorBase::InArgs<Scalar> & inArgs,
497  const Teuchos::RCP<const Thyra::VectorBase<Scalar> > & delta_p,
498  const Teuchos::RCP<Thyra::LinearOpBase<Scalar> > & D2fDxDp) const;
499 
500 protected:
501 
504 
506  Thyra::ModelEvaluatorBase::OutArgs<Scalar> createOutArgsImpl() const;
507 
509  virtual void evalModelImpl(const Thyra::ModelEvaluatorBase::InArgs<Scalar> &inArgs,
510  const Thyra::ModelEvaluatorBase::OutArgs<Scalar> &outArgs) const;
511 
513 
515  virtual void evalModelImpl_basic(const Thyra::ModelEvaluatorBase::InArgs<Scalar> &inArgs,
516  const Thyra::ModelEvaluatorBase::OutArgs<Scalar> &outArgs) const;
517 
519  virtual void evalModelImpl_basic_g(const Thyra::ModelEvaluatorBase::InArgs<Scalar> &inArgs,
520  const Thyra::ModelEvaluatorBase::OutArgs<Scalar> &outArgs) const;
521 
527  virtual void evalModelImpl_basic_dgdx(const Thyra::ModelEvaluatorBase::InArgs<Scalar> &inArgs,
528  const Thyra::ModelEvaluatorBase::OutArgs<Scalar> &outArgs) const;
529 
535  virtual void evalModelImpl_basic_dgdp_scalar(const Thyra::ModelEvaluatorBase::InArgs<Scalar> &inArgs,
536  const Thyra::ModelEvaluatorBase::OutArgs<Scalar> &outArgs) const;
537 
543  virtual void evalModelImpl_basic_dgdp_distro(const Thyra::ModelEvaluatorBase::InArgs<Scalar> &inArgs,
544  const Thyra::ModelEvaluatorBase::OutArgs<Scalar> &outArgs) const;
545 
551  virtual void evalModelImpl_basic_dfdp_scalar(const Thyra::ModelEvaluatorBase::InArgs<Scalar> &inArgs,
552  const Thyra::ModelEvaluatorBase::OutArgs<Scalar> &outArgs) const;
553 
559  virtual void evalModelImpl_basic_dfdp_scalar_fd(const Thyra::ModelEvaluatorBase::InArgs<Scalar> &inArgs,
560  const Thyra::ModelEvaluatorBase::OutArgs<Scalar> &outArgs) const;
561 
567  virtual void evalModelImpl_basic_dfdp_distro(const Thyra::ModelEvaluatorBase::InArgs<Scalar> &inArgs,
568  const Thyra::ModelEvaluatorBase::OutArgs<Scalar> &outArgs) const;
569 
571  bool required_basic_g(const Thyra::ModelEvaluatorBase::OutArgs<Scalar> &outArgs) const;
572 
574  bool required_basic_dgdx(const Thyra::ModelEvaluatorBase::OutArgs<Scalar> &outArgs) const;
575 
577  bool required_basic_dgdp_scalar(const Thyra::ModelEvaluatorBase::OutArgs<Scalar> &outArgs) const;
578 
580  bool required_basic_dgdp_distro(const Thyra::ModelEvaluatorBase::OutArgs<Scalar> &outArgs) const;
581 
583  bool required_basic_dfdp_scalar(const Thyra::ModelEvaluatorBase::OutArgs<Scalar> &outArgs) const;
584 
586  bool required_basic_dfdp_distro(const Thyra::ModelEvaluatorBase::OutArgs<Scalar> &outArgs) const;
587 
589  void initializeNominalValues() const;
590 
591  // Set parameters in model supplied by inArgs
592  void setParameters(const Thyra::ModelEvaluatorBase::InArgs<Scalar> &inArgs) const;
593 
594  // Reset parameters back to nominal values
595  void resetParameters() const;
596 
597 private: // data members
598 
600  bool is_distributed; // or (is scalar?)
604 
605  // for distributed parameters
608  // for residual sensitivities with respect to a distributed parameter
610  // for response sensitivities with respect to a distributed parameter
611 
612  // for scalar parameters
614  };
615 
616  struct ResponseObject {
617  std::string name;
619 
620  // for distributed parameter sensitivities
622  // used for delayed construction of dgdp (distributed parameter) responses
623  std::vector<WorksetDescriptor> wkst_desc;
624  // used for delayed construction of dgdp (distributed parameter) responses
625 
626  struct SearchName {
627  std::string name;
628  SearchName(const std::string & n) : name(n) {}
629  bool operator()(const Teuchos::RCP<ResponseObject> & ro) { return name==ro->name; }
630  };
631  };
632 
634  const Teuchos::Array<Scalar> & in_values) const;
636  const Teuchos::RCP<const Thyra::VectorSpaceBase<Scalar> > & vs,
637  const Teuchos::RCP<const Thyra::VectorBase<Scalar> > & initial,
639 
640  double t_init_;
641 
644 
645  mutable Thyra::ModelEvaluatorBase::InArgs<Scalar> prototypeInArgs_;
646  mutable Thyra::ModelEvaluatorBase::OutArgs<Scalar> prototypeOutArgs_;
647 
648  mutable Thyra::ModelEvaluatorBase::InArgs<Scalar> nominalValues_;
649 
650  mutable panzer::AssemblyEngine_TemplateManager<panzer::Traits> ae_tm_; // they control and provide access to evaluate
651 
652  std::vector<Teuchos::RCP<ParameterObject> > parameters_;
653  std::vector<Teuchos::RCP<Thyra::VectorSpaceBase<double> > > tangent_space_;
657 
660 
661  // responses
663  std::vector<Teuchos::RCP<ResponseObject> > responses_;
664 
667 
668  // basic specific linear object objects
673 
674 
676 
679 
681  mutable Scalar oneTimeDirichletBeta_;
682 };
683 
684 // Inline definition of the add response (its template on the builder type)
685 template<typename Scalar>
686 template <typename ResponseEvaluatorFactory_BuilderT>
688 addResponse(const std::string & responseName,
689  const std::vector<WorksetDescriptor> & wkst_desc,
690  const ResponseEvaluatorFactory_BuilderT & builder)
691 {
692  using Teuchos::RCP;
693  using Teuchos::rcp;
694 
695  // see if the response evaluators have been constucted yet
696  TEUCHOS_TEST_FOR_EXCEPTION(responseLibrary_->responseEvaluatorsBuilt(),std::logic_error,
697  "panzer::ModelEvaluator::addResponse: Response with name \"" << responseName << "\" "
698  "cannot be added to the model evaluator because evalModel has already been called!");
699 
700  // add the response
701  responseLibrary_->addResponse(responseName,wkst_desc,builder);
702 
703  // check that the response can be found
704  TEUCHOS_TEST_FOR_EXCEPTION(std::find_if(responses_.begin(),responses_.end(),typename ResponseObject::SearchName(responseName))!=responses_.end(),
705  std::logic_error,
706  "panzer::ModelEvaluator::addResponse: Response with name \"" << responseName << "\" "
707  "has already been added to the model evaluator!");
708 
709  // allocate response object
710  RCP<ResponseObject> respObject = rcp(new ResponseObject);
711 
712  // handle panzer::Traits::Residual
713  {
714  // check that at least there is a response value
715  Teuchos::RCP<panzer::ResponseBase> respBase = responseLibrary_->getResponse<panzer::Traits::Residual>(responseName);
716  TEUCHOS_TEST_FOR_EXCEPTION(respBase==Teuchos::null,std::logic_error,
717  "panzer::ModelEvaluator::addResponse: Response with name \"" << responseName << "\" "
718  "has no residual type! Not sure what is going on!");
719 
720  // check that the response supports interactions with the model evaluator
722  Teuchos::rcp_dynamic_cast<panzer::ResponseMESupportBase<panzer::Traits::Residual> >(respBase);
723  TEUCHOS_TEST_FOR_EXCEPTION(resp==Teuchos::null,std::logic_error,
724  "panzer::ModelEvaluator::addResponse: Response with name \"" << responseName << "\" "
725  "resulted in bad cast to panzer::ResponseMESupportBase, the type of the response is incompatible!");
726 
727  // set the response in the model evaluator
728  Teuchos::RCP<const Thyra::VectorSpaceBase<double> > vs = resp->getVectorSpace();
729  respObject->space = vs;
730 
731  // lets be cautious and set a vector on the response
732  resp->setVector(Thyra::createMember(vs));
733  }
734 
735  // handle panzer::Traits::Jacobian (do a quick safety check, response is null or appropriate for jacobian)
736  Teuchos::RCP<panzer::ResponseBase> respJacBase = responseLibrary_->getResponse<panzer::Traits::Jacobian>(responseName);
737  if(respJacBase!=Teuchos::null) {
738  typedef panzer::Traits::Jacobian RespEvalT;
739 
740  // check that the response supports interactions with the model evaluator
742  Teuchos::rcp_dynamic_cast<panzer::ResponseMESupportBase<RespEvalT> >(respJacBase);
743  TEUCHOS_TEST_FOR_EXCEPTION(resp==Teuchos::null,std::logic_error,
744  "panzer::ModelEvaluator::addResponse: Response with name \"" << responseName <<
745  "\" resulted in bad cast to panzer::ResponseMESupportBase<Jacobian>, the type "
746  "of the response is incompatible!");
747 
748  // setup the vector (register response as epetra)
749  if(resp->supportsDerivative())
750  resp->setDerivative(resp->buildDerivative());
751  }
752 
753 #ifdef Panzer_BUILD_HESSIAN_SUPPORT
754  // handle panzer::Traits::Hessian (do a quick safety check, response is null or appropriate for jacobian)
755  Teuchos::RCP<panzer::ResponseBase> respHesBase = responseLibrary_->getResponse<panzer::Traits::Hessian>(responseName);
756  if(respHesBase!=Teuchos::null) {
757  typedef panzer::Traits::Hessian RespEvalT;
758 
759  // check that the response supports interactions with the model evaluator
761  Teuchos::rcp_dynamic_cast<panzer::ResponseMESupportBase<RespEvalT> >(respHesBase);
762  TEUCHOS_TEST_FOR_EXCEPTION(resp==Teuchos::null,std::logic_error,
763  "panzer::ModelEvaluator::addResponse: Response with name \"" << responseName <<
764  "\" resulted in bad cast to panzer::ResponseMESupportBase<Hessian>, the type "
765  "of the response is incompatible!");
766 
767  // setup the vector (register response as epetra)
768  if(resp->supportsDerivative())
769  resp->setDerivative(resp->buildDerivative());
770  }
771 #endif
772 
773  respObject->name = responseName;
774  respObject->wkst_desc = wkst_desc;
775 
776  responses_.push_back(respObject);
777 
778  require_in_args_refresh_ = true;
779  require_out_args_refresh_ = true;
780 
781  return responses_.size()-1;
782 }
783 
784 
785 }
786 
787 // #include "Panzer_ModelEvaluator_impl.hpp"
788 
789 #endif
Interface for constructing a BCStrategy_TemplateManager.
std::vector< WorksetDescriptor > wkst_desc
virtual void evalModelImpl_basic_dfdp_scalar(const Thyra::ModelEvaluatorBase::InArgs< Scalar > &inArgs, const Thyra::ModelEvaluatorBase::OutArgs< Scalar > &outArgs) const
Allocates and initializes an equation set template manager.
Thyra::ModelEvaluatorBase::InArgs< Scalar > getNominalValues() const
Teuchos::RCP< panzer::ResponseLibrary< panzer::Traits > > dfdp_rl
std::vector< Teuchos::RCP< ParameterObject > > parameters_
void setParameters(const Thyra::ModelEvaluatorBase::InArgs< Scalar > &inArgs) const
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_g_space(int i) const
GlobalEvaluationDataContainer nonParamGlobalEvaluationData_
int addFlexibleResponse(const std::string &responseName, const std::vector< WorksetDescriptor > &wkst_desc, const Teuchos::RCP< ResponseMESupportBuilderBase > &builder)
int getXDotTangentVectorIndex(const int index) const
Teuchos::RCP< const Thyra::VectorBase< Scalar > > initial_value
const std::string & get_g_name(int i) const
Thyra::ModelEvaluatorBase::InArgs< Scalar > prototypeInArgs_
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > space
#define TEUCHOS_TEST_FOR_EXCEPTION(throw_exception_test, Exception, msg)
panzer::AssemblyEngine_TemplateManager< panzer::Traits > ae_tm_
void buildDistroParamDgDp_RL(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 > &cm_factory, const Teuchos::ParameterList &closure_models, const Teuchos::ParameterList &user_data, const bool write_graphviz_file=false, const std::string &graphviz_file_prefix="")
std::vector< Teuchos::RCP< ResponseObject > > responses_
void evalModel_D2gDpDx(int rIndex, int pIndex, const Thyra::ModelEvaluatorBase::InArgs< Scalar > &inArgs, const Teuchos::RCP< const Thyra::VectorBase< Scalar > > &delta_x, const Teuchos::RCP< Thyra::VectorBase< Scalar > > &D2gDpDx) const
void evalModel_D2gDx2(int rIndex, const Thyra::ModelEvaluatorBase::InArgs< Scalar > &inArgs, const Teuchos::RCP< const Thyra::VectorBase< Scalar > > &delta_x, const Teuchos::RCP< Thyra::VectorBase< Scalar > > &D2gDx2) const
Teuchos::RCP< Teuchos::Array< std::string > > names
Teuchos::RCP< panzer::GlobalData > global_data_
void setupAssemblyInArgs(const Thyra::ModelEvaluatorBase::InArgs< Scalar > &inArgs, panzer::AssemblyEngineInArgs &ae_inargs) const
virtual void evalModelImpl_basic_dgdp_distro(const Thyra::ModelEvaluatorBase::InArgs< Scalar > &inArgs, const Thyra::ModelEvaluatorBase::OutArgs< Scalar > &outArgs) const
void evalModel_D2fDp2(int pIndex, const Thyra::ModelEvaluatorBase::InArgs< Scalar > &inArgs, const Teuchos::RCP< const Thyra::VectorBase< Scalar > > &delta_x, const Teuchos::RCP< Thyra::LinearOpBase< Scalar > > &D2fDp2) const
void initializeNominalValues() const
Initialize the nominal values with good starting conditions.
Thyra::ModelEvaluatorBase::OutArgs< Scalar > createOutArgsImpl() const
virtual void evalModelImpl_basic_dgdp_scalar(const Thyra::ModelEvaluatorBase::InArgs< Scalar > &inArgs, const Thyra::ModelEvaluatorBase::OutArgs< Scalar > &outArgs) const
Sacado::ScalarParameterVector< panzer::EvaluationTraits > ParamVec
Teuchos::RCP< panzer::ResponseLibrary< panzer::Traits > > dgdp_rl
Teuchos::RCP< const Thyra::LinearOpWithSolveFactoryBase< Scalar > > get_W_factory() const
bool operator()(const Teuchos::RCP< ResponseObject > &ro)
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_x_space() const
int addResponse(const std::string &responseName, const std::vector< WorksetDescriptor > &wkst_desc, const ResponseEvaluatorFactory_BuilderT &builder)
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > x_space_
void evalModel_D2gDp2(int rIndex, int pIndex, const Thyra::ModelEvaluatorBase::InArgs< Scalar > &inArgs, const Teuchos::RCP< const Thyra::VectorBase< Scalar > > &delta_x, const Teuchos::RCP< Thyra::VectorBase< Scalar > > &D2gDp2) const
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > f_space_
bool required_basic_dfdp_scalar(const Thyra::ModelEvaluatorBase::OutArgs< Scalar > &outArgs) const
Are derivatives of the residual with respect to the scalar parameters in the out args? DfDp.
Teuchos::ArrayView< const std::string > get_g_names(int i) const override
Teuchos::RCP< const UniqueGlobalIndexerBase > global_indexer
panzer::AssemblyEngine_TemplateManager< panzer::Traits > getAssemblyEngineTemplateManager() const
return a copy of the model evaluators template manager, this is shallow class so pass by value ...
bool required_basic_dfdp_distro(const Thyra::ModelEvaluatorBase::OutArgs< Scalar > &outArgs) const
Are derivatives of the residual with respect to the distributed parameters in the out args...
TEUCHOS_DEPRECATED RCP< T > rcp(T *p, Dealloc_T dealloc, bool owns_mem)
virtual void evalModelImpl(const Thyra::ModelEvaluatorBase::InArgs< Scalar > &inArgs, const Thyra::ModelEvaluatorBase::OutArgs< Scalar > &outArgs) const
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_p_space(int i) const
GlobalEvaluationDataContainer distrParamGlobalEvaluationData_
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > space
Teuchos::RCP< Thyra::LinearOpBase< Scalar > > create_DfDp_op(int i) const
Teuchos::RCP< const panzer::LinearObjFactory< panzer::Traits > > lof_
void buildResponses(const std::vector< Teuchos::RCP< panzer::PhysicsBlock > > &physicsBlocks, const panzer::ClosureModelFactory_TemplateManager< panzer::Traits > &cm_factory, const Teuchos::ParameterList &closure_models, const Teuchos::ParameterList &user_data, const bool write_graphviz_file=false, const std::string &graphviz_file_prefix="")
void applyDirichletBCs(const Teuchos::RCP< Thyra::VectorBase< Scalar > > &x, const Teuchos::RCP< Thyra::VectorBase< Scalar > > &f) const
Thyra::ModelEvaluatorBase::OutArgs< Scalar > prototypeOutArgs_
std::vector< Teuchos::RCP< Thyra::VectorSpaceBase< double > > > tangent_space_
Teuchos::RCP< const Thyra::LinearOpWithSolveFactoryBase< Scalar > > solverFactory_
int getXTangentVectorIndex(const int index) const
Teuchos::RCP< ParameterObject > createDistributedParameter(const std::string &key, const Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > &vs, const Teuchos::RCP< const Thyra::VectorBase< Scalar > > &initial, const Teuchos::RCP< const UniqueGlobalIndexerBase > &ugi) const
Teuchos::RCP< const Teuchos::Array< std::string > > get_p_names(int i) const
bool required_basic_g(const Thyra::ModelEvaluatorBase::OutArgs< Scalar > &outArgs) const
Does this set of out args require a simple response?
bool required_basic_dgdx(const Thyra::ModelEvaluatorBase::OutArgs< Scalar > &outArgs) const
Are their required responses in the out args? DgDx.
Teuchos::RCP< ParameterObject > createScalarParameter(const Teuchos::Array< std::string > &names, const Teuchos::Array< Scalar > &in_values) const
Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > get_f_space() const
void buildDistroParamDfDp_RL(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 > &cm_factory, const Teuchos::ParameterList &closure_models, const Teuchos::ParameterList &user_data, const bool write_graphviz_file=false, const std::string &graphviz_file_prefix="")
Thyra::ModelEvaluatorBase::InArgs< Scalar > createInArgs() const
void setupModel(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 Teuchos::ParameterList &user_data, bool writeGraph=false, const std::string &graphPrefix="", const Teuchos::ParameterList &me_params=Teuchos::ParameterList())
void evalModel_D2fDpDx(int pIndex, const Thyra::ModelEvaluatorBase::InArgs< Scalar > &inArgs, const Teuchos::RCP< const Thyra::VectorBase< Scalar > > &delta_x, const Teuchos::RCP< Thyra::LinearOpBase< Scalar > > &D2fDpDx) const
void evalModel_D2gDxDp(int rIndex, int pIndex, const Thyra::ModelEvaluatorBase::InArgs< Scalar > &inArgs, const Teuchos::RCP< const Thyra::VectorBase< Scalar > > &delta_p, const Teuchos::RCP< Thyra::VectorBase< Scalar > > &D2gDxDp) const
void evalModel_D2fDx2(const Thyra::ModelEvaluatorBase::InArgs< Scalar > &inArgs, const Teuchos::RCP< const Thyra::VectorBase< Scalar > > &delta_x, const Teuchos::RCP< Thyra::LinearOpBase< Scalar > > &D2fDx2) const
Teuchos::RCP< panzer::ResponseLibrary< panzer::Traits > > getResponseLibrary() const
virtual void evalModelImpl_basic_g(const Thyra::ModelEvaluatorBase::InArgs< Scalar > &inArgs, const Thyra::ModelEvaluatorBase::OutArgs< Scalar > &outArgs) const
Construct a simple response dicatated by this set of out args.
Teuchos::RCP< Thyra::LinearOpBase< Scalar > > create_W_op() const
int addDistributedParameter(const std::string &name, const Teuchos::RCP< const Thyra::VectorSpaceBase< Scalar > > &vs, const Teuchos::RCP< GlobalEvaluationData > &ged, const Teuchos::RCP< const Thyra::VectorBase< Scalar > > &initial, const Teuchos::RCP< const UniqueGlobalIndexerBase > &ugi=Teuchos::null)
void setOneTimeDirichletBeta(const Scalar &beta) const
Teuchos::RCP< ResponseMESupportBuilderBase > builder
Teuchos::RCP< panzer::LinearObjContainer > ghostedContainer_
void buildResponses(const std::vector< Teuchos::RCP< panzer::PhysicsBlock > > &physicsBlocks, const panzer::EquationSetFactory &eqset_factory, const panzer::ClosureModelFactory_TemplateManager< panzer::Traits > &cm_factory, const Teuchos::ParameterList &closure_models, const Teuchos::ParameterList &user_data, const bool write_graphviz_file=false, const std::string &graphviz_file_prefix="")
Teuchos::RCP< panzer::ResponseLibrary< panzer::Traits > > responseLibrary_
void addNonParameterGlobalEvaluationData(const std::string &name, const Teuchos::RCP< GlobalEvaluationData > &ged)
virtual void evalModelImpl_basic(const Thyra::ModelEvaluatorBase::InArgs< Scalar > &inArgs, const Thyra::ModelEvaluatorBase::OutArgs< Scalar > &outArgs) const
Evaluate a simple model, meaning a residual and a jacobian, no fancy stochastic galerkin or multipoin...
Teuchos::RCP< ReadOnlyVector_GlobalEvaluationData > xdotContainer_
Teuchos::RCP< const Thyra::VectorBase< Scalar > > get_parameter_vector(int index) const
Thyra::ModelEvaluatorBase::InArgs< Scalar > nominalValues_
bool required_basic_dgdp_scalar(const Thyra::ModelEvaluatorBase::OutArgs< Scalar > &outArgs) const
Are their required responses in the out args? DgDp.
virtual void evalModelImpl_basic_dfdp_distro(const Thyra::ModelEvaluatorBase::InArgs< Scalar > &inArgs, const Thyra::ModelEvaluatorBase::OutArgs< Scalar > &outArgs) const
virtual void evalModelImpl_basic_dgdx(const Thyra::ModelEvaluatorBase::InArgs< Scalar > &inArgs, const Thyra::ModelEvaluatorBase::OutArgs< Scalar > &outArgs) const
bool required_basic_dgdp_distro(const Thyra::ModelEvaluatorBase::OutArgs< Scalar > &outArgs) const
Are their required responses in the out args? DgDp.
Teuchos::RCP< ReadOnlyVector_GlobalEvaluationData > xContainer_
int addParameter(const std::string &name, const Scalar &initial)
void evalModel_D2fDxDp(int pIndex, const Thyra::ModelEvaluatorBase::InArgs< Scalar > &inArgs, const Teuchos::RCP< const Thyra::VectorBase< Scalar > > &delta_p, const Teuchos::RCP< Thyra::LinearOpBase< Scalar > > &D2fDxDp) const
virtual void evalModelImpl_basic_dfdp_scalar_fd(const Thyra::ModelEvaluatorBase::InArgs< Scalar > &inArgs, const Thyra::ModelEvaluatorBase::OutArgs< Scalar > &outArgs) const