Tpetra parallel linear algebra  Version of the Day
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Tpetra_BlockVector_def.hpp
1 // @HEADER
2 // *****************************************************************************
3 // Tpetra: Templated Linear Algebra Services Package
4 //
5 // Copyright 2008 NTESS and the Tpetra contributors.
6 // SPDX-License-Identifier: BSD-3-Clause
7 // *****************************************************************************
8 // @HEADER
9 
10 #ifndef TPETRA_BLOCKVECTOR_DEF_HPP
11 #define TPETRA_BLOCKVECTOR_DEF_HPP
12 
13 namespace Tpetra {
14 
15  template<class Scalar, class LO, class GO, class Node>
16  BlockVector<Scalar, LO, GO, Node>::
17  BlockVector () :
18  base_type ()
19  {}
20 
21  template<class Scalar, class LO, class GO, class Node>
22  BlockVector<Scalar, LO, GO, Node>::
23  BlockVector (const BlockVector<Scalar, LO, GO, Node>& in,
24  const Teuchos::DataAccess copyOrView) :
25  base_type (in, copyOrView)
26  {}
27 
28  template<class Scalar, class LO, class GO, class Node>
29  BlockVector<Scalar, LO, GO, Node>::
30  BlockVector (const map_type& meshMap, const LO blockSize) :
31  base_type (meshMap, blockSize, 1)
32  {}
33 
34  template<class Scalar, class LO, class GO, class Node>
35  BlockVector<Scalar, LO, GO, Node>::
36  BlockVector (const map_type& meshMap,
37  const map_type& pointMap,
38  const LO blockSize) :
39  base_type (meshMap, pointMap, blockSize, 1)
40  {}
41 
42  template<class Scalar, class LO, class GO, class Node>
43  BlockVector<Scalar, LO, GO, Node>::
44  BlockVector (const mv_type& X_mv,
45  const map_type& meshMap,
46  const LO blockSize) :
47  base_type (X_mv, meshMap, blockSize)
48  {
49  TEUCHOS_TEST_FOR_EXCEPTION(
50  X_mv.getNumVectors () != 1, std::invalid_argument,
51  "Tpetra::BlockVector: Input MultiVector has "
52  << X_mv.getNumVectors () << " != 1 columns.");
53  }
54 
55  template<class Scalar, class LO, class GO, class Node>
56  BlockVector<Scalar, LO, GO, Node>::
57  BlockVector (const vec_type& X_vec,
58  const map_type& meshMap,
59  const LO blockSize) :
60  base_type (X_vec, meshMap, blockSize)
61  {}
62 
63  template<class Scalar, class LO, class GO, class Node>
64  BlockVector<Scalar, LO, GO, Node>::
65  BlockVector (const BlockVector<Scalar, LO, GO, Node>& X,
66  const map_type& newMeshMap,
67  const map_type& newPointMap,
68  const size_t offset) :
69  base_type (X, newMeshMap, newPointMap, offset)
70  {}
71 
72  template<class Scalar, class LO, class GO, class Node>
73  BlockVector<Scalar, LO, GO, Node>::
74  BlockVector (const BlockVector<Scalar, LO, GO, Node>& X,
75  const map_type& newMeshMap,
76  const size_t offset) :
77  base_type (X, newMeshMap, offset)
78  {}
79 
80  template<class Scalar, class LO, class GO, class Node>
81  typename BlockVector<Scalar, LO, GO, Node>::vec_type
82  BlockVector<Scalar, LO, GO, Node>::getVectorView () {
83  Teuchos::RCP<vec_type> vPtr = this->mv_.getVectorNonConst (0);
84  return *vPtr; // shallow copy
85  }
86 
87  template<class Scalar, class LO, class GO, class Node>
88  bool
89  BlockVector<Scalar, LO, GO, Node>::
90  replaceLocalValues (const LO localRowIndex, const Scalar vals[]) {
91  return ((base_type*) this)->replaceLocalValues (localRowIndex, 0, vals);
92  }
93 
94  template<class Scalar, class LO, class GO, class Node>
95  bool
96  BlockVector<Scalar, LO, GO, Node>::
97  replaceGlobalValues (const GO globalRowIndex, const Scalar vals[]) {
98  return ((base_type*) this)->replaceGlobalValues (globalRowIndex, 0, vals);
99  }
100 
101  template<class Scalar, class LO, class GO, class Node>
102  bool
103  BlockVector<Scalar, LO, GO, Node>::
104  sumIntoLocalValues (const LO localRowIndex, const Scalar vals[]) {
105  return ((base_type*) this)->sumIntoLocalValues (localRowIndex, 0, vals);
106  }
107 
108  template<class Scalar, class LO, class GO, class Node>
109  bool
110  BlockVector<Scalar, LO, GO, Node>::
111  sumIntoGlobalValues (const GO globalRowIndex, const Scalar vals[]) {
112  return ((base_type*) this)->sumIntoLocalValues (globalRowIndex, 0, vals);
113  }
114 
115  template<class Scalar, class LO, class GO, class Node>
116  typename BlockVector<Scalar, LO, GO, Node>::const_little_host_vec_type
117  BlockVector<Scalar, LO, GO, Node>::
118  getLocalBlockHost (const LO localRowIndex, Access::ReadOnlyStruct) const
119  {
120  return ((const base_type*) this)->getLocalBlockHost(localRowIndex, 0,
121  Access::ReadOnly);
122  }
123 
124  template<class Scalar, class LO, class GO, class Node>
125  typename BlockVector<Scalar, LO, GO, Node>::little_host_vec_type
126  BlockVector<Scalar, LO, GO, Node>::
127  getLocalBlockHost (const LO localRowIndex, Access::ReadWriteStruct)
128  {
129  return ((base_type*) this)->getLocalBlockHost(localRowIndex, 0,
130  Access::ReadWrite);
131  }
132 
133  template<class Scalar, class LO, class GO, class Node>
134  typename BlockVector<Scalar, LO, GO, Node>::little_host_vec_type
135  BlockVector<Scalar, LO, GO, Node>::
136  getLocalBlockHost (const LO localRowIndex, Access::OverwriteAllStruct)
137  {
138  return ((base_type*) this)->getLocalBlockHost(localRowIndex, 0,
139  Access::OverwriteAll);
140  }
141 
142 } // namespace Tpetra
143 
144 //
145 // Explicit instantiation macro
146 //
147 // Must be expanded from within the Tpetra namespace!
148 //
149 #define TPETRA_BLOCKVECTOR_INSTANT(S,LO,GO,NODE) \
150  template class BlockVector< S, LO, GO, NODE >;
151 
152 #endif // TPETRA_BLOCKVECTOR_DEF_HPP
Tpetra::BlockVector::sumIntoGlobalValues
bool sumIntoGlobalValues(const GO globalRowIndex, const Scalar vals[])
Sum into all values at the given mesh point, using a global index.
Definition: Tpetra_BlockVector_def.hpp:111
Tpetra::MultiVector::getNumVectors
size_t getNumVectors() const
Number of columns in the multivector.
Definition: Tpetra_MultiVector_def.hpp:2206
Tpetra::BlockVector::sumIntoLocalValues
bool sumIntoLocalValues(const LO localRowIndex, const Scalar vals[])
Sum into all values at the given mesh point, using a local index.
Definition: Tpetra_BlockVector_def.hpp:104
Tpetra::BlockMultiVector
MultiVector for multiple degrees of freedom per mesh point.
Definition: Tpetra_BlockMultiVector_decl.hpp:110
Tpetra::BlockVector::getVectorView
vec_type getVectorView()
Get a Tpetra::Vector that views this BlockVector&#39;s data.
Definition: Tpetra_BlockVector_def.hpp:82
Tpetra::MultiVector::getVectorNonConst
Teuchos::RCP< Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > getVectorNonConst(const size_t j)
Return a Vector which is a nonconst view of column j.
Definition: Tpetra_MultiVector_def.hpp:3749
Tpetra::BlockVector::getLocalBlockHost
const_little_host_vec_type getLocalBlockHost(const LO localRowIndex, Access::ReadOnlyStruct) const
Get a view of the degrees of freedom at the given mesh point, using a local index.
Definition: Tpetra_BlockVector_def.hpp:118
Tpetra::BlockVector::replaceLocalValues
bool replaceLocalValues(const LO localRowIndex, const Scalar vals[])
Replace all values at the given mesh point, using a local index.
Definition: Tpetra_BlockVector_def.hpp:90
Tpetra::Vector
A distributed dense vector.
Definition: Tpetra_Vector_decl.hpp:44
Tpetra::BlockVector::replaceGlobalValues
bool replaceGlobalValues(const GO globalRowIndex, const Scalar vals[])
Replace all values at the given mesh point, using a global index.
Definition: Tpetra_BlockVector_def.hpp:97
Tpetra::BlockVector::BlockVector
BlockVector()
Default constructor.
Definition: Tpetra_BlockVector_def.hpp:17
Tpetra::BlockVector
Vector for multiple degrees of freedom per mesh point.
Definition: Tpetra_BlockVector_decl.hpp:48
Generated on Thu Nov 21 2024 09:19:24 for Tpetra parallel linear algebra by   doxygen 1.8.5