Tpetra parallel linear algebra Version of the Day
Loading...
Searching...
No Matches
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
13namespace Tpetra {
14
15 template<class Scalar, class LO, class GO, class Node>
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::BlockMultiVector::mv_
mv_type mv_
The Tpetra::MultiVector used to represent the data.
Definition Tpetra_BlockMultiVector_decl.hpp:606
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::mv_type
Tpetra::MultiVector< Scalar, LO, GO, Node > mv_type
The specialization of Tpetra::MultiVector that this class uses.
Definition Tpetra_BlockVector_decl.hpp:73
Tpetra::BlockVector::map_type
Tpetra::Map< LO, GO, Node > map_type
The specialization of Tpetra::Map that this class uses.
Definition Tpetra_BlockVector_decl.hpp:71
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::BlockVector::vec_type
Tpetra::Vector< Scalar, LO, GO, Node > vec_type
The specialization of Tpetra::Vector that this class uses.
Definition Tpetra_BlockVector_decl.hpp:75
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::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::BlockVector::getVectorView
vec_type getVectorView()
Get a Tpetra::Vector that views this BlockVector's data.
Definition Tpetra_BlockVector_def.hpp:82
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::MultiVector::getNumVectors
size_t getNumVectors() const
Number of columns in the multivector.
Definition Tpetra_MultiVector_def.hpp:2206
Tpetra
Namespace Tpetra contains the class and methods constituting the Tpetra library.
Generated on for Tpetra parallel linear algebra by doxygen 1.15.0