Xpetra Version of the Day
Xpetra_TpetraMultiVector_def.hpp
Go to the documentation of this file.
1// @HEADER
2//
3// ***********************************************************************
4//
5// Xpetra: A linear algebra interface package
6// Copyright 2012 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
39// Jonathan Hu (jhu@sandia.gov)
40// Andrey Prokopenko (aprokop@sandia.gov)
41// Ray Tuminaro (rstumin@sandia.gov)
42//
43// ***********************************************************************
44//
45// @HEADER
46#ifndef XPETRA_TPETRAMULTIVECTOR_DEF_HPP
47#define XPETRA_TPETRAMULTIVECTOR_DEF_HPP
48#include "Xpetra_TpetraConfigDefs.hpp"
49
50#include "Xpetra_TpetraMap.hpp" //TMP
51#include "Xpetra_Utils.hpp"
52#include "Xpetra_TpetraImport.hpp"
53#include "Xpetra_TpetraExport.hpp"
54
55#include "Xpetra_TpetraMultiVector_decl.hpp"
56#include "Tpetra_MultiVector.hpp"
57#include "Tpetra_Vector.hpp"
58
59namespace Xpetra {
60
61
63 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
64 TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::TpetraMultiVector(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map, size_t NumVectors, bool zeroOut)
65 : vec_(Teuchos::rcp(new Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node >(toTpetra(map), NumVectors, zeroOut))) {
66 // TAW 1/30/2016: even though Tpetra allows numVecs == 0, Epetra does not. Introduce exception to keep behavior of Epetra and Tpetra consistent.
67 TEUCHOS_TEST_FOR_EXCEPTION(NumVectors < 1, std::invalid_argument, "Xpetra::TpetraMultiVector(map,numVecs,zeroOut): numVecs = " << NumVectors << " < 1.");
68 }
69
71 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
72 TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
73 TpetraMultiVector(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Teuchos::DataAccess copyOrView)
74 : vec_(Teuchos::rcp(new Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node >(toTpetra(source), copyOrView))) { }
75
77 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
78 TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
79 TpetraMultiVector(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map, const Teuchos::ArrayView< const Scalar > &A, size_t LDA, size_t NumVectors)
80 : vec_(Teuchos::rcp(new Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node >(toTpetra(map), A, LDA, NumVectors))) {
81 // TAW 1/30/2016: even though Tpetra allows numVecs == 0, Epetra does not. Introduce exception to keep behavior of Epetra and Tpetra consistent.
82 TEUCHOS_TEST_FOR_EXCEPTION(NumVectors < 1, std::invalid_argument, "Xpetra::TpetraMultiVector(map,A,LDA,numVecs): numVecs = " << NumVectors << " < 1.");
83 }
84
86 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
87 TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
88 TpetraMultiVector(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map, const Teuchos::ArrayView< const Teuchos::ArrayView< const Scalar > > &ArrayOfPtrs, size_t NumVectors)
89 : vec_(Teuchos::rcp(new Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node >(toTpetra(map), ArrayOfPtrs, NumVectors))) {
90 // TAW 1/30/2016: even though Tpetra allows numVecs == 0, Epetra does not. Introduce exception to keep behavior of Epetra and Tpetra consistent.
91 TEUCHOS_TEST_FOR_EXCEPTION(NumVectors < 1, std::invalid_argument, "Xpetra::TpetraMultiVector(map,ArrayOfPtrs,numVecs): numVecs = " << NumVectors << " < 1.");
92 }
93
94
96 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
97 TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
98 ~TpetraMultiVector() { }
99
101 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
102 void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
103 replaceGlobalValue(GlobalOrdinal globalRow, size_t vectorIndex, const Scalar &value) { XPETRA_MONITOR("TpetraMultiVector::replaceGlobalValue"); vec_->replaceGlobalValue(globalRow, vectorIndex, value); }
104
106 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
107 void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
108 sumIntoGlobalValue(GlobalOrdinal globalRow, size_t vectorIndex, const Scalar &value) { XPETRA_MONITOR("TpetraMultiVector::sumIntoGlobalValue"); vec_->sumIntoGlobalValue(globalRow, vectorIndex, value); }
109
111 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
112 void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
113 replaceLocalValue(LocalOrdinal myRow, size_t vectorIndex, const Scalar &value) { XPETRA_MONITOR("TpetraMultiVector::replaceLocalValue"); vec_->replaceLocalValue(myRow, vectorIndex, value); }
114
116 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
117 void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
118 sumIntoLocalValue(LocalOrdinal myRow, size_t vectorIndex, const Scalar &value) { XPETRA_MONITOR("TpetraMultiVector::sumIntoLocalValue"); vec_->sumIntoLocalValue(myRow, vectorIndex, value); }
119
121 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
122 void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
123 putScalar(const Scalar &value) { XPETRA_MONITOR("TpetraMultiVector::putScalar"); vec_->putScalar(value); }
124
126 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
127 void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
128 reduce() { XPETRA_MONITOR("TpetraMultiVector::reduce"); vec_->reduce(); }
129
130
131 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
132 Teuchos::RCP< const Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > >
133 TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
134 getVector(size_t j) const { XPETRA_MONITOR("TpetraMultiVector::getVector"); return toXpetra(vec_->getVector(j)); }
135
137 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
138 Teuchos::RCP< Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > >
139 TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
140 getVectorNonConst(size_t j) { XPETRA_MONITOR("TpetraMultiVector::getVectorNonConst"); return toXpetra(vec_->getVectorNonConst(j)); }
141
143 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
144 Teuchos::ArrayRCP< const Scalar >
145 TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
146 getData(size_t j) const { XPETRA_MONITOR("TpetraMultiVector::getData"); return vec_->getData(j); }
147
149 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
150 Teuchos::ArrayRCP< Scalar >
151 TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
152 getDataNonConst(size_t j) { XPETRA_MONITOR("TpetraMultiVector::getDataNonConst"); return vec_->getDataNonConst(j); }
153
155 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
156 void
157 TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
158 get1dCopy(Teuchos::ArrayView< Scalar > A, size_t LDA) const { XPETRA_MONITOR("TpetraMultiVector::get1dCopy"); vec_->get1dCopy(A, LDA); }
159
161 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
162 void
163 TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
164 get2dCopy(Teuchos::ArrayView< const Teuchos::ArrayView< Scalar > > ArrayOfPtrs) const { XPETRA_MONITOR("TpetraMultiVector::get2dCopy"); vec_->get2dCopy(ArrayOfPtrs); }
165
167 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
168 Teuchos::ArrayRCP< const Scalar >
169 TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
170 get1dView() const { XPETRA_MONITOR("TpetraMultiVector::get1dView"); return vec_->get1dView(); }
171
173 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
174 Teuchos::ArrayRCP< Teuchos::ArrayRCP< const Scalar > >
175 TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
176 get2dView() const { XPETRA_MONITOR("TpetraMultiVector::get2dView"); return vec_->get2dView(); }
177
179 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
180 Teuchos::ArrayRCP< Scalar >
181 TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
182 get1dViewNonConst() { XPETRA_MONITOR("TpetraMultiVector::get1dViewNonConst"); return vec_->get1dViewNonConst(); }
183
185 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
186 Teuchos::ArrayRCP< Teuchos::ArrayRCP< Scalar > >
187 TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
188 get2dViewNonConst() { XPETRA_MONITOR("TpetraMultiVector::get2dViewNonConst"); return vec_->get2dViewNonConst(); }
189
191 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
192 void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
193 dot(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const Teuchos::ArrayView< Scalar > &dots) const { XPETRA_MONITOR("TpetraMultiVector::dot"); vec_->dot(toTpetra(A), dots); }
194
196 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
197 void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
198 abs(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A) { XPETRA_MONITOR("TpetraMultiVector::abs"); vec_->abs(toTpetra(A)); }
199
201 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
202 void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
203 reciprocal(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A) { XPETRA_MONITOR("TpetraMultiVector::reciprocal"); vec_->reciprocal(toTpetra(A)); }
204
206 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
207 void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
208 scale(const Scalar &alpha) { XPETRA_MONITOR("TpetraMultiVector::scale"); vec_->scale(alpha); }
209
211 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
212 void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
213 scale(Teuchos::ArrayView< const Scalar > alpha) { XPETRA_MONITOR("TpetraMultiVector::scale"); vec_->scale(alpha); }
214
216 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
217 void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
218 scale(const Scalar &alpha, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A) { XPETRA_MONITOR("TpetraMultiVector::scale"); vec_->scale(alpha, toTpetra(A)); }
219
221 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
222 void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
223 update(const Scalar &alpha, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const Scalar &beta) { XPETRA_MONITOR("TpetraMultiVector::update"); vec_->update(alpha, toTpetra(A), beta); }
224
226 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
227 void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
228 update(const Scalar &alpha, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const Scalar &beta, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &B, const Scalar &gamma) { XPETRA_MONITOR("TpetraMultiVector::update"); vec_->update(alpha, toTpetra(A), beta, toTpetra(B), gamma); }
229
231 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
232 void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
233 norm1(const Teuchos::ArrayView< typename Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const { XPETRA_MONITOR("TpetraMultiVector::norm1"); vec_->norm1(norms); }
234
236 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
237 void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
238 norm2(const Teuchos::ArrayView< typename Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const { XPETRA_MONITOR("TpetraMultiVector::norm2"); vec_->norm2(norms); }
239
241 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
242 void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
243 normInf(const Teuchos::ArrayView< typename Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const { XPETRA_MONITOR("TpetraMultiVector::normInf"); vec_->normInf(norms); }
244
246 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
247 void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
248 meanValue(const Teuchos::ArrayView< Scalar > &means) const { XPETRA_MONITOR("TpetraMultiVector::meanValue"); vec_->meanValue(means); }
249
251 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
252 void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
253 multiply(Teuchos::ETransp transA, Teuchos::ETransp transB, const Scalar &alpha, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &B, const Scalar &beta) { XPETRA_MONITOR("TpetraMultiVector::multiply"); vec_->multiply(transA, transB, alpha, toTpetra(A), toTpetra(B), beta); }
254
255
257 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
258 size_t TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
259 getNumVectors() const { XPETRA_MONITOR("TpetraMultiVector::getNumVectors"); return vec_->getNumVectors(); }
260
262 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
263 size_t TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
264 getLocalLength() const { XPETRA_MONITOR("TpetraMultiVector::getLocalLength"); return vec_->getLocalLength(); }
265
267 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
268 global_size_t TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
269 getGlobalLength() const { XPETRA_MONITOR("TpetraMultiVector::getGlobalLength"); return vec_->getGlobalLength(); }
270
271 // \brief Checks to see if the local length, number of vectors and size of Scalar type match
272 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
273 bool TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
274 isSameSize(const MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> & vec) const { XPETRA_MONITOR("TpetraMultiVector::isSameSize"); return vec_->isSameSize(toTpetra(vec));}
275
277 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
278 std::string TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
279 description() const { XPETRA_MONITOR("TpetraMultiVector::description"); return vec_->description(); }
280
282 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
283 void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
284 describe(Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel) const { XPETRA_MONITOR("TpetraMultiVector::describe"); vec_->describe(out, verbLevel); }
285
287 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
288 void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
289 randomize(bool bUseXpetraImplementation) {
290 XPETRA_MONITOR("TpetraMultiVector::randomize");
291
292 if(bUseXpetraImplementation)
293 MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node >::Xpetra_randomize();
294 else
295 vec_->randomize();
296 }
297
298 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
299 Teuchos::RCP< const Map<LocalOrdinal,GlobalOrdinal,Node> >
300 TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
301 getMap() const { XPETRA_MONITOR("TpetraMultiVector::getMap"); return toXpetra(vec_->getMap()); }
302
303 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
304 void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
305 doImport(const DistObject< Scalar, LocalOrdinal,GlobalOrdinal,Node> &source, const Import<LocalOrdinal,GlobalOrdinal,Node> &importer, CombineMode CM) {
306 XPETRA_MONITOR("TpetraMultiVector::doImport");
307
308 XPETRA_DYNAMIC_CAST(const TpetraMultiVectorClass, source, tSource, "Xpetra::TpetraMultiVector::doImport only accept Xpetra::TpetraMultiVector as input arguments."); //TODO: remove and use toTpetra()
309 RCP< const Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal,Node> > v = tSource.getTpetra_MultiVector();
310 this->getTpetra_MultiVector()->doImport(*v, toTpetra(importer), toTpetra(CM));
311 }
312
313 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
314 void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
315 beginImport(const DistObject< Scalar, LocalOrdinal,GlobalOrdinal,Node> &source, const Import<LocalOrdinal,GlobalOrdinal,Node> &importer, CombineMode CM) {
316 XPETRA_MONITOR("TpetraMultiVector::beginImport");
317
318 XPETRA_DYNAMIC_CAST(const TpetraMultiVectorClass, source, tSource, "Xpetra::TpetraMultiVector::doImport only accept Xpetra::TpetraMultiVector as input arguments."); //TODO: remove and use toTpetra()
319 RCP< const Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal,Node> > v = tSource.getTpetra_MultiVector();
320 this->getTpetra_MultiVector()->beginImport(*v, toTpetra(importer), toTpetra(CM));
321 }
322
323 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
324 void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
325 endImport(const DistObject< Scalar, LocalOrdinal,GlobalOrdinal,Node> &source, const Import<LocalOrdinal,GlobalOrdinal,Node> &importer, CombineMode CM) {
326 XPETRA_MONITOR("TpetraMultiVector::endImport");
327
328 XPETRA_DYNAMIC_CAST(const TpetraMultiVectorClass, source, tSource, "Xpetra::TpetraMultiVector::doImport only accept Xpetra::TpetraMultiVector as input arguments."); //TODO: remove and use toTpetra()
329 RCP< const Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal,Node> > v = tSource.getTpetra_MultiVector();
330 this->getTpetra_MultiVector()->endImport(*v, toTpetra(importer), toTpetra(CM));
331 }
332
333
334 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
335 void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
336 doExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Import<LocalOrdinal,GlobalOrdinal,Node>& importer, CombineMode CM) {
337 XPETRA_MONITOR("TpetraMultiVector::doExport");
338
339 XPETRA_DYNAMIC_CAST(const TpetraMultiVectorClass, dest, tDest, "Xpetra::TpetraMultiVector::doImport only accept Xpetra::TpetraMultiVector as input arguments."); //TODO: remove and use toTpetra()
340 RCP< const Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal,Node> > v = tDest.getTpetra_MultiVector();
341 this->getTpetra_MultiVector()->doExport(*v, toTpetra(importer), toTpetra(CM));
342
343 }
344
345 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
346 void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
347 beginExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Import<LocalOrdinal,GlobalOrdinal,Node>& importer, CombineMode CM) {
348 XPETRA_MONITOR("TpetraMultiVector::beginExport");
349
350 XPETRA_DYNAMIC_CAST(const TpetraMultiVectorClass, dest, tDest, "Xpetra::TpetraMultiVector::doImport only accept Xpetra::TpetraMultiVector as input arguments."); //TODO: remove and use toTpetra()
351 RCP< const Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal,Node> > v = tDest.getTpetra_MultiVector();
352 this->getTpetra_MultiVector()->beginExport(*v, toTpetra(importer), toTpetra(CM));
353
354 }
355
356 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
357 void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
358 endExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Import<LocalOrdinal,GlobalOrdinal,Node>& importer, CombineMode CM) {
359 XPETRA_MONITOR("TpetraMultiVector::endExport");
360
361 XPETRA_DYNAMIC_CAST(const TpetraMultiVectorClass, dest, tDest, "Xpetra::TpetraMultiVector::doImport only accept Xpetra::TpetraMultiVector as input arguments."); //TODO: remove and use toTpetra()
362 RCP< const Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal,Node> > v = tDest.getTpetra_MultiVector();
363 this->getTpetra_MultiVector()->endExport(*v, toTpetra(importer), toTpetra(CM));
364
365 }
366
367 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
368 void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
369 doImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Export<LocalOrdinal,GlobalOrdinal,Node>& exporter, CombineMode CM) {
370 XPETRA_MONITOR("TpetraMultiVector::doImport");
371
372 XPETRA_DYNAMIC_CAST(const TpetraMultiVectorClass, source, tSource, "Xpetra::TpetraMultiVector::doImport only accept Xpetra::TpetraMultiVector as input arguments."); //TODO: remove and use toTpetra()
373 RCP< const Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal,Node> > v = tSource.getTpetra_MultiVector();
374 this->getTpetra_MultiVector()->doImport(*v, toTpetra(exporter), toTpetra(CM));
375
376 }
377
378 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
379 void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
380 beginImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Export<LocalOrdinal,GlobalOrdinal,Node>& exporter, CombineMode CM) {
381 XPETRA_MONITOR("TpetraMultiVector::beginImport");
382
383 XPETRA_DYNAMIC_CAST(const TpetraMultiVectorClass, source, tSource, "Xpetra::TpetraMultiVector::doImport only accept Xpetra::TpetraMultiVector as input arguments."); //TODO: remove and use toTpetra()
384 RCP< const Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal,Node> > v = tSource.getTpetra_MultiVector();
385 this->getTpetra_MultiVector()->beginImport(*v, toTpetra(exporter), toTpetra(CM));
386
387 }
388
389
390 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
391 void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
392 endImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Export<LocalOrdinal,GlobalOrdinal,Node>& exporter, CombineMode CM) {
393 XPETRA_MONITOR("TpetraMultiVector::endImport");
394
395 XPETRA_DYNAMIC_CAST(const TpetraMultiVectorClass, source, tSource, "Xpetra::TpetraMultiVector::doImport only accept Xpetra::TpetraMultiVector as input arguments."); //TODO: remove and use toTpetra()
396 RCP< const Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal,Node> > v = tSource.getTpetra_MultiVector();
397 this->getTpetra_MultiVector()->endImport(*v, toTpetra(exporter), toTpetra(CM));
398
399 }
400
401 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
402 void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
403 doExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Export<LocalOrdinal,GlobalOrdinal,Node>& exporter, CombineMode CM) {
404 XPETRA_MONITOR("TpetraMultiVector::doExport");
405
406 XPETRA_DYNAMIC_CAST(const TpetraMultiVectorClass, dest, tDest, "Xpetra::TpetraMultiVector::doImport only accept Xpetra::TpetraMultiVector as input arguments."); //TODO: remove and use toTpetra()
407 RCP< const Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal,Node> > v = tDest.getTpetra_MultiVector();
408 this->getTpetra_MultiVector()->doExport(*v, toTpetra(exporter), toTpetra(CM));
409
410 }
411
412 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
413 void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
414 beginExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Export<LocalOrdinal,GlobalOrdinal,Node>& exporter, CombineMode CM) {
415 XPETRA_MONITOR("TpetraMultiVector::beginExport");
416
417 XPETRA_DYNAMIC_CAST(const TpetraMultiVectorClass, dest, tDest, "Xpetra::TpetraMultiVector::doImport only accept Xpetra::TpetraMultiVector as input arguments."); //TODO: remove and use toTpetra()
418 RCP< const Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal,Node> > v = tDest.getTpetra_MultiVector();
419 this->getTpetra_MultiVector()->beginExport(*v, toTpetra(exporter), toTpetra(CM));
420
421 }
422
423 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
424 void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
425 endExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Export<LocalOrdinal,GlobalOrdinal,Node>& exporter, CombineMode CM) {
426 XPETRA_MONITOR("TpetraMultiVector::endExport");
427
428 XPETRA_DYNAMIC_CAST(const TpetraMultiVectorClass, dest, tDest, "Xpetra::TpetraMultiVector::doImport only accept Xpetra::TpetraMultiVector as input arguments."); //TODO: remove and use toTpetra()
429 RCP< const Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal,Node> > v = tDest.getTpetra_MultiVector();
430 this->getTpetra_MultiVector()->endExport(*v, toTpetra(exporter), toTpetra(CM));
431
432 }
433
434 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
435 void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
436 replaceMap(const RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> >& map) {
437 XPETRA_MONITOR("TpetraMultiVector::replaceMap");
438 this->getTpetra_MultiVector()->replaceMap(toTpetra(map));
439 }
440
442 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
443 TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
444 TpetraMultiVector(const Teuchos::RCP<Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node> > &vec) : vec_(vec) { } //TODO removed const
445
447 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
448 RCP< Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node> >
449 TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
450 getTpetra_MultiVector() const { return vec_; }
451
453 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
454 void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
455 setSeed(unsigned int seed) { XPETRA_MONITOR("TpetraMultiVector::seedrandom"); Teuchos::ScalarTraits< Scalar >::seedrandom(seed); }
456
457
458#ifdef HAVE_XPETRA_KOKKOS_REFACTOR
459
460 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
461 typename TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::dual_view_type::t_host_const_um
462 TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
463 getHostLocalView (Access::ReadOnlyStruct) const {
464 return subview(vec_->getLocalViewHost(Tpetra::Access::ReadOnly),Kokkos::ALL(), Kokkos::ALL());
465
466 }
467
468 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
469 typename TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::dual_view_type::t_dev_const_um
470 TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
471 getDeviceLocalView(Access::ReadOnlyStruct) const {
472 return subview(vec_->getLocalViewDevice(Tpetra::Access::ReadOnly),Kokkos::ALL(), Kokkos::ALL());
473 }
474
475 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
476 typename TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::dual_view_type::t_host_um
477 TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
478 getHostLocalView (Access::OverwriteAllStruct) const {
479 return subview(vec_->getLocalViewHost(Tpetra::Access::OverwriteAll),Kokkos::ALL(), Kokkos::ALL());
480
481 }
482
483 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
484 typename TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::dual_view_type::t_dev_um
485 TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
486 getDeviceLocalView(Access::OverwriteAllStruct) const {
487 return subview(vec_->getLocalViewDevice(Tpetra::Access::OverwriteAll),Kokkos::ALL(), Kokkos::ALL());
488 }
489
490 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
491 typename TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::dual_view_type::t_host_um
492 TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
493 getHostLocalView (Access::ReadWriteStruct) const {
494 return subview(vec_->getLocalViewHost(Tpetra::Access::ReadWrite),Kokkos::ALL(), Kokkos::ALL());
495
496 }
497
498 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
499 typename TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::dual_view_type::t_dev_um
500 TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
501 getDeviceLocalView(Access::ReadWriteStruct) const {
502 return subview(vec_->getLocalViewDevice(Tpetra::Access::ReadWrite),Kokkos::ALL(), Kokkos::ALL());
503 }
504
505#endif
506
509 template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
510 void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
511 assign (const MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>& rhs)
512 {
513 typedef TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> this_type;
514 const this_type* rhsPtr = dynamic_cast<const this_type*> (&rhs);
515 TEUCHOS_TEST_FOR_EXCEPTION(
516 rhsPtr == NULL, std::invalid_argument, "Xpetra::MultiVector::operator=:"
517 " The left-hand side (LHS) of the assignment has a different type than "
518 "the right-hand side (RHS). The LHS has type Xpetra::TpetraMultiVector"
519 " (which means it wraps a Tpetra::MultiVector), but the RHS has some "
520 "other type. This probably means that the RHS wraps an "
521 "Epetra_MultiVector. Xpetra::MultiVector does not currently implement "
522 "assignment from an Epetra object to a Tpetra object, though this could"
523 " be added with sufficient interest.");
524
525 typedef Tpetra::MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> TMV;
526 RCP<const TMV> rhsImpl = rhsPtr->getTpetra_MultiVector ();
527 RCP<TMV> lhsImpl = this->getTpetra_MultiVector ();
528
529 TEUCHOS_TEST_FOR_EXCEPTION(
530 rhsImpl.is_null (), std::logic_error, "Xpetra::MultiVector::operator= "
531 "(in Xpetra::TpetraMultiVector::assign): *this (the right-hand side of "
532 "the assignment) has a null RCP<Tpetra::MultiVector> inside. Please "
533 "report this bug to the Xpetra developers.");
534 TEUCHOS_TEST_FOR_EXCEPTION(
535 lhsImpl.is_null (), std::logic_error, "Xpetra::MultiVector::operator= "
536 "(in Xpetra::TpetraMultiVector::assign): The left-hand side of the "
537 "assignment has a null RCP<Tpetra::MultiVector> inside. Please report "
538 "this bug to the Xpetra developers.");
539
540 Tpetra::deep_copy (*lhsImpl, *rhsImpl);
541 }
542
543
544#ifdef HAVE_XPETRA_EPETRA
545
546#if ((defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_OPENMP) || !defined(HAVE_TPETRA_INST_INT_INT))) || \
547 (!defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_SERIAL) || !defined(HAVE_TPETRA_INST_INT_INT))))
548
549 // specialization for TpetraMultiVector on EpetraNode and GO=int
550 template <class Scalar>
551 class TpetraMultiVector<Scalar,int,int,EpetraNode>
552 : public virtual MultiVector< Scalar, int, int, EpetraNode >
553 {
554 typedef int LocalOrdinal;
555 typedef int GlobalOrdinal;
556 typedef EpetraNode Node;
557
558 // The following typedef are used by the XPETRA_DYNAMIC_CAST() macro.
559 typedef TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> TpetraMultiVectorClass;
560
561 public:
562
564
565
567 TpetraMultiVector () {
568 XPETRA_TPETRA_ETI_EXCEPTION( typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name() , typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name(), "int", typeid(EpetraNode).name() );
569 }
570
572 TpetraMultiVector(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map, size_t NumVectors, bool zeroOut=true) {
573 XPETRA_TPETRA_ETI_EXCEPTION( typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name() , typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name(), "int", typeid(EpetraNode).name() );
574 }
575
577 TpetraMultiVector(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Teuchos::DataAccess copyOrView) {
578 XPETRA_TPETRA_ETI_EXCEPTION( typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name() , typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name(), "int", typeid(EpetraNode).name() );
579 }
580
582 TpetraMultiVector(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map, const Teuchos::ArrayView< const Scalar > &A, size_t LDA, size_t NumVectors) {
583 XPETRA_TPETRA_ETI_EXCEPTION( typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name() , typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name(), "int", typeid(EpetraNode).name() );
584 }
585
587 TpetraMultiVector(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map, const Teuchos::ArrayView< const Teuchos::ArrayView< const Scalar > > &ArrayOfPtrs, size_t NumVectors) {
588 XPETRA_TPETRA_ETI_EXCEPTION( typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name() , typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name(), "int", typeid(EpetraNode).name() );
589 }
590
591
593 virtual ~TpetraMultiVector() { }
594
596
598
599
601 void replaceGlobalValue(GlobalOrdinal globalRow, size_t vectorIndex, const Scalar &value) { }
602
604 void sumIntoGlobalValue(GlobalOrdinal globalRow, size_t vectorIndex, const Scalar &value) { }
605
607 void replaceLocalValue(LocalOrdinal myRow, size_t vectorIndex, const Scalar &value) { }
608
610 void sumIntoLocalValue(LocalOrdinal myRow, size_t vectorIndex, const Scalar &value) { }
611
613 void putScalar(const Scalar &value) { }
614
616 void reduce() { }
617
619
621
622
624 Teuchos::RCP< const Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > getVector(size_t j) const { return Teuchos::null; }
625
627 Teuchos::RCP< Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > getVectorNonConst(size_t j) { return Teuchos::null; }
628
630 Teuchos::ArrayRCP< const Scalar > getData(size_t j) const { return Teuchos::ArrayRCP< const Scalar >(); }
631
633 Teuchos::ArrayRCP< Scalar > getDataNonConst(size_t j) { return Teuchos::ArrayRCP< Scalar >(); }
634
636 void get1dCopy(Teuchos::ArrayView< Scalar > A, size_t LDA) const { }
637
639 void get2dCopy(Teuchos::ArrayView< const Teuchos::ArrayView< Scalar > > ArrayOfPtrs) const { }
640
642 Teuchos::ArrayRCP< const Scalar > get1dView() const { return Teuchos::ArrayRCP< const Scalar >(); }
643
645 Teuchos::ArrayRCP< Teuchos::ArrayRCP< const Scalar > > get2dView() const { return Teuchos::ArrayRCP< Teuchos::ArrayRCP< const Scalar > >(); }
646
648 Teuchos::ArrayRCP< Scalar > get1dViewNonConst() { return Teuchos::ArrayRCP< Scalar >(); }
649
651 Teuchos::ArrayRCP< Teuchos::ArrayRCP< Scalar > > get2dViewNonConst() { return Teuchos::ArrayRCP< Teuchos::ArrayRCP< Scalar > >(); }
652
654
656
657
659 void dot(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const Teuchos::ArrayView< Scalar > &dots) const { }
660
662 void abs(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A) { }
663
665 void reciprocal(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A) { }
666
668 void scale(const Scalar &alpha) { }
669
671 void scale(Teuchos::ArrayView< const Scalar > alpha) { }
672
674 void scale(const Scalar &alpha, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A) { }
675
677 void update(const Scalar &alpha, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const Scalar &beta) { }
678
680 void update(const Scalar &alpha, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const Scalar &beta, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &B, const Scalar &gamma) { }
681
683 void norm1(const Teuchos::ArrayView< typename Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const { }
684
686 void norm2(const Teuchos::ArrayView< typename Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const { }
687
689 void normInf(const Teuchos::ArrayView< typename Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const { }
690
692 void meanValue(const Teuchos::ArrayView< Scalar > &means) const { }
693
695 void multiply(Teuchos::ETransp transA, Teuchos::ETransp transB, const Scalar &alpha, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &B, const Scalar &beta) { }
696
698
700
701
703 size_t getNumVectors() const { return 0; }
704
706 size_t getLocalLength() const { return 0; }
707
709 global_size_t getGlobalLength() const { return 0; }
710
711 // \! Checks to see if the local length, number of vectors and size of Scalar type match
712 bool isSameSize(const MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> & vec) const { return false; }
713
715
717
718
720 std::string description() const { return std::string(""); }
721
723 void describe(Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel=Teuchos::Describable::verbLevel_default) const { }
724
726
728 void elementWiseMultiply(Scalar scalarAB, const Vector<Scalar,LocalOrdinal,GlobalOrdinal,Node> &A, const MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> &B, Scalar scalarThis) {};
729
731 void randomize(bool bUseXpetraImplementation = false) { }
732
733 //{@
734 // Implements DistObject interface
735
736 Teuchos::RCP< const Map<LocalOrdinal,GlobalOrdinal,Node> > getMap() const { return Teuchos::null; }
737
738 void doImport(const DistObject< Scalar, LocalOrdinal,GlobalOrdinal,Node> &source, const Import<LocalOrdinal,GlobalOrdinal,Node> &importer, CombineMode CM) { }
739
740 void beginImport(const DistObject< Scalar, LocalOrdinal,GlobalOrdinal,Node> &source, const Import<LocalOrdinal,GlobalOrdinal,Node> &importer, CombineMode CM) { }
741
742 void endImport(const DistObject< Scalar, LocalOrdinal,GlobalOrdinal,Node> &source, const Import<LocalOrdinal,GlobalOrdinal,Node> &importer, CombineMode CM) { }
743
744 void doExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Import<LocalOrdinal,GlobalOrdinal,Node>& importer, CombineMode CM) { }
745
746 void beginExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Import<LocalOrdinal,GlobalOrdinal,Node>& importer, CombineMode CM) { }
747
748 void endExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Import<LocalOrdinal,GlobalOrdinal,Node>& importer, CombineMode CM) { }
749
750 void doImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Export<LocalOrdinal,GlobalOrdinal,Node>& exporter, CombineMode CM) { }
751
752 void beginImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Export<LocalOrdinal,GlobalOrdinal,Node>& exporter, CombineMode CM) { }
753
754 void endImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Export<LocalOrdinal,GlobalOrdinal,Node>& exporter, CombineMode CM) { }
755
756 void doExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Export<LocalOrdinal,GlobalOrdinal,Node>& exporter, CombineMode CM) { }
757
758 void beginExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Export<LocalOrdinal,GlobalOrdinal,Node>& exporter, CombineMode CM) { }
759
760 void endExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Export<LocalOrdinal,GlobalOrdinal,Node>& exporter, CombineMode CM) { }
761
762 void replaceMap(const RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> >& map) { }
763
765
767
768
770 TpetraMultiVector(const Teuchos::RCP<Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node> > &vec) {
771 XPETRA_TPETRA_ETI_EXCEPTION( typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name() , typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name(), "int", typeid(EpetraNode).name() );
772 }
773
775 RCP< Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node> > getTpetra_MultiVector() const { return Teuchos::null; }
776
778 void setSeed(unsigned int seed) { }
779
781
782 protected:
785 virtual void
786 assign (const MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>& rhs)
787 { }
788 }; // TpetraMultiVector class (specialization GO=int, NO=EpetraNode)
789#endif
790
791#if ((defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_OPENMP) || !defined(HAVE_TPETRA_INST_INT_LONG_LONG))) || \
792 (!defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_SERIAL) || !defined(HAVE_TPETRA_INST_INT_LONG_LONG))))
793
794 // specialization for TpetraMultiVector on EpetraNode and GO=long long
795 template <class Scalar>
796 class TpetraMultiVector<Scalar,int,long long,EpetraNode>
797 : public virtual MultiVector< Scalar, int, long long, EpetraNode >
798 {
799 typedef int LocalOrdinal;
800 typedef long long GlobalOrdinal;
801 typedef EpetraNode Node;
802
803 // The following typedef are used by the XPETRA_DYNAMIC_CAST() macro.
804 typedef TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> TpetraMultiVectorClass;
805
806 public:
807
809
810
812 TpetraMultiVector(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map, size_t NumVectors, bool zeroOut=true) {
813 XPETRA_TPETRA_ETI_EXCEPTION( typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name() , typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name(), "long long", typeid(EpetraNode).name() );
814 }
815
817 TpetraMultiVector(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Teuchos::DataAccess copyOrView) {
818 XPETRA_TPETRA_ETI_EXCEPTION( typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name() , typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name(), "long long", typeid(EpetraNode).name() );
819 }
820
822 TpetraMultiVector(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map, const Teuchos::ArrayView< const Scalar > &A, size_t LDA, size_t NumVectors) {
823 XPETRA_TPETRA_ETI_EXCEPTION( typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name() , typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name(), "long long", typeid(EpetraNode).name() );
824 }
825
827 TpetraMultiVector(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map, const Teuchos::ArrayView< const Teuchos::ArrayView< const Scalar > > &ArrayOfPtrs, size_t NumVectors) {
828 XPETRA_TPETRA_ETI_EXCEPTION( typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name() , typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name(), "long long", typeid(EpetraNode).name() );
829 }
830
831
833 virtual ~TpetraMultiVector() { }
834
836
838
839
841 void replaceGlobalValue(GlobalOrdinal globalRow, size_t vectorIndex, const Scalar &value) { }
842
844 void sumIntoGlobalValue(GlobalOrdinal globalRow, size_t vectorIndex, const Scalar &value) { }
845
847 void replaceLocalValue(LocalOrdinal myRow, size_t vectorIndex, const Scalar &value) { }
848
850 void sumIntoLocalValue(LocalOrdinal myRow, size_t vectorIndex, const Scalar &value) { }
851
853 void putScalar(const Scalar &value) { }
854
856 void reduce() { }
857
859
861
862
864 Teuchos::RCP< const Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > getVector(size_t j) const { return Teuchos::null; }
865
867 Teuchos::RCP< Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > getVectorNonConst(size_t j) { return Teuchos::null; }
868
870 Teuchos::ArrayRCP< const Scalar > getData(size_t j) const { return Teuchos::ArrayRCP< const Scalar >(); }
871
873 Teuchos::ArrayRCP< Scalar > getDataNonConst(size_t j) { return Teuchos::ArrayRCP< Scalar >(); }
874
876 void get1dCopy(Teuchos::ArrayView< Scalar > A, size_t LDA) const { }
877
879 void get2dCopy(Teuchos::ArrayView< const Teuchos::ArrayView< Scalar > > ArrayOfPtrs) const { }
880
882 Teuchos::ArrayRCP< const Scalar > get1dView() const { return Teuchos::ArrayRCP< const Scalar >(); }
883
885 Teuchos::ArrayRCP< Teuchos::ArrayRCP< const Scalar > > get2dView() const { return Teuchos::ArrayRCP< Teuchos::ArrayRCP< const Scalar > >(); }
886
888 Teuchos::ArrayRCP< Scalar > get1dViewNonConst() { return Teuchos::ArrayRCP< Scalar >(); }
889
891 Teuchos::ArrayRCP< Teuchos::ArrayRCP< Scalar > > get2dViewNonConst() { return Teuchos::ArrayRCP< Teuchos::ArrayRCP< Scalar > >(); }
892
894
896
897
899 void dot(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const Teuchos::ArrayView< Scalar > &dots) const { }
900
902 void abs(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A) { }
903
905 void reciprocal(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A) { }
906
908 void scale(const Scalar &alpha) { }
909
911 void scale(Teuchos::ArrayView< const Scalar > alpha) { }
912
914 void scale(const Scalar &alpha, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A) { }
915
917 void update(const Scalar &alpha, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const Scalar &beta) { }
918
920 void update(const Scalar &alpha, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const Scalar &beta, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &B, const Scalar &gamma) { }
921
923 void norm1(const Teuchos::ArrayView< typename Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const { }
924
926 void norm2(const Teuchos::ArrayView< typename Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const { }
927
929 void normInf(const Teuchos::ArrayView< typename Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const { }
930
932 void meanValue(const Teuchos::ArrayView< Scalar > &means) const { }
933
935 void multiply(Teuchos::ETransp transA, Teuchos::ETransp transB, const Scalar &alpha, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &B, const Scalar &beta) { }
936
938
940
941
943 size_t getNumVectors() const { return 0; }
944
946 size_t getLocalLength() const { return 0; }
947
949 global_size_t getGlobalLength() const { return 0; }
950
951 // \! Checks to see if the local length, number of vectors and size of Scalar type match
952 bool isSameSize(const MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> & vec) const { return false; }
953
955
957
958
960 std::string description() const { return std::string(""); }
961
963 void describe(Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel=Teuchos::Describable::verbLevel_default) const { }
964
966
968 void elementWiseMultiply(Scalar scalarAB, const Vector<Scalar,LocalOrdinal,GlobalOrdinal,Node> &A, const MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> &B, Scalar scalarThis) {};
969
971 void randomize(bool bUseXpetraImplementation = false) { }
972
973 //{@
974 // Implements DistObject interface
975
976 Teuchos::RCP< const Map<LocalOrdinal,GlobalOrdinal,Node> > getMap() const { return Teuchos::null; }
977
978 void doImport(const DistObject< Scalar, LocalOrdinal,GlobalOrdinal,Node> &source, const Import<LocalOrdinal,GlobalOrdinal,Node> &importer, CombineMode CM) { }
979
980 void beginImport(const DistObject< Scalar, LocalOrdinal,GlobalOrdinal,Node> &source, const Import<LocalOrdinal,GlobalOrdinal,Node> &importer, CombineMode CM) { }
981
982 void endImport(const DistObject< Scalar, LocalOrdinal,GlobalOrdinal,Node> &source, const Import<LocalOrdinal,GlobalOrdinal,Node> &importer, CombineMode CM) { }
983
984 void doExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Import<LocalOrdinal,GlobalOrdinal,Node>& importer, CombineMode CM) { }
985
986 void beginExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Import<LocalOrdinal,GlobalOrdinal,Node>& importer, CombineMode CM) { }
987
988 void endExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Import<LocalOrdinal,GlobalOrdinal,Node>& importer, CombineMode CM) { }
989
990 void doImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Export<LocalOrdinal,GlobalOrdinal,Node>& exporter, CombineMode CM) { }
991
992 void beginImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Export<LocalOrdinal,GlobalOrdinal,Node>& exporter, CombineMode CM) { }
993
994 void endImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Export<LocalOrdinal,GlobalOrdinal,Node>& exporter, CombineMode CM) { }
995
996 void doExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Export<LocalOrdinal,GlobalOrdinal,Node>& exporter, CombineMode CM) { }
997
998 void beginExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Export<LocalOrdinal,GlobalOrdinal,Node>& exporter, CombineMode CM) { }
999
1000 void endExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Export<LocalOrdinal,GlobalOrdinal,Node>& exporter, CombineMode CM) { }
1001
1002 void replaceMap(const RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> >& map) { }
1003
1005
1007
1008
1010 TpetraMultiVector(const Teuchos::RCP<Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node> > &vec) {
1011 XPETRA_TPETRA_ETI_EXCEPTION( typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name() , typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name(), "long long", typeid(EpetraNode).name() );
1012 }
1013
1015 RCP< Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node> > getTpetra_MultiVector() const { return Teuchos::null; }
1016
1018 void setSeed(unsigned int seed) { }
1019
1021
1022 protected:
1025 virtual void
1026 assign (const MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>& rhs)
1027 { }
1028 }; // TpetraMultiVector class (specialization GO=int, NO=EpetraNode)
1029
1030#endif // TpetraMultiVector class (specialization GO=long long, NO=EpetraNode)
1031
1032#endif // HAVE_XPETRA_EPETRA
1033
1034} // Xpetra namespace
1035
1036// Following header file inculsion is needed for the dynamic_cast to TpetraVector in
1037// elementWiseMultiply (because we cannot dynamic_cast if target is not a complete type)
1038// It is included here to avoid circular dependency between Vector and MultiVector
1039// TODO: there is certainly a more elegant solution...
1040#include "Xpetra_TpetraVector.hpp"
1041
1042namespace Xpetra {
1043
1044 template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
1045 void
1046 TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
1047 elementWiseMultiply(Scalar scalarAB,
1048 const Vector<Scalar,LocalOrdinal,GlobalOrdinal,Node> &A,
1049 const MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> &B,
1050 Scalar scalarThis)
1051 {
1052 XPETRA_MONITOR("TpetraMultiVector::elementWiseMultiply");
1053
1054 // XPETRA_DYNAMIC_CAST won't take TpetraVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>
1055 // as an argument, hence the following typedef.
1056 typedef TpetraVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> tpv;
1057 XPETRA_DYNAMIC_CAST(const tpv, A, tA, "Xpetra::TpetraMultiVectorMatrix->multiply() only accept Xpetra::TpetraMultiVector as input arguments.");
1058 XPETRA_DYNAMIC_CAST(const TpetraMultiVector, B, tB, "Xpetra::TpetraMultiVectorMatrix->multiply() only accept Xpetra::TpetraMultiVector as input arguments.");
1059 vec_->elementWiseMultiply(scalarAB, *tA.getTpetra_Vector(), *tB.getTpetra_MultiVector(), scalarThis);
1060 }
1061
1062} // Xpetra namespace
1063
1064#define XPETRA_TPETRAMULTIVECTOR_SHORT
1065#endif // XPETRA_TPETRAMULTIVECTOR_HPP
XPETRA_MONITOR
#define XPETRA_MONITOR(funcName)
Definition: Xpetra_ConfigDefs.hpp:132
XPETRA_DYNAMIC_CAST
#define XPETRA_DYNAMIC_CAST(type, obj, newObj, exceptionMsg)
Definition: Xpetra_Exceptions.hpp:58
XPETRA_TPETRA_ETI_EXCEPTION
#define XPETRA_TPETRA_ETI_EXCEPTION(cl, obj, go, node)
Definition: Xpetra_Exceptions.hpp:79
Xpetra_TpetraConfigDefs.hpp
Xpetra_TpetraMultiVector_decl.hpp
Xpetra_Utils.hpp
Teuchos::Describable::verbLevel_default
static const EVerbosityLevel verbLevel_default
Teuchos::RCP::is_null
bool is_null() const
Teuchos::basic_FancyOStream
Xpetra::MultiVector::Xpetra_randomize
virtual void Xpetra_randomize()
Set multi-vector values to random numbers. XPetra implementation.
Definition: Xpetra_MultiVector_def.hpp:75
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::replaceMap
void replaceMap(const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map)
Definition: Xpetra_TpetraMultiVector_def.hpp:762
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::endExport
void endExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Export< LocalOrdinal, GlobalOrdinal, Node > &exporter, CombineMode CM)
Export data into this object using an Export object ("forward mode").
Definition: Xpetra_TpetraMultiVector_def.hpp:760
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::beginImport
void beginImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Export< LocalOrdinal, GlobalOrdinal, Node > &exporter, CombineMode CM)
Import data into this object using an Export object ("reverse mode").
Definition: Xpetra_TpetraMultiVector_def.hpp:752
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::endImport
void endImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Import data into this object using an Import object ("forward mode").
Definition: Xpetra_TpetraMultiVector_def.hpp:742
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::elementWiseMultiply
void elementWiseMultiply(Scalar scalarAB, const Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &B, Scalar scalarThis)
Element-wise multiply of a Vector A with a TpetraMultiVector B.
Definition: Xpetra_TpetraMultiVector_def.hpp:728
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::get2dView
Teuchos::ArrayRCP< Teuchos::ArrayRCP< const Scalar > > get2dView() const
Return const persisting pointers to values.
Definition: Xpetra_TpetraMultiVector_def.hpp:645
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::getData
Teuchos::ArrayRCP< const Scalar > getData(size_t j) const
Const view of the local values in a particular vector of this multivector.
Definition: Xpetra_TpetraMultiVector_def.hpp:630
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::TpetraMultiVector
TpetraMultiVector(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Teuchos::DataAccess copyOrView)
Copy constructor (performs a deep copy).
Definition: Xpetra_TpetraMultiVector_def.hpp:577
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::doImport
void doImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Import data into this object using an Import object ("forward mode").
Definition: Xpetra_TpetraMultiVector_def.hpp:738
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::sumIntoGlobalValue
void sumIntoGlobalValue(GlobalOrdinal globalRow, size_t vectorIndex, const Scalar &value)
Add value to existing value, using global (row) index.
Definition: Xpetra_TpetraMultiVector_def.hpp:604
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::TpetraMultiVector
TpetraMultiVector(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map, const Teuchos::ArrayView< const Teuchos::ArrayView< const Scalar > > &ArrayOfPtrs, size_t NumVectors)
Create multivector by copying array of views of local data.
Definition: Xpetra_TpetraMultiVector_def.hpp:587
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::TpetraMultiVector
TpetraMultiVector(const Teuchos::RCP< Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &vec)
TpetraMultiVector constructor to wrap a Tpetra::MultiVector object.
Definition: Xpetra_TpetraMultiVector_def.hpp:770
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::beginExport
void beginExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Export< LocalOrdinal, GlobalOrdinal, Node > &exporter, CombineMode CM)
Export data into this object using an Export object ("forward mode").
Definition: Xpetra_TpetraMultiVector_def.hpp:758
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::getMap
Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > getMap() const
The Map describing the parallel distribution of this object.
Definition: Xpetra_TpetraMultiVector_def.hpp:736
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::scale
void scale(Teuchos::ArrayView< const Scalar > alpha)
Scale the current values of a multi-vector, this[j] = alpha[j]*this[j].
Definition: Xpetra_TpetraMultiVector_def.hpp:671
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::getGlobalLength
global_size_t getGlobalLength() const
Global number of rows in the multivector.
Definition: Xpetra_TpetraMultiVector_def.hpp:709
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::update
void update(const Scalar &alpha, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const Scalar &beta)
Update multi-vector values with scaled values of A, this = beta*this + alpha*A.
Definition: Xpetra_TpetraMultiVector_def.hpp:677
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::reduce
void reduce()
Sum values of a locally replicated multivector across all processes.
Definition: Xpetra_TpetraMultiVector_def.hpp:616
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::assign
virtual void assign(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &rhs)
Implementation of the assignment operator (operator=); does a deep copy.
Definition: Xpetra_TpetraMultiVector_def.hpp:786
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::endImport
void endImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Export< LocalOrdinal, GlobalOrdinal, Node > &exporter, CombineMode CM)
Import data into this object using an Export object ("reverse mode").
Definition: Xpetra_TpetraMultiVector_def.hpp:754
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::describe
void describe(Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel=Teuchos::Describable::verbLevel_default) const
Print the object with the given verbosity level to a FancyOStream.
Definition: Xpetra_TpetraMultiVector_def.hpp:723
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::TpetraMultiVectorClass
TpetraMultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > TpetraMultiVectorClass
Definition: Xpetra_TpetraMultiVector_def.hpp:559
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::norm2
void norm2(const Teuchos::ArrayView< typename Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const
Definition: Xpetra_TpetraMultiVector_def.hpp:686
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::scale
void scale(const Scalar &alpha)
Scale the current values of a multi-vector, this = alpha*this.
Definition: Xpetra_TpetraMultiVector_def.hpp:668
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::replaceLocalValue
void replaceLocalValue(LocalOrdinal myRow, size_t vectorIndex, const Scalar &value)
Replace value, using local (row) index.
Definition: Xpetra_TpetraMultiVector_def.hpp:607
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::doExport
void doExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Export< LocalOrdinal, GlobalOrdinal, Node > &exporter, CombineMode CM)
Export data into this object using an Export object ("forward mode").
Definition: Xpetra_TpetraMultiVector_def.hpp:756
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::doExport
void doExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Export data into this object using an Import object ("reverse mode").
Definition: Xpetra_TpetraMultiVector_def.hpp:744
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::TpetraMultiVector
TpetraMultiVector(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map, size_t NumVectors, bool zeroOut=true)
Basic constuctor.
Definition: Xpetra_TpetraMultiVector_def.hpp:572
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::dot
void dot(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const Teuchos::ArrayView< Scalar > &dots) const
Compute dot product of each corresponding pair of vectors, dots[i] = this[i].dot(A[i]).
Definition: Xpetra_TpetraMultiVector_def.hpp:659
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::getVectorNonConst
Teuchos::RCP< Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > getVectorNonConst(size_t j)
Return a Vector which is a nonconst view of column j.
Definition: Xpetra_TpetraMultiVector_def.hpp:627
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::~TpetraMultiVector
virtual ~TpetraMultiVector()
Destructor (virtual for memory safety of derived classes).
Definition: Xpetra_TpetraMultiVector_def.hpp:593
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::doImport
void doImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Export< LocalOrdinal, GlobalOrdinal, Node > &exporter, CombineMode CM)
Import data into this object using an Export object ("reverse mode").
Definition: Xpetra_TpetraMultiVector_def.hpp:750
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::getTpetra_MultiVector
RCP< Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > getTpetra_MultiVector() const
Get the underlying Tpetra multivector.
Definition: Xpetra_TpetraMultiVector_def.hpp:775
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::sumIntoLocalValue
void sumIntoLocalValue(LocalOrdinal myRow, size_t vectorIndex, const Scalar &value)
Add value to existing value, using local (row) index.
Definition: Xpetra_TpetraMultiVector_def.hpp:610
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::multiply
void multiply(Teuchos::ETransp transA, Teuchos::ETransp transB, const Scalar &alpha, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &B, const Scalar &beta)
Matrix-matrix multiplication: this = beta*this + alpha*op(A)*op(B).
Definition: Xpetra_TpetraMultiVector_def.hpp:695
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::getDataNonConst
Teuchos::ArrayRCP< Scalar > getDataNonConst(size_t j)
View of the local values in a particular vector of this multivector.
Definition: Xpetra_TpetraMultiVector_def.hpp:633
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::norm1
void norm1(const Teuchos::ArrayView< typename Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const
Compute 1-norm of each vector in multi-vector.
Definition: Xpetra_TpetraMultiVector_def.hpp:683
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::get1dView
Teuchos::ArrayRCP< const Scalar > get1dView() const
Const persisting (1-D) view of this multivector's local values.
Definition: Xpetra_TpetraMultiVector_def.hpp:642
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::putScalar
void putScalar(const Scalar &value)
Set all values in the multivector with the given value.
Definition: Xpetra_TpetraMultiVector_def.hpp:613
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::getLocalLength
size_t getLocalLength() const
Local number of rows on the calling process.
Definition: Xpetra_TpetraMultiVector_def.hpp:706
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::isSameSize
bool isSameSize(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &vec) const
Definition: Xpetra_TpetraMultiVector_def.hpp:712
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::beginExport
void beginExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Export data into this object using an Import object ("reverse mode").
Definition: Xpetra_TpetraMultiVector_def.hpp:746
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::randomize
void randomize(bool bUseXpetraImplementation=false)
Set multi-vector values to random numbers.
Definition: Xpetra_TpetraMultiVector_def.hpp:731
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::setSeed
void setSeed(unsigned int seed)
Set seed for Random function.
Definition: Xpetra_TpetraMultiVector_def.hpp:778
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::description
std::string description() const
A simple one-line description of this object.
Definition: Xpetra_TpetraMultiVector_def.hpp:720
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::scale
void scale(const Scalar &alpha, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A)
Replace multi-vector values with scaled values of A, this = alpha*A.
Definition: Xpetra_TpetraMultiVector_def.hpp:674
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::update
void update(const Scalar &alpha, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const Scalar &beta, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &B, const Scalar &gamma)
Update multi-vector with scaled values of A and B, this = gamma*this + alpha*A + beta*B.
Definition: Xpetra_TpetraMultiVector_def.hpp:680
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::normInf
void normInf(const Teuchos::ArrayView< typename Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const
Compute Inf-norm of each vector in multi-vector.
Definition: Xpetra_TpetraMultiVector_def.hpp:689
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::get1dViewNonConst
Teuchos::ArrayRCP< Scalar > get1dViewNonConst()
Nonconst persisting (1-D) view of this multivector's local values.
Definition: Xpetra_TpetraMultiVector_def.hpp:648
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::getVector
Teuchos::RCP< const Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > getVector(size_t j) const
Return a Vector which is a const view of column j.
Definition: Xpetra_TpetraMultiVector_def.hpp:624
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::meanValue
void meanValue(const Teuchos::ArrayView< Scalar > &means) const
Compute mean (average) value of each vector in multi-vector. The outcome of this routine is undefined...
Definition: Xpetra_TpetraMultiVector_def.hpp:692
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::getNumVectors
size_t getNumVectors() const
Number of columns in the multivector.
Definition: Xpetra_TpetraMultiVector_def.hpp:703
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::get1dCopy
void get1dCopy(Teuchos::ArrayView< Scalar > A, size_t LDA) const
Fill the given array with a copy of this multivector's local values.
Definition: Xpetra_TpetraMultiVector_def.hpp:636
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::get2dViewNonConst
Teuchos::ArrayRCP< Teuchos::ArrayRCP< Scalar > > get2dViewNonConst()
Return non-const persisting pointers to values.
Definition: Xpetra_TpetraMultiVector_def.hpp:651
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::endExport
void endExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Export data into this object using an Import object ("reverse mode").
Definition: Xpetra_TpetraMultiVector_def.hpp:748
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::abs
void abs(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A)
Put element-wise absolute values of input Multi-vector in target: A = abs(this).
Definition: Xpetra_TpetraMultiVector_def.hpp:662
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::replaceGlobalValue
void replaceGlobalValue(GlobalOrdinal globalRow, size_t vectorIndex, const Scalar &value)
Replace value, using global (row) index.
Definition: Xpetra_TpetraMultiVector_def.hpp:601
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::beginImport
void beginImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Import data into this object using an Import object ("forward mode").
Definition: Xpetra_TpetraMultiVector_def.hpp:740
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::get2dCopy
void get2dCopy(Teuchos::ArrayView< const Teuchos::ArrayView< Scalar > > ArrayOfPtrs) const
Fill the given array with a copy of this multivector's local values.
Definition: Xpetra_TpetraMultiVector_def.hpp:639
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::TpetraMultiVector
TpetraMultiVector(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map, const Teuchos::ArrayView< const Scalar > &A, size_t LDA, size_t NumVectors)
Create multivector by copying two-dimensional array of local data.
Definition: Xpetra_TpetraMultiVector_def.hpp:582
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::reciprocal
void reciprocal(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A)
Put element-wise reciprocal values of input Multi-vector in target, this(i,j) = 1/A(i,...
Definition: Xpetra_TpetraMultiVector_def.hpp:665
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::beginImport
void beginImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Export< LocalOrdinal, GlobalOrdinal, Node > &exporter, CombineMode CM)
Import data into this object using an Export object ("reverse mode").
Definition: Xpetra_TpetraMultiVector_def.hpp:992
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::get2dViewNonConst
Teuchos::ArrayRCP< Teuchos::ArrayRCP< Scalar > > get2dViewNonConst()
Return non-const persisting pointers to values.
Definition: Xpetra_TpetraMultiVector_def.hpp:891
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::doExport
void doExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Export data into this object using an Import object ("reverse mode").
Definition: Xpetra_TpetraMultiVector_def.hpp:984
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::scale
void scale(Teuchos::ArrayView< const Scalar > alpha)
Scale the current values of a multi-vector, this[j] = alpha[j]*this[j].
Definition: Xpetra_TpetraMultiVector_def.hpp:911
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::norm2
void norm2(const Teuchos::ArrayView< typename Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const
Definition: Xpetra_TpetraMultiVector_def.hpp:926
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::doImport
void doImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Import data into this object using an Import object ("forward mode").
Definition: Xpetra_TpetraMultiVector_def.hpp:978
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::getLocalLength
size_t getLocalLength() const
Local number of rows on the calling process.
Definition: Xpetra_TpetraMultiVector_def.hpp:946
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::get2dCopy
void get2dCopy(Teuchos::ArrayView< const Teuchos::ArrayView< Scalar > > ArrayOfPtrs) const
Fill the given array with a copy of this multivector's local values.
Definition: Xpetra_TpetraMultiVector_def.hpp:879
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::description
std::string description() const
A simple one-line description of this object.
Definition: Xpetra_TpetraMultiVector_def.hpp:960
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::getVector
Teuchos::RCP< const Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > getVector(size_t j) const
Return a Vector which is a const view of column j.
Definition: Xpetra_TpetraMultiVector_def.hpp:864
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::getGlobalLength
global_size_t getGlobalLength() const
Global number of rows in the multivector.
Definition: Xpetra_TpetraMultiVector_def.hpp:949
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::multiply
void multiply(Teuchos::ETransp transA, Teuchos::ETransp transB, const Scalar &alpha, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &B, const Scalar &beta)
Matrix-matrix multiplication: this = beta*this + alpha*op(A)*op(B).
Definition: Xpetra_TpetraMultiVector_def.hpp:935
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::getVectorNonConst
Teuchos::RCP< Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > getVectorNonConst(size_t j)
Return a Vector which is a nonconst view of column j.
Definition: Xpetra_TpetraMultiVector_def.hpp:867
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::update
void update(const Scalar &alpha, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const Scalar &beta, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &B, const Scalar &gamma)
Update multi-vector with scaled values of A and B, this = gamma*this + alpha*A + beta*B.
Definition: Xpetra_TpetraMultiVector_def.hpp:920
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::meanValue
void meanValue(const Teuchos::ArrayView< Scalar > &means) const
Compute mean (average) value of each vector in multi-vector. The outcome of this routine is undefined...
Definition: Xpetra_TpetraMultiVector_def.hpp:932
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::describe
void describe(Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel=Teuchos::Describable::verbLevel_default) const
Print the object with the given verbosity level to a FancyOStream.
Definition: Xpetra_TpetraMultiVector_def.hpp:963
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::endExport
void endExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Export< LocalOrdinal, GlobalOrdinal, Node > &exporter, CombineMode CM)
Export data into this object using an Export object ("forward mode").
Definition: Xpetra_TpetraMultiVector_def.hpp:1000
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::sumIntoGlobalValue
void sumIntoGlobalValue(GlobalOrdinal globalRow, size_t vectorIndex, const Scalar &value)
Add value to existing value, using global (row) index.
Definition: Xpetra_TpetraMultiVector_def.hpp:844
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::TpetraMultiVector
TpetraMultiVector(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Teuchos::DataAccess copyOrView)
Copy constructor (performs a deep copy).
Definition: Xpetra_TpetraMultiVector_def.hpp:817
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::endImport
void endImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Import data into this object using an Import object ("forward mode").
Definition: Xpetra_TpetraMultiVector_def.hpp:982
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::get2dView
Teuchos::ArrayRCP< Teuchos::ArrayRCP< const Scalar > > get2dView() const
Return const persisting pointers to values.
Definition: Xpetra_TpetraMultiVector_def.hpp:885
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::getMap
Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > getMap() const
The Map describing the parallel distribution of this object.
Definition: Xpetra_TpetraMultiVector_def.hpp:976
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::beginImport
void beginImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Import data into this object using an Import object ("forward mode").
Definition: Xpetra_TpetraMultiVector_def.hpp:980
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::normInf
void normInf(const Teuchos::ArrayView< typename Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const
Compute Inf-norm of each vector in multi-vector.
Definition: Xpetra_TpetraMultiVector_def.hpp:929
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::abs
void abs(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A)
Put element-wise absolute values of input Multi-vector in target: A = abs(this).
Definition: Xpetra_TpetraMultiVector_def.hpp:902
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::update
void update(const Scalar &alpha, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const Scalar &beta)
Update multi-vector values with scaled values of A, this = beta*this + alpha*A.
Definition: Xpetra_TpetraMultiVector_def.hpp:917
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::sumIntoLocalValue
void sumIntoLocalValue(LocalOrdinal myRow, size_t vectorIndex, const Scalar &value)
Add value to existing value, using local (row) index.
Definition: Xpetra_TpetraMultiVector_def.hpp:850
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::reduce
void reduce()
Sum values of a locally replicated multivector across all processes.
Definition: Xpetra_TpetraMultiVector_def.hpp:856
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::TpetraMultiVector
TpetraMultiVector(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map, const Teuchos::ArrayView< const Teuchos::ArrayView< const Scalar > > &ArrayOfPtrs, size_t NumVectors)
Create multivector by copying array of views of local data.
Definition: Xpetra_TpetraMultiVector_def.hpp:827
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::get1dView
Teuchos::ArrayRCP< const Scalar > get1dView() const
Const persisting (1-D) view of this multivector's local values.
Definition: Xpetra_TpetraMultiVector_def.hpp:882
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::assign
virtual void assign(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &rhs)
Implementation of the assignment operator (operator=); does a deep copy.
Definition: Xpetra_TpetraMultiVector_def.hpp:1026
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::getDataNonConst
Teuchos::ArrayRCP< Scalar > getDataNonConst(size_t j)
View of the local values in a particular vector of this multivector.
Definition: Xpetra_TpetraMultiVector_def.hpp:873
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::getNumVectors
size_t getNumVectors() const
Number of columns in the multivector.
Definition: Xpetra_TpetraMultiVector_def.hpp:943
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::replaceLocalValue
void replaceLocalValue(LocalOrdinal myRow, size_t vectorIndex, const Scalar &value)
Replace value, using local (row) index.
Definition: Xpetra_TpetraMultiVector_def.hpp:847
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::TpetraMultiVector
TpetraMultiVector(const Teuchos::RCP< Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &vec)
TpetraMultiVector constructor to wrap a Tpetra::MultiVector object.
Definition: Xpetra_TpetraMultiVector_def.hpp:1010
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::TpetraMultiVectorClass
TpetraMultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > TpetraMultiVectorClass
Definition: Xpetra_TpetraMultiVector_def.hpp:804
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::isSameSize
bool isSameSize(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &vec) const
Definition: Xpetra_TpetraMultiVector_def.hpp:952
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::getData
Teuchos::ArrayRCP< const Scalar > getData(size_t j) const
Const view of the local values in a particular vector of this multivector.
Definition: Xpetra_TpetraMultiVector_def.hpp:870
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::beginExport
void beginExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Export< LocalOrdinal, GlobalOrdinal, Node > &exporter, CombineMode CM)
Export data into this object using an Export object ("forward mode").
Definition: Xpetra_TpetraMultiVector_def.hpp:998
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::TpetraMultiVector
TpetraMultiVector(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map, const Teuchos::ArrayView< const Scalar > &A, size_t LDA, size_t NumVectors)
Create multivector by copying two-dimensional array of local data.
Definition: Xpetra_TpetraMultiVector_def.hpp:822
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::randomize
void randomize(bool bUseXpetraImplementation=false)
Set multi-vector values to random numbers.
Definition: Xpetra_TpetraMultiVector_def.hpp:971
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::putScalar
void putScalar(const Scalar &value)
Set all values in the multivector with the given value.
Definition: Xpetra_TpetraMultiVector_def.hpp:853
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::reciprocal
void reciprocal(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A)
Put element-wise reciprocal values of input Multi-vector in target, this(i,j) = 1/A(i,...
Definition: Xpetra_TpetraMultiVector_def.hpp:905
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::doImport
void doImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Export< LocalOrdinal, GlobalOrdinal, Node > &exporter, CombineMode CM)
Import data into this object using an Export object ("reverse mode").
Definition: Xpetra_TpetraMultiVector_def.hpp:990
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::replaceGlobalValue
void replaceGlobalValue(GlobalOrdinal globalRow, size_t vectorIndex, const Scalar &value)
Replace value, using global (row) index.
Definition: Xpetra_TpetraMultiVector_def.hpp:841
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::beginExport
void beginExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Export data into this object using an Import object ("reverse mode").
Definition: Xpetra_TpetraMultiVector_def.hpp:986
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::getTpetra_MultiVector
RCP< Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > getTpetra_MultiVector() const
Get the underlying Tpetra multivector.
Definition: Xpetra_TpetraMultiVector_def.hpp:1015
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::elementWiseMultiply
void elementWiseMultiply(Scalar scalarAB, const Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &B, Scalar scalarThis)
Element-wise multiply of a Vector A with a TpetraMultiVector B.
Definition: Xpetra_TpetraMultiVector_def.hpp:968
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::scale
void scale(const Scalar &alpha)
Scale the current values of a multi-vector, this = alpha*this.
Definition: Xpetra_TpetraMultiVector_def.hpp:908
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::TpetraMultiVector
TpetraMultiVector(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map, size_t NumVectors, bool zeroOut=true)
Basic constuctor.
Definition: Xpetra_TpetraMultiVector_def.hpp:812
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::norm1
void norm1(const Teuchos::ArrayView< typename Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const
Compute 1-norm of each vector in multi-vector.
Definition: Xpetra_TpetraMultiVector_def.hpp:923
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::~TpetraMultiVector
virtual ~TpetraMultiVector()
Destructor (virtual for memory safety of derived classes).
Definition: Xpetra_TpetraMultiVector_def.hpp:833
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::endImport
void endImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Export< LocalOrdinal, GlobalOrdinal, Node > &exporter, CombineMode CM)
Import data into this object using an Export object ("reverse mode").
Definition: Xpetra_TpetraMultiVector_def.hpp:994
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::scale
void scale(const Scalar &alpha, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A)
Replace multi-vector values with scaled values of A, this = alpha*A.
Definition: Xpetra_TpetraMultiVector_def.hpp:914
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::dot
void dot(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const Teuchos::ArrayView< Scalar > &dots) const
Compute dot product of each corresponding pair of vectors, dots[i] = this[i].dot(A[i]).
Definition: Xpetra_TpetraMultiVector_def.hpp:899
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::endExport
void endExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Export data into this object using an Import object ("reverse mode").
Definition: Xpetra_TpetraMultiVector_def.hpp:988
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::replaceMap
void replaceMap(const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map)
Definition: Xpetra_TpetraMultiVector_def.hpp:1002
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::doExport
void doExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Export< LocalOrdinal, GlobalOrdinal, Node > &exporter, CombineMode CM)
Export data into this object using an Export object ("forward mode").
Definition: Xpetra_TpetraMultiVector_def.hpp:996
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::get1dViewNonConst
Teuchos::ArrayRCP< Scalar > get1dViewNonConst()
Nonconst persisting (1-D) view of this multivector's local values.
Definition: Xpetra_TpetraMultiVector_def.hpp:888
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::get1dCopy
void get1dCopy(Teuchos::ArrayView< Scalar > A, size_t LDA) const
Fill the given array with a copy of this multivector's local values.
Definition: Xpetra_TpetraMultiVector_def.hpp:876
Xpetra::TpetraMultiVector::replaceMap
void replaceMap(const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map)
Definition: Xpetra_TpetraMultiVector_def.hpp:436
Xpetra::TpetraMultiVector::replaceGlobalValue
void replaceGlobalValue(GlobalOrdinal globalRow, size_t vectorIndex, const Scalar &value)
Replace value, using global (row) index.
Definition: Xpetra_TpetraMultiVector_def.hpp:103
Xpetra::TpetraMultiVector::doExport
void doExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Export data into this object using an Import object ("reverse mode").
Definition: Xpetra_TpetraMultiVector_def.hpp:336
Xpetra::TpetraMultiVector::isSameSize
bool isSameSize(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &vec) const
Definition: Xpetra_TpetraMultiVector_def.hpp:274
Xpetra::TpetraMultiVector::describe
void describe(Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel=Teuchos::Describable::verbLevel_default) const
Print the object with the given verbosity level to a FancyOStream.
Definition: Xpetra_TpetraMultiVector_def.hpp:284
Xpetra::TpetraMultiVector::multiply
void multiply(Teuchos::ETransp transA, Teuchos::ETransp transB, const Scalar &alpha, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &B, const Scalar &beta)
Matrix-matrix multiplication: this = beta*this + alpha*op(A)*op(B).
Definition: Xpetra_TpetraMultiVector_def.hpp:253
Xpetra::TpetraMultiVector::getVectorNonConst
Teuchos::RCP< Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > getVectorNonConst(size_t j)
Return a Vector which is a nonconst view of column j.
Definition: Xpetra_TpetraMultiVector_def.hpp:140
Xpetra::TpetraMultiVector::doImport
void doImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Import data into this object using an Import object ("forward mode").
Definition: Xpetra_TpetraMultiVector_def.hpp:305
Xpetra::TpetraMultiVector::assign
virtual void assign(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &rhs)
Implementation of the assignment operator (operator=); does a deep copy.
Definition: Xpetra_TpetraMultiVector_def.hpp:511
Xpetra::TpetraMultiVector::beginExport
void beginExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Export data into this object using an Import object ("reverse mode").
Definition: Xpetra_TpetraMultiVector_def.hpp:347
Xpetra::TpetraMultiVector::meanValue
void meanValue(const Teuchos::ArrayView< Scalar > &means) const
Compute mean (average) value of each vector in multi-vector. The outcome of this routine is undefined...
Definition: Xpetra_TpetraMultiVector_def.hpp:248
Xpetra::TpetraMultiVector::dot
void dot(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const Teuchos::ArrayView< Scalar > &dots) const
Compute dot product of each corresponding pair of vectors, dots[i] = this[i].dot(A[i]).
Definition: Xpetra_TpetraMultiVector_def.hpp:193
Xpetra::TpetraMultiVector::putScalar
void putScalar(const Scalar &value)
Set all values in the multivector with the given value.
Definition: Xpetra_TpetraMultiVector_def.hpp:123
Xpetra::TpetraMultiVector::get1dViewNonConst
Teuchos::ArrayRCP< Scalar > get1dViewNonConst()
Nonconst persisting (1-D) view of this multivector's local values.
Definition: Xpetra_TpetraMultiVector_def.hpp:182
Xpetra::TpetraMultiVector::getLocalLength
size_t getLocalLength() const
Local number of rows on the calling process.
Definition: Xpetra_TpetraMultiVector_def.hpp:264
Xpetra::TpetraMultiVector::getNumVectors
size_t getNumVectors() const
Number of columns in the multivector.
Definition: Xpetra_TpetraMultiVector_def.hpp:259
Xpetra::TpetraMultiVector::sumIntoGlobalValue
void sumIntoGlobalValue(GlobalOrdinal globalRow, size_t vectorIndex, const Scalar &value)
Add value to existing value, using global (row) index.
Definition: Xpetra_TpetraMultiVector_def.hpp:108
Xpetra::TpetraMultiVector::setSeed
void setSeed(unsigned int seed)
Set seed for Random function.
Definition: Xpetra_TpetraMultiVector_def.hpp:455
Xpetra::TpetraMultiVector::reduce
void reduce()
Sum values of a locally replicated multivector across all processes.
Definition: Xpetra_TpetraMultiVector_def.hpp:128
Xpetra::TpetraMultiVector::replaceLocalValue
void replaceLocalValue(LocalOrdinal myRow, size_t vectorIndex, const Scalar &value)
Replace value, using local (row) index.
Definition: Xpetra_TpetraMultiVector_def.hpp:113
Xpetra::TpetraMultiVector::norm2
void norm2(const Teuchos::ArrayView< typename Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const
Definition: Xpetra_TpetraMultiVector_def.hpp:238
Xpetra::TpetraMultiVector::getVector
Teuchos::RCP< const Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > getVector(size_t j) const
Return a Vector which is a const view of column j.
Definition: Xpetra_TpetraMultiVector_def.hpp:134
Xpetra::TpetraMultiVector::endExport
void endExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Export data into this object using an Import object ("reverse mode").
Definition: Xpetra_TpetraMultiVector_def.hpp:358
Xpetra::TpetraMultiVector::description
std::string description() const
A simple one-line description of this object.
Definition: Xpetra_TpetraMultiVector_def.hpp:279
Xpetra::TpetraMultiVector::scale
void scale(const Scalar &alpha)
Scale the current values of a multi-vector, this = alpha*this.
Definition: Xpetra_TpetraMultiVector_def.hpp:208
Xpetra::TpetraMultiVector::TpetraMultiVector
TpetraMultiVector(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map, size_t NumVectors, bool zeroOut=true)
Basic constuctor.
Definition: Xpetra_TpetraMultiVector_def.hpp:64
Xpetra::TpetraMultiVector::update
void update(const Scalar &alpha, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const Scalar &beta)
Update multi-vector values with scaled values of A, this = beta*this + alpha*A.
Definition: Xpetra_TpetraMultiVector_def.hpp:223
Xpetra::TpetraMultiVector::getGlobalLength
global_size_t getGlobalLength() const
Global number of rows in the multivector.
Definition: Xpetra_TpetraMultiVector_def.hpp:269
Xpetra::TpetraMultiVector::reciprocal
void reciprocal(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A)
Put element-wise reciprocal values of input Multi-vector in target, this(i,j) = 1/A(i,...
Definition: Xpetra_TpetraMultiVector_def.hpp:203
Xpetra::TpetraMultiVector::~TpetraMultiVector
virtual ~TpetraMultiVector()
Destructor (virtual for memory safety of derived classes).
Definition: Xpetra_TpetraMultiVector_def.hpp:98
Xpetra::TpetraMultiVector::endImport
void endImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Import data into this object using an Import object ("forward mode").
Definition: Xpetra_TpetraMultiVector_def.hpp:325
Xpetra::TpetraMultiVector::abs
void abs(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A)
Put element-wise absolute values of input Multi-vector in target: A = abs(this).
Definition: Xpetra_TpetraMultiVector_def.hpp:198
Xpetra::TpetraMultiVector::elementWiseMultiply
void elementWiseMultiply(Scalar scalarAB, const Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &B, Scalar scalarThis)
Element-wise multiply of a Vector A with a TpetraMultiVector B.
Definition: Xpetra_TpetraMultiVector_def.hpp:1047
Xpetra::TpetraMultiVector::sumIntoLocalValue
void sumIntoLocalValue(LocalOrdinal myRow, size_t vectorIndex, const Scalar &value)
Add value to existing value, using local (row) index.
Definition: Xpetra_TpetraMultiVector_def.hpp:118
Xpetra::TpetraMultiVector::get1dView
Teuchos::ArrayRCP< const Scalar > get1dView() const
Const persisting (1-D) view of this multivector's local values.
Definition: Xpetra_TpetraMultiVector_def.hpp:170
Xpetra::TpetraMultiVector::getData
Teuchos::ArrayRCP< const Scalar > getData(size_t j) const
Const view of the local values in a particular vector of this multivector.
Definition: Xpetra_TpetraMultiVector_def.hpp:146
Xpetra::TpetraMultiVector::getTpetra_MultiVector
RCP< Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > getTpetra_MultiVector() const
Get the underlying Tpetra multivector.
Definition: Xpetra_TpetraMultiVector_def.hpp:450
Xpetra::TpetraMultiVector::getDataNonConst
Teuchos::ArrayRCP< Scalar > getDataNonConst(size_t j)
View of the local values in a particular vector of this multivector.
Definition: Xpetra_TpetraMultiVector_def.hpp:152
Xpetra::TpetraMultiVector::get2dCopy
void get2dCopy(Teuchos::ArrayView< const Teuchos::ArrayView< Scalar > > ArrayOfPtrs) const
Fill the given array with a copy of this multivector's local values.
Definition: Xpetra_TpetraMultiVector_def.hpp:164
Xpetra::TpetraMultiVector::beginImport
void beginImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Import data into this object using an Import object ("forward mode").
Definition: Xpetra_TpetraMultiVector_def.hpp:315
Xpetra::TpetraMultiVector::get2dViewNonConst
Teuchos::ArrayRCP< Teuchos::ArrayRCP< Scalar > > get2dViewNonConst()
Return non-const persisting pointers to values.
Definition: Xpetra_TpetraMultiVector_def.hpp:188
Xpetra::TpetraMultiVector::getMap
Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > getMap() const
The Map describing the parallel distribution of this object.
Definition: Xpetra_TpetraMultiVector_def.hpp:301
Xpetra::TpetraMultiVector::norm1
void norm1(const Teuchos::ArrayView< typename Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const
Compute 1-norm of each vector in multi-vector.
Definition: Xpetra_TpetraMultiVector_def.hpp:233
Xpetra::TpetraMultiVector::normInf
void normInf(const Teuchos::ArrayView< typename Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const
Compute Inf-norm of each vector in multi-vector.
Definition: Xpetra_TpetraMultiVector_def.hpp:243
Xpetra::TpetraMultiVector::randomize
void randomize(bool bUseXpetraImplementation=false)
Set multi-vector values to random numbers.
Definition: Xpetra_TpetraMultiVector_def.hpp:289
Xpetra::TpetraMultiVector::get1dCopy
void get1dCopy(Teuchos::ArrayView< Scalar > A, size_t LDA) const
Fill the given array with a copy of this multivector's local values.
Definition: Xpetra_TpetraMultiVector_def.hpp:158
Xpetra::TpetraMultiVector::get2dView
Teuchos::ArrayRCP< Teuchos::ArrayRCP< const Scalar > > get2dView() const
Return const persisting pointers to values.
Definition: Xpetra_TpetraMultiVector_def.hpp:176
TEUCHOS_TEST_FOR_EXCEPTION
#define TEUCHOS_TEST_FOR_EXCEPTION(throw_exception_test, Exception, msg)
Teuchos::EVerbosityLevel
EVerbosityLevel
rcp
TEUCHOS_DEPRECATED RCP< T > rcp(T *p, Dealloc_T dealloc, bool owns_mem)
Xpetra::Access::ReadWrite
constexpr struct ReadWriteStruct ReadWrite
Definition: Xpetra_Access.hpp:65
Xpetra::Access::ReadOnly
constexpr struct ReadOnlyStruct ReadOnly
Definition: Xpetra_Access.hpp:61
Xpetra::Access::OverwriteAll
constexpr struct OverwriteAllStruct OverwriteAll
Definition: Xpetra_Access.hpp:63
Xpetra
Xpetra namespace
Definition: Xpetra_BlockedCrsMatrix.hpp:89
Xpetra::global_size_t
size_t global_size_t
Global size_t object.
Definition: Xpetra_ConfigDefs.hpp:174
Xpetra::toXpetra
RCP< const CrsGraph< int, GlobalOrdinal, Node > > toXpetra(const Epetra_CrsGraph &g)
Definition: Xpetra_EpetraCrsGraph.cpp:65
Xpetra::toTpetra
RCP< const Tpetra::CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > toTpetra(const RCP< const CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > &graph)
Definition: Xpetra_TpetraCrsGraph_decl.hpp:425
Xpetra::CombineMode
CombineMode
Xpetra::Combine Mode enumerable type.
Definition: Xpetra_ConfigDefs.hpp:212
Teuchos::ScalarTraits::seedrandom
static void seedrandom(unsigned int s)