Panzer  Version of the Day
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
Panzer_TpetraVector_ReadOnly_GlobalEvaluationData.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_TpetraVector_ReadOnly_GlobalEvaluationData_hpp__
44 #define __Panzer_TpetraVector_ReadOnly_GlobalEvaluationData_hpp__
45 
46 #include "Tpetra_Import.hpp"
47 #include "Tpetra_Vector.hpp"
48 #include "Tpetra_Map.hpp"
49 
50 #include "Teuchos_RCP.hpp"
51 
52 #include "Thyra_VectorSpaceBase.hpp"
53 #include "Thyra_VectorBase.hpp"
54 
55 #include "Panzer_NodeType.hpp"
57 
58 namespace panzer {
59 
63 template <typename ScalarT,typename LocalOrdinalT,typename GlobalOrdinalT,
64  typename NodeT=panzer::TpetraNodeType>
66 public:
67  typedef Tpetra::Vector<ScalarT,LocalOrdinalT,GlobalOrdinalT,NodeT> VectorType;
68  typedef Tpetra::Map<LocalOrdinalT,GlobalOrdinalT,NodeT> MapType;
69  typedef Tpetra::Import<LocalOrdinalT,GlobalOrdinalT,NodeT> ImportType;
70 
73  : isInitialized_(false) { }
74 
76  : isInitialized_(false) { initialize(src.importer_, src.ghostedMap_, src.ownedMap_); }
77 
87  const Teuchos::RCP<const MapType>& ghostedMap,
88  const Teuchos::RCP<const MapType>& ownedMap)
89  : isInitialized_(false) { initialize(importer, ghostedMap, ownedMap); }
90 
98  void useConstantValues(const std::vector<GlobalOrdinalT> & indices,double value);
99 
108  void initialize(const Teuchos::RCP<const ImportType>& importer,
109  const Teuchos::RCP<const MapType>& ghostedMap,
110  const Teuchos::RCP<const MapType>& ownedMap);
111 
115  virtual void globalToGhost(int mem);
116 
118  virtual void initializeData();
119 
121  virtual void ghostToGlobal(int /* mem */) {}
122 
124  virtual bool requiresDirichletAdjustment() const { return false; }
125 
127  void setOwnedVector_Tpetra(const Teuchos::RCP<const VectorType>& ownedVector);
128 
131 
134 
136  void setOwnedVector(const Teuchos::RCP<const Thyra::VectorBase<double> >& ownedVector);
137 
140 
143 
145  bool isInitialized() const { return isInitialized_; }
146 
148  void print(std::ostream & os) const;
149 
150 private:
152 
155 
158 
162 
163  typedef std::pair<std::vector<GlobalOrdinalT>,double> FilteredGlobalPair;
164  typedef std::pair<std::vector<LocalOrdinalT>,double> FilteredLocalPair;
165  std::vector<FilteredGlobalPair> globalFilteredPairs_;
166  std::vector<FilteredLocalPair> filteredPairs_;
167 };
168 
169 }
170 
171 #endif
virtual void ghostToGlobal(int)
For this class this method does nothing.
Tpetra::Vector< ScalarT, LocalOrdinalT, GlobalOrdinalT, NodeT > VectorType
void setOwnedVector(const Teuchos::RCP< const Thyra::VectorBase< double > > &ownedVector)
Set the owned vector (Thyra version)
virtual bool requiresDirichletAdjustment() const
Nothing to do (its read only)
TpetraVector_ReadOnly_GlobalEvaluationData(const TpetraVector_ReadOnly_GlobalEvaluationData &src)
void useConstantValues(const std::vector< GlobalOrdinalT > &indices, double value)
Teuchos::RCP< const Thyra::VectorBase< double > > getOwnedVector() const
Get the owned vector (Thyra version)
Kokkos::Compat::KokkosDeviceWrapperNode< PHX::Device > TpetraNodeType
TpetraVector_ReadOnly_GlobalEvaluationData(const Teuchos::RCP< const ImportType > &importer, const Teuchos::RCP< const MapType > &ghostedMap, const Teuchos::RCP< const MapType > &ownedMap)
void setOwnedVector_Tpetra(const Teuchos::RCP< const VectorType > &ownedVector)
Set the owned vector (Tpetra version)
Teuchos::RCP< const VectorType > getOwnedVector_Tpetra() const
Get the owned vector (Tpetra version)
Teuchos::RCP< VectorType > getGhostedVector_Tpetra() const
Get the ghosted vector (Tpetra version)
Teuchos::RCP< Thyra::VectorBase< double > > getGhostedVector() const
Get the ghosted vector (Thyra version)
void initialize(const Teuchos::RCP< const ImportType > &importer, const Teuchos::RCP< const MapType > &ghostedMap, const Teuchos::RCP< const MapType > &ownedMap)