All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
Xpetra_TpetraMultiVector_def.hpp
Go to the documentation of this file.
1 // @HEADER
2 //
3 // ***********************************************************************
4 //
5 // Xpetra: A linear algebra interface package
6 // Copyright 2012 Sandia Corporation
7 //
8 // Under the terms of Contract DE-AC04-94AL85000 with Sandia Corporation,
9 // the U.S. Government retains certain rights in this software.
10 //
11 // Redistribution and use in source and binary forms, with or without
12 // modification, are permitted provided that the following conditions are
13 // met:
14 //
15 // 1. Redistributions of source code must retain the above copyright
16 // notice, this list of conditions and the following disclaimer.
17 //
18 // 2. Redistributions in binary form must reproduce the above copyright
19 // notice, this list of conditions and the following disclaimer in the
20 // documentation and/or other materials provided with the distribution.
21 //
22 // 3. Neither the name of the Corporation nor the names of the
23 // contributors may be used to endorse or promote products derived from
24 // this software without specific prior written permission.
25 //
26 // THIS SOFTWARE IS PROVIDED BY SANDIA CORPORATION "AS IS" AND ANY
27 // EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
28 // IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
29 // PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL SANDIA CORPORATION OR THE
30 // CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
31 // EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
32 // PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
33 // PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
34 // LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
35 // NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
36 // SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
37 //
38 // Questions? Contact
39 // Jonathan Hu (jhu@sandia.gov)
40 // Andrey Prokopenko (aprokop@sandia.gov)
41 // Ray Tuminaro (rstumin@sandia.gov)
42 //
43 // ***********************************************************************
44 //
45 // @HEADER
46 #ifndef XPETRA_TPETRAMULTIVECTOR_DEF_HPP
47 #define XPETRA_TPETRAMULTIVECTOR_DEF_HPP
48 #include "Xpetra_TpetraConfigDefs.hpp"
49 
50 #include "Xpetra_TpetraMap.hpp" //TMP
51 #include "Xpetra_Utils.hpp"
52 #include "Xpetra_TpetraImport.hpp"
53 #include "Xpetra_TpetraExport.hpp"
54 
55 #include "Xpetra_TpetraMultiVector_decl.hpp"
56 #include "Tpetra_MultiVector.hpp"
57 #include "Tpetra_Vector.hpp"
58 
59 namespace Xpetra {
60 
61 
63  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
64  TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::TpetraMultiVector(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map, size_t NumVectors, bool zeroOut)
65  : vec_(Teuchos::rcp(new Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node >(toTpetra(map), NumVectors, zeroOut))) {
66  // TAW 1/30/2016: even though Tpetra allows numVecs == 0, Epetra does not. Introduce exception to keep behavior of Epetra and Tpetra consistent.
67  TEUCHOS_TEST_FOR_EXCEPTION(NumVectors < 1, std::invalid_argument, "Xpetra::TpetraMultiVector(map,numVecs,zeroOut): numVecs = " << NumVectors << " < 1.");
68  }
69 
71  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
72  TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
73  TpetraMultiVector(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source)
74  : vec_(Teuchos::rcp(new Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node >(Tpetra::createCopy(toTpetra(source))))) { }
75 
77  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
78  TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
79  TpetraMultiVector(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map, const Teuchos::ArrayView< const Scalar > &A, size_t LDA, size_t NumVectors)
80  : vec_(Teuchos::rcp(new Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node >(toTpetra(map), A, LDA, NumVectors))) {
81  // TAW 1/30/2016: even though Tpetra allows numVecs == 0, Epetra does not. Introduce exception to keep behavior of Epetra and Tpetra consistent.
82  TEUCHOS_TEST_FOR_EXCEPTION(NumVectors < 1, std::invalid_argument, "Xpetra::TpetraMultiVector(map,A,LDA,numVecs): numVecs = " << NumVectors << " < 1.");
83  }
84 
86  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
87  TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
88  TpetraMultiVector(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map, const Teuchos::ArrayView< const Teuchos::ArrayView< const Scalar > > &ArrayOfPtrs, size_t NumVectors)
89  : vec_(Teuchos::rcp(new Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node >(toTpetra(map), ArrayOfPtrs, NumVectors))) {
90  // TAW 1/30/2016: even though Tpetra allows numVecs == 0, Epetra does not. Introduce exception to keep behavior of Epetra and Tpetra consistent.
91  TEUCHOS_TEST_FOR_EXCEPTION(NumVectors < 1, std::invalid_argument, "Xpetra::TpetraMultiVector(map,ArrayOfPtrs,numVecs): numVecs = " << NumVectors << " < 1.");
92  }
93 
94 
96  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
97  TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
98  ~TpetraMultiVector() { }
99 
101  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
102  void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
103  replaceGlobalValue(GlobalOrdinal globalRow, size_t vectorIndex, const Scalar &value) { XPETRA_MONITOR("TpetraMultiVector::replaceGlobalValue"); vec_->replaceGlobalValue(globalRow, vectorIndex, value); }
104 
106  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
107  void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
108  sumIntoGlobalValue(GlobalOrdinal globalRow, size_t vectorIndex, const Scalar &value) { XPETRA_MONITOR("TpetraMultiVector::sumIntoGlobalValue"); vec_->sumIntoGlobalValue(globalRow, vectorIndex, value); }
109 
111  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
112  void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
113  replaceLocalValue(LocalOrdinal myRow, size_t vectorIndex, const Scalar &value) { XPETRA_MONITOR("TpetraMultiVector::replaceLocalValue"); vec_->replaceLocalValue(myRow, vectorIndex, value); }
114 
116  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
117  void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
118  sumIntoLocalValue(LocalOrdinal myRow, size_t vectorIndex, const Scalar &value) { XPETRA_MONITOR("TpetraMultiVector::sumIntoLocalValue"); vec_->sumIntoLocalValue(myRow, vectorIndex, value); }
119 
121  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
122  void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
123  putScalar(const Scalar &value) { XPETRA_MONITOR("TpetraMultiVector::putScalar"); vec_->putScalar(value); }
124 
126  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
127  void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
128  reduce() { XPETRA_MONITOR("TpetraMultiVector::reduce"); vec_->reduce(); }
129 
130 
131  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
132  Teuchos::RCP< const Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > >
133  TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
134  getVector(size_t j) const { XPETRA_MONITOR("TpetraMultiVector::getVector"); return toXpetra(vec_->getVector(j)); }
135 
137  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
138  Teuchos::RCP< Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > >
139  TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
140  getVectorNonConst(size_t j) { XPETRA_MONITOR("TpetraMultiVector::getVectorNonConst"); return toXpetra(vec_->getVectorNonConst(j)); }
141 
143  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
144  Teuchos::ArrayRCP< const Scalar >
145  TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
146  getData(size_t j) const { XPETRA_MONITOR("TpetraMultiVector::getData"); return vec_->getData(j); }
147 
149  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
150  Teuchos::ArrayRCP< Scalar >
151  TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
152  getDataNonConst(size_t j) { XPETRA_MONITOR("TpetraMultiVector::getDataNonConst"); return vec_->getDataNonConst(j); }
153 
155  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
156  void
157  TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
158  get1dCopy(Teuchos::ArrayView< Scalar > A, size_t LDA) const { XPETRA_MONITOR("TpetraMultiVector::get1dCopy"); vec_->get1dCopy(A, LDA); }
159 
161  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
162  void
163  TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
164  get2dCopy(Teuchos::ArrayView< const Teuchos::ArrayView< Scalar > > ArrayOfPtrs) const { XPETRA_MONITOR("TpetraMultiVector::get2dCopy"); vec_->get2dCopy(ArrayOfPtrs); }
165 
167  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
168  Teuchos::ArrayRCP< const Scalar >
169  TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
170  get1dView() const { XPETRA_MONITOR("TpetraMultiVector::get1dView"); return vec_->get1dView(); }
171 
173  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
174  Teuchos::ArrayRCP< Teuchos::ArrayRCP< const Scalar > >
175  TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
176  get2dView() const { XPETRA_MONITOR("TpetraMultiVector::get2dView"); return vec_->get2dView(); }
177 
179  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
180  Teuchos::ArrayRCP< Scalar >
181  TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
182  get1dViewNonConst() { XPETRA_MONITOR("TpetraMultiVector::get1dViewNonConst"); return vec_->get1dViewNonConst(); }
183 
185  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
186  Teuchos::ArrayRCP< Teuchos::ArrayRCP< Scalar > >
187  TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
188  get2dViewNonConst() { XPETRA_MONITOR("TpetraMultiVector::get2dViewNonConst"); return vec_->get2dViewNonConst(); }
189 
191  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
192  void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
193  dot(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const Teuchos::ArrayView< Scalar > &dots) const { XPETRA_MONITOR("TpetraMultiVector::dot"); vec_->dot(toTpetra(A), dots); }
194 
196  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
197  void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
198  abs(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A) { XPETRA_MONITOR("TpetraMultiVector::abs"); vec_->abs(toTpetra(A)); }
199 
201  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
202  void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
203  reciprocal(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A) { XPETRA_MONITOR("TpetraMultiVector::reciprocal"); vec_->reciprocal(toTpetra(A)); }
204 
206  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
207  void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
208  scale(const Scalar &alpha) { XPETRA_MONITOR("TpetraMultiVector::scale"); vec_->scale(alpha); }
209 
211  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
212  void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
213  scale(Teuchos::ArrayView< const Scalar > alpha) { XPETRA_MONITOR("TpetraMultiVector::scale"); vec_->scale(alpha); }
214 
216  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
217  void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
218  scale(const Scalar &alpha, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A) { XPETRA_MONITOR("TpetraMultiVector::scale"); vec_->scale(alpha, toTpetra(A)); }
219 
221  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
222  void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
223  update(const Scalar &alpha, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const Scalar &beta) { XPETRA_MONITOR("TpetraMultiVector::update"); vec_->update(alpha, toTpetra(A), beta); }
224 
226  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
227  void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
228  update(const Scalar &alpha, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const Scalar &beta, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &B, const Scalar &gamma) { XPETRA_MONITOR("TpetraMultiVector::update"); vec_->update(alpha, toTpetra(A), beta, toTpetra(B), gamma); }
229 
231  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
232  void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
233  norm1(const Teuchos::ArrayView< typename Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const { XPETRA_MONITOR("TpetraMultiVector::norm1"); vec_->norm1(norms); }
234 
236  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
237  void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
238  norm2(const Teuchos::ArrayView< typename Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const { XPETRA_MONITOR("TpetraMultiVector::norm2"); vec_->norm2(norms); }
239 
241  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
242  void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
243  normInf(const Teuchos::ArrayView< typename Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const { XPETRA_MONITOR("TpetraMultiVector::normInf"); vec_->normInf(norms); }
244 
246  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
247  void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
248  meanValue(const Teuchos::ArrayView< Scalar > &means) const { XPETRA_MONITOR("TpetraMultiVector::meanValue"); vec_->meanValue(means); }
249 
251  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
252  void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
253  multiply(Teuchos::ETransp transA, Teuchos::ETransp transB, const Scalar &alpha, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &B, const Scalar &beta) { XPETRA_MONITOR("TpetraMultiVector::multiply"); vec_->multiply(transA, transB, alpha, toTpetra(A), toTpetra(B), beta); }
254 
255 
257  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
258  size_t TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
259  getNumVectors() const { XPETRA_MONITOR("TpetraMultiVector::getNumVectors"); return vec_->getNumVectors(); }
260 
262  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
263  size_t TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
264  getLocalLength() const { XPETRA_MONITOR("TpetraMultiVector::getLocalLength"); return vec_->getLocalLength(); }
265 
267  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
268  global_size_t TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
269  getGlobalLength() const { XPETRA_MONITOR("TpetraMultiVector::getGlobalLength"); return vec_->getGlobalLength(); }
270 
271  // \brief Checks to see if the local length, number of vectors and size of Scalar type match
272  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
273  bool TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
274  isSameSize(const MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> & vec) const { XPETRA_MONITOR("TpetraMultiVector::isSameSize"); return vec_->isSameSize(toTpetra(vec));}
275 
277  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
278  std::string TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
279  description() const { XPETRA_MONITOR("TpetraMultiVector::description"); return vec_->description(); }
280 
282  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
283  void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
284  describe(Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel) const { XPETRA_MONITOR("TpetraMultiVector::describe"); vec_->describe(out, verbLevel); }
285 
287  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
288  void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
289  randomize(bool bUseXpetraImplementation) {
290  XPETRA_MONITOR("TpetraMultiVector::randomize");
291 
292  if(bUseXpetraImplementation)
293  MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node >::Xpetra_randomize();
294  else
295  vec_->randomize();
296  }
297 
298  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
299  Teuchos::RCP< const Map<LocalOrdinal,GlobalOrdinal,Node> >
300  TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
301  getMap() const { XPETRA_MONITOR("TpetraMultiVector::getMap"); return toXpetra(vec_->getMap()); }
302 
303  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
304  void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
305  doImport(const DistObject< Scalar, LocalOrdinal,GlobalOrdinal,Node> &source, const Import<LocalOrdinal,GlobalOrdinal,Node> &importer, CombineMode CM) {
306  XPETRA_MONITOR("TpetraMultiVector::doImport");
307 
308  XPETRA_DYNAMIC_CAST(const TpetraMultiVectorClass, source, tSource, "Xpetra::TpetraMultiVector::doImport only accept Xpetra::TpetraMultiVector as input arguments."); //TODO: remove and use toTpetra()
309  RCP< const Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal,Node> > v = tSource.getTpetra_MultiVector();
310  this->getTpetra_MultiVector()->doImport(*v, toTpetra(importer), toTpetra(CM));
311  }
312 
313  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
314  void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
315  doExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Import<LocalOrdinal,GlobalOrdinal,Node>& importer, CombineMode CM) {
316  XPETRA_MONITOR("TpetraMultiVector::doExport");
317 
318  XPETRA_DYNAMIC_CAST(const TpetraMultiVectorClass, dest, tDest, "Xpetra::TpetraMultiVector::doImport only accept Xpetra::TpetraMultiVector as input arguments."); //TODO: remove and use toTpetra()
319  RCP< const Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal,Node> > v = tDest.getTpetra_MultiVector();
320  this->getTpetra_MultiVector()->doExport(*v, toTpetra(importer), toTpetra(CM));
321 
322  }
323 
324  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
325  void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
326  doImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Export<LocalOrdinal,GlobalOrdinal,Node>& exporter, CombineMode CM) {
327  XPETRA_MONITOR("TpetraMultiVector::doImport");
328 
329  XPETRA_DYNAMIC_CAST(const TpetraMultiVectorClass, source, tSource, "Xpetra::TpetraMultiVector::doImport only accept Xpetra::TpetraMultiVector as input arguments."); //TODO: remove and use toTpetra()
330  RCP< const Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal,Node> > v = tSource.getTpetra_MultiVector();
331  this->getTpetra_MultiVector()->doImport(*v, toTpetra(exporter), toTpetra(CM));
332 
333  }
334 
335  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
336  void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
337  doExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Export<LocalOrdinal,GlobalOrdinal,Node>& exporter, CombineMode CM) {
338  XPETRA_MONITOR("TpetraMultiVector::doExport");
339 
340  XPETRA_DYNAMIC_CAST(const TpetraMultiVectorClass, dest, tDest, "Xpetra::TpetraMultiVector::doImport only accept Xpetra::TpetraMultiVector as input arguments."); //TODO: remove and use toTpetra()
341  RCP< const Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal,Node> > v = tDest.getTpetra_MultiVector();
342  this->getTpetra_MultiVector()->doExport(*v, toTpetra(exporter), toTpetra(CM));
343 
344  }
345 
346  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
347  void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
348  replaceMap(const RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> >& map) {
349  XPETRA_MONITOR("TpetraMultiVector::replaceMap");
350  this->getTpetra_MultiVector()->replaceMap(toTpetra(map));
351  }
352 
353 #ifdef XPETRA_ENABLE_DEPRECATED_CODE
354  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
355  template<class Node2>
356  RCP<MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node2> > XPETRA_DEPRECATED
357  TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
358  clone(const RCP<Node2> &node2) const {
359  XPETRA_MONITOR("TpetraMultiVector::clone");
360  return RCP<MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node2> >(new TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node2>(vec_->clone(node2)));
361  }
362 #endif
363 
364 
366  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
367  TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
368  TpetraMultiVector(const Teuchos::RCP<Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node> > &vec) : vec_(vec) { } //TODO removed const
369 
371  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
372  RCP< Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node> >
373  TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
374  getTpetra_MultiVector() const { return vec_; }
375 
377  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
378  void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
379  setSeed(unsigned int seed) { XPETRA_MONITOR("TpetraMultiVector::seedrandom"); Teuchos::ScalarTraits< Scalar >::seedrandom(seed); }
380 
381 
382 #ifdef HAVE_XPETRA_KOKKOS_REFACTOR
383 #if 0
394  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
395  template<class TargetDeviceType>
396  typename Kokkos::Impl::if_c<
397  Kokkos::Impl::is_same<
398  typename dual_view_type::t_dev_um::execution_space::memory_space,
399  typename TargetDeviceType::memory_space>::value,
400  typename dual_view_type::t_dev_um,
401  typename dual_view_type::t_host_um>::type
402  TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
403  getLocalView () const {
404  return this->MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node >::template getLocalView<TargetDeviceType>();
405  }
406 #endif
407 
408  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
409  typename TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::dual_view_type::t_host_um
410  TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
411  getHostLocalView () const {
412  return subview(vec_->template getLocalView<typename TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::dual_view_type::host_mirror_space> (),
413  Kokkos::ALL(), Kokkos::ALL());
414  }
415 
416  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
417  typename TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::dual_view_type::t_dev_um
418  TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
419  getDeviceLocalView() const {
420  return subview(vec_->template getLocalView<typename TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::dual_view_type::t_dev_um::execution_space> (),
421  Kokkos::ALL(), Kokkos::ALL());
422  }
423 
424 #endif
425 
428  template<class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
429  void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::
430  assign (const MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>& rhs)
431  {
432  typedef TpetraMultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> this_type;
433  const this_type* rhsPtr = dynamic_cast<const this_type*> (&rhs);
434  TEUCHOS_TEST_FOR_EXCEPTION(
435  rhsPtr == NULL, std::invalid_argument, "Xpetra::MultiVector::operator=:"
436  " The left-hand side (LHS) of the assignment has a different type than "
437  "the right-hand side (RHS). The LHS has type Xpetra::TpetraMultiVector"
438  " (which means it wraps a Tpetra::MultiVector), but the RHS has some "
439  "other type. This probably means that the RHS wraps an "
440  "Epetra_MultiVector. Xpetra::MultiVector does not currently implement "
441  "assignment from an Epetra object to a Tpetra object, though this could"
442  " be added with sufficient interest.");
443 
444  typedef Tpetra::MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node> TMV;
445  RCP<const TMV> rhsImpl = rhsPtr->getTpetra_MultiVector ();
446  RCP<TMV> lhsImpl = this->getTpetra_MultiVector ();
447 
448  TEUCHOS_TEST_FOR_EXCEPTION(
449  rhsImpl.is_null (), std::logic_error, "Xpetra::MultiVector::operator= "
450  "(in Xpetra::TpetraMultiVector::assign): *this (the right-hand side of "
451  "the assignment) has a null RCP<Tpetra::MultiVector> inside. Please "
452  "report this bug to the Xpetra developers.");
453  TEUCHOS_TEST_FOR_EXCEPTION(
454  lhsImpl.is_null (), std::logic_error, "Xpetra::MultiVector::operator= "
455  "(in Xpetra::TpetraMultiVector::assign): The left-hand side of the "
456  "assignment has a null RCP<Tpetra::MultiVector> inside. Please report "
457  "this bug to the Xpetra developers.");
458 
459  Tpetra::deep_copy (*lhsImpl, *rhsImpl);
460  }
461 
462 
463 #ifdef HAVE_XPETRA_EPETRA
464 
465 #if ((defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_OPENMP) || !defined(HAVE_TPETRA_INST_INT_INT))) || \
466  (!defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_SERIAL) || !defined(HAVE_TPETRA_INST_INT_INT))))
467 
468  // specialization for TpetraMultiVector on EpetraNode and GO=int
469  template <class Scalar>
470  class TpetraMultiVector<Scalar,int,int,EpetraNode>
471  : public virtual MultiVector< Scalar, int, int, EpetraNode >
472  {
473  typedef int LocalOrdinal;
474  typedef int GlobalOrdinal;
475  typedef EpetraNode Node;
476 
477  // The following typedef are used by the XPETRA_DYNAMIC_CAST() macro.
478  typedef TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> TpetraMultiVectorClass;
479 
480  public:
481 
483 
484 
486  TpetraMultiVector () {
487  XPETRA_TPETRA_ETI_EXCEPTION( typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name() , typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name(), "int", typeid(EpetraNode).name() );
488  }
489 
491  TpetraMultiVector(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map, size_t NumVectors, bool zeroOut=true) {
492  XPETRA_TPETRA_ETI_EXCEPTION( typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name() , typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name(), "int", typeid(EpetraNode).name() );
493  }
494 
496  TpetraMultiVector(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source) {
497  XPETRA_TPETRA_ETI_EXCEPTION( typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name() , typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name(), "int", typeid(EpetraNode).name() );
498  }
499 
501  TpetraMultiVector(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map, const Teuchos::ArrayView< const Scalar > &A, size_t LDA, size_t NumVectors) {
502  XPETRA_TPETRA_ETI_EXCEPTION( typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name() , typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name(), "int", typeid(EpetraNode).name() );
503  }
504 
506  TpetraMultiVector(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map, const Teuchos::ArrayView< const Teuchos::ArrayView< const Scalar > > &ArrayOfPtrs, size_t NumVectors) {
507  XPETRA_TPETRA_ETI_EXCEPTION( typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name() , typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name(), "int", typeid(EpetraNode).name() );
508  }
509 
510 
512  virtual ~TpetraMultiVector() { }
513 
515 
517 
518 
520  void replaceGlobalValue(GlobalOrdinal globalRow, size_t vectorIndex, const Scalar &value) { }
521 
523  void sumIntoGlobalValue(GlobalOrdinal globalRow, size_t vectorIndex, const Scalar &value) { }
524 
526  void replaceLocalValue(LocalOrdinal myRow, size_t vectorIndex, const Scalar &value) { }
527 
529  void sumIntoLocalValue(LocalOrdinal myRow, size_t vectorIndex, const Scalar &value) { }
530 
532  void putScalar(const Scalar &value) { }
533 
535  void reduce() { }
536 
538 
540 
541 
543  Teuchos::RCP< const Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > getVector(size_t j) const { return Teuchos::null; }
544 
546  Teuchos::RCP< Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > getVectorNonConst(size_t j) { return Teuchos::null; }
547 
549  Teuchos::ArrayRCP< const Scalar > getData(size_t j) const { return Teuchos::ArrayRCP< const Scalar >(); }
550 
552  Teuchos::ArrayRCP< Scalar > getDataNonConst(size_t j) { return Teuchos::ArrayRCP< Scalar >(); }
553 
555  void get1dCopy(Teuchos::ArrayView< Scalar > A, size_t LDA) const { }
556 
558  void get2dCopy(Teuchos::ArrayView< const Teuchos::ArrayView< Scalar > > ArrayOfPtrs) const { }
559 
561  Teuchos::ArrayRCP< const Scalar > get1dView() const { return Teuchos::ArrayRCP< const Scalar >(); }
562 
564  Teuchos::ArrayRCP< Teuchos::ArrayRCP< const Scalar > > get2dView() const { return Teuchos::ArrayRCP< Teuchos::ArrayRCP< const Scalar > >(); }
565 
567  Teuchos::ArrayRCP< Scalar > get1dViewNonConst() { return Teuchos::ArrayRCP< Scalar >(); }
568 
570  Teuchos::ArrayRCP< Teuchos::ArrayRCP< Scalar > > get2dViewNonConst() { return Teuchos::ArrayRCP< Teuchos::ArrayRCP< Scalar > >(); }
571 
573 
575 
576 
578  void dot(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const Teuchos::ArrayView< Scalar > &dots) const { }
579 
581  void abs(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A) { }
582 
584  void reciprocal(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A) { }
585 
587  void scale(const Scalar &alpha) { }
588 
590  void scale(Teuchos::ArrayView< const Scalar > alpha) { }
591 
593  void scale(const Scalar &alpha, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A) { }
594 
596  void update(const Scalar &alpha, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const Scalar &beta) { }
597 
599  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) { }
600 
602  void norm1(const Teuchos::ArrayView< typename Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const { }
603 
605  void norm2(const Teuchos::ArrayView< typename Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const { }
606 
608  void normInf(const Teuchos::ArrayView< typename Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const { }
609 
611  void meanValue(const Teuchos::ArrayView< Scalar > &means) const { }
612 
614  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) { }
615 
617 
619 
620 
622  size_t getNumVectors() const { return 0; }
623 
625  size_t getLocalLength() const { return 0; }
626 
628  global_size_t getGlobalLength() const { return 0; }
629 
630  // \! Checks to see if the local length, number of vectors and size of Scalar type match
631  bool isSameSize(const MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> & vec) const { return false; }
632 
634 
636 
637 
639  std::string description() const { return std::string(""); }
640 
642  void describe(Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel=Teuchos::Describable::verbLevel_default) const { }
643 
645 
647  void elementWiseMultiply(Scalar scalarAB, const Vector<Scalar,LocalOrdinal,GlobalOrdinal,Node> &A, const MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> &B, Scalar scalarThis) {};
648 
650  void randomize(bool bUseXpetraImplementation = false) { }
651 
652  //{@
653  // Implements DistObject interface
654 
655  Teuchos::RCP< const Map<LocalOrdinal,GlobalOrdinal,Node> > getMap() const { return Teuchos::null; }
656 
657  void doImport(const DistObject< Scalar, LocalOrdinal,GlobalOrdinal,Node> &source, const Import<LocalOrdinal,GlobalOrdinal,Node> &importer, CombineMode CM) { }
658 
659  void doExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Import<LocalOrdinal,GlobalOrdinal,Node>& importer, CombineMode CM) { }
660 
661  void doImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Export<LocalOrdinal,GlobalOrdinal,Node>& exporter, CombineMode CM) { }
662 
663  void doExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Export<LocalOrdinal,GlobalOrdinal,Node>& exporter, CombineMode CM) { }
664 
665  void replaceMap(const RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> >& map) { }
666 
667 #ifdef XPETRA_ENABLE_DEPRECATED_CODE
668  template<class Node2>
669  RCP<MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node2> > clone(const RCP<Node2> &node2) const { return Teuchos::null; }
670 #endif
671 
672 
674 
675 
677  TpetraMultiVector(const Teuchos::RCP<Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node> > &vec) {
678  XPETRA_TPETRA_ETI_EXCEPTION( typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name() , typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name(), "int", typeid(EpetraNode).name() );
679  }
680 
682  RCP< Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node> > getTpetra_MultiVector() const { return Teuchos::null; }
683 
685  void setSeed(unsigned int seed) { }
686 
687 
688 #ifdef HAVE_XPETRA_KOKKOS_REFACTOR
689  typedef typename Xpetra::MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::dual_view_type dual_view_type;
690  //typedef typename Xpetra::MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::host_execution_space host_execution_space;
691  //typedef typename Xpetra::MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::dev_execution_space dev_execution_space;
692 
703  template<class TargetDeviceType>
704  typename Kokkos::Impl::if_c<
705  Kokkos::Impl::is_same<
706  typename dual_view_type::t_dev_um::execution_space::memory_space,
707  typename TargetDeviceType::memory_space>::value,
708  typename dual_view_type::t_dev_um,
709  typename dual_view_type::t_host_um>::type
710  getLocalView () const {
711  typename Kokkos::Impl::if_c<
712  Kokkos::Impl::is_same<
713  typename dual_view_type::t_dev_um::execution_space::memory_space,
714  typename TargetDeviceType::memory_space>::value,
715  typename dual_view_type::t_dev_um,
716  typename dual_view_type::t_host_um>::type dummy;
717  return dummy;
718  }
719 
720  typename dual_view_type::t_host_um getHostLocalView () const {
721  //return subview(vec_->template getLocalView<typename dual_view_type::host_mirror_space> (),
722  // Kokkos::ALL(), Kokkos::ALL());
723  return typename dual_view_type::t_host_um();
724  }
725 
726  typename dual_view_type::t_dev_um getDeviceLocalView() const {
727  //return subview(vec_->template getLocalView<typename dual_view_type::t_dev_um::execution_space> (),
728  // Kokkos::ALL(), Kokkos::ALL());
729  return typename dual_view_type::t_dev_um();
730  }
731 
732 #endif
733 
735 
736  protected:
739  virtual void
740  assign (const MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>& rhs)
741  { }
742  }; // TpetraMultiVector class (specialization GO=int, NO=EpetraNode)
743 #endif
744 
745 #if ((defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_OPENMP) || !defined(HAVE_TPETRA_INST_INT_LONG_LONG))) || \
746  (!defined(EPETRA_HAVE_OMP) && (!defined(HAVE_TPETRA_INST_SERIAL) || !defined(HAVE_TPETRA_INST_INT_LONG_LONG))))
747 
748  // specialization for TpetraMultiVector on EpetraNode and GO=long long
749  template <class Scalar>
750  class TpetraMultiVector<Scalar,int,long long,EpetraNode>
751  : public virtual MultiVector< Scalar, int, long long, EpetraNode >
752  {
753  typedef int LocalOrdinal;
754  typedef long long GlobalOrdinal;
755  typedef EpetraNode Node;
756 
757  // The following typedef are used by the XPETRA_DYNAMIC_CAST() macro.
758  typedef TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> TpetraMultiVectorClass;
759 
760  public:
761 
763 
764 
766  TpetraMultiVector(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map, size_t NumVectors, bool zeroOut=true) {
767  XPETRA_TPETRA_ETI_EXCEPTION( typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name() , typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name(), "long long", typeid(EpetraNode).name() );
768  }
769 
771  TpetraMultiVector(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source) {
772  XPETRA_TPETRA_ETI_EXCEPTION( typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name() , typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name(), "long long", typeid(EpetraNode).name() );
773  }
774 
776  TpetraMultiVector(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map, const Teuchos::ArrayView< const Scalar > &A, size_t LDA, size_t NumVectors) {
777  XPETRA_TPETRA_ETI_EXCEPTION( typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name() , typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name(), "long long", typeid(EpetraNode).name() );
778  }
779 
781  TpetraMultiVector(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map, const Teuchos::ArrayView< const Teuchos::ArrayView< const Scalar > > &ArrayOfPtrs, size_t NumVectors) {
782  XPETRA_TPETRA_ETI_EXCEPTION( typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name() , typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name(), "long long", typeid(EpetraNode).name() );
783  }
784 
785 
787  virtual ~TpetraMultiVector() { }
788 
790 
792 
793 
795  void replaceGlobalValue(GlobalOrdinal globalRow, size_t vectorIndex, const Scalar &value) { }
796 
798  void sumIntoGlobalValue(GlobalOrdinal globalRow, size_t vectorIndex, const Scalar &value) { }
799 
801  void replaceLocalValue(LocalOrdinal myRow, size_t vectorIndex, const Scalar &value) { }
802 
804  void sumIntoLocalValue(LocalOrdinal myRow, size_t vectorIndex, const Scalar &value) { }
805 
807  void putScalar(const Scalar &value) { }
808 
810  void reduce() { }
811 
813 
815 
816 
818  Teuchos::RCP< const Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > getVector(size_t j) const { return Teuchos::null; }
819 
821  Teuchos::RCP< Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > getVectorNonConst(size_t j) { return Teuchos::null; }
822 
824  Teuchos::ArrayRCP< const Scalar > getData(size_t j) const { return Teuchos::ArrayRCP< const Scalar >(); }
825 
827  Teuchos::ArrayRCP< Scalar > getDataNonConst(size_t j) { return Teuchos::ArrayRCP< Scalar >(); }
828 
830  void get1dCopy(Teuchos::ArrayView< Scalar > A, size_t LDA) const { }
831 
833  void get2dCopy(Teuchos::ArrayView< const Teuchos::ArrayView< Scalar > > ArrayOfPtrs) const { }
834 
836  Teuchos::ArrayRCP< const Scalar > get1dView() const { return Teuchos::ArrayRCP< const Scalar >(); }
837 
839  Teuchos::ArrayRCP< Teuchos::ArrayRCP< const Scalar > > get2dView() const { return Teuchos::ArrayRCP< Teuchos::ArrayRCP< const Scalar > >(); }
840 
842  Teuchos::ArrayRCP< Scalar > get1dViewNonConst() { return Teuchos::ArrayRCP< Scalar >(); }
843 
845  Teuchos::ArrayRCP< Teuchos::ArrayRCP< Scalar > > get2dViewNonConst() { return Teuchos::ArrayRCP< Teuchos::ArrayRCP< Scalar > >(); }
846 
848 
850 
851 
853  void dot(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const Teuchos::ArrayView< Scalar > &dots) const { }
854 
856  void abs(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A) { }
857 
859  void reciprocal(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A) { }
860 
862  void scale(const Scalar &alpha) { }
863 
865  void scale(Teuchos::ArrayView< const Scalar > alpha) { }
866 
868  void scale(const Scalar &alpha, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A) { }
869 
871  void update(const Scalar &alpha, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const Scalar &beta) { }
872 
874  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) { }
875 
877  void norm1(const Teuchos::ArrayView< typename Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const { }
878 
880  void norm2(const Teuchos::ArrayView< typename Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const { }
881 
883  void normInf(const Teuchos::ArrayView< typename Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const { }
884 
886  void meanValue(const Teuchos::ArrayView< Scalar > &means) const { }
887 
889  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) { }
890 
892 
894 
895 
897  size_t getNumVectors() const { return 0; }
898 
900  size_t getLocalLength() const { return 0; }
901 
903  global_size_t getGlobalLength() const { return 0; }
904 
905  // \! Checks to see if the local length, number of vectors and size of Scalar type match
906  bool isSameSize(const MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> & vec) const { return false; }
907 
909 
911 
912 
914  std::string description() const { return std::string(""); }
915 
917  void describe(Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel=Teuchos::Describable::verbLevel_default) const { }
918 
920 
922  void elementWiseMultiply(Scalar scalarAB, const Vector<Scalar,LocalOrdinal,GlobalOrdinal,Node> &A, const MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> &B, Scalar scalarThis) {};
923 
925  void randomize(bool bUseXpetraImplementation = false) { }
926 
927  //{@
928  // Implements DistObject interface
929 
930  Teuchos::RCP< const Map<LocalOrdinal,GlobalOrdinal,Node> > getMap() const { return Teuchos::null; }
931 
932  void doImport(const DistObject< Scalar, LocalOrdinal,GlobalOrdinal,Node> &source, const Import<LocalOrdinal,GlobalOrdinal,Node> &importer, CombineMode CM) { }
933 
934  void doExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Import<LocalOrdinal,GlobalOrdinal,Node>& importer, CombineMode CM) { }
935 
936  void doImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Export<LocalOrdinal,GlobalOrdinal,Node>& exporter, CombineMode CM) { }
937 
938  void doExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Export<LocalOrdinal,GlobalOrdinal,Node>& exporter, CombineMode CM) { }
939 
940  void replaceMap(const RCP<const Map<LocalOrdinal,GlobalOrdinal,Node> >& map) { }
941 
942 #ifdef XPETRA_ENABLE_DEPRECATED_CODE
943  template<class Node2>
944  RCP<MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node2> > clone(const RCP<Node2> &node2) const { return Teuchos::null; }
945 #endif
946 
947 
949 
950 
952  TpetraMultiVector(const Teuchos::RCP<Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node> > &vec) {
953  XPETRA_TPETRA_ETI_EXCEPTION( typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name() , typeid(TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,EpetraNode>).name(), "long long", typeid(EpetraNode).name() );
954  }
955 
957  RCP< Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node> > getTpetra_MultiVector() const { return Teuchos::null; }
958 
960  void setSeed(unsigned int seed) { }
961 
962 
963 #ifdef HAVE_XPETRA_KOKKOS_REFACTOR
964  typedef typename Xpetra::MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::dual_view_type dual_view_type;
965  //typedef typename Xpetra::MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::host_execution_space host_execution_space;
966  //typedef typename Xpetra::MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>::dev_execution_space dev_execution_space;
967 
978  template<class TargetDeviceType>
979  typename Kokkos::Impl::if_c<
980  Kokkos::Impl::is_same<
981  typename dual_view_type::t_dev_um::execution_space::memory_space,
982  typename TargetDeviceType::memory_space>::value,
983  typename dual_view_type::t_dev_um,
984  typename dual_view_type::t_host_um>::type
985  getLocalView () const {
986  typename Kokkos::Impl::if_c<
987  Kokkos::Impl::is_same<
988  typename dual_view_type::t_dev_um::execution_space::memory_space,
989  typename TargetDeviceType::memory_space>::value,
990  typename dual_view_type::t_dev_um,
991  typename dual_view_type::t_host_um>::type dummy;
992  return dummy;
993  }
994 
995  typename dual_view_type::t_host_um getHostLocalView () const {
996  //return subview(vec_->template getLocalView<typename dual_view_type::host_mirror_space> (),
997  // Kokkos::ALL(), Kokkos::ALL());
998  return typename dual_view_type::t_host_um();
999  }
1000 
1001  typename dual_view_type::t_dev_um getDeviceLocalView() const {
1002  //return subview(vec_->template getLocalView<typename dual_view_type::t_dev_um::execution_space> (),
1003  // Kokkos::ALL(), Kokkos::ALL());
1004  return typename dual_view_type::t_dev_um();
1005  }
1006 
1007 #endif
1008 
1010 
1011  protected:
1014  virtual void
1015  assign (const MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>& rhs)
1016  { }
1017  }; // TpetraMultiVector class (specialization GO=int, NO=EpetraNode)
1018 
1019 #endif // TpetraMultiVector class (specialization GO=long long, NO=EpetraNode)
1020 
1021 #endif // HAVE_XPETRA_EPETRA
1022 
1023 } // Xpetra namespace
1024 
1025 // Following header file inculsion is needed for the dynamic_cast to TpetraVector in elementWiseMultiply (because we cannot dynamic_cast if target is not a complete type)
1026 // It is included here to avoid circular dependency between Vector and MultiVector
1027 // TODO: there is certainly a more elegant solution...
1028 #include "Xpetra_TpetraVector.hpp"
1029 
1030 namespace Xpetra {
1031  template <class Scalar, class LocalOrdinal, class GlobalOrdinal, class Node>
1032  void TpetraMultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>::elementWiseMultiply(Scalar scalarAB, const Vector<Scalar,LocalOrdinal,GlobalOrdinal,Node> &A, const MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> &B, Scalar scalarThis) {
1033  XPETRA_MONITOR("TpetraMultiVector::elementWiseMultiply");
1034 
1035  // XPETRA_DYNAMIC_CAST won't take TpetraVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>
1036  // as an argument, hence the following typedef.
1037  typedef TpetraVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> tpv;
1038  XPETRA_DYNAMIC_CAST(const tpv, A, tA, "Xpetra::TpetraMultiVectorMatrix->multiply() only accept Xpetra::TpetraMultiVector as input arguments.");
1039  XPETRA_DYNAMIC_CAST(const TpetraMultiVector, B, tB, "Xpetra::TpetraMultiVectorMatrix->multiply() only accept Xpetra::TpetraMultiVector as input arguments.");
1040  vec_->elementWiseMultiply(scalarAB, *tA.getTpetra_Vector(), *tB.getTpetra_MultiVector(), scalarThis);
1041  }
1042 
1043 } // Xpetra namespace
1044 
1045 #define XPETRA_TPETRAMULTIVECTOR_SHORT
1046 #endif // XPETRA_TPETRAMULTIVECTOR_HPP
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::replaceMap
void replaceMap(const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map)
Definition: Xpetra_TpetraMultiVector_def.hpp:940
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::getLocalLength
size_t getLocalLength() const
Local number of rows on the calling process.
Definition: Xpetra_TpetraMultiVector_def.hpp:625
Xpetra::clone
RCP< Map< LocalOrdinal, GlobalOrdinal, Node2 > > clone(const Map< LocalOrdinal, GlobalOrdinal, Node1 > &map, const RCP< Node2 > &node2)
Definition: Xpetra_Cloner.hpp:72
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:596
Xpetra::TpetraMultiVector::norm1
void norm1(const Teuchos::ArrayView< typename Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const
Compute 1-norm of each vector in multi-vector.
Definition: Xpetra_TpetraMultiVector_def.hpp:233
Xpetra::TpetraMultiVector::get1dView
Teuchos::ArrayRCP< const Scalar > get1dView() const
Const persisting (1-D) view of this multivector&#39;s local values.
Definition: Xpetra_TpetraMultiVector_def.hpp:170
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:682
Xpetra::TpetraMultiVector::get2dView
Teuchos::ArrayRCP< Teuchos::ArrayRCP< const Scalar > > get2dView() const
Return const persisting pointers to values.
Definition: Xpetra_TpetraMultiVector_def.hpp:176
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:934
Xpetra::TpetraMultiVector::dot
void dot(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const Teuchos::ArrayView< Scalar > &dots) const
Compute dot product of each corresponding pair of vectors, dots[i] = this[i].dot(A[i]).
Definition: Xpetra_TpetraMultiVector_def.hpp:193
Xpetra::TpetraMultiVector< 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:587
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:491
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:520
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::setSeed
void setSeed(unsigned int seed)
Set seed for Random function.
Definition: Xpetra_TpetraMultiVector_def.hpp:685
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:917
Xpetra::TpetraMultiVector::get2dViewNonConst
Teuchos::ArrayRCP< Teuchos::ArrayRCP< Scalar > > get2dViewNonConst()
Return non-const persisting pointers to values.
Definition: Xpetra_TpetraMultiVector_def.hpp:188
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:810
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:578
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:608
Xpetra::TpetraMultiVector::description
std::string description() const
A simple one-line description of this object.
Definition: Xpetra_TpetraMultiVector_def.hpp:279
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:781
Xpetra_TpetraMap.hpp
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,j).
Definition: Xpetra_TpetraMultiVector_def.hpp:859
Xpetra::TpetraMultiVector::getDataNonConst
Teuchos::ArrayRCP< Scalar > getDataNonConst(size_t j)
View of the local values in a particular vector of this multivector.
Definition: Xpetra_TpetraMultiVector_def.hpp:152
Xpetra::TpetraMultiVector< 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:938
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:821
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::~TpetraMultiVector
virtual ~TpetraMultiVector()
Destructor (virtual for memory safety of derived classes).
Definition: Xpetra_TpetraMultiVector_def.hpp:512
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::setSeed
void setSeed(unsigned int seed)
Set seed for Random function.
Definition: Xpetra_TpetraMultiVector_def.hpp:960
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:1032
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:766
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:663
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::getNumVectors
size_t getNumVectors() const
Number of columns in the multivector.
Definition: Xpetra_TpetraMultiVector_def.hpp:897
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::norm2
void norm2(const Teuchos::ArrayView< typename Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const
Definition: Xpetra_TpetraMultiVector_def.hpp:605
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:922
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:546
TEUCHOS_TEST_FOR_EXCEPTION
#define TEUCHOS_TEST_FOR_EXCEPTION(throw_exception_test, Exception, msg)
Xpetra::MultiVector::Xpetra_randomize
virtual void Xpetra_randomize()
Set multi-vector values to random numbers. XPetra implementation.
Definition: Xpetra_MultiVector.hpp:256
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:914
Xpetra::TpetraMultiVector::normInf
void normInf(const Teuchos::ArrayView< typename Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const
Compute Inf-norm of each vector in multi-vector.
Definition: Xpetra_TpetraMultiVector_def.hpp:243
Xpetra::TpetraMultiVector::getNumVectors
size_t getNumVectors() const
Number of columns in the multivector.
Definition: Xpetra_TpetraMultiVector_def.hpp:259
Xpetra_TpetraImport.hpp
Xpetra::TpetraMultiVector::setSeed
void setSeed(unsigned int seed)
Set seed for Random function.
Definition: Xpetra_TpetraMultiVector_def.hpp:379
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:552
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:862
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::TpetraMultiVectorClass
TpetraMultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > TpetraMultiVectorClass
Definition: Xpetra_TpetraMultiVector_def.hpp:758
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:657
Xpetra::TpetraMultiVector::randomize
void randomize(bool bUseXpetraImplementation=false)
Set multi-vector values to random numbers.
Definition: Xpetra_TpetraMultiVector_def.hpp:289
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:593
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::replaceMap
void replaceMap(const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map)
Definition: Xpetra_TpetraMultiVector_def.hpp:665
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:564
Xpetra::TpetraMultiVector::reduce
void reduce()
Sum values of a locally replicated multivector across all processes.
Definition: Xpetra_TpetraMultiVector_def.hpp:128
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:795
Teuchos::EVerbosityLevel
EVerbosityLevel
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:824
XPETRA_TPETRA_ETI_EXCEPTION
#define XPETRA_TPETRA_ETI_EXCEPTION(cl, obj, go, node)
Definition: Xpetra_Exceptions.hpp:79
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:526
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&#39;s local values.
Definition: Xpetra_TpetraMultiVector_def.hpp:833
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:818
Xpetra::TpetraMultiVector::replaceGlobalValue
void replaceGlobalValue(GlobalOrdinal globalRow, size_t vectorIndex, const Scalar &value)
Replace value, using global (row) index.
Definition: Xpetra_TpetraMultiVector_def.hpp:103
Xpetra::TpetraMultiVector::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 (&quot;forward mode&quot;).
Definition: Xpetra_TpetraMultiVector_def.hpp:305
Xpetra::TpetraMultiVector::meanValue
void meanValue(const Teuchos::ArrayView< Scalar > &means) const
Compute mean (average) value of each vector in multi-vector. The outcome of this routine is undefined...
Definition: Xpetra_TpetraMultiVector_def.hpp:248
Xpetra::TpetraMultiVector< 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:952
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&#39;s local values.
Definition: Xpetra_TpetraMultiVector_def.hpp:555
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:647
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:880
Xpetra::TpetraMultiVector::multiply
void multiply(Teuchos::ETransp transA, Teuchos::ETransp transB, const Scalar &alpha, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &B, const Scalar &beta)
Matrix-matrix multiplication: this = beta*this + alpha*op(A)*op(B).
Definition: Xpetra_TpetraMultiVector_def.hpp:253
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::get1dViewNonConst
Teuchos::ArrayRCP< Scalar > get1dViewNonConst()
Nonconst persisting (1-D) view of this multivector&#39;s local values.
Definition: Xpetra_TpetraMultiVector_def.hpp:567
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:677
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:570
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:599
Xpetra::TpetraMultiVector::abs
void abs(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A)
Put element-wise absolute values of input Multi-vector in target: A = abs(this).
Definition: Xpetra_TpetraMultiVector_def.hpp:198
Xpetra::TpetraMultiVector::putScalar
void putScalar(const Scalar &value)
Set all values in the multivector with the given value.
Definition: Xpetra_TpetraMultiVector_def.hpp:123
Xpetra::TpetraMultiVector::getLocalLength
size_t getLocalLength() const
Local number of rows on the calling process.
Definition: Xpetra_TpetraMultiVector_def.hpp:264
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:614
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:807
Xpetra_TpetraConfigDefs.hpp
rcp
TEUCHOS_DEPRECATED RCP< T > rcp(T *p, Dealloc_T dealloc, bool owns_mem)
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:523
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:655
Xpetra::TpetraMultiVector::get2dCopy
void get2dCopy(Teuchos::ArrayView< const Teuchos::ArrayView< Scalar > > ArrayOfPtrs) const
Fill the given array with a copy of this multivector&#39;s local values.
Definition: Xpetra_TpetraMultiVector_def.hpp:164
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:889
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:877
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::description
std::string description() const
A simple one-line description of this object.
Definition: Xpetra_TpetraMultiVector_def.hpp:639
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:590
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:853
Teuchos::basic_FancyOStream
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::getNumVectors
size_t getNumVectors() const
Number of columns in the multivector.
Definition: Xpetra_TpetraMultiVector_def.hpp:622
Xpetra::TpetraMultiVector::get1dViewNonConst
Teuchos::ArrayRCP< Scalar > get1dViewNonConst()
Nonconst persisting (1-D) view of this multivector&#39;s local values.
Definition: Xpetra_TpetraMultiVector_def.hpp:182
Xpetra::TpetraMultiVector::replaceLocalValue
void replaceLocalValue(LocalOrdinal myRow, size_t vectorIndex, const Scalar &value)
Replace value, using local (row) index.
Definition: Xpetra_TpetraMultiVector_def.hpp:113
Xpetra::TpetraMultiVector< 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:740
Xpetra::TpetraMultiVector::update
void update(const Scalar &alpha, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const Scalar &beta)
Update multi-vector values with scaled values of A, this = beta*this + alpha*A.
Definition: Xpetra_TpetraMultiVector_def.hpp:223
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::reduce
void reduce()
Sum values of a locally replicated multivector across all processes.
Definition: Xpetra_TpetraMultiVector_def.hpp:535
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::isSameSize
bool isSameSize(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &vec) const
Definition: Xpetra_TpetraMultiVector_def.hpp:906
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:827
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&#39;s local values.
Definition: Xpetra_TpetraMultiVector_def.hpp:558
Xpetra::TpetraMultiVector::sumIntoLocalValue
void sumIntoLocalValue(LocalOrdinal myRow, size_t vectorIndex, const Scalar &value)
Add value to existing value, using local (row) index.
Definition: Xpetra_TpetraMultiVector_def.hpp:118
Xpetra::TpetraMultiVector< 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:661
Xpetra::TpetraMultiVector::sumIntoGlobalValue
void sumIntoGlobalValue(GlobalOrdinal globalRow, size_t vectorIndex, const Scalar &value)
Add value to existing value, using global (row) index.
Definition: Xpetra_TpetraMultiVector_def.hpp:108
Xpetra::TpetraMultiVector::scale
void scale(const Scalar &alpha)
Scale the current values of a multi-vector, this = alpha*this.
Definition: Xpetra_TpetraMultiVector_def.hpp:208
Xpetra::TpetraMultiVector::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,j).
Definition: Xpetra_TpetraMultiVector_def.hpp:203
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::~TpetraMultiVector
virtual ~TpetraMultiVector()
Destructor (virtual for memory safety of derived classes).
Definition: Xpetra_TpetraMultiVector_def.hpp:787
Xpetra::TpetraMultiVector::isSameSize
bool isSameSize(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &vec) const
Definition: Xpetra_TpetraMultiVector_def.hpp:274
Teuchos::ScalarTraits::seedrandom
static void seedrandom(unsigned int s)
XPETRA_DYNAMIC_CAST
#define XPETRA_DYNAMIC_CAST(type, obj, newObj, exceptionMsg)
Definition: Xpetra_Exceptions.hpp:58
Xpetra::global_size_t
size_t global_size_t
Global size_t object.
Definition: Xpetra_ConfigDefs.hpp:170
Xpetra::TpetraMultiVector::TpetraMultiVector
TpetraMultiVector(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map, size_t NumVectors, bool zeroOut=true)
Basic constuctor.
Definition: Xpetra_TpetraMultiVector_def.hpp:64
Xpetra::TpetraMultiVector< 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:874
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:659
Teuchos::Describable::verbLevel_default
static const EVerbosityLevel verbLevel_default
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:856
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:430
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:930
Xpetra::TpetraMultiVector::getVector
Teuchos::RCP< const Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > getVector(size_t j) const
Return a Vector which is a const view of column j.
Definition: Xpetra_TpetraMultiVector_def.hpp:134
Xpetra::TpetraMultiVector< 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:936
Xpetra::TpetraMultiVector::get1dCopy
void get1dCopy(Teuchos::ArrayView< Scalar > A, size_t LDA) const
Fill the given array with a copy of this multivector&#39;s local values.
Definition: Xpetra_TpetraMultiVector_def.hpp:158
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&#39;s local values.
Definition: Xpetra_TpetraMultiVector_def.hpp:830
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,j).
Definition: Xpetra_TpetraMultiVector_def.hpp:584
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:903
Xpetra::toTpetra
RCP< const Tpetra::CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > toTpetra(const RCP< const CrsGraph< LocalOrdinal, GlobalOrdinal, Node > > &graph)
Definition: Xpetra_TpetraCrsGraph.hpp:430
Xpetra::TpetraMultiVector::~TpetraMultiVector
virtual ~TpetraMultiVector()
Destructor (virtual for memory safety of derived classes).
Definition: Xpetra_TpetraMultiVector_def.hpp:98
Xpetra::TpetraMultiVector::norm2
void norm2(const Teuchos::ArrayView< typename Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const
Definition: Xpetra_TpetraMultiVector_def.hpp:238
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:506
Xpetra::toXpetra
RCP< const CrsGraph< int, GlobalOrdinal, Node > > toXpetra(const Epetra_CrsGraph &g)
Definition: Xpetra_EpetraCrsGraph.cpp:168
Xpetra::TpetraMultiVector::getVectorNonConst
Teuchos::RCP< Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > getVectorNonConst(size_t j)
Return a Vector which is a nonconst view of column j.
Definition: Xpetra_TpetraMultiVector_def.hpp:140
Xpetra::TpetraMultiVector< 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:543
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:532
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::isSameSize
bool isSameSize(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &vec) const
Definition: Xpetra_TpetraMultiVector_def.hpp:631
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:900
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:868
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:776
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::TpetraMultiVector
TpetraMultiVector(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source)
Copy constructor (performs a deep copy).
Definition: Xpetra_TpetraMultiVector_def.hpp:771
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:798
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:501
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:886
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::get1dView
Teuchos::ArrayRCP< const Scalar > get1dView() const
Const persisting (1-D) view of this multivector&#39;s local values.
Definition: Xpetra_TpetraMultiVector_def.hpp:836
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::get1dView
Teuchos::ArrayRCP< const Scalar > get1dView() const
Const persisting (1-D) view of this multivector&#39;s local values.
Definition: Xpetra_TpetraMultiVector_def.hpp:561
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:883
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:1015
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:865
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 (&quot;reverse mode&quot;).
Definition: Xpetra_TpetraMultiVector_def.hpp:315
Xpetra::CombineMode
CombineMode
Xpetra::Combine Mode enumerable type.
Definition: Xpetra_ConfigDefs.hpp:214
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:932
Xpetra::TpetraMultiVector::getTpetra_MultiVector
RCP< Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > getTpetra_MultiVector() const
Get the underlying Tpetra multivector.
Definition: Xpetra_TpetraMultiVector_def.hpp:374
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::TpetraMultiVectorClass
TpetraMultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > TpetraMultiVectorClass
Definition: Xpetra_TpetraMultiVector_def.hpp:478
XPETRA_MONITOR
#define XPETRA_MONITOR(funcName)
Definition: Xpetra_ConfigDefs.hpp:128
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:611
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:549
Xpetra::TpetraMultiVector< Scalar, int, long long, EpetraNode >::get1dViewNonConst
Teuchos::ArrayRCP< Scalar > get1dViewNonConst()
Nonconst persisting (1-D) view of this multivector&#39;s local values.
Definition: Xpetra_TpetraMultiVector_def.hpp:842
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::getGlobalLength
global_size_t getGlobalLength() const
Global number of rows in the multivector.
Definition: Xpetra_TpetraMultiVector_def.hpp:628
Xpetra_Utils.hpp
Xpetra_TpetraMultiVector_decl.hpp
Xpetra::TpetraMultiVector::describe
void describe(Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel=Teuchos::Describable::verbLevel_default) const
Print the object with the given verbosity level to a FancyOStream.
Definition: Xpetra_TpetraMultiVector_def.hpp:284
Xpetra_TpetraExport.hpp
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:581
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:925
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:845
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:804
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::TpetraMultiVector
TpetraMultiVector(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source)
Copy constructor (performs a deep copy).
Definition: Xpetra_TpetraMultiVector_def.hpp:496
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:602
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:839
Xpetra::TpetraMultiVector::getMap
Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > getMap() const
The Map describing the parallel distribution of this object.
Definition: Xpetra_TpetraMultiVector_def.hpp:301
Xpetra::TpetraMultiVector< Scalar, int, int, EpetraNode >::randomize
void randomize(bool bUseXpetraImplementation=false)
Set multi-vector values to random numbers.
Definition: Xpetra_TpetraMultiVector_def.hpp:650
Xpetra::TpetraMultiVector::getData
Teuchos::ArrayRCP< const Scalar > getData(size_t j) const
Const view of the local values in a particular vector of this multivector.
Definition: Xpetra_TpetraMultiVector_def.hpp:146
Xpetra::TpetraMultiVector< 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:871
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:801
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:529
Xpetra::TpetraMultiVector::getGlobalLength
global_size_t getGlobalLength() const
Global number of rows in the multivector.
Definition: Xpetra_TpetraMultiVector_def.hpp:269
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:642
Xpetra::TpetraMultiVector::replaceMap
void replaceMap(const RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map)
Definition: Xpetra_TpetraMultiVector_def.hpp:348
Teuchos::RCP::is_null
bool is_null() const
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:957
Generated by   doxygen 1.8.5