Panzer Version of the Day
Loading...
Searching...
No Matches
Panzer_GatherSolution_Epetra_Hessian.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_GatherSolution_Epetra_Hessian_hpp__
44#define __Panzer_GatherSolution_Epetra_Hessian_hpp__
45
46// Only do this if required by the user.
47#ifdef Panzer_BUILD_HESSIAN_SUPPORT
48
49// the includes for this file come in as a result of the includes in the main
50// Epetra gather solution file
51
52namespace panzer
53{
64 template<typename TRAITS, typename LO, typename GO>
66 :
67 public panzer::EvaluatorWithBaseImpl<TRAITS>,
68 public PHX::EvaluatorDerived<panzer::Traits::Hessian, TRAITS>,
70 {
71 public:
72
82 const Teuchos::RCP<const panzer::GlobalIndexer>& indexer)
83 :
84 globalIndexer_(indexer)
85 {
86 } // end of Constructor
87
100 const Teuchos::RCP<const panzer::GlobalIndexer>& indexer,
101 const Teuchos::ParameterList& p);
102
111 void
112 postRegistrationSetup(
113 typename TRAITS::SetupData d,
115
132 void
133 preEvaluate(
134 typename TRAITS::PreEvalData d);
135
145 void
146 evaluateFields(
147 typename TRAITS::EvalData d);
148
161 virtual Teuchos::RCP<CloneableEvaluator>
163 const Teuchos::ParameterList& pl) const
164 {
165 using panzer::Traits;
166 using Teuchos::rcp;
168 (globalIndexer_, pl));
169 } // end of clone()
170
171 private:
172
177
182
187 Teuchos::RCP<const panzer::GlobalIndexer> globalIndexer_;
188
192 std::vector<std::string> indexerNames_;
193
197 std::vector<int> fieldIds_;
198
202 std::vector< PHX::MDField<ScalarT, Cell, NODE>> gatherFields_;
203
207 std::vector< PHX::MDField<ScalarT, Cell, NODE>> sensFields_;
208
214
218 std::string globalDataKey_;
219
224
231
237
243
248
254
260
264 Teuchos::RCP<Epetra_Vector> x_;
265
270 Teuchos::RCP<panzer::EpetraVector_ReadOnly_GlobalEvaluationData>
272
277 Teuchos::RCP<panzer::EpetraVector_ReadOnly_GlobalEvaluationData>
279
284
285 }; // end of class GatherSolution_Epetra (Hessian Specialization)
286
287} // end of namespace panzer
288
289#endif // Panzer_BUILD_HESSIAN_SUPPORT
290
291#endif // __Panzer_GatherSolution_Epetra_Hessian_hpp__
Non-templated empty base class for template managers.
Wrapper to PHX::EvaluatorWithBaseImpl that implements Panzer-specific helpers.
bool secondSensitivitiesAvailable_
A flag indicating whether or not we're to be working with the second derivative sensitivities.
std::vector< PHX::MDField< ScalarT, Cell, NODE > > gatherFields_
The fields to be gathered.
virtual Teuchos::RCP< CloneableEvaluator > clone(const Teuchos::ParameterList &pl) const
Create a copy.
bool secondApplySensitivities_
Used by evaluateFields() to turn on/off the second derivative sensitivities.
bool firstApplySensitivities_
Used by evaluateFields() to turn on/off the first derivative sensitivities.
std::vector< PHX::MDField< ScalarT, Cell, NODE > > sensFields_
The sensitivity fields.
bool firstSensitivitiesAvailable_
A flag indicating whether or not we're to be working with the first derivative sensitivities.
bool useTimeDerivativeSolutionVector_
A flag indicating whether we're to be working with or .
Teuchos::RCP< const panzer::GlobalIndexer > globalIndexer_
Maps the local (field, element, basis) triplet to a global ID for scattering.
std::string sensitivities2ndPrefix_
The prefix for the field containing the second sensitivities.
GatherSolution_Epetra(const Teuchos::RCP< const panzer::GlobalIndexer > &indexer)
Constructor.
Teuchos::RCP< panzer::EpetraVector_ReadOnly_GlobalEvaluationData > xEvRoGed_
The GlobalEvaluationData containing both the owned and ghosted solution vectors.
std::vector< std::string > indexerNames_
A list of the names of the fields to be gathered.
Teuchos::RCP< panzer::EpetraVector_ReadOnly_GlobalEvaluationData > dxEvRoGed_
The GlobalEvaluationData containing both the owned and ghosted derivative vectors.
Gathers solution values from the Newton solution vector into the nodal fields of the field manager.