Panzer  Version of the Day
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
Panzer_ResponseEvaluatorFactory_Probe_impl.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_ResponseEvaluatorFactory_Probe_impl_hpp__
44 #define __Panzer_ResponseEvaluatorFactory_Probe_impl_hpp__
45 
46 #include <string>
47 
48 #include "PanzerDiscFE_config.hpp"
49 
51 #include "Panzer_PhysicsBlock.hpp"
52 #include "Panzer_CellExtreme.hpp"
55 
56 namespace panzer {
57 
58 template <typename EvalT,typename LO,typename GO>
60 buildResponseObject(const std::string & responseName) const
61 {
62  Teuchos::RCP<ResponseBase> response = Teuchos::rcp(new Response_Probe<EvalT>(responseName,comm_,linearObjFactory_));
63  response->setRequiresDirichletAdjustment(applyDirichletToDerivative_);
64 
65  return response;
66 }
67 
68 template <typename EvalT,typename LO,typename GO>
70 buildAndRegisterEvaluators(const std::string & responseName,
72  const panzer::PhysicsBlock & physicsBlock,
73  const Teuchos::ParameterList & user_data) const
74 {
75  using Teuchos::RCP;
76  using Teuchos::rcp;
77 
78  // build scatter evaluator
79  {
81  (globalIndexer_!=Teuchos::null) ? Teuchos::rcp(new ProbeScatter<LO,GO>(globalIndexer_)) : Teuchos::null;
82  std::string field = (fieldName_=="" ? responseName : fieldName_);
83 
84  // Get basis and integration rule associated with field
85  std::vector<panzer::StrPureBasisPair> blockFields = physicsBlock.getProvidedDOFs();
87  for (auto&& v : blockFields) {
88  if (v.first == field) {
89  basis = v.second;
90  break;
91  }
92  }
93  RCP<panzer::IntegrationRule> ir = physicsBlock.getIntegrationRules().at(cubatureDegree_);
94 
95  // build useful evaluator
98  field,
99  fieldComponent_,
100  point_,
101  *ir,
102  basis,
103  globalIndexer_,
104  scatterObj));
105 
106  this->template registerEvaluator<EvalT>(fm, eval);
107 
108  // require last field
109  fm.template requireField<EvalT>(*eval->evaluatedFields()[0]);
110  }
111 }
112 
113 template <typename EvalT,typename LO,typename GO>
116 {
117  if(PHX::typeAsString<EvalT>()==PHX::typeAsString<panzer::Traits::Residual>() ||
118  PHX::typeAsString<EvalT>()==PHX::typeAsString<panzer::Traits::Tangent>() ||
119  PHX::typeAsString<EvalT>()==PHX::typeAsString<panzer::Traits::Jacobian>()
120  )
121  return true;
122 
123  return false;
124 }
125 
126 }
127 
128 #endif
Object that contains information on the physics and discretization of a block of elements with the SA...
virtual Teuchos::RCP< ResponseBase > buildResponseObject(const std::string &responseName) const
TEUCHOS_DEPRECATED RCP< T > rcp(T *p, Dealloc_T dealloc, bool owns_mem)
const std::map< int, Teuchos::RCP< panzer::IntegrationRule > > & getIntegrationRules() const
Returns the unique set of point rules, key is the unique panzer::PointRule::name() ...
PHX::MDField< ScalarT, panzer::Cell, panzer::BASIS > field
A field to which we&#39;ll contribute, or in which we&#39;ll store, the result of computing this integral...
virtual void buildAndRegisterEvaluators(const std::string &responseName, PHX::FieldManager< panzer::Traits > &fm, const panzer::PhysicsBlock &physicsBlock, const Teuchos::ParameterList &user_data) const
const std::vector< StrPureBasisPair > & getProvidedDOFs() const