Xpetra Version of the Day
Loading...
Searching...
No Matches
Xpetra_TpetraMultiVector_def.hpp
Go to the documentation of this file.
1// @HEADER
2// *****************************************************************************
3// Xpetra: A linear algebra interface package
4//
5// Copyright 2012 NTESS and the Xpetra contributors.
6// SPDX-License-Identifier: BSD-3-Clause
7// *****************************************************************************
8// @HEADER
9
10#ifndef XPETRA_TPETRAMULTIVECTOR_DEF_HPP
11#define XPETRA_TPETRAMULTIVECTOR_DEF_HPP
12#include "Xpetra_TpetraConfigDefs.hpp"
13
14#include "Xpetra_TpetraMap.hpp" //TMP
15#include "Xpetra_Utils.hpp"
16#include "Xpetra_TpetraImport.hpp"
17#include "Xpetra_TpetraExport.hpp"
18
19#include "Xpetra_TpetraMultiVector_decl.hpp"
20#include "Tpetra_MultiVector.hpp"
21#include "Tpetra_Vector.hpp"
22#include "Tpetra_Details_Random.hpp"
23
24namespace Xpetra {
25
27template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
28TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::TpetraMultiVector(const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &map, size_t NumVectors, bool zeroOut)
29 : vec_(Teuchos::rcp(new Tpetra::MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>(toTpetra(map), NumVectors, zeroOut))) {
30 // TAW 1/30/2016: even though Tpetra allows numVecs == 0, Epetra does not. Introduce exception to keep behavior of Epetra and Tpetra consistent.
31 TEUCHOS_TEST_FOR_EXCEPTION(NumVectors < 1, std::invalid_argument, "Xpetra::TpetraMultiVector(map,numVecs,zeroOut): numVecs = " << NumVectors << " < 1.");
32}
33
35template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
36TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
37 TpetraMultiVector(const MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &source, const Teuchos::DataAccess copyOrView)
38 : vec_(Teuchos::rcp(new Tpetra::MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>(toTpetra(source), copyOrView))) {}
39
41template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
42TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
43 TpetraMultiVector(const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &map, const Teuchos::ArrayView<const Scalar> &A, size_t LDA, size_t NumVectors)
44 : vec_(Teuchos::rcp(new Tpetra::MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>(toTpetra(map), A, LDA, NumVectors))) {
45 // TAW 1/30/2016: even though Tpetra allows numVecs == 0, Epetra does not. Introduce exception to keep behavior of Epetra and Tpetra consistent.
46 TEUCHOS_TEST_FOR_EXCEPTION(NumVectors < 1, std::invalid_argument, "Xpetra::TpetraMultiVector(map,A,LDA,numVecs): numVecs = " << NumVectors << " < 1.");
47}
48
50template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
51TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
52 TpetraMultiVector(const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &map, const Teuchos::ArrayView<const Teuchos::ArrayView<const Scalar> > &ArrayOfPtrs, size_t NumVectors)
53 : vec_(Teuchos::rcp(new Tpetra::MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>(toTpetra(map), ArrayOfPtrs, NumVectors))) {
54 // TAW 1/30/2016: even though Tpetra allows numVecs == 0, Epetra does not. Introduce exception to keep behavior of Epetra and Tpetra consistent.
55 TEUCHOS_TEST_FOR_EXCEPTION(NumVectors < 1, std::invalid_argument, "Xpetra::TpetraMultiVector(map,ArrayOfPtrs,numVecs): numVecs = " << NumVectors << " < 1.");
56}
57
59template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
62
64template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
65void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
66 replaceGlobalValue(GlobalOrdinal globalRow, size_t vectorIndex, const Scalar &value) {
67 XPETRA_MONITOR("TpetraMultiVector::replaceGlobalValue");
68 vec_->replaceGlobalValue(globalRow, vectorIndex, value);
69}
70
72template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
73void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
74 sumIntoGlobalValue(GlobalOrdinal globalRow, size_t vectorIndex, const Scalar &value) {
75 XPETRA_MONITOR("TpetraMultiVector::sumIntoGlobalValue");
76 vec_->sumIntoGlobalValue(globalRow, vectorIndex, value);
77}
78
80template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
81void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
82 replaceLocalValue(LocalOrdinal myRow, size_t vectorIndex, const Scalar &value) {
83 XPETRA_MONITOR("TpetraMultiVector::replaceLocalValue");
84 vec_->replaceLocalValue(myRow, vectorIndex, value);
85}
86
88template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
89void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
90 sumIntoLocalValue(LocalOrdinal myRow, size_t vectorIndex, const Scalar &value) {
91 XPETRA_MONITOR("TpetraMultiVector::sumIntoLocalValue");
92 vec_->sumIntoLocalValue(myRow, vectorIndex, value);
93}
94
96template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
97void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
98 putScalar(const Scalar &value) {
99 XPETRA_MONITOR("TpetraMultiVector::putScalar");
100 vec_->putScalar(value);
101}
102
104template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
105void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
106 reduce() {
107 XPETRA_MONITOR("TpetraMultiVector::reduce");
108 vec_->reduce();
109}
110
111template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
112Teuchos::RCP<const Vector<Scalar, LocalOrdinal, GlobalOrdinal, Node> >
113TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
114 getVector(size_t j) const {
115 XPETRA_MONITOR("TpetraMultiVector::getVector");
116 return toXpetra(vec_->getVector(j));
117}
118
120template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
121Teuchos::RCP<Vector<Scalar, LocalOrdinal, GlobalOrdinal, Node> >
122TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
123 getVectorNonConst(size_t j) {
124 XPETRA_MONITOR("TpetraMultiVector::getVectorNonConst");
125 return toXpetra(vec_->getVectorNonConst(j));
126}
127
129template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
130Teuchos::ArrayRCP<const Scalar>
131TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
132 getData(size_t j) const {
133 XPETRA_MONITOR("TpetraMultiVector::getData");
134 return vec_->getData(j);
135}
136
138template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
139Teuchos::ArrayRCP<Scalar>
140TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
141 getDataNonConst(size_t j) {
142 XPETRA_MONITOR("TpetraMultiVector::getDataNonConst");
143 return vec_->getDataNonConst(j);
144}
145
147template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
148void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
149 get1dCopy(Teuchos::ArrayView<Scalar> A, size_t LDA) const {
150 XPETRA_MONITOR("TpetraMultiVector::get1dCopy");
151 vec_->get1dCopy(A, LDA);
152}
153
155template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
156void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
157 get2dCopy(Teuchos::ArrayView<const Teuchos::ArrayView<Scalar> > ArrayOfPtrs) const {
158 XPETRA_MONITOR("TpetraMultiVector::get2dCopy");
159 vec_->get2dCopy(ArrayOfPtrs);
160}
161
163template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
164Teuchos::ArrayRCP<const Scalar>
165TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
166 get1dView() const {
167 XPETRA_MONITOR("TpetraMultiVector::get1dView");
168 return vec_->get1dView();
169}
170
172template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
173Teuchos::ArrayRCP<Teuchos::ArrayRCP<const Scalar> >
174TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
175 get2dView() const {
176 XPETRA_MONITOR("TpetraMultiVector::get2dView");
177 return vec_->get2dView();
178}
179
181template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
182Teuchos::ArrayRCP<Scalar>
183TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
184 get1dViewNonConst() {
185 XPETRA_MONITOR("TpetraMultiVector::get1dViewNonConst");
186 return vec_->get1dViewNonConst();
187}
188
190template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
191Teuchos::ArrayRCP<Teuchos::ArrayRCP<Scalar> >
192TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
193 get2dViewNonConst() {
194 XPETRA_MONITOR("TpetraMultiVector::get2dViewNonConst");
195 return vec_->get2dViewNonConst();
196}
197
199template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
200void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
201 dot(const MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &A, const Teuchos::ArrayView<Scalar> &dots) const {
202 XPETRA_MONITOR("TpetraMultiVector::dot");
203 vec_->dot(toTpetra(A), dots);
204}
205
207template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
213
215template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
216void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
217 reciprocal(const MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &A) {
218 XPETRA_MONITOR("TpetraMultiVector::reciprocal");
219 vec_->reciprocal(toTpetra(A));
220}
221
223template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
224void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
225 scale(const Scalar &alpha) {
226 XPETRA_MONITOR("TpetraMultiVector::scale");
227 vec_->scale(alpha);
228}
229
231template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
232void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
233 scale(Teuchos::ArrayView<const Scalar> alpha) {
234 XPETRA_MONITOR("TpetraMultiVector::scale");
235 vec_->scale(alpha);
236}
237
239template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
240void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
241 scale(const Scalar &alpha, const MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &A) {
242 XPETRA_MONITOR("TpetraMultiVector::scale");
243 vec_->scale(alpha, toTpetra(A));
244}
245
247template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
248void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
249 update(const Scalar &alpha, const MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &A, const Scalar &beta) {
250 XPETRA_MONITOR("TpetraMultiVector::update");
251 vec_->update(alpha, toTpetra(A), beta);
252}
253
255template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
256void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
257 update(const Scalar &alpha, const MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &A, const Scalar &beta, const MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &B, const Scalar &gamma) {
258 XPETRA_MONITOR("TpetraMultiVector::update");
259 vec_->update(alpha, toTpetra(A), beta, toTpetra(B), gamma);
260}
261
263template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
264void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
265 norm1(const Teuchos::ArrayView<typename Teuchos::ScalarTraits<Scalar>::magnitudeType> &norms) const {
266 XPETRA_MONITOR("TpetraMultiVector::norm1");
267 vec_->norm1(norms);
268}
269
271template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
272void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
273 norm2(const Teuchos::ArrayView<typename Teuchos::ScalarTraits<Scalar>::magnitudeType> &norms) const {
274 XPETRA_MONITOR("TpetraMultiVector::norm2");
275 vec_->norm2(norms);
276}
277
279template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
280void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
281 normInf(const Teuchos::ArrayView<typename Teuchos::ScalarTraits<Scalar>::magnitudeType> &norms) const {
282 XPETRA_MONITOR("TpetraMultiVector::normInf");
283 vec_->normInf(norms);
284}
285
287template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
288void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
289 meanValue(const Teuchos::ArrayView<Scalar> &means) const {
290 XPETRA_MONITOR("TpetraMultiVector::meanValue");
291 vec_->meanValue(means);
292}
293
295template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
296void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
297 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) {
298 XPETRA_MONITOR("TpetraMultiVector::multiply");
299 vec_->multiply(transA, transB, alpha, toTpetra(A), toTpetra(B), beta);
300}
301
303template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
304size_t TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
305 getNumVectors() const {
306 XPETRA_MONITOR("TpetraMultiVector::getNumVectors");
307 return vec_->getNumVectors();
308}
309
311template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
312size_t TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
313 getLocalLength() const {
314 XPETRA_MONITOR("TpetraMultiVector::getLocalLength");
315 return vec_->getLocalLength();
316}
317
319template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
320global_size_t TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
321 getGlobalLength() const {
322 XPETRA_MONITOR("TpetraMultiVector::getGlobalLength");
323 return vec_->getGlobalLength();
324}
325
326// \brief Checks to see if the local length, number of vectors and size of Scalar type match
327template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
328bool TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
329 isSameSize(const MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &vec) const {
330 XPETRA_MONITOR("TpetraMultiVector::isSameSize");
331 return vec_->isSameSize(toTpetra(vec));
332}
333
335template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
336std::string TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
337 description() const {
338 XPETRA_MONITOR("TpetraMultiVector::description");
339 return vec_->description();
340}
341
343template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
344void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
345 describe(Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel) const {
346 XPETRA_MONITOR("TpetraMultiVector::describe");
347 vec_->describe(out, verbLevel);
348}
349
351template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
352void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
353 randomize(bool bUseXpetraImplementation) {
354 XPETRA_MONITOR("TpetraMultiVector::randomize");
355
356 if (bUseXpetraImplementation)
357 MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::Xpetra_randomize();
358 else
359 vec_->randomize();
360}
361
363template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
364void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
365 randomize(const Scalar &minVal, const Scalar &maxVal, bool bUseXpetraImplementation) {
366 XPETRA_MONITOR("TpetraMultiVector::randomize");
367
368 if (bUseXpetraImplementation)
369 MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::Xpetra_randomize(minVal, maxVal);
370 else
371 vec_->randomize(minVal, maxVal);
372}
373
374template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
375Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> >
376TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
377 getMap() const {
378 XPETRA_MONITOR("TpetraMultiVector::getMap");
379 return toXpetra(vec_->getMap());
380}
381
382template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
383void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
384 doImport(const DistObject<Scalar, LocalOrdinal, GlobalOrdinal, Node> &source, const Import<LocalOrdinal, GlobalOrdinal, Node> &importer, CombineMode CM) {
385 XPETRA_MONITOR("TpetraMultiVector::doImport");
386
387 XPETRA_DYNAMIC_CAST(const TpetraMultiVectorClass, source, tSource, "Xpetra::TpetraMultiVector::doImport only accept Xpetra::TpetraMultiVector as input arguments."); // TODO: remove and use toTpetra()
388 RCP<const Tpetra::MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> > v = tSource.getTpetra_MultiVector();
389 this->getTpetra_MultiVector()->doImport(*v, toTpetra(importer), toTpetra(CM));
390}
391
392template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
393void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
394 beginImport(const DistObject<Scalar, LocalOrdinal, GlobalOrdinal, Node> &source, const Import<LocalOrdinal, GlobalOrdinal, Node> &importer, CombineMode CM) {
395 XPETRA_MONITOR("TpetraMultiVector::beginImport");
396
397 XPETRA_DYNAMIC_CAST(const TpetraMultiVectorClass, source, tSource, "Xpetra::TpetraMultiVector::doImport only accept Xpetra::TpetraMultiVector as input arguments."); // TODO: remove and use toTpetra()
398 RCP<const Tpetra::MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> > v = tSource.getTpetra_MultiVector();
399 this->getTpetra_MultiVector()->beginImport(*v, toTpetra(importer), toTpetra(CM));
400}
401
402template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
403void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
404 endImport(const DistObject<Scalar, LocalOrdinal, GlobalOrdinal, Node> &source, const Import<LocalOrdinal, GlobalOrdinal, Node> &importer, CombineMode CM) {
405 XPETRA_MONITOR("TpetraMultiVector::endImport");
406
407 XPETRA_DYNAMIC_CAST(const TpetraMultiVectorClass, source, tSource, "Xpetra::TpetraMultiVector::doImport only accept Xpetra::TpetraMultiVector as input arguments."); // TODO: remove and use toTpetra()
408 RCP<const Tpetra::MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> > v = tSource.getTpetra_MultiVector();
409 this->getTpetra_MultiVector()->endImport(*v, toTpetra(importer), toTpetra(CM));
410}
411
412template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
413void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
414 doExport(const DistObject<Scalar, LocalOrdinal, GlobalOrdinal, Node> &dest, const Import<LocalOrdinal, GlobalOrdinal, Node> &importer, CombineMode CM) {
415 XPETRA_MONITOR("TpetraMultiVector::doExport");
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()->doExport(*v, toTpetra(importer), toTpetra(CM));
420}
421
422template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
423void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
424 beginExport(const DistObject<Scalar, LocalOrdinal, GlobalOrdinal, Node> &dest, const Import<LocalOrdinal, GlobalOrdinal, Node> &importer, CombineMode CM) {
425 XPETRA_MONITOR("TpetraMultiVector::beginExport");
426
427 XPETRA_DYNAMIC_CAST(const TpetraMultiVectorClass, dest, tDest, "Xpetra::TpetraMultiVector::doImport only accept Xpetra::TpetraMultiVector as input arguments."); // TODO: remove and use toTpetra()
428 RCP<const Tpetra::MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> > v = tDest.getTpetra_MultiVector();
429 this->getTpetra_MultiVector()->beginExport(*v, toTpetra(importer), toTpetra(CM));
430}
431
432template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
433void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
434 endExport(const DistObject<Scalar, LocalOrdinal, GlobalOrdinal, Node> &dest, const Import<LocalOrdinal, GlobalOrdinal, Node> &importer, CombineMode CM) {
435 XPETRA_MONITOR("TpetraMultiVector::endExport");
436
437 XPETRA_DYNAMIC_CAST(const TpetraMultiVectorClass, dest, tDest, "Xpetra::TpetraMultiVector::doImport only accept Xpetra::TpetraMultiVector as input arguments."); // TODO: remove and use toTpetra()
438 RCP<const Tpetra::MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> > v = tDest.getTpetra_MultiVector();
439 this->getTpetra_MultiVector()->endExport(*v, toTpetra(importer), toTpetra(CM));
440}
441
442template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
443void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
444 doImport(const DistObject<Scalar, LocalOrdinal, GlobalOrdinal, Node> &source, const Export<LocalOrdinal, GlobalOrdinal, Node> &exporter, CombineMode CM) {
445 XPETRA_MONITOR("TpetraMultiVector::doImport");
446
447 XPETRA_DYNAMIC_CAST(const TpetraMultiVectorClass, source, tSource, "Xpetra::TpetraMultiVector::doImport only accept Xpetra::TpetraMultiVector as input arguments."); // TODO: remove and use toTpetra()
448 RCP<const Tpetra::MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> > v = tSource.getTpetra_MultiVector();
449 this->getTpetra_MultiVector()->doImport(*v, toTpetra(exporter), toTpetra(CM));
450}
451
452template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
453void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
454 beginImport(const DistObject<Scalar, LocalOrdinal, GlobalOrdinal, Node> &source, const Export<LocalOrdinal, GlobalOrdinal, Node> &exporter, CombineMode CM) {
455 XPETRA_MONITOR("TpetraMultiVector::beginImport");
456
457 XPETRA_DYNAMIC_CAST(const TpetraMultiVectorClass, source, tSource, "Xpetra::TpetraMultiVector::doImport only accept Xpetra::TpetraMultiVector as input arguments."); // TODO: remove and use toTpetra()
458 RCP<const Tpetra::MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> > v = tSource.getTpetra_MultiVector();
459 this->getTpetra_MultiVector()->beginImport(*v, toTpetra(exporter), toTpetra(CM));
460}
461
462template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
463void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
464 endImport(const DistObject<Scalar, LocalOrdinal, GlobalOrdinal, Node> &source, const Export<LocalOrdinal, GlobalOrdinal, Node> &exporter, CombineMode CM) {
465 XPETRA_MONITOR("TpetraMultiVector::endImport");
466
467 XPETRA_DYNAMIC_CAST(const TpetraMultiVectorClass, source, tSource, "Xpetra::TpetraMultiVector::doImport only accept Xpetra::TpetraMultiVector as input arguments."); // TODO: remove and use toTpetra()
468 RCP<const Tpetra::MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> > v = tSource.getTpetra_MultiVector();
469 this->getTpetra_MultiVector()->endImport(*v, toTpetra(exporter), toTpetra(CM));
470}
471
472template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
473void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
474 doExport(const DistObject<Scalar, LocalOrdinal, GlobalOrdinal, Node> &dest, const Export<LocalOrdinal, GlobalOrdinal, Node> &exporter, CombineMode CM) {
475 XPETRA_MONITOR("TpetraMultiVector::doExport");
476
477 XPETRA_DYNAMIC_CAST(const TpetraMultiVectorClass, dest, tDest, "Xpetra::TpetraMultiVector::doImport only accept Xpetra::TpetraMultiVector as input arguments."); // TODO: remove and use toTpetra()
478 RCP<const Tpetra::MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> > v = tDest.getTpetra_MultiVector();
479 this->getTpetra_MultiVector()->doExport(*v, toTpetra(exporter), toTpetra(CM));
480}
481
482template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
483void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
484 beginExport(const DistObject<Scalar, LocalOrdinal, GlobalOrdinal, Node> &dest, const Export<LocalOrdinal, GlobalOrdinal, Node> &exporter, CombineMode CM) {
485 XPETRA_MONITOR("TpetraMultiVector::beginExport");
486
487 XPETRA_DYNAMIC_CAST(const TpetraMultiVectorClass, dest, tDest, "Xpetra::TpetraMultiVector::doImport only accept Xpetra::TpetraMultiVector as input arguments."); // TODO: remove and use toTpetra()
488 RCP<const Tpetra::MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> > v = tDest.getTpetra_MultiVector();
489 this->getTpetra_MultiVector()->beginExport(*v, toTpetra(exporter), toTpetra(CM));
490}
491
492template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
493void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
494 endExport(const DistObject<Scalar, LocalOrdinal, GlobalOrdinal, Node> &dest, const Export<LocalOrdinal, GlobalOrdinal, Node> &exporter, CombineMode CM) {
495 XPETRA_MONITOR("TpetraMultiVector::endExport");
496
497 XPETRA_DYNAMIC_CAST(const TpetraMultiVectorClass, dest, tDest, "Xpetra::TpetraMultiVector::doImport only accept Xpetra::TpetraMultiVector as input arguments."); // TODO: remove and use toTpetra()
498 RCP<const Tpetra::MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> > v = tDest.getTpetra_MultiVector();
499 this->getTpetra_MultiVector()->endExport(*v, toTpetra(exporter), toTpetra(CM));
500}
501
502template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
503void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
504 replaceMap(const RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &map) {
505 XPETRA_MONITOR("TpetraMultiVector::replaceMap");
506 this->getTpetra_MultiVector()->replaceMap(toTpetra(map));
507}
508
510template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
514
516template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
517RCP<Tpetra::MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> >
520
522template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
523void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
524 setSeed(unsigned int seed) {
525 XPETRA_MONITOR("TpetraMultiVector::seedrandom");
526 Teuchos::ScalarTraits<Scalar>::seedrandom(seed);
527 // Tell Tpetra to update its RNG pool for the new random seed
528 Tpetra::Details::Static_Random_XorShift64_Pool<typename Node::device_type::execution_space>::resetPool(getMap()->getComm()->getRank());
529}
530
531template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
532typename TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::dual_view_type::t_host_const_um
533TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
534 getLocalViewHost(Tpetra::Access::ReadOnlyStruct) const {
535 return subview(vec_->getLocalViewHost(Tpetra::Access::ReadOnly), Kokkos::ALL(), Kokkos::ALL());
536}
537
538template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
539typename TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::dual_view_type::t_dev_const_um
540TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
541 getLocalViewDevice(Tpetra::Access::ReadOnlyStruct) const {
542 return subview(vec_->getLocalViewDevice(Tpetra::Access::ReadOnly), Kokkos::ALL(), Kokkos::ALL());
543}
544
545template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
546typename TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::dual_view_type::t_host_um
547TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
548 getLocalViewHost(Tpetra::Access::OverwriteAllStruct) const {
549 return subview(vec_->getLocalViewHost(Tpetra::Access::OverwriteAll), Kokkos::ALL(), Kokkos::ALL());
550}
551
552template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
553typename TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::dual_view_type::t_dev_um
554TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
555 getLocalViewDevice(Tpetra::Access::OverwriteAllStruct) const {
556 return subview(vec_->getLocalViewDevice(Tpetra::Access::OverwriteAll), Kokkos::ALL(), Kokkos::ALL());
557}
558
559template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
560typename TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::dual_view_type::t_host_um
561TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
562 getLocalViewHost(Tpetra::Access::ReadWriteStruct) const {
563 return subview(vec_->getLocalViewHost(Tpetra::Access::ReadWrite), Kokkos::ALL(), Kokkos::ALL());
564}
565
566template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
567typename TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::dual_view_type::t_dev_um
568TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
569 getLocalViewDevice(Tpetra::Access::ReadWriteStruct) const {
570 return subview(vec_->getLocalViewDevice(Tpetra::Access::ReadWrite), Kokkos::ALL(), Kokkos::ALL());
571}
572
575template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
576void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
577 assign(const MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &rhs) {
578 typedef TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> this_type;
579 const this_type *rhsPtr = dynamic_cast<const this_type *>(&rhs);
580 TEUCHOS_TEST_FOR_EXCEPTION(
581 rhsPtr == NULL, std::invalid_argument,
582 "Xpetra::MultiVector::operator=:"
583 " The left-hand side (LHS) of the assignment has a different type than "
584 "the right-hand side (RHS). The LHS has type Xpetra::TpetraMultiVector"
585 " (which means it wraps a Tpetra::MultiVector), but the RHS has some "
586 "other type. This probably means that the RHS wraps an "
587 "Epetra_MultiVector. Xpetra::MultiVector does not currently implement "
588 "assignment from an Epetra object to a Tpetra object, though this could"
589 " be added with sufficient interest.");
590
591 typedef Tpetra::MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> TMV;
592 RCP<const TMV> rhsImpl = rhsPtr->getTpetra_MultiVector();
593 RCP<TMV> lhsImpl = this->getTpetra_MultiVector();
594
595 TEUCHOS_TEST_FOR_EXCEPTION(
596 rhsImpl.is_null(), std::logic_error,
597 "Xpetra::MultiVector::operator= "
598 "(in Xpetra::TpetraMultiVector::assign): *this (the right-hand side of "
599 "the assignment) has a null RCP<Tpetra::MultiVector> inside. Please "
600 "report this bug to the Xpetra developers.");
601 TEUCHOS_TEST_FOR_EXCEPTION(
602 lhsImpl.is_null(), std::logic_error,
603 "Xpetra::MultiVector::operator= "
604 "(in Xpetra::TpetraMultiVector::assign): The left-hand side of the "
605 "assignment has a null RCP<Tpetra::MultiVector> inside. Please report "
606 "this bug to the Xpetra developers.");
607
608 Tpetra::deep_copy(*lhsImpl, *rhsImpl);
609}
610
611} // namespace Xpetra
612
613// Following header file inculsion is needed for the dynamic_cast to TpetraVector in
614// elementWiseMultiply (because we cannot dynamic_cast if target is not a complete type)
615// It is included here to avoid circular dependency between Vector and MultiVector
616// TODO: there is certainly a more elegant solution...
617#include "Xpetra_TpetraVector.hpp"
618
619namespace Xpetra {
620
621template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
622void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
623 elementWiseMultiply(Scalar scalarAB,
624 const Vector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &A,
625 const MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &B,
626 Scalar scalarThis) {
627 XPETRA_MONITOR("TpetraMultiVector::elementWiseMultiply");
628
629 // XPETRA_DYNAMIC_CAST won't take TpetraVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>
630 // as an argument, hence the following typedef.
631 typedef TpetraVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> tpv;
632 XPETRA_DYNAMIC_CAST(const tpv, A, tA, "Xpetra::TpetraMultiVectorMatrix->multiply() only accept Xpetra::TpetraMultiVector as input arguments.");
633 XPETRA_DYNAMIC_CAST(const TpetraMultiVector, B, tB, "Xpetra::TpetraMultiVectorMatrix->multiply() only accept Xpetra::TpetraMultiVector as input arguments.");
634 vec_->elementWiseMultiply(scalarAB, *tA.getTpetra_Vector(), *tB.getTpetra_MultiVector(), scalarThis);
635}
636
637} // namespace Xpetra
638
639#define XPETRA_TPETRAMULTIVECTOR_SHORT
640#endif // XPETRA_TPETRAMULTIVECTOR_HPP
XPETRA_MONITOR
#define XPETRA_MONITOR(funcName)
Definition Xpetra_ConfigDefs.hpp:94
XPETRA_DYNAMIC_CAST
#define XPETRA_DYNAMIC_CAST(type, obj, newObj, exceptionMsg)
Definition Xpetra_Exceptions.hpp:22
Xpetra_TpetraConfigDefs.hpp
Xpetra_TpetraMultiVector_decl.hpp
Xpetra_Utils.hpp
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:32
Xpetra::TpetraMultiVector::replaceMap
void replaceMap(const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map)
Definition Xpetra_TpetraMultiVector_def.hpp:504
Xpetra::TpetraMultiVector::replaceGlobalValue
void replaceGlobalValue(GlobalOrdinal globalRow, size_t vectorIndex, const Scalar &value)
Replace value, using global (row) index.
Definition Xpetra_TpetraMultiVector_def.hpp:66
Xpetra::TpetraMultiVector::doExport
void doExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Definition Xpetra_TpetraMultiVector_def.hpp:414
Xpetra::TpetraMultiVector::isSameSize
bool isSameSize(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &vec) const
Definition Xpetra_TpetraMultiVector_def.hpp:329
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:345
Xpetra::TpetraMultiVector::getLocalViewDevice
virtual dual_view_type::t_dev_const_um getLocalViewDevice(Tpetra::Access::ReadOnlyStruct) const
Definition Xpetra_TpetraMultiVector_def.hpp:541
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:297
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:123
Xpetra::TpetraMultiVector::doImport
void doImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Definition Xpetra_TpetraMultiVector_def.hpp:384
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:577
Xpetra::TpetraMultiVector::beginExport
void beginExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Definition Xpetra_TpetraMultiVector_def.hpp:424
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:289
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:201
Xpetra::TpetraMultiVector::putScalar
void putScalar(const Scalar &value)
Set all values in the multivector with the given value.
Definition Xpetra_TpetraMultiVector_def.hpp:98
Xpetra::TpetraMultiVector::get1dViewNonConst
Teuchos::ArrayRCP< Scalar > get1dViewNonConst()
Nonconst persisting (1-D) view of this multivector's local values.
Definition Xpetra_TpetraMultiVector_def.hpp:184
Xpetra::TpetraMultiVector::getLocalLength
size_t getLocalLength() const
Local number of rows on the calling process.
Definition Xpetra_TpetraMultiVector_def.hpp:313
Xpetra::TpetraMultiVector::getNumVectors
size_t getNumVectors() const
Number of columns in the multivector.
Definition Xpetra_TpetraMultiVector_def.hpp:305
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:74
Xpetra::TpetraMultiVector::setSeed
void setSeed(unsigned int seed)
Set seed for Random function.
Definition Xpetra_TpetraMultiVector_def.hpp:524
Xpetra::TpetraMultiVector::reduce
void reduce()
Sum values of a locally replicated multivector across all processes.
Definition Xpetra_TpetraMultiVector_def.hpp:106
Xpetra::TpetraMultiVector::replaceLocalValue
void replaceLocalValue(LocalOrdinal myRow, size_t vectorIndex, const Scalar &value)
Replace value, using local (row) index.
Definition Xpetra_TpetraMultiVector_def.hpp:82
Xpetra::TpetraMultiVector::norm2
void norm2(const Teuchos::ArrayView< typename Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const
Definition Xpetra_TpetraMultiVector_def.hpp:273
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:114
Xpetra::TpetraMultiVector::endExport
void endExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Definition Xpetra_TpetraMultiVector_def.hpp:434
Xpetra::TpetraMultiVector::description
std::string description() const
A simple one-line description of this object.
Definition Xpetra_TpetraMultiVector_def.hpp:337
Xpetra::TpetraMultiVector::scale
void scale(const Scalar &alpha)
Scale the current values of a multi-vector, this = alpha*this.
Definition Xpetra_TpetraMultiVector_def.hpp:225
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:28
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:249
Xpetra::TpetraMultiVector::getLocalViewHost
virtual dual_view_type::t_host_const_um getLocalViewHost(Tpetra::Access::ReadOnlyStruct) const
Definition Xpetra_TpetraMultiVector_def.hpp:534
Xpetra::TpetraMultiVector::getGlobalLength
global_size_t getGlobalLength() const
Global number of rows in the multivector.
Definition Xpetra_TpetraMultiVector_def.hpp:321
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:217
Xpetra::TpetraMultiVector::~TpetraMultiVector
virtual ~TpetraMultiVector()
Destructor (virtual for memory safety of derived classes).
Definition Xpetra_TpetraMultiVector_def.hpp:61
Xpetra::TpetraMultiVector::endImport
void endImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Definition Xpetra_TpetraMultiVector_def.hpp:404
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:209
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:623
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:90
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:166
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:132
Xpetra::TpetraMultiVector::getTpetra_MultiVector
RCP< Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > getTpetra_MultiVector() const
Get the underlying Tpetra multivector.
Definition Xpetra_TpetraMultiVector_def.hpp:519
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:141
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:157
Xpetra::TpetraMultiVector::beginImport
void beginImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Definition Xpetra_TpetraMultiVector_def.hpp:394
Xpetra::TpetraMultiVector::get2dViewNonConst
Teuchos::ArrayRCP< Teuchos::ArrayRCP< Scalar > > get2dViewNonConst()
Return non-const persisting pointers to values.
Definition Xpetra_TpetraMultiVector_def.hpp:193
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:377
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:265
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:281
Xpetra::TpetraMultiVector::randomize
void randomize(bool bUseXpetraImplementation=false)
Set multi-vector values to random numbers.
Definition Xpetra_TpetraMultiVector_def.hpp:353
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:149
Xpetra::TpetraMultiVector::get2dView
Teuchos::ArrayRCP< Teuchos::ArrayRCP< const Scalar > > get2dView() const
Return const persisting pointers to values.
Definition Xpetra_TpetraMultiVector_def.hpp:175
TEUCHOS_TEST_FOR_EXCEPTION
#define TEUCHOS_TEST_FOR_EXCEPTION(throw_exception_test, Exception, msg)
Teuchos::EVerbosityLevel
EVerbosityLevel
Tpetra::deep_copy
void deep_copy(MultiVector< DS, DL, DG, DN > &dst, const MultiVector< SS, SL, SG, SN > &src)
Xpetra::toXpetra
RCP< const CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > toXpetra(RCP< const Tpetra::CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > graph)
Definition Xpetra_TpetraCrsGraph_decl.hpp:361
Xpetra::global_size_t
size_t global_size_t
Global size_t object.
Definition Xpetra_ConfigDefs.hpp:136
Xpetra::toTpetra
RCP< const Tpetra::CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > toTpetra(const RCP< const CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > &graph)
Definition Xpetra_TpetraCrsGraph_decl.hpp:374
Xpetra::CombineMode
CombineMode
Xpetra::Combine Mode enumerable type.
Definition Xpetra_ConfigDefs.hpp:174
Teuchos::ScalarTraits::seedrandom
static void seedrandom(unsigned int s)