Panzer  Version of the Day
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
Panzer_ResponseEvaluatorFactory_ExtremeValue_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_ExtremeValue_impl_hpp__
44 #define __Panzer_ResponseEvaluatorFactory_ExtremeValue_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_ExtremeValue<EvalT>(responseName,comm_,useMax_,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 
79  // build integration evaluator (integrate over element)
80  if(requiresCellExtreme_) {
81  std::string field = (quadPointField_=="" ? responseName : quadPointField_);
82 
83  // build integration rule to use in cell integral
84  RCP<IntegrationRule> ir = rcp(new IntegrationRule(cubatureDegree_,physicsBlock.cellData()));
85 
87  // add prefix_ to help identify
88  pl.set("Extreme Name",prefix_+field);
89  pl.set("Field Name",field);
90  pl.set("IR",ir);
91  pl.set("Use Max",useMax_);
92 
95 
96  this->template registerEvaluator<EvalT>(fm, eval);
97  }
98 
99 
100  // build scatter evaluator
101  {
103  (globalIndexer_!=Teuchos::null) ? Teuchos::rcp(new ExtremeValueScatter<LO,GO>(globalIndexer_)) : Teuchos::null;
104  std::string field = (quadPointField_=="" ? responseName : quadPointField_);
105  field = prefix_+field; // add prefix to help identify
106 
107  // build useful evaluator
110  responseName,
111  physicsBlock.cellData(),
112  useMax_,
113  scatterObj));
114 
115  this->template registerEvaluator<EvalT>(fm, eval);
116 
117  // require last field
118  fm.template requireField<EvalT>(*eval->evaluatedFields()[0]);
119  }
120 }
121 
122 template <typename EvalT,typename LO,typename GO>
125 {
126  if( PHX::typeAsString<EvalT>()==PHX::typeAsString<panzer::Traits::Residual>() ||
127  PHX::typeAsString<EvalT>()==PHX::typeAsString<panzer::Traits::Tangent>()
128  )
129  return true;
130 
131  if(PHX::typeAsString<EvalT>()==PHX::typeAsString<panzer::Traits::Jacobian>())
132  return false;
133 
134  return false;
135 }
136 
137 }
138 
139 #endif
Object that contains information on the physics and discretization of a block of elements with the SA...
const panzer::CellData & cellData() const
ParameterList & set(std::string const &name, T const &value, std::string const &docString="", RCP< const ParameterEntryValidator > const &validator=null)
TEUCHOS_DEPRECATED RCP< T > rcp(T *p, Dealloc_T dealloc, bool owns_mem)
virtual Teuchos::RCP< ResponseBase > buildResponseObject(const std::string &responseName) const
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