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>
511TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
512 TpetraMultiVector(const Teuchos::RCP<Tpetra::MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> > &vec)
513 : vec_(vec) {} // TODO removed const
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(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(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(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(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(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(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#ifdef HAVE_XPETRA_EPETRA
612
613#if ((defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_OPENMP) || !defined(HAVE_TPETRA_INST_INT_INT))) || \
614 (!defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_SERIAL) || !defined(HAVE_TPETRA_INST_INT_INT))))
615
616// specialization for TpetraMultiVector on EpetraNode and GO=int
617template <class Scalar>
618class TpetraMultiVector<Scalar, int, int, EpetraNode>
619 : public virtual MultiVector<Scalar, int, int, EpetraNode> {
620 typedef int LocalOrdinal;
621 typedef int GlobalOrdinal;
622 typedef EpetraNode Node;
623
624 // The following typedef are used by the XPETRA_DYNAMIC_CAST() macro.
625 typedef TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> TpetraMultiVectorClass;
626
627 public:
629
630
635
640
645
650
655
657 virtual ~TpetraMultiVector() {}
658
660
662
663
665 void replaceGlobalValue(GlobalOrdinal globalRow, size_t vectorIndex, const Scalar &value) {}
666
668 void sumIntoGlobalValue(GlobalOrdinal globalRow, size_t vectorIndex, const Scalar &value) {}
669
671 void replaceLocalValue(LocalOrdinal myRow, size_t vectorIndex, const Scalar &value) {}
672
674 void sumIntoLocalValue(LocalOrdinal myRow, size_t vectorIndex, const Scalar &value) {}
675
677 void putScalar(const Scalar &value) {}
678
680 void reduce() {}
681
683
685
686
688 Teuchos::RCP<const Vector<Scalar, LocalOrdinal, GlobalOrdinal, Node> > getVector(size_t j) const { return Teuchos::null; }
689
691 Teuchos::RCP<Vector<Scalar, LocalOrdinal, GlobalOrdinal, Node> > getVectorNonConst(size_t j) { return Teuchos::null; }
692
694 Teuchos::ArrayRCP<const Scalar> getData(size_t j) const { return Teuchos::ArrayRCP<const Scalar>(); }
695
697 Teuchos::ArrayRCP<Scalar> getDataNonConst(size_t j) { return Teuchos::ArrayRCP<Scalar>(); }
698
700 void get1dCopy(Teuchos::ArrayView<Scalar> A, size_t LDA) const {}
701
703 void get2dCopy(Teuchos::ArrayView<const Teuchos::ArrayView<Scalar> > ArrayOfPtrs) const {}
704
706 Teuchos::ArrayRCP<const Scalar> get1dView() const { return Teuchos::ArrayRCP<const Scalar>(); }
707
709 Teuchos::ArrayRCP<Teuchos::ArrayRCP<const Scalar> > get2dView() const { return Teuchos::ArrayRCP<Teuchos::ArrayRCP<const Scalar> >(); }
710
712 Teuchos::ArrayRCP<Scalar> get1dViewNonConst() { return Teuchos::ArrayRCP<Scalar>(); }
713
715 Teuchos::ArrayRCP<Teuchos::ArrayRCP<Scalar> > get2dViewNonConst() { return Teuchos::ArrayRCP<Teuchos::ArrayRCP<Scalar> >(); }
716
718
720
721
723 void dot(const MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &A, const Teuchos::ArrayView<Scalar> &dots) const {}
724
726 void abs(const MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &A) {}
727
729 void reciprocal(const MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &A) {}
730
732 void scale(const Scalar &alpha) {}
733
735 void scale(Teuchos::ArrayView<const Scalar> alpha) {}
736
738 void scale(const Scalar &alpha, const MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &A) {}
739
741 void update(const Scalar &alpha, const MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &A, const Scalar &beta) {}
742
744 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) {}
745
747 void norm1(const Teuchos::ArrayView<typename Teuchos::ScalarTraits<Scalar>::magnitudeType> &norms) const {}
748
750 void norm2(const Teuchos::ArrayView<typename Teuchos::ScalarTraits<Scalar>::magnitudeType> &norms) const {}
751
753 void normInf(const Teuchos::ArrayView<typename Teuchos::ScalarTraits<Scalar>::magnitudeType> &norms) const {}
754
756 void meanValue(const Teuchos::ArrayView<Scalar> &means) const {}
757
759 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) {}
760
762
764
765
767 size_t getNumVectors() const { return 0; }
768
770 size_t getLocalLength() const { return 0; }
771
773 global_size_t getGlobalLength() const { return 0; }
774
775 // \! Checks to see if the local length, number of vectors and size of Scalar type match
776 bool isSameSize(const MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &vec) const { return false; }
777
779
781
782
784 std::string description() const { return std::string(""); }
785
787 void describe(Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel = Teuchos::Describable::verbLevel_default) const {}
788
790
792 void elementWiseMultiply(Scalar scalarAB, const Vector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &A, const MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &B, Scalar scalarThis){};
793
795 void randomize(bool bUseXpetraImplementation = false) {}
796
798 void randomize(const Scalar &minVal, const Scalar &maxVal, bool bUseXpetraImplementation = false) {}
799
800 //{@
801 // Implements DistObject interface
802
803 Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > getMap() const { return Teuchos::null; }
804
805 void doImport(const DistObject<Scalar, LocalOrdinal, GlobalOrdinal, Node> &source, const Import<LocalOrdinal, GlobalOrdinal, Node> &importer, CombineMode CM) {}
806
807 void beginImport(const DistObject<Scalar, LocalOrdinal, GlobalOrdinal, Node> &source, const Import<LocalOrdinal, GlobalOrdinal, Node> &importer, CombineMode CM) {}
808
809 void endImport(const DistObject<Scalar, LocalOrdinal, GlobalOrdinal, Node> &source, const Import<LocalOrdinal, GlobalOrdinal, Node> &importer, CombineMode CM) {}
810
811 void doExport(const DistObject<Scalar, LocalOrdinal, GlobalOrdinal, Node> &dest, const Import<LocalOrdinal, GlobalOrdinal, Node> &importer, CombineMode CM) {}
812
813 void beginExport(const DistObject<Scalar, LocalOrdinal, GlobalOrdinal, Node> &dest, const Import<LocalOrdinal, GlobalOrdinal, Node> &importer, CombineMode CM) {}
814
815 void endExport(const DistObject<Scalar, LocalOrdinal, GlobalOrdinal, Node> &dest, const Import<LocalOrdinal, GlobalOrdinal, Node> &importer, CombineMode CM) {}
816
817 void doImport(const DistObject<Scalar, LocalOrdinal, GlobalOrdinal, Node> &source, const Export<LocalOrdinal, GlobalOrdinal, Node> &exporter, CombineMode CM) {}
818
819 void beginImport(const DistObject<Scalar, LocalOrdinal, GlobalOrdinal, Node> &source, const Export<LocalOrdinal, GlobalOrdinal, Node> &exporter, CombineMode CM) {}
820
821 void endImport(const DistObject<Scalar, LocalOrdinal, GlobalOrdinal, Node> &source, const Export<LocalOrdinal, GlobalOrdinal, Node> &exporter, CombineMode CM) {}
822
823 void doExport(const DistObject<Scalar, LocalOrdinal, GlobalOrdinal, Node> &dest, const Export<LocalOrdinal, GlobalOrdinal, Node> &exporter, CombineMode CM) {}
824
825 void beginExport(const DistObject<Scalar, LocalOrdinal, GlobalOrdinal, Node> &dest, const Export<LocalOrdinal, GlobalOrdinal, Node> &exporter, CombineMode CM) {}
826
827 void endExport(const DistObject<Scalar, LocalOrdinal, GlobalOrdinal, Node> &dest, const Export<LocalOrdinal, GlobalOrdinal, Node> &exporter, CombineMode CM) {}
828
829 void replaceMap(const RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &map) {}
830
832
834
835
837 TpetraMultiVector(const Teuchos::RCP<Tpetra::MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> > &vec) {
838 XPETRA_TPETRA_ETI_EXCEPTION(typeid(TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), typeid(TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), "int", typeid(EpetraNode).name());
839 }
840
842 RCP<Tpetra::MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> > getTpetra_MultiVector() const { return Teuchos::null; }
843
845 void setSeed(unsigned int seed) {}
846
848
849 protected:
852 virtual void
853 assign(const MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &rhs) {}
854}; // TpetraMultiVector class (specialization GO=int, NO=EpetraNode)
855#endif
856
857#if ((defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_OPENMP) || !defined(HAVE_TPETRA_INST_INT_LONG_LONG))) || \
858 (!defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_SERIAL) || !defined(HAVE_TPETRA_INST_INT_LONG_LONG))))
859
860// specialization for TpetraMultiVector on EpetraNode and GO=long long
861template <class Scalar>
862class TpetraMultiVector<Scalar, int, long long, EpetraNode>
863 : public virtual MultiVector<Scalar, int, long long, EpetraNode> {
864 typedef int LocalOrdinal;
865 typedef long long GlobalOrdinal;
866 typedef EpetraNode Node;
867
868 // The following typedef are used by the XPETRA_DYNAMIC_CAST() macro.
869 typedef TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> TpetraMultiVectorClass;
870
871 public:
873
874
879
884
889
894
896 virtual ~TpetraMultiVector() {}
897
899
901
902
904 void replaceGlobalValue(GlobalOrdinal globalRow, size_t vectorIndex, const Scalar &value) {}
905
907 void sumIntoGlobalValue(GlobalOrdinal globalRow, size_t vectorIndex, const Scalar &value) {}
908
910 void replaceLocalValue(LocalOrdinal myRow, size_t vectorIndex, const Scalar &value) {}
911
913 void sumIntoLocalValue(LocalOrdinal myRow, size_t vectorIndex, const Scalar &value) {}
914
916 void putScalar(const Scalar &value) {}
917
919 void reduce() {}
920
922
924
925
927 Teuchos::RCP<const Vector<Scalar, LocalOrdinal, GlobalOrdinal, Node> > getVector(size_t j) const { return Teuchos::null; }
928
930 Teuchos::RCP<Vector<Scalar, LocalOrdinal, GlobalOrdinal, Node> > getVectorNonConst(size_t j) { return Teuchos::null; }
931
933 Teuchos::ArrayRCP<const Scalar> getData(size_t j) const { return Teuchos::ArrayRCP<const Scalar>(); }
934
936 Teuchos::ArrayRCP<Scalar> getDataNonConst(size_t j) { return Teuchos::ArrayRCP<Scalar>(); }
937
939 void get1dCopy(Teuchos::ArrayView<Scalar> A, size_t LDA) const {}
940
942 void get2dCopy(Teuchos::ArrayView<const Teuchos::ArrayView<Scalar> > ArrayOfPtrs) const {}
943
945 Teuchos::ArrayRCP<const Scalar> get1dView() const { return Teuchos::ArrayRCP<const Scalar>(); }
946
948 Teuchos::ArrayRCP<Teuchos::ArrayRCP<const Scalar> > get2dView() const { return Teuchos::ArrayRCP<Teuchos::ArrayRCP<const Scalar> >(); }
949
951 Teuchos::ArrayRCP<Scalar> get1dViewNonConst() { return Teuchos::ArrayRCP<Scalar>(); }
952
954 Teuchos::ArrayRCP<Teuchos::ArrayRCP<Scalar> > get2dViewNonConst() { return Teuchos::ArrayRCP<Teuchos::ArrayRCP<Scalar> >(); }
955
957
959
960
962 void dot(const MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &A, const Teuchos::ArrayView<Scalar> &dots) const {}
963
965 void abs(const MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &A) {}
966
968 void reciprocal(const MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &A) {}
969
971 void scale(const Scalar &alpha) {}
972
974 void scale(Teuchos::ArrayView<const Scalar> alpha) {}
975
977 void scale(const Scalar &alpha, const MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &A) {}
978
980 void update(const Scalar &alpha, const MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &A, const Scalar &beta) {}
981
983 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) {}
984
986 void norm1(const Teuchos::ArrayView<typename Teuchos::ScalarTraits<Scalar>::magnitudeType> &norms) const {}
987
989 void norm2(const Teuchos::ArrayView<typename Teuchos::ScalarTraits<Scalar>::magnitudeType> &norms) const {}
990
992 void normInf(const Teuchos::ArrayView<typename Teuchos::ScalarTraits<Scalar>::magnitudeType> &norms) const {}
993
995 void meanValue(const Teuchos::ArrayView<Scalar> &means) const {}
996
998 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) {}
999
1001
1003
1004
1006 size_t getNumVectors() const { return 0; }
1007
1009 size_t getLocalLength() const { return 0; }
1010
1012 global_size_t getGlobalLength() const { return 0; }
1013
1014 // \! Checks to see if the local length, number of vectors and size of Scalar type match
1015 bool isSameSize(const MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &vec) const { return false; }
1016
1018
1020
1021
1023 std::string description() const { return std::string(""); }
1024
1026 void describe(Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel = Teuchos::Describable::verbLevel_default) const {}
1027
1029
1031 void elementWiseMultiply(Scalar scalarAB, const Vector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &A, const MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &B, Scalar scalarThis){};
1032
1034 void randomize(bool bUseXpetraImplementation = false) {}
1035
1037 void randomize(const Scalar &minVal, const Scalar &maxVal, bool bUseXpetraImplementation = false) {}
1038
1039 //{@
1040 // Implements DistObject interface
1041
1042 Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > getMap() const { return Teuchos::null; }
1043
1044 void doImport(const DistObject<Scalar, LocalOrdinal, GlobalOrdinal, Node> &source, const Import<LocalOrdinal, GlobalOrdinal, Node> &importer, CombineMode CM) {}
1045
1046 void beginImport(const DistObject<Scalar, LocalOrdinal, GlobalOrdinal, Node> &source, const Import<LocalOrdinal, GlobalOrdinal, Node> &importer, CombineMode CM) {}
1047
1048 void endImport(const DistObject<Scalar, LocalOrdinal, GlobalOrdinal, Node> &source, const Import<LocalOrdinal, GlobalOrdinal, Node> &importer, CombineMode CM) {}
1049
1050 void doExport(const DistObject<Scalar, LocalOrdinal, GlobalOrdinal, Node> &dest, const Import<LocalOrdinal, GlobalOrdinal, Node> &importer, CombineMode CM) {}
1051
1052 void beginExport(const DistObject<Scalar, LocalOrdinal, GlobalOrdinal, Node> &dest, const Import<LocalOrdinal, GlobalOrdinal, Node> &importer, CombineMode CM) {}
1053
1054 void endExport(const DistObject<Scalar, LocalOrdinal, GlobalOrdinal, Node> &dest, const Import<LocalOrdinal, GlobalOrdinal, Node> &importer, CombineMode CM) {}
1055
1056 void doImport(const DistObject<Scalar, LocalOrdinal, GlobalOrdinal, Node> &source, const Export<LocalOrdinal, GlobalOrdinal, Node> &exporter, CombineMode CM) {}
1057
1058 void beginImport(const DistObject<Scalar, LocalOrdinal, GlobalOrdinal, Node> &source, const Export<LocalOrdinal, GlobalOrdinal, Node> &exporter, CombineMode CM) {}
1059
1060 void endImport(const DistObject<Scalar, LocalOrdinal, GlobalOrdinal, Node> &source, const Export<LocalOrdinal, GlobalOrdinal, Node> &exporter, CombineMode CM) {}
1061
1062 void doExport(const DistObject<Scalar, LocalOrdinal, GlobalOrdinal, Node> &dest, const Export<LocalOrdinal, GlobalOrdinal, Node> &exporter, CombineMode CM) {}
1063
1064 void beginExport(const DistObject<Scalar, LocalOrdinal, GlobalOrdinal, Node> &dest, const Export<LocalOrdinal, GlobalOrdinal, Node> &exporter, CombineMode CM) {}
1065
1066 void endExport(const DistObject<Scalar, LocalOrdinal, GlobalOrdinal, Node> &dest, const Export<LocalOrdinal, GlobalOrdinal, Node> &exporter, CombineMode CM) {}
1067
1068 void replaceMap(const RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &map) {}
1069
1071
1073
1074
1076 TpetraMultiVector(const Teuchos::RCP<Tpetra::MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> > &vec) {
1077 XPETRA_TPETRA_ETI_EXCEPTION(typeid(TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), typeid(TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, EpetraNode>).name(), "long long", typeid(EpetraNode).name());
1078 }
1079
1081 RCP<Tpetra::MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> > getTpetra_MultiVector() const { return Teuchos::null; }
1082
1084 void setSeed(unsigned int seed) {}
1085
1087
1088 protected:
1091 virtual void
1092 assign(const MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &rhs) {}
1093}; // TpetraMultiVector class (specialization GO=int, NO=EpetraNode)
1094
1095#endif // TpetraMultiVector class (specialization GO=long long, NO=EpetraNode)
1096
1097#endif // HAVE_XPETRA_EPETRA
1098
1099} // namespace Xpetra
1100
1101// Following header file inculsion is needed for the dynamic_cast to TpetraVector in
1102// elementWiseMultiply (because we cannot dynamic_cast if target is not a complete type)
1103// It is included here to avoid circular dependency between Vector and MultiVector
1104// TODO: there is certainly a more elegant solution...
1105#include "Xpetra_TpetraVector.hpp"
1106
1107namespace Xpetra {
1108
1109template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
1110void TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::
1111 elementWiseMultiply(Scalar scalarAB,
1112 const Vector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &A,
1113 const MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> &B,
1114 Scalar scalarThis) {
1115 XPETRA_MONITOR("TpetraMultiVector::elementWiseMultiply");
1116
1117 // XPETRA_DYNAMIC_CAST won't take TpetraVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>
1118 // as an argument, hence the following typedef.
1119 typedef TpetraVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> tpv;
1120 XPETRA_DYNAMIC_CAST(const tpv, A, tA, "Xpetra::TpetraMultiVectorMatrix->multiply() only accept Xpetra::TpetraMultiVector as input arguments.");
1121 XPETRA_DYNAMIC_CAST(const TpetraMultiVector, B, tB, "Xpetra::TpetraMultiVectorMatrix->multiply() only accept Xpetra::TpetraMultiVector as input arguments.");
1122 vec_->elementWiseMultiply(scalarAB, *tA.getTpetra_Vector(), *tB.getTpetra_MultiVector(), scalarThis);
1123}
1124
1125} // namespace Xpetra
1126
1127#define XPETRA_TPETRAMULTIVECTOR_SHORT
1128#endif // XPETRA_TPETRAMULTIVECTOR_HPP
XPETRA_MONITOR
#define XPETRA_MONITOR(funcName)
Definition Xpetra_ConfigDefs.hpp:96
XPETRA_DYNAMIC_CAST
#define XPETRA_DYNAMIC_CAST(type, obj, newObj, exceptionMsg)
Definition Xpetra_Exceptions.hpp:22
XPETRA_TPETRA_ETI_EXCEPTION
#define XPETRA_TPETRA_ETI_EXCEPTION(cl, obj, go, node)
Definition Xpetra_Exceptions.hpp:43
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:32
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::TpetraMultiVectorClass
TpetraMultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > TpetraMultiVectorClass
Definition Xpetra_TpetraMultiVector_def.hpp:625
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::replaceMap
void replaceMap(const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map)
Definition Xpetra_TpetraMultiVector_def.hpp:829
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::endExport
void endExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Export< LocalOrdinal, GlobalOrdinal, Node > &exporter, CombineMode CM)
Definition Xpetra_TpetraMultiVector_def.hpp:827
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::beginImport
void beginImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Export< LocalOrdinal, GlobalOrdinal, Node > &exporter, CombineMode CM)
Definition Xpetra_TpetraMultiVector_def.hpp:819
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::endImport
void endImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Definition Xpetra_TpetraMultiVector_def.hpp:809
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:792
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:709
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:694
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:642
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::doImport
void doImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Definition Xpetra_TpetraMultiVector_def.hpp:805
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:668
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:652
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:837
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::beginExport
void beginExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Export< LocalOrdinal, GlobalOrdinal, Node > &exporter, CombineMode CM)
Definition Xpetra_TpetraMultiVector_def.hpp:825
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:803
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:735
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::getGlobalLength
global_size_t getGlobalLength() const
Global number of rows in the multivector.
Definition Xpetra_TpetraMultiVector_def.hpp:773
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:741
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::reduce
void reduce()
Sum values of a locally replicated multivector across all processes.
Definition Xpetra_TpetraMultiVector_def.hpp:680
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:853
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::endImport
void endImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Export< LocalOrdinal, GlobalOrdinal, Node > &exporter, CombineMode CM)
Definition Xpetra_TpetraMultiVector_def.hpp:821
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:787
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::norm2
void norm2(const Teuchos::ArrayView< typename Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const
Definition Xpetra_TpetraMultiVector_def.hpp:750
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:732
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:671
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::doExport
void doExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Export< LocalOrdinal, GlobalOrdinal, Node > &exporter, CombineMode CM)
Definition Xpetra_TpetraMultiVector_def.hpp:823
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::doExport
void doExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Definition Xpetra_TpetraMultiVector_def.hpp:811
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:637
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:723
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:691
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::~TpetraMultiVector
virtual ~TpetraMultiVector()
Destructor (virtual for memory safety of derived classes).
Definition Xpetra_TpetraMultiVector_def.hpp:657
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::doImport
void doImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Export< LocalOrdinal, GlobalOrdinal, Node > &exporter, CombineMode CM)
Definition Xpetra_TpetraMultiVector_def.hpp:817
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:842
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:674
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:759
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:697
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:747
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:706
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:677
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::getLocalLength
size_t getLocalLength() const
Local number of rows on the calling process.
Definition Xpetra_TpetraMultiVector_def.hpp:770
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::randomize
void randomize(const Scalar &minVal, const Scalar &maxVal, bool bUseXpetraImplementation=false)
Set multi-vector values to random numbers.
Definition Xpetra_TpetraMultiVector_def.hpp:798
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::isSameSize
bool isSameSize(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &vec) const
Definition Xpetra_TpetraMultiVector_def.hpp:776
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::beginExport
void beginExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Definition Xpetra_TpetraMultiVector_def.hpp:813
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::randomize
void randomize(bool bUseXpetraImplementation=false)
Set multi-vector values to random numbers.
Definition Xpetra_TpetraMultiVector_def.hpp:795
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::setSeed
void setSeed(unsigned int seed)
Set seed for Random function.
Definition Xpetra_TpetraMultiVector_def.hpp:845
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::description
std::string description() const
A simple one-line description of this object.
Definition Xpetra_TpetraMultiVector_def.hpp:784
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:738
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:744
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:753
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:712
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:688
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:756
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::getNumVectors
size_t getNumVectors() const
Number of columns in the multivector.
Definition Xpetra_TpetraMultiVector_def.hpp:767
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:700
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:715
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::endExport
void endExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Definition Xpetra_TpetraMultiVector_def.hpp:815
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:726
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:665
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::beginImport
void beginImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Definition Xpetra_TpetraMultiVector_def.hpp:807
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:703
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:647
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:729
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)
Definition Xpetra_TpetraMultiVector_def.hpp:1058
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:954
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)
Definition Xpetra_TpetraMultiVector_def.hpp:1050
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:974
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:989
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)
Definition Xpetra_TpetraMultiVector_def.hpp:1044
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:1009
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:942
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:1023
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::randomize
void randomize(const Scalar &minVal, const Scalar &maxVal, bool bUseXpetraImplementation=false)
Set multi-vector values to random numbers.
Definition Xpetra_TpetraMultiVector_def.hpp:1037
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:927
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:1012
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:998
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:930
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:983
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:995
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:1026
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)
Definition Xpetra_TpetraMultiVector_def.hpp:1066
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:907
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:881
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)
Definition Xpetra_TpetraMultiVector_def.hpp:1048
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:948
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:1042
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)
Definition Xpetra_TpetraMultiVector_def.hpp:1046
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:992
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:965
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:980
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:913
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:919
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:891
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:945
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:1092
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:936
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:910
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:1076
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::isSameSize
bool isSameSize(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &vec) const
Definition Xpetra_TpetraMultiVector_def.hpp:1015
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:933
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)
Definition Xpetra_TpetraMultiVector_def.hpp:1064
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:886
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:1034
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:916
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:968
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)
Definition Xpetra_TpetraMultiVector_def.hpp:1056
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:904
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)
Definition Xpetra_TpetraMultiVector_def.hpp:1052
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:1081
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::TpetraMultiVectorClass
TpetraMultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > TpetraMultiVectorClass
Definition Xpetra_TpetraMultiVector_def.hpp:869
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:1031
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:971
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:876
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:986
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::~TpetraMultiVector
virtual ~TpetraMultiVector()
Destructor (virtual for memory safety of derived classes).
Definition Xpetra_TpetraMultiVector_def.hpp:896
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)
Definition Xpetra_TpetraMultiVector_def.hpp:1060
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:977
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:962
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)
Definition Xpetra_TpetraMultiVector_def.hpp:1054
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::replaceMap
void replaceMap(const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map)
Definition Xpetra_TpetraMultiVector_def.hpp:1068
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)
Definition Xpetra_TpetraMultiVector_def.hpp:1062
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:951
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:939
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)
Export data into this object using an Import object ("reverse mode").
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::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)
Import data into this object using an Import object ("forward mode").
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)
Export data into this object using an Import object ("reverse mode").
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::getLocalViewHost
virtual dual_view_type::t_host_const_um getLocalViewHost(Access::ReadOnlyStruct) const
Definition Xpetra_TpetraMultiVector_def.hpp:534
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)
Export data into this object using an Import object ("reverse mode").
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::getLocalViewDevice
virtual dual_view_type::t_dev_const_um getLocalViewDevice(Access::ReadOnlyStruct) const
Definition Xpetra_TpetraMultiVector_def.hpp:541
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)
Import data into this object using an Import object ("forward mode").
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:1111
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)
Import data into this object using an Import object ("forward mode").
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
Xpetra::EpetraNode
Tpetra::KokkosCompat::KokkosSerialWrapperNode EpetraNode
Definition Xpetra_Map_decl.hpp:35
Xpetra::global_size_t
size_t global_size_t
Global size_t object.
Definition Xpetra_ConfigDefs.hpp:138
Xpetra::toXpetra
RCP< const CrsGraph< int, GlobalOrdinal, Node > > toXpetra(const Epetra_CrsGraph &g)
Definition Xpetra_EpetraCrsGraph.cpp:29
Xpetra::toTpetra
RCP< const Tpetra::CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > toTpetra(const RCP< const CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > &graph)
Definition Xpetra_TpetraCrsGraph_decl.hpp:376
Xpetra::CombineMode
CombineMode
Xpetra::Combine Mode enumerable type.
Definition Xpetra_ConfigDefs.hpp:176
Teuchos::ScalarTraits::seedrandom
static void seedrandom(unsigned int s)