Xpetra Version of the Day
Loading...
Searching...
No Matches
Xpetra_EpetraIntVector.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_EPETRAINTVECTOR_HPP
47#define XPETRA_EPETRAINTVECTOR_HPP
48
49#include "Xpetra_EpetraConfigDefs.hpp"
50
51#include "Xpetra_ConfigDefs.hpp"
52#include "Xpetra_MultiVector.hpp"
53#include "Xpetra_Vector.hpp"
54#include "Xpetra_Exceptions.hpp"
55
56#include "Xpetra_EpetraMap.hpp"
57#include "Xpetra_EpetraMultiVector.hpp"
58#include "Epetra_IntVector.h"
59
60namespace Xpetra {
61
62// TODO: move that elsewhere
63template<class GlobalOrdinal, class Node>
64Epetra_IntVector & toEpetra(Vector<int, int, GlobalOrdinal, Node> &);
65
66template<class GlobalOrdinal, class Node>
67const Epetra_IntVector & toEpetra(const Vector<int, int, GlobalOrdinal, Node> &);
68//
69
70 // stub implementation for EpetraIntVectorT
71 template<class EpetraGlobalOrdinal, class Node>
72 class EpetraIntVectorT
73 : public Vector<int,int,EpetraGlobalOrdinal, Node>
74 {
75 typedef int Scalar;
76 typedef int LocalOrdinal;
77 typedef EpetraGlobalOrdinal GlobalOrdinal;
78
79 public:
80
82
83
85 explicit EpetraIntVectorT(const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &map, bool zeroOut=true) { }
86
88 ~EpetraIntVectorT() { };
89
91
93
94
96 int dot(const Vector<Scalar,LocalOrdinal,GlobalOrdinal,Node> &a) const { XPETRA_MONITOR("EpetraIntVectorT::dot"); TEUCHOS_TEST_FOR_EXCEPTION(-1, Xpetra::Exceptions::NotImplemented, "TODO"); TEUCHOS_UNREACHABLE_RETURN(-1); }
97
98
100 Teuchos::ScalarTraits<int>::magnitudeType norm1() const { XPETRA_MONITOR("EpetraIntVectorT::norm1"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); TEUCHOS_UNREACHABLE_RETURN(-1); }
101
102
104 Teuchos::ScalarTraits<int>::magnitudeType norm2() const { XPETRA_MONITOR("EpetraIntVectorT::norm2"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); TEUCHOS_UNREACHABLE_RETURN(-1); }
105
107 Teuchos::ScalarTraits<int>::magnitudeType normInf() const { XPETRA_MONITOR("EpetraIntVectorT::normInf"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); TEUCHOS_UNREACHABLE_RETURN(-1); }
108
110 int meanValue() const { XPETRA_MONITOR("EpetraIntVectorT::meanValue"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); TEUCHOS_UNREACHABLE_RETURN(-1); }
111
113 int maxValue() const { XPETRA_MONITOR("EpetraIntVectorT::maxValue"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); TEUCHOS_UNREACHABLE_RETURN(-1); }
114
115
117
119
120
122 void replaceGlobalValue(GlobalOrdinal globalRow, const Scalar &value) { XPETRA_MONITOR("EpetraIntVectorT::replaceGlobalValue"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); }
123
125 void sumIntoGlobalValue(GlobalOrdinal globalRow, const Scalar &value) { XPETRA_MONITOR("EpetraIntVectorT::sumIntoGlobalValue"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); }
126
128 void replaceLocalValue(LocalOrdinal myRow, const Scalar &value) { XPETRA_MONITOR("EpetraIntVectorT::replaceLocalValue"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); }
129
131 void sumIntoLocalValue(LocalOrdinal myRow, const Scalar &value) { XPETRA_MONITOR("EpetraIntVectorT::sumIntoLocalValue"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); }
132
134 void putScalar(const int &value) { }
135
137 void randomize(bool bUseXpetraImplementation = true) { XPETRA_MONITOR("EpetraIntVectorT::randomize"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "Xpetra::EpetraIntVectorT::randomize(): Functionnality not available in Epetra"); }
138
140 void randomize(const Scalar& /*minVal*/, const Scalar& /*maxVal*/, bool bUseXpetraImplementation = true) { XPETRA_MONITOR("EpetraIntVectorT::randomize"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "Xpetra::EpetraIntVectorT::randomize(): Functionnality not available in Epetra"); }
141
142
144
145 void setSeed(unsigned int seed) { XPETRA_MONITOR("EpetraIntVectorT::setSeed"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "Xpetra::EpetraIntVectorT::setSeed(): Functionnality not available in Epetra"); }
146
147
149
151
152
157
162
165 Teuchos::ArrayRCP<const int> getData(size_t j) const { return Teuchos::ArrayRCP<const int>(); }
166
169 Teuchos::ArrayRCP<int> getDataNonConst(size_t j) { return Teuchos::ArrayRCP<int>(); }
170
172
174
175
176 void dot(const MultiVector<int,int,GlobalOrdinal,Node> &A, const Teuchos::ArrayView<int> &dots) const { }
177
179 void abs(const MultiVector<int,int,GlobalOrdinal,Node> &A) { }
180
182 void reciprocal(const MultiVector<int,int,GlobalOrdinal,Node> &A) { }
183
185 void scale(const int &alpha) { }
186
188 void scale (Teuchos::ArrayView< const int > alpha) {
189 XPETRA_MONITOR("EpetraIntVectorT::scale");
190 TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO");
191 }
192
194 void update(const int &alpha, const MultiVector<int,int,GlobalOrdinal,Node> &A, const int &beta) {
195 XPETRA_MONITOR("EpetraIntVectorT::update");
196
197 // XPETRA_DYNAMIC_CAST(const EpetraMultiVectorT, A, eA, "This Xpetra::EpetraMultiVectorT method only accept Xpetra::EpetraMultiVectorT as input arguments.");
198 TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO");
199 }
200
202 void update(const int &alpha, const MultiVector<int,int,GlobalOrdinal,Node> &A, const int &beta, const MultiVector<int,int,GlobalOrdinal,Node> &B, const int &gamma) {
203 XPETRA_MONITOR("EpetraIntVectorT::update");
204
205 //XPETRA_DYNAMIC_CAST(const EpetraMultiVectorT, A, eA, "This Xpetra::EpetraMultiVectorT method only accept Xpetra::EpetraMultiVectorT as input arguments.");
206 //XPETRA_DYNAMIC_CAST(const EpetraMultiVectorT, B, eB, "This Xpetra::EpetraMultiVectorT method only accept Xpetra::EpetraMultiVectorT as input arguments.");
207 TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO");
208 }
209
211 void norm1(const Teuchos::ArrayView<Teuchos::ScalarTraits<int>::magnitudeType> &norms) const { XPETRA_MONITOR("EpetraIntVectorT::norm1"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); }
212
214 void norm2(const Teuchos::ArrayView<Teuchos::ScalarTraits<int>::magnitudeType> &norms) const { XPETRA_MONITOR("EpetraIntVectorT::norm2"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); }
215
217 void normInf(const Teuchos::ArrayView<Teuchos::ScalarTraits<int>::magnitudeType> &norms) const { XPETRA_MONITOR("EpetraIntVectorT::normInf"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); }
218
220 void meanValue(const Teuchos::ArrayView<int> &means) const { XPETRA_MONITOR("EpetraIntVectorT::meanValue"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); }
221
223 void maxValue(const Teuchos::ArrayView<int> &maxs) const { XPETRA_MONITOR("EpetraIntVectorT::maxValue"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); }
224
226 void multiply(Teuchos::ETransp transA, Teuchos::ETransp transB, const int &alpha, const MultiVector<int,int,GlobalOrdinal,Node> &A, const MultiVector<int,int,GlobalOrdinal,Node> &B, const int &beta) { XPETRA_MONITOR("EpetraIntVectorT::multiply"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "Not available in Epetra"); }
227
229 void elementWiseMultiply(int scalarAB, const Vector<int,int,GlobalOrdinal,Node> &A, const MultiVector<int,int,GlobalOrdinal,Node> &B, int scalarThis) {
230 XPETRA_MONITOR("EpetraIntVectorT::elementWiseMultiply");
231 TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "Xpetra_EpetraIntVector: elementWiseMultiply not implemented because Epetra_IntVector does not support this operation");
232 }
233
235
237
238
240 void replaceGlobalValue(GlobalOrdinal globalRow, size_t vectorIndex, const Scalar &value) { XPETRA_MONITOR("EpetraIntVectorT::replaceGlobalValue"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); }
241
243 void sumIntoGlobalValue(GlobalOrdinal globalRow, size_t vectorIndex, const Scalar &value) { XPETRA_MONITOR("EpetraIntVectorT::sumIntoGlobalValue"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); }
244
246 void replaceLocalValue(LocalOrdinal myRow, size_t vectorIndex, const Scalar &value) { XPETRA_MONITOR("EpetraIntVectorT::replaceLocalValue"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); }
247
249 void sumIntoLocalValue(LocalOrdinal myRow, size_t vectorIndex, const Scalar &value) { XPETRA_MONITOR("EpetraIntVectorT::sumIntoLocalValue"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); }
250
252
254
255
257 size_t getNumVectors() const { XPETRA_MONITOR("EpetraIntVectorT::getNumVectors"); return 1; }
258
259
261 size_t getLocalLength() const { return 0; }
262
264 global_size_t getGlobalLength() const { return 0; }
265
267 bool isSameSize(const MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> & vec) const { return false; }
268
270
272
273
275 std::string description() const {
276 return std::string("");
277 }
278
280 void describe(Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel=Teuchos::Describable::verbLevel_default) const { }
281
283
284 RCP< Epetra_IntVector > getEpetra_IntVector() const { return Teuchos::null; }
285
286 const RCP<const Comm<int> > getComm() const {
287 TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO getComm Epetra MultiVector not implemented");
288 }
289
290 // Implementing DistObject
291 Teuchos::RCP<const Map<int, GlobalOrdinal, Node> > getMap () const {
292 return Teuchos::null;
293 }
294
297
300
303
306
307 void replaceMap(const RCP<const Map<int, GlobalOrdinal, Node> >& map) {
308 // do nothing
309 }
310
311 protected:
314 virtual void
317
318
319 private:
321 //RCP< Epetra_IntVector > vec_;
322
323 }; // class EpetraIntVectorT
324
325 // specialization on GO=int and Node=Serial
326#ifndef XPETRA_EPETRA_NO_32BIT_GLOBAL_INDICES
327 template<>
328 class EpetraIntVectorT<int, EpetraNode>
329 : public virtual Vector<int,int,int,EpetraNode>
330 {
331 typedef int Scalar;
332 typedef int LocalOrdinal;
333 typedef int GlobalOrdinal;
334 typedef EpetraNode Node;
335
336 public:
337
339
340
342 explicit EpetraIntVectorT(const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &map, bool zeroOut=true)
343 {
344 vec_ = rcp(new Epetra_IntVector(toEpetra<GlobalOrdinal,Node>(map), zeroOut));
345 }
346
348 ~EpetraIntVectorT() { };
349
351
353
354
356 int dot(const Vector<Scalar,LocalOrdinal,GlobalOrdinal,Node> &/* a */) const { XPETRA_MONITOR("EpetraIntVectorT::dot"); TEUCHOS_TEST_FOR_EXCEPTION(-1, Xpetra::Exceptions::NotImplemented, "TODO"); TEUCHOS_UNREACHABLE_RETURN(-1); }
357
358
360 Teuchos::ScalarTraits<int>::magnitudeType norm1() const { XPETRA_MONITOR("EpetraIntVectorT::norm1"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); TEUCHOS_UNREACHABLE_RETURN(-1); }
361
362
364 Teuchos::ScalarTraits<int>::magnitudeType norm2() const { XPETRA_MONITOR("EpetraIntVectorT::norm2"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); TEUCHOS_UNREACHABLE_RETURN(-1); }
365
367 Teuchos::ScalarTraits<int>::magnitudeType normInf() const { XPETRA_MONITOR("EpetraIntVectorT::normInf"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); TEUCHOS_UNREACHABLE_RETURN(-1); }
368
370 int meanValue() const { XPETRA_MONITOR("EpetraIntVectorT::meanValue"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); TEUCHOS_UNREACHABLE_RETURN(-1); }
371
373 int maxValue() const { XPETRA_MONITOR("EpetraIntVectorT::maxValue"); return vec_->MaxValue(); }
374
375
377
379
380
382 void replaceGlobalValue(GlobalOrdinal /* globalRow */, const Scalar &/* value */) { XPETRA_MONITOR("EpetraIntVectorT::replaceGlobalValue"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); }
383
385 void sumIntoGlobalValue(GlobalOrdinal /* globalRow */, const Scalar &/* value */) { XPETRA_MONITOR("EpetraIntVectorT::sumIntoGlobalValue"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); }
386
388 void replaceLocalValue(LocalOrdinal myRow, const Scalar &value) { XPETRA_MONITOR("EpetraIntVectorT::replaceLocalValue");(*vec_)[myRow] = value; }
389
391 void sumIntoLocalValue(LocalOrdinal myRow, const Scalar &value) { XPETRA_MONITOR("EpetraIntVectorT::sumIntoLocalValue"); (*vec_)[myRow] += value;}
392
394 void putScalar(const int &value) { vec_->PutValue(value); }
395
397 void randomize(bool /* bUseXpetraImplementation */ = true) { XPETRA_MONITOR("EpetraIntVectorT::randomize"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "Xpetra::EpetraIntVectorT::randomize(): Functionnality not available in Epetra"); }
398
400 void randomize(const Scalar& /*minVal*/, const Scalar& /*maxVal*/, bool /* bUseXpetraImplementation */ = true) { XPETRA_MONITOR("EpetraIntVectorT::randomize"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "Xpetra::EpetraIntVectorT::randomize(): Functionnality not available in Epetra"); }
401
403
404 void setSeed(unsigned int /* seed */) { XPETRA_MONITOR("EpetraIntVectorT::setSeed"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "Xpetra::EpetraIntVectorT::setSeed(): Functionnality not available in Epetra"); }
405
406
408
410
411
413 Teuchos::RCP< const Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > getVector(size_t /* j */) const {
414 TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO");
415 }
416
418 Teuchos::RCP< Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > getVectorNonConst(size_t /* j */) {
419 TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO");
420 }
421
424 Teuchos::ArrayRCP<const int> getData(size_t /* j */) const {
425 XPETRA_MONITOR("EpetraIntVectorT::getData");
426
427 int * data = vec_->Values();
428 int localLength = vec_->MyLength();
429
430 return ArrayRCP<int>(data, 0, localLength, false); // not ownership
431 }
432
435 Teuchos::ArrayRCP<int> getDataNonConst(size_t /* j */) {
436 XPETRA_MONITOR("EpetraIntVectorT::getDataNonConst");
437
438 int * data = vec_->Values();
439 int localLength = vec_->MyLength();
440
441 return ArrayRCP<int>(data, 0, localLength, false); // not ownership
442 }
443
445
447
448
449 void dot(const MultiVector<int,int,GlobalOrdinal,Node> &/* A */, const Teuchos::ArrayView<int> &/* dots */) const {
450 XPETRA_MONITOR("EpetraIntVectorT::dot");
451
452 //XPETRA_DYNAMIC_CAST(const EpetraMultiVectorT, A, eA, "This Xpetra::EpetraMultiVectorT method only accept Xpetra::EpetraMultiVectorT as input arguments.");
453 TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO");
454 }
455
457 void abs(const MultiVector<int,int,GlobalOrdinal,Node> &/* A */) {
458 XPETRA_MONITOR("EpetraIntVectorT::abs");
459
460 //XPETRA_DYNAMIC_CAST(const EpetraMultiVectorT, A, eA, "This Xpetra::EpetraMultiVectorT method only accept Xpetra::EpetraMultiVectorT as input arguments.");
461 TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO");
462 }
463
465 void reciprocal(const MultiVector<int,int,GlobalOrdinal,Node> &/* A */) {
466 XPETRA_MONITOR("EpetraIntVectorT::reciprocal");
467
468 //XPETRA_DYNAMIC_CAST(const EpetraMultiVectorT, A, eA, "This Xpetra::EpetraMultiVectorT method only accept Xpetra::EpetraMultiVectorT as input arguments.");
469 TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO");
470 }
471
473 void scale(const int &/* alpha */) {
474 XPETRA_MONITOR("EpetraIntVectorT::scale");
475 TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO");
476 }
477
479 void scale (Teuchos::ArrayView< const int > /* alpha */) {
480 XPETRA_MONITOR("EpetraIntVectorT::scale");
481 TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO");
482 }
483
485 void update(const int &/* alpha */, const MultiVector<int,int,GlobalOrdinal,Node> &/* A */, const int &/* beta */) {
486 XPETRA_MONITOR("EpetraIntVectorT::update");
487
488 // XPETRA_DYNAMIC_CAST(const EpetraMultiVectorT, A, eA, "This Xpetra::EpetraMultiVectorT method only accept Xpetra::EpetraMultiVectorT as input arguments.");
489 TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO");
490 }
491
493 void update(const int &/* alpha */, const MultiVector<int,int,GlobalOrdinal,Node> &/* A */, const int &/* beta */, const MultiVector<int,int,GlobalOrdinal,Node> &/* B */, const int &/* gamma */) {
494 XPETRA_MONITOR("EpetraIntVectorT::update");
495
496 //XPETRA_DYNAMIC_CAST(const EpetraMultiVectorT, A, eA, "This Xpetra::EpetraMultiVectorT method only accept Xpetra::EpetraMultiVectorT as input arguments.");
497 //XPETRA_DYNAMIC_CAST(const EpetraMultiVectorT, B, eB, "This Xpetra::EpetraMultiVectorT method only accept Xpetra::EpetraMultiVectorT as input arguments.");
498 TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO");
499 }
500
502 void norm1(const Teuchos::ArrayView<Teuchos::ScalarTraits<int>::magnitudeType> &/* norms */) const { XPETRA_MONITOR("EpetraIntVectorT::norm1"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); }
503
505 void norm2(const Teuchos::ArrayView<Teuchos::ScalarTraits<int>::magnitudeType> &/* norms */) const { XPETRA_MONITOR("EpetraIntVectorT::norm2"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); }
506
508 void normInf(const Teuchos::ArrayView<Teuchos::ScalarTraits<int>::magnitudeType> &/* norms */) const { XPETRA_MONITOR("EpetraIntVectorT::normInf"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); }
509
511 void meanValue(const Teuchos::ArrayView<int> &/* means */) const { XPETRA_MONITOR("EpetraIntVectorT::meanValue"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); }
512
514 void maxValue(const Teuchos::ArrayView<int> &/* maxs */) const { XPETRA_MONITOR("EpetraIntVectorT::maxValue"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); }
515
517 void multiply(Teuchos::ETransp /* transA */, Teuchos::ETransp /* transB */, const int &/* alpha */, const MultiVector<int,int,GlobalOrdinal,Node> &/* A */, const MultiVector<int,int,GlobalOrdinal,Node> &/* B */, const int &/* beta */) { XPETRA_MONITOR("EpetraIntVectorT::multiply"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "Not available in Epetra"); }
518
520 void elementWiseMultiply(int /* scalarAB */, const Vector<int,int,GlobalOrdinal,Node> &/* A */, const MultiVector<int,int,GlobalOrdinal,Node> &/* B */, int /* scalarThis */) {
521 XPETRA_MONITOR("EpetraIntVectorT::elementWiseMultiply");
522 TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "Xpetra_EpetraIntVector: elementWiseMultiply not implemented because Epetra_IntVector does not support this operation");
523 }
524
526
528
529
531 void replaceGlobalValue(GlobalOrdinal /* globalRow */, size_t /* vectorIndex */, const Scalar &/* value */) { XPETRA_MONITOR("EpetraIntVectorT::replaceGlobalValue"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); }
532
534 void sumIntoGlobalValue(GlobalOrdinal /* globalRow */, size_t /* vectorIndex */, const Scalar &/* value */) { XPETRA_MONITOR("EpetraIntVectorT::sumIntoGlobalValue"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); }
535
537 void replaceLocalValue(LocalOrdinal /* myRow */, size_t /* vectorIndex */, const Scalar &/* value */) { XPETRA_MONITOR("EpetraIntVectorT::replaceLocalValue"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); }
538
540 void sumIntoLocalValue(LocalOrdinal /* myRow */, size_t /* vectorIndex */, const Scalar &/* value */) { XPETRA_MONITOR("EpetraIntVectorT::sumIntoLocalValue"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); }
541
543
545
546
548 size_t getNumVectors() const { XPETRA_MONITOR("EpetraIntVectorT::getNumVectors"); return 1; }
549
550
552 size_t getLocalLength() const { return vec_->MyLength(); }
553
555 global_size_t getGlobalLength() const { return vec_->GlobalLength64(); }
556
558 bool isSameSize(const MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> & vec) const {
559 XPETRA_MONITOR("EpetraIntVectorT::isSameSize");
560 const Vector<Scalar,LocalOrdinal,GlobalOrdinal,Node> *asvec = dynamic_cast<const Vector<Scalar,LocalOrdinal,GlobalOrdinal,Node> *>(&vec);
561 if(!asvec) return false;
562 auto vv = toEpetra(*asvec);
563 return ( (vec_->MyLength() == vv.MyLength()) && (getNumVectors() == vec.getNumVectors()));
564 }
565
567
569
570
572 std::string description() const {
573 XPETRA_MONITOR("EpetraIntVectorT::description");
574
575 // This implementation come from Epetra_Vector_def.hpp (without modification)
576 std::ostringstream oss;
577 oss << Teuchos::Describable::description();
578 oss << "{length="<<this->getGlobalLength()
579 << "}";
580 return oss.str();
581 }
582
584 void describe(Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel=Teuchos::Describable::verbLevel_default) const {
585 XPETRA_MONITOR("EpetraIntVectorT::describe");
586
587 // This implementation come from Tpetra_Vector_def.hpp (without modification) // JG: true?
588 using std::endl;
589 using std::setw;
590 using Teuchos::VERB_DEFAULT;
591 using Teuchos::VERB_NONE;
592 using Teuchos::VERB_LOW;
593 using Teuchos::VERB_MEDIUM;
594 using Teuchos::VERB_HIGH;
595 using Teuchos::VERB_EXTREME;
596
597 if (verbLevel > Teuchos::VERB_NONE)
598 vec_->Print(out);
599 }
600
602
603 RCP< Epetra_IntVector > getEpetra_IntVector() const { return vec_; }
604
605 const RCP<const Comm<int> > getComm() const {
606 TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO getComm Epetra MultiVector not implemented");
607 }
608
609 // Implementing DistObject
610 Teuchos::RCP<const Map<int, GlobalOrdinal, Node> > getMap () const {
611 RCP<const Epetra_BlockMap> map = rcp(new Epetra_BlockMap(vec_->Map()));
612 return rcp (new Xpetra::EpetraMapT<GlobalOrdinal, Node>(map));
613 }
614
615 void doImport(const DistObject<int, int, GlobalOrdinal, Node> &source,
616 const Import<int, GlobalOrdinal, Node> &importer, CombineMode CM) {
617 XPETRA_MONITOR("EpetraIntVectorT::doImport");
618
619 XPETRA_DYNAMIC_CAST(const EpetraIntVectorT<GlobalOrdinal XPETRA_COMMA Node>, source, tSource, "Xpetra::EpetraIntVectorT::doImport only accept Xpetra::EpetraIntVectorT as input arguments.");
620 XPETRA_DYNAMIC_CAST(const EpetraImportT<GlobalOrdinal XPETRA_COMMA Node>, importer, tImporter, "Xpetra::EpetraIntVectorT::doImport only accept Xpetra::EpetraImportT as input arguments.");
621
622 const Epetra_IntVector & v = *tSource.getEpetra_IntVector();
623 int err = vec_->Import(v, *tImporter.getEpetra_Import(), toEpetra(CM));
624 TEUCHOS_TEST_FOR_EXCEPTION(err != 0, std::runtime_error, "Catch error code returned by Epetra.");
625 }
626
627 void doExport(const DistObject<int, LocalOrdinal, GlobalOrdinal, Node> &dest,
628 const Import<int, GlobalOrdinal, Node>& importer, CombineMode CM) {
629 XPETRA_MONITOR("EpetraIntVectorT::doExport");
630
631 XPETRA_DYNAMIC_CAST(const EpetraIntVectorT<GlobalOrdinal XPETRA_COMMA Node>, dest, tDest, "Xpetra::EpetraIntVectorT::doImport only accept Xpetra::EpetraIntVectorT as input arguments.");
632 XPETRA_DYNAMIC_CAST(const EpetraImportT<GlobalOrdinal XPETRA_COMMA Node>, importer, tImporter, "Xpetra::EpetraIntVectorT::doImport only accept Xpetra::EpetraImportT as input arguments.");
633
634 const Epetra_IntVector & v = *tDest.getEpetra_IntVector();
635 int err = vec_->Import(v, *tImporter.getEpetra_Import(), toEpetra(CM));
636 TEUCHOS_TEST_FOR_EXCEPTION(err != 0, std::runtime_error, "Catch error code returned by Epetra.");
637 }
638
639 void doImport(const DistObject<int, LocalOrdinal, GlobalOrdinal, Node> &source,
640 const Export<int, GlobalOrdinal, Node>& exporter, CombineMode CM) {
641 XPETRA_MONITOR("EpetraIntVectorT::doImport");
642
643 XPETRA_DYNAMIC_CAST(const EpetraIntVectorT<GlobalOrdinal XPETRA_COMMA Node>, source, tSource, "Xpetra::EpetraIntVectorT::doImport only accept Xpetra::EpetraIntVectorT as input arguments.");
644 XPETRA_DYNAMIC_CAST(const EpetraExportT<GlobalOrdinal XPETRA_COMMA Node>, exporter, tExporter, "Xpetra::EpetraIntVectorT::doImport only accept Xpetra::EpetraImportT as input arguments.");
645
646 const Epetra_IntVector & v = *tSource.getEpetra_IntVector();
647 int err = vec_->Import(v, *tExporter.getEpetra_Export(), toEpetra(CM));
648 TEUCHOS_TEST_FOR_EXCEPTION(err != 0, std::runtime_error, "Catch error code returned by Epetra.");
649 }
650
651 void doExport(const DistObject<int, LocalOrdinal, GlobalOrdinal, Node> &dest,
652 const Export<int, GlobalOrdinal, Node>& exporter, CombineMode CM) {
653 XPETRA_MONITOR("EpetraIntVectorT::doExport");
654
655 XPETRA_DYNAMIC_CAST(const EpetraIntVectorT<GlobalOrdinal XPETRA_COMMA Node>, dest, tDest, "Xpetra::EpetraIntVectorT::doImport only accept Xpetra::EpetraIntVectorT as input arguments.");
656 XPETRA_DYNAMIC_CAST(const EpetraExportT<GlobalOrdinal XPETRA_COMMA Node>, exporter, tExporter, "Xpetra::EpetraIntVectorT::doImport only accept Xpetra::EpetraImportT as input arguments.");
657
658 const Epetra_IntVector & v = *tDest.getEpetra_IntVector();
659 int err = vec_->Export(v, *tExporter.getEpetra_Export(), toEpetra(CM));
660 TEUCHOS_TEST_FOR_EXCEPTION(err != 0, std::runtime_error, "Catch error code returned by Epetra.");
661 }
662
663 void replaceMap(const RCP<const Map<int, GlobalOrdinal, Node> >& /* map */) {
664 // do nothing
665 }
666
667 protected:
670 virtual void
671 assign (const MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>& rhs)
672 {
673 typedef EpetraIntVectorT<GlobalOrdinal, Node> this_type;
674 const this_type* rhsPtr = dynamic_cast<const this_type*> (&rhs);
675 TEUCHOS_TEST_FOR_EXCEPTION(
676 rhsPtr == NULL, std::invalid_argument, "Xpetra::MultiVector::operator=: "
677 "The left-hand side (LHS) of the assignment has a different type than "
678 "the right-hand side (RHS). The LHS has type Xpetra::EpetraIntVectorT "
679 "(which means it wraps an Epetra_IntVector), but the RHS has some "
680 "other type. This probably means that the RHS wraps either an "
681 "Tpetra::MultiVector, or an Epetra_MultiVector. Xpetra::MultiVector "
682 "does not currently implement assignment from a Tpetra object to an "
683 "Epetra object, though this could be added with sufficient interest.");
684
685 RCP<const Epetra_IntVector> rhsImpl = rhsPtr->getEpetra_IntVector ();
686 RCP<Epetra_IntVector> lhsImpl = this->getEpetra_IntVector ();
687
688 TEUCHOS_TEST_FOR_EXCEPTION(
689 rhsImpl.is_null (), std::logic_error, "Xpetra::MultiVector::operator= "
690 "(in Xpetra::EpetraIntVectorT::assign): *this (the right-hand side of "
691 "the assignment) has a null RCP<Epetra_IntVector> inside. Please "
692 "report this bug to the Xpetra developers.");
693 TEUCHOS_TEST_FOR_EXCEPTION(
694 lhsImpl.is_null (), std::logic_error, "Xpetra::MultiVector::operator= "
695 "(in Xpetra::EpetraIntVectorT::assign): The left-hand side of the "
696 "assignment has a null RCP<Epetra_IntVector> inside. Please report "
697 "this bug to the Xpetra developers.");
698
699 // Epetra_IntVector's assignment operator does a deep copy.
700 *lhsImpl = *rhsImpl;
701 }
702
703
704 private:
706 RCP< Epetra_IntVector > vec_;
707 };
708#endif
709
710 // specialization on GO=long long and Node=Serial
711#ifndef XPETRA_EPETRA_NO_64BIT_GLOBAL_INDICES
712 template<>
713 class EpetraIntVectorT<long long, EpetraNode>
714 : public virtual Vector<int,int,long long,EpetraNode>
715 {
716 typedef int Scalar;
717 typedef int LocalOrdinal;
718 typedef long long GlobalOrdinal;
719 typedef EpetraNode Node;
720
721 public:
722
724
725
727 explicit EpetraIntVectorT(const Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> > &map, bool zeroOut=true)
728 {
729 vec_ = rcp(new Epetra_IntVector(toEpetra<GlobalOrdinal,Node>(map), zeroOut));
730 }
731
733 ~EpetraIntVectorT() { };
734
736
738
739
741 int dot(const Vector<Scalar,LocalOrdinal,GlobalOrdinal,Node> &/* a */) const { XPETRA_MONITOR("EpetraIntVectorT::dot"); TEUCHOS_TEST_FOR_EXCEPTION(true, Xpetra::Exceptions::NotImplemented, "TODO"); /* return -1; */ }
742
743
745 Teuchos::ScalarTraits<int>::magnitudeType norm1() const { XPETRA_MONITOR("EpetraIntVectorT::norm1"); TEUCHOS_TEST_FOR_EXCEPTION(true, Xpetra::Exceptions::NotImplemented, "TODO"); /* return -1; */ }
746
747
749 Teuchos::ScalarTraits<int>::magnitudeType norm2() const { XPETRA_MONITOR("EpetraIntVectorT::norm2"); TEUCHOS_TEST_FOR_EXCEPTION(true, Xpetra::Exceptions::NotImplemented, "TODO"); /* return -1; */ }
750
752 Teuchos::ScalarTraits<int>::magnitudeType normInf() const { XPETRA_MONITOR("EpetraIntVectorT::normInf"); TEUCHOS_TEST_FOR_EXCEPTION(true, Xpetra::Exceptions::NotImplemented, "TODO"); /* return -1; */ }
753
755 int meanValue() const { XPETRA_MONITOR("EpetraIntVectorT::meanValue"); TEUCHOS_TEST_FOR_EXCEPTION(true, Xpetra::Exceptions::NotImplemented, "TODO"); /* return -1; */ }
756
758 int maxValue() const { XPETRA_MONITOR("EpetraIntVectorT::maxValue"); return Teuchos::as<int>(vec_->MaxValue()); }
759
761
763
764
766 void replaceGlobalValue(GlobalOrdinal /* globalRow */, const Scalar &/* value */) { XPETRA_MONITOR("EpetraIntVectorT::replaceGlobalValue"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); }
767
769 void sumIntoGlobalValue(GlobalOrdinal /* globalRow */, const Scalar &/* value */) { XPETRA_MONITOR("EpetraIntVectorT::sumIntoGlobalValue"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); }
770
772 void replaceLocalValue(LocalOrdinal myRow, const Scalar &value) { XPETRA_MONITOR("EpetraIntVectorT::replaceLocalValue");(*vec_)[myRow] = value;}
773
775 void sumIntoLocalValue(LocalOrdinal myRow, const Scalar &value) { XPETRA_MONITOR("EpetraIntVectorT::sumIntoLocalValue"); (*vec_)[myRow] += value;}
776
778 void putScalar(const int &value) { vec_->PutValue(value); }
779
781 void randomize(bool /* bUseXpetraImplementation */ = true) { XPETRA_MONITOR("EpetraIntVectorT::randomize"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "Xpetra::EpetraIntVectorT::randomize(): Functionnality not available in Epetra"); }
782
784 void randomize(const Scalar& /*minVal*/, const Scalar& /*maxVal*/, bool /* bUseXpetraImplementation */ = true) { XPETRA_MONITOR("EpetraIntVectorT::randomize"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "Xpetra::EpetraIntVectorT::randomize(): Functionnality not available in Epetra"); }
785
787
788 void setSeed(unsigned int /* seed */) { XPETRA_MONITOR("EpetraIntVectorT::setSeed"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "Xpetra::EpetraIntVectorT::setSeed(): Functionnality not available in Epetra"); }
789
790
792
794
795
800
805
808 Teuchos::ArrayRCP<const int> getData(size_t /* j */) const {
809 XPETRA_MONITOR("EpetraIntVectorT::getData");
810
811 int * data = vec_->Values();
812 int localLength = vec_->MyLength();
813
814 return ArrayRCP<int>(data, 0, localLength, false); // not ownership
815 }
816
819 Teuchos::ArrayRCP<int> getDataNonConst(size_t /* j */) {
820 XPETRA_MONITOR("EpetraIntVectorT::getDataNonConst");
821
822 int * data = vec_->Values();
823 int localLength = vec_->MyLength();
824
825 return ArrayRCP<int>(data, 0, localLength, false); // not ownership
826 }
827
829
831
832
833 void dot(const MultiVector<int,int,GlobalOrdinal,Node> &/* A */, const Teuchos::ArrayView<int> &/* dots */) const {
834 XPETRA_MONITOR("EpetraIntVectorT::dot");
835
836 //XPETRA_DYNAMIC_CAST(const EpetraMultiVectorT, A, eA, "This Xpetra::EpetraMultiVectorT method only accept Xpetra::EpetraMultiVectorT as input arguments.");
837 TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO");
838 }
839
841 void abs(const MultiVector<int,int,GlobalOrdinal,Node> &/* A */) {
842 XPETRA_MONITOR("EpetraIntVectorT::abs");
843
844 //XPETRA_DYNAMIC_CAST(const EpetraMultiVectorT, A, eA, "This Xpetra::EpetraMultiVectorT method only accept Xpetra::EpetraMultiVectorT as input arguments.");
845 TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO");
846 }
847
849 void reciprocal(const MultiVector<int,int,GlobalOrdinal,Node> &/* A */) {
850 XPETRA_MONITOR("EpetraIntVectorT::reciprocal");
851
852 //XPETRA_DYNAMIC_CAST(const EpetraMultiVectorT, A, eA, "This Xpetra::EpetraMultiVectorT method only accept Xpetra::EpetraMultiVectorT as input arguments.");
853 TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO");
854 }
855
857 void scale(const int &/* alpha */) {
858 XPETRA_MONITOR("EpetraIntVectorT::scale");
859 TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO");
860 }
861
863 void scale (Teuchos::ArrayView< const int > /* alpha */) {
864 XPETRA_MONITOR("EpetraIntVectorT::scale");
865 TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO");
866 }
867
869 void update(const int &/* alpha */, const MultiVector<int,int,GlobalOrdinal,Node> &/* A */, const int &/* beta */) {
870 XPETRA_MONITOR("EpetraIntVectorT::update");
871
872 // XPETRA_DYNAMIC_CAST(const EpetraMultiVectorT, A, eA, "This Xpetra::EpetraMultiVectorT method only accept Xpetra::EpetraMultiVectorT as input arguments.");
873 TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO");
874 }
875
877 void update(const int &/* alpha */, const MultiVector<int,int,GlobalOrdinal,Node> &/* A */, const int &/* beta */, const MultiVector<int,int,GlobalOrdinal,Node> &/* B */, const int &/* gamma */) {
878 XPETRA_MONITOR("EpetraIntVectorT::update");
879
880 //XPETRA_DYNAMIC_CAST(const EpetraMultiVectorT, A, eA, "This Xpetra::EpetraMultiVectorT method only accept Xpetra::EpetraMultiVectorT as input arguments.");
881 //XPETRA_DYNAMIC_CAST(const EpetraMultiVectorT, B, eB, "This Xpetra::EpetraMultiVectorT method only accept Xpetra::EpetraMultiVectorT as input arguments.");
882 TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO");
883 }
884
886 void norm1(const Teuchos::ArrayView<Teuchos::ScalarTraits<int>::magnitudeType> &/* norms */) const { XPETRA_MONITOR("EpetraIntVectorT::norm1"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); }
887
889 void norm2(const Teuchos::ArrayView<Teuchos::ScalarTraits<int>::magnitudeType> &/* norms */) const { XPETRA_MONITOR("EpetraIntVectorT::norm2"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); }
890
892 void normInf(const Teuchos::ArrayView<Teuchos::ScalarTraits<int>::magnitudeType> &/* norms */) const { XPETRA_MONITOR("EpetraIntVectorT::normInf"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); }
893
895 void meanValue(const Teuchos::ArrayView<int> &/* means */) const { XPETRA_MONITOR("EpetraIntVectorT::meanValue"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); }
896
898 void maxValue(const Teuchos::ArrayView<int> &/* maxs */) const { XPETRA_MONITOR("EpetraIntVectorT::maxValue"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); }
899
901 void multiply(Teuchos::ETransp /* transA */, Teuchos::ETransp /* transB */, const int &/* alpha */, const MultiVector<int,int,GlobalOrdinal,Node> &/* A */, const MultiVector<int,int,GlobalOrdinal,Node> &/* B */, const int &/* beta */) { XPETRA_MONITOR("EpetraIntVectorT::multiply"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "Not available in Epetra"); }
902
904 void elementWiseMultiply(int /* scalarAB */, const Vector<int,int,GlobalOrdinal,Node> &/* A */, const MultiVector<int,int,GlobalOrdinal,Node> &/* B */, int /* scalarThis */) {
905 XPETRA_MONITOR("EpetraIntVectorT::elementWiseMultiply");
906 TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "Xpetra_EpetraIntVector: elementWiseMultiply not implemented because Epetra_IntVector does not support this operation");
907 }
908
910
912
913
915 void replaceGlobalValue(GlobalOrdinal /* globalRow */, size_t /* vectorIndex */, const Scalar &/* value */) { XPETRA_MONITOR("EpetraIntVectorT::replaceGlobalValue"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); }
916
918 void sumIntoGlobalValue(GlobalOrdinal /* globalRow */, size_t /* vectorIndex */, const Scalar &/* value */) { XPETRA_MONITOR("EpetraIntVectorT::sumIntoGlobalValue"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); }
919
921 void replaceLocalValue(LocalOrdinal /* myRow */, size_t /* vectorIndex */, const Scalar &/* value */) { XPETRA_MONITOR("EpetraIntVectorT::replaceLocalValue"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); }
922
924 void sumIntoLocalValue(LocalOrdinal /* myRow */, size_t /* vectorIndex */, const Scalar &/* value */) { XPETRA_MONITOR("EpetraIntVectorT::sumIntoLocalValue"); TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO"); }
925
927
929
930
932 size_t getNumVectors() const { XPETRA_MONITOR("EpetraIntVectorT::getNumVectors"); return 1; }
933
934
936 size_t getLocalLength() const { return vec_->MyLength(); }
937
939 global_size_t getGlobalLength() const { return vec_->GlobalLength64(); }
940
941
943 bool isSameSize(const MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node> & vec) const {
944 XPETRA_MONITOR("EpetraIntVectorT::isSameSize");
945 const Vector<Scalar,LocalOrdinal,GlobalOrdinal,Node> *asvec = dynamic_cast<const Vector<Scalar,LocalOrdinal,GlobalOrdinal,Node>* >(&vec);
946 if(!asvec) return false;
947 auto vv = toEpetra(*asvec);
948 return ( (vec_->MyLength() == vv.MyLength()) && (getNumVectors() == vec.getNumVectors()));
949 }
950
951
953
954
956 std::string description() const {
957 XPETRA_MONITOR("EpetraIntVectorT::description");
958
959 // This implementation come from Epetra_Vector_def.hpp (without modification)
960 std::ostringstream oss;
961 oss << Teuchos::Describable::description();
962 oss << "{length="<<this->getGlobalLength()
963 << "}";
964 return oss.str();
965 }
966
968 void describe(Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel=Teuchos::Describable::verbLevel_default) const {
969 XPETRA_MONITOR("EpetraIntVectorT::describe");
970
971 // This implementation come from Tpetra_Vector_def.hpp (without modification) // JG: true?
972 using std::endl;
973 using std::setw;
974 using Teuchos::VERB_DEFAULT;
975 using Teuchos::VERB_NONE;
976 using Teuchos::VERB_LOW;
977 using Teuchos::VERB_MEDIUM;
978 using Teuchos::VERB_HIGH;
979 using Teuchos::VERB_EXTREME;
980
981 if (verbLevel > Teuchos::VERB_NONE)
982 vec_->Print(out);
983 }
984
986
987 RCP< Epetra_IntVector > getEpetra_IntVector() const { return vec_; }
988
989 const RCP<const Comm<int> > getComm() const {
990 TEUCHOS_TEST_FOR_EXCEPTION(1, Xpetra::Exceptions::NotImplemented, "TODO getComm Epetra MultiVector not implemented");
991 }
992
993 // Implementing DistObject
998
999 void doImport(const DistObject<int, int, GlobalOrdinal, Node> &source,
1000 const Import<int, GlobalOrdinal, Node> &importer, CombineMode CM) {
1001 XPETRA_MONITOR("EpetraIntVectorT::doImport");
1002
1003 XPETRA_DYNAMIC_CAST(const EpetraIntVectorT<GlobalOrdinal XPETRA_COMMA Node>, source, tSource, "Xpetra::EpetraIntVectorT::doImport only accept Xpetra::EpetraIntVectorT as input arguments.");
1004 XPETRA_DYNAMIC_CAST(const EpetraImportT<GlobalOrdinal XPETRA_COMMA Node>, importer, tImporter, "Xpetra::EpetraIntVectorT::doImport only accept Xpetra::EpetraImportT as input arguments.");
1005
1006 const Epetra_IntVector & v = *tSource.getEpetra_IntVector();
1007 int err = vec_->Import(v, *tImporter.getEpetra_Import(), toEpetra(CM));
1008 TEUCHOS_TEST_FOR_EXCEPTION(err != 0, std::runtime_error, "Catch error code returned by Epetra.");
1009 }
1010
1011 void doExport(const DistObject<int, LocalOrdinal, GlobalOrdinal, Node> &dest,
1012 const Import<int, GlobalOrdinal, Node>& importer, CombineMode CM) {
1013 XPETRA_MONITOR("EpetraIntVectorT::doExport");
1014
1015 XPETRA_DYNAMIC_CAST(const EpetraIntVectorT<GlobalOrdinal XPETRA_COMMA Node>, dest, tDest, "Xpetra::EpetraIntVectorT::doImport only accept Xpetra::EpetraIntVectorT as input arguments.");
1016 XPETRA_DYNAMIC_CAST(const EpetraImportT<GlobalOrdinal XPETRA_COMMA Node>, importer, tImporter, "Xpetra::EpetraIntVectorT::doImport only accept Xpetra::EpetraImportT as input arguments.");
1017
1018 const Epetra_IntVector & v = *tDest.getEpetra_IntVector();
1019 int err = vec_->Import(v, *tImporter.getEpetra_Import(), toEpetra(CM));
1020 TEUCHOS_TEST_FOR_EXCEPTION(err != 0, std::runtime_error, "Catch error code returned by Epetra.");
1021 }
1022
1023 void doImport(const DistObject<int, LocalOrdinal, GlobalOrdinal, Node> &source,
1024 const Export<int, GlobalOrdinal, Node>& exporter, CombineMode CM) {
1025 XPETRA_MONITOR("EpetraIntVectorT::doImport");
1026
1027 XPETRA_DYNAMIC_CAST(const EpetraIntVectorT<GlobalOrdinal XPETRA_COMMA Node>, source, tSource, "Xpetra::EpetraIntVectorT::doImport only accept Xpetra::EpetraIntVectorT as input arguments.");
1028 XPETRA_DYNAMIC_CAST(const EpetraExportT<GlobalOrdinal XPETRA_COMMA Node>, exporter, tExporter, "Xpetra::EpetraIntVectorT::doImport only accept Xpetra::EpetraImportT as input arguments.");
1029
1030 const Epetra_IntVector & v = *tSource.getEpetra_IntVector();
1031 int err = vec_->Import(v, *tExporter.getEpetra_Export(), toEpetra(CM));
1032 TEUCHOS_TEST_FOR_EXCEPTION(err != 0, std::runtime_error, "Catch error code returned by Epetra.");
1033 }
1034
1035 void doExport(const DistObject<int, LocalOrdinal, GlobalOrdinal, Node> &dest,
1036 const Export<int, GlobalOrdinal, Node>& exporter, CombineMode CM) {
1037 XPETRA_MONITOR("EpetraIntVectorT::doExport");
1038
1039 XPETRA_DYNAMIC_CAST(const EpetraIntVectorT<GlobalOrdinal XPETRA_COMMA Node>, dest, tDest, "Xpetra::EpetraIntVectorT::doImport only accept Xpetra::EpetraIntVectorT as input arguments.");
1040 XPETRA_DYNAMIC_CAST(const EpetraExportT<GlobalOrdinal XPETRA_COMMA Node>, exporter, tExporter, "Xpetra::EpetraIntVectorT::doImport only accept Xpetra::EpetraImportT as input arguments.");
1041
1042 const Epetra_IntVector & v = *tDest.getEpetra_IntVector();
1043 int err = vec_->Export(v, *tExporter.getEpetra_Export(), toEpetra(CM));
1044 TEUCHOS_TEST_FOR_EXCEPTION(err != 0, std::runtime_error, "Catch error code returned by Epetra.");
1045 }
1046
1047 void replaceMap(const RCP<const Map<int, GlobalOrdinal, Node> >& /* map */) {
1048 // do nothing
1049 }
1050
1051 protected:
1054 virtual void
1055 assign (const MultiVector<Scalar, LocalOrdinal, GlobalOrdinal, Node>& rhs)
1056 {
1057 typedef EpetraIntVectorT<GlobalOrdinal, Node> this_type;
1058 const this_type* rhsPtr = dynamic_cast<const this_type*> (&rhs);
1059 TEUCHOS_TEST_FOR_EXCEPTION(
1060 rhsPtr == NULL, std::invalid_argument, "Xpetra::MultiVector::operator=: "
1061 "The left-hand side (LHS) of the assignment has a different type than "
1062 "the right-hand side (RHS). The LHS has type Xpetra::EpetraIntVectorT "
1063 "(which means it wraps an Epetra_IntVector), but the RHS has some "
1064 "other type. This probably means that the RHS wraps either an "
1065 "Tpetra::MultiVector, or an Epetra_MultiVector. Xpetra::MultiVector "
1066 "does not currently implement assignment from a Tpetra object to an "
1067 "Epetra object, though this could be added with sufficient interest.");
1068
1069 RCP<const Epetra_IntVector> rhsImpl = rhsPtr->getEpetra_IntVector ();
1070 RCP<Epetra_IntVector> lhsImpl = this->getEpetra_IntVector ();
1071
1072 TEUCHOS_TEST_FOR_EXCEPTION(
1073 rhsImpl.is_null (), std::logic_error, "Xpetra::MultiVector::operator= "
1074 "(in Xpetra::EpetraIntVectorT::assign): *this (the right-hand side of "
1075 "the assignment) has a null RCP<Epetra_IntVector> inside. Please "
1076 "report this bug to the Xpetra developers.");
1077 TEUCHOS_TEST_FOR_EXCEPTION(
1078 lhsImpl.is_null (), std::logic_error, "Xpetra::MultiVector::operator= "
1079 "(in Xpetra::EpetraIntVectorT::assign): The left-hand side of the "
1080 "assignment has a null RCP<Epetra_IntVector> inside. Please report "
1081 "this bug to the Xpetra developers.");
1082
1083 // Epetra_IntVector's assignment operator does a deep copy.
1084 *lhsImpl = *rhsImpl;
1085 }
1086
1087
1088 private:
1090 RCP< Epetra_IntVector > vec_;
1091 };
1092#endif
1093
1094
1095} // namespace Xpetra
1096
1097#endif // XPETRA_EPETRAINTVECTOR_HPP
XPETRA_MONITOR
#define XPETRA_MONITOR(funcName)
Xpetra_ConfigDefs.hpp
Xpetra_EpetraConfigDefs.hpp
Xpetra_EpetraMap.hpp
Xpetra_EpetraMultiVector.hpp
Xpetra_Exceptions.hpp
XPETRA_DYNAMIC_CAST
#define XPETRA_DYNAMIC_CAST(type, obj, newObj, exceptionMsg)
Definition Xpetra_Exceptions.hpp:58
Xpetra_Vector.hpp
Epetra_BlockMap
Epetra_DistObject::Import
int Import(const Epetra_SrcDistObject &A, const Epetra_Import &Importer, Epetra_CombineMode CombineMode, const Epetra_OffsetIndex *Indexor=0)
Epetra_DistObject::Export
int Export(const Epetra_SrcDistObject &A, const Epetra_Import &Importer, Epetra_CombineMode CombineMode, const Epetra_OffsetIndex *Indexor=0)
Epetra_IntVector
Teuchos::Describable::verbLevel_default
static const EVerbosityLevel verbLevel_default
Teuchos::Describable::description
virtual std::string description() const
Teuchos::RCP::is_null
bool is_null() const
Xpetra::EpetraIntVectorT< long long, EpetraNode >::doImport
void doImport(const DistObject< int, LocalOrdinal, GlobalOrdinal, Node > &source, const Export< int, GlobalOrdinal, Node > &exporter, CombineMode CM)
Import data into this object using an Export object ("reverse mode").
Definition Xpetra_EpetraIntVector.hpp:1023
Xpetra::EpetraIntVectorT< long long, EpetraNode >::scale
void scale(const int &)
Scale the current values of a multi-vector, this = alpha*this.
Definition Xpetra_EpetraIntVector.hpp:857
Xpetra::EpetraIntVectorT< long long, EpetraNode >::dot
int dot(const Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &) const
TODO missing comment.
Definition Xpetra_EpetraIntVector.hpp:741
Xpetra::EpetraIntVectorT< long long, EpetraNode >::normInf
void normInf(const Teuchos::ArrayView< Teuchos::ScalarTraits< int >::magnitudeType > &) const
Compute Inf-norm of each vector in multi-vector.
Definition Xpetra_EpetraIntVector.hpp:892
Xpetra::EpetraIntVectorT< long long, EpetraNode >::replaceLocalValue
void replaceLocalValue(LocalOrdinal myRow, const Scalar &value)
Replace current value at the specified location with specified values.
Definition Xpetra_EpetraIntVector.hpp:772
Xpetra::EpetraIntVectorT< long long, EpetraNode >::getData
Teuchos::ArrayRCP< const int > getData(size_t) const
Definition Xpetra_EpetraIntVector.hpp:808
Xpetra::EpetraIntVectorT< long long, EpetraNode >::sumIntoLocalValue
void sumIntoLocalValue(LocalOrdinal myRow, const Scalar &value)
Adds specified value to existing value at the specified location.
Definition Xpetra_EpetraIntVector.hpp:775
Xpetra::EpetraIntVectorT< long long, EpetraNode >::reciprocal
void reciprocal(const MultiVector< int, int, GlobalOrdinal, Node > &)
Puts element-wise reciprocal values of input Multi-vector in target, this(i,j) = 1/A(i,...
Definition Xpetra_EpetraIntVector.hpp:849
Xpetra::EpetraIntVectorT< long long, EpetraNode >::abs
void abs(const MultiVector< int, int, GlobalOrdinal, Node > &)
Puts element-wise absolute values of input Multi-vector in target: A = abs(this)
Definition Xpetra_EpetraIntVector.hpp:841
Xpetra::EpetraIntVectorT< long long, EpetraNode >::replaceGlobalValue
void replaceGlobalValue(GlobalOrdinal, size_t, const Scalar &)
Replace value, using global (row) index.
Definition Xpetra_EpetraIntVector.hpp:915
Xpetra::EpetraIntVectorT< long long, EpetraNode >::vec_
RCP< Epetra_IntVector > vec_
The Epetra_IntVector which this class wraps.
Definition Xpetra_EpetraIntVector.hpp:1090
Xpetra::EpetraIntVectorT< long long, EpetraNode >::norm1
void norm1(const Teuchos::ArrayView< Teuchos::ScalarTraits< int >::magnitudeType > &) const
Compute 1-norm of each vector in multi-vector.
Definition Xpetra_EpetraIntVector.hpp:886
Xpetra::EpetraIntVectorT< long long, EpetraNode >::maxValue
void maxValue(const Teuchos::ArrayView< int > &) const
Compute max value of each vector in multi-vector.
Definition Xpetra_EpetraIntVector.hpp:898
Xpetra::EpetraIntVectorT< long long, EpetraNode >::describe
void describe(Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel=Teuchos::Describable::verbLevel_default) const
Print the object with some verbosity level to an FancyOStream object.
Definition Xpetra_EpetraIntVector.hpp:968
Xpetra::EpetraIntVectorT< long long, EpetraNode >::sumIntoLocalValue
void sumIntoLocalValue(LocalOrdinal, size_t, const Scalar &)
Add value to existing value, using local (row) index.
Definition Xpetra_EpetraIntVector.hpp:924
Xpetra::EpetraIntVectorT< long long, EpetraNode >::doExport
void doExport(const DistObject< int, LocalOrdinal, GlobalOrdinal, Node > &dest, const Export< int, GlobalOrdinal, Node > &exporter, CombineMode CM)
Export data into this object using an Export object ("forward mode").
Definition Xpetra_EpetraIntVector.hpp:1035
Xpetra::EpetraIntVectorT< long long, EpetraNode >::getVector
Teuchos::RCP< const Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > getVector(size_t) const
Return a Vector which is a const view of column j.
Definition Xpetra_EpetraIntVector.hpp:797
Xpetra::EpetraIntVectorT< long long, EpetraNode >::randomize
void randomize(const Scalar &, const Scalar &, bool=true)
Set multi-vector values to random numbers.
Definition Xpetra_EpetraIntVector.hpp:784
Xpetra::EpetraIntVectorT< long long, EpetraNode >::elementWiseMultiply
void elementWiseMultiply(int, const Vector< int, int, GlobalOrdinal, Node > &, const MultiVector< int, int, GlobalOrdinal, Node > &, int)
Element-wise multiply of a Vector A with a EpetraMultiVector B.
Definition Xpetra_EpetraIntVector.hpp:904
Xpetra::EpetraIntVectorT< long long, EpetraNode >::meanValue
void meanValue(const Teuchos::ArrayView< int > &) const
Compute mean (average) value of each vector in multi-vector.
Definition Xpetra_EpetraIntVector.hpp:895
Xpetra::EpetraIntVectorT< long long, EpetraNode >::meanValue
int meanValue() const
Compute mean (average) value of this Vector.
Definition Xpetra_EpetraIntVector.hpp:755
Xpetra::EpetraIntVectorT< long long, EpetraNode >::norm2
Teuchos::ScalarTraits< int >::magnitudeType norm2() const
Compute 2-norm of this Vector.
Definition Xpetra_EpetraIntVector.hpp:749
Xpetra::EpetraIntVectorT< 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_EpetraIntVector.hpp:1055
Xpetra::EpetraIntVectorT< long long, EpetraNode >::putScalar
void putScalar(const int &value)
Initialize all values in a multi-vector with specified value.
Definition Xpetra_EpetraIntVector.hpp:778
Xpetra::EpetraIntVectorT< long long, EpetraNode >::normInf
Teuchos::ScalarTraits< int >::magnitudeType normInf() const
Compute Inf-norm of this Vector.
Definition Xpetra_EpetraIntVector.hpp:752
Xpetra::EpetraIntVectorT< long long, EpetraNode >::getLocalLength
size_t getLocalLength() const
Returns the local vector length on the calling processor of vectors in the multi-vector.
Definition Xpetra_EpetraIntVector.hpp:936
Xpetra::EpetraIntVectorT< long long, EpetraNode >::sumIntoGlobalValue
void sumIntoGlobalValue(GlobalOrdinal, size_t, const Scalar &)
Add value to existing value, using global (row) index.
Definition Xpetra_EpetraIntVector.hpp:918
Xpetra::EpetraIntVectorT< long long, EpetraNode >::norm1
Teuchos::ScalarTraits< int >::magnitudeType norm1() const
Return 1-norm of this Vector.
Definition Xpetra_EpetraIntVector.hpp:745
Xpetra::EpetraIntVectorT< long long, EpetraNode >::update
void update(const int &, const MultiVector< int, int, GlobalOrdinal, Node > &, const int &, const MultiVector< int, int, GlobalOrdinal, Node > &, const int &)
Update multi-vector with scaled values of A and B, this = gamma*this + alpha*A + beta*B.
Definition Xpetra_EpetraIntVector.hpp:877
Xpetra::EpetraIntVectorT< long long, EpetraNode >::getMap
Teuchos::RCP< const Map< int, GlobalOrdinal, Node > > getMap() const
The Map describing the parallel distribution of this object.
Definition Xpetra_EpetraIntVector.hpp:994
Xpetra::EpetraIntVectorT< long long, EpetraNode >::norm2
void norm2(const Teuchos::ArrayView< Teuchos::ScalarTraits< int >::magnitudeType > &) const
Compute 2-norm of each vector in multi-vector.
Definition Xpetra_EpetraIntVector.hpp:889
Xpetra::EpetraIntVectorT< long long, EpetraNode >::sumIntoGlobalValue
void sumIntoGlobalValue(GlobalOrdinal, const Scalar &)
Adds specified value to existing value at the specified location.
Definition Xpetra_EpetraIntVector.hpp:769
Xpetra::EpetraIntVectorT< long long, EpetraNode >::replaceLocalValue
void replaceLocalValue(LocalOrdinal, size_t, const Scalar &)
Replace value, using local (row) index.
Definition Xpetra_EpetraIntVector.hpp:921
Xpetra::EpetraIntVectorT< long long, EpetraNode >::getVectorNonConst
Teuchos::RCP< Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > > getVectorNonConst(size_t)
Return a Vector which is a nonconst view of column j.
Definition Xpetra_EpetraIntVector.hpp:802
Xpetra::EpetraIntVectorT< long long, EpetraNode >::doImport
void doImport(const DistObject< int, int, GlobalOrdinal, Node > &source, const Import< int, GlobalOrdinal, Node > &importer, CombineMode CM)
Import data into this object using an Import object ("forward mode").
Definition Xpetra_EpetraIntVector.hpp:999
Xpetra::EpetraIntVectorT< long long, EpetraNode >::getGlobalLength
global_size_t getGlobalLength() const
Returns the global vector length of vectors in the multi-vector.
Definition Xpetra_EpetraIntVector.hpp:939
Xpetra::EpetraIntVectorT< long long, EpetraNode >::randomize
void randomize(bool=true)
Set multi-vector values to random numbers.
Definition Xpetra_EpetraIntVector.hpp:781
Xpetra::EpetraIntVectorT< long long, EpetraNode >::doExport
void doExport(const DistObject< int, LocalOrdinal, GlobalOrdinal, Node > &dest, const Import< int, GlobalOrdinal, Node > &importer, CombineMode CM)
Export data into this object using an Import object ("reverse mode").
Definition Xpetra_EpetraIntVector.hpp:1011
Xpetra::EpetraIntVectorT< long long, EpetraNode >::setSeed
void setSeed(unsigned int)
Set seed for Random function.
Definition Xpetra_EpetraIntVector.hpp:788
Xpetra::EpetraIntVectorT< long long, EpetraNode >::replaceGlobalValue
void replaceGlobalValue(GlobalOrdinal, const Scalar &)
Replace current value at the specified location with specified value.
Definition Xpetra_EpetraIntVector.hpp:766
Xpetra::EpetraIntVectorT< long long, EpetraNode >::description
std::string description() const
Return a simple one-line description of this object.
Definition Xpetra_EpetraIntVector.hpp:956
Xpetra::EpetraIntVectorT< long long, EpetraNode >::replaceMap
void replaceMap(const RCP< const Map< int, GlobalOrdinal, Node > > &)
Definition Xpetra_EpetraIntVector.hpp:1047
Xpetra::EpetraIntVectorT< long long, EpetraNode >::maxValue
int maxValue() const
Compute max value of this Vector.
Definition Xpetra_EpetraIntVector.hpp:758
Xpetra::EpetraIntVectorT< long long, EpetraNode >::update
void update(const int &, const MultiVector< int, int, GlobalOrdinal, Node > &, const int &)
Update multi-vector values with scaled values of A, this = beta*this + alpha*A.
Definition Xpetra_EpetraIntVector.hpp:869
Xpetra::EpetraIntVectorT< long long, EpetraNode >::isSameSize
bool isSameSize(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &vec) const
Checks to see if the local length, number of vectors and size of Scalar type match.
Definition Xpetra_EpetraIntVector.hpp:943
Xpetra::EpetraIntVectorT< long long, EpetraNode >::dot
void dot(const MultiVector< int, int, GlobalOrdinal, Node > &, const Teuchos::ArrayView< int > &) const
Computes dot product of each corresponding pair of vectors, dots[i] = this[i].dot(A[i])
Definition Xpetra_EpetraIntVector.hpp:833
Xpetra::EpetraIntVectorT< long long, EpetraNode >::EpetraIntVectorT
EpetraIntVectorT(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map, bool zeroOut=true)
Sets all vector entries to zero.
Definition Xpetra_EpetraIntVector.hpp:727
Xpetra::EpetraIntVectorT< long long, EpetraNode >::getNumVectors
size_t getNumVectors() const
Returns the number of vectors in the multi-vector.
Definition Xpetra_EpetraIntVector.hpp:932
Xpetra::EpetraIntVectorT< long long, EpetraNode >::multiply
void multiply(Teuchos::ETransp, Teuchos::ETransp, const int &, const MultiVector< int, int, GlobalOrdinal, Node > &, const MultiVector< int, int, GlobalOrdinal, Node > &, const int &)
Matrix-Matrix multiplication, this = beta*this + alpha*op(A)*op(B).
Definition Xpetra_EpetraIntVector.hpp:901
Xpetra::EpetraIntVectorT< long long, EpetraNode >::scale
void scale(Teuchos::ArrayView< const int >)
Scale the current values of a multi-vector, this[j] = alpha[j]*this[j].
Definition Xpetra_EpetraIntVector.hpp:863
Xpetra::EpetraIntVectorT::replaceMap
void replaceMap(const RCP< const Map< int, GlobalOrdinal, Node > > &map)
Definition Xpetra_EpetraIntVector.hpp:307
Xpetra::EpetraIntVectorT::abs
void abs(const MultiVector< int, int, GlobalOrdinal, Node > &A)
Puts element-wise absolute values of input Multi-vector in target: A = abs(this)
Definition Xpetra_EpetraIntVector.hpp:179
Xpetra::EpetraIntVectorT::multiply
void multiply(Teuchos::ETransp transA, Teuchos::ETransp transB, const int &alpha, const MultiVector< int, int, GlobalOrdinal, Node > &A, const MultiVector< int, int, GlobalOrdinal, Node > &B, const int &beta)
Matrix-Matrix multiplication, this = beta*this + alpha*op(A)*op(B).
Definition Xpetra_EpetraIntVector.hpp:226
Xpetra::EpetraIntVectorT::normInf
void normInf(const Teuchos::ArrayView< Teuchos::ScalarTraits< int >::magnitudeType > &norms) const
Compute Inf-norm of each vector in multi-vector.
Definition Xpetra_EpetraIntVector.hpp:217
Xpetra::EpetraIntVectorT::getMap
Teuchos::RCP< const Map< int, GlobalOrdinal, Node > > getMap() const
The Map describing the parallel distribution of this object.
Definition Xpetra_EpetraIntVector.hpp:291
Xpetra::EpetraIntVectorT::norm1
void norm1(const Teuchos::ArrayView< Teuchos::ScalarTraits< int >::magnitudeType > &norms) const
Compute 1-norm of each vector in multi-vector.
Definition Xpetra_EpetraIntVector.hpp:211
Xpetra::EpetraIntVectorT::doExport
void doExport(const DistObject< int, LocalOrdinal, GlobalOrdinal, Node > &dest, const Import< int, GlobalOrdinal, Node > &importer, CombineMode CM)
Export data into this object using an Import object ("reverse mode").
Definition Xpetra_EpetraIntVector.hpp:298
Xpetra::EpetraIntVectorT::EpetraIntVectorT
EpetraIntVectorT(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map, bool zeroOut=true)
Sets all vector entries to zero.
Definition Xpetra_EpetraIntVector.hpp:85
Xpetra::EpetraIntVectorT::getEpetra_IntVector
RCP< Epetra_IntVector > getEpetra_IntVector() const
Definition Xpetra_EpetraIntVector.hpp:284
Xpetra::EpetraIntVectorT::replaceLocalValue
void replaceLocalValue(LocalOrdinal myRow, size_t vectorIndex, const Scalar &value)
Replace value, using local (row) index.
Definition Xpetra_EpetraIntVector.hpp:246
Xpetra::EpetraIntVectorT::description
std::string description() const
Return a simple one-line description of this object.
Definition Xpetra_EpetraIntVector.hpp:275
Xpetra::EpetraIntVectorT::putScalar
void putScalar(const int &value)
Initialize all values in a multi-vector with specified value.
Definition Xpetra_EpetraIntVector.hpp:134
Xpetra::EpetraIntVectorT::maxValue
int maxValue() const
Compute max value of this Vector.
Definition Xpetra_EpetraIntVector.hpp:113
Xpetra::EpetraIntVectorT::randomize
void randomize(const Scalar &, const Scalar &, bool bUseXpetraImplementation=true)
Set multi-vector values to random numbers.
Definition Xpetra_EpetraIntVector.hpp:140
Xpetra::EpetraIntVectorT::doImport
void doImport(const DistObject< int, LocalOrdinal, GlobalOrdinal, Node > &source, const Export< int, GlobalOrdinal, Node > &exporter, CombineMode CM)
Import data into this object using an Export object ("reverse mode").
Definition Xpetra_EpetraIntVector.hpp:301
Xpetra::EpetraIntVectorT::doImport
void doImport(const DistObject< int, int, GlobalOrdinal, Node > &source, const Import< int, GlobalOrdinal, Node > &importer, CombineMode CM)
Import data into this object using an Import object ("forward mode").
Definition Xpetra_EpetraIntVector.hpp:295
Xpetra::EpetraIntVectorT::getDataNonConst
Teuchos::ArrayRCP< int > getDataNonConst(size_t j)
Definition Xpetra_EpetraIntVector.hpp:169
Xpetra::EpetraIntVectorT::getGlobalLength
global_size_t getGlobalLength() const
Returns the global vector length of vectors in the multi-vector.
Definition Xpetra_EpetraIntVector.hpp:264
Xpetra::EpetraIntVectorT::dot
void dot(const MultiVector< int, int, GlobalOrdinal, Node > &A, const Teuchos::ArrayView< int > &dots) const
Computes dot product of each corresponding pair of vectors, dots[i] = this[i].dot(A[i])
Definition Xpetra_EpetraIntVector.hpp:176
Xpetra::EpetraIntVectorT::getData
Teuchos::ArrayRCP< const int > getData(size_t j) const
Definition Xpetra_EpetraIntVector.hpp:165
Xpetra::EpetraIntVectorT::meanValue
void meanValue(const Teuchos::ArrayView< int > &means) const
Compute mean (average) value of each vector in multi-vector.
Definition Xpetra_EpetraIntVector.hpp:220
Xpetra::EpetraIntVectorT::norm1
Teuchos::ScalarTraits< int >::magnitudeType norm1() const
Return 1-norm of this Vector.
Definition Xpetra_EpetraIntVector.hpp:100
Xpetra::EpetraIntVectorT::randomize
void randomize(bool bUseXpetraImplementation=true)
Set multi-vector values to random numbers.
Definition Xpetra_EpetraIntVector.hpp:137
Xpetra::EpetraIntVectorT::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_EpetraIntVector.hpp:159
Xpetra::EpetraIntVectorT::getNumVectors
size_t getNumVectors() const
Returns the number of vectors in the multi-vector.
Definition Xpetra_EpetraIntVector.hpp:257
Xpetra::EpetraIntVectorT::getComm
const RCP< const Comm< int > > getComm() const
Definition Xpetra_EpetraIntVector.hpp:286
Xpetra::EpetraIntVectorT::sumIntoGlobalValue
void sumIntoGlobalValue(GlobalOrdinal globalRow, const Scalar &value)
Adds specified value to existing value at the specified location.
Definition Xpetra_EpetraIntVector.hpp:125
Xpetra::EpetraIntVectorT::reciprocal
void reciprocal(const MultiVector< int, int, GlobalOrdinal, Node > &A)
Puts element-wise reciprocal values of input Multi-vector in target, this(i,j) = 1/A(i,...
Definition Xpetra_EpetraIntVector.hpp:182
Xpetra::EpetraIntVectorT::assign
virtual void assign(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &rhs)
Implementation of the assignment operator (operator=); does a deep copy.
Definition Xpetra_EpetraIntVector.hpp:315
Xpetra::EpetraIntVectorT::sumIntoGlobalValue
void sumIntoGlobalValue(GlobalOrdinal globalRow, size_t vectorIndex, const Scalar &value)
Add value to existing value, using global (row) index.
Definition Xpetra_EpetraIntVector.hpp:243
Xpetra::EpetraIntVectorT::sumIntoLocalValue
void sumIntoLocalValue(LocalOrdinal myRow, const Scalar &value)
Adds specified value to existing value at the specified location.
Definition Xpetra_EpetraIntVector.hpp:131
Xpetra::EpetraIntVectorT::replaceGlobalValue
void replaceGlobalValue(GlobalOrdinal globalRow, size_t vectorIndex, const Scalar &value)
Replace value, using global (row) index.
Definition Xpetra_EpetraIntVector.hpp:240
Xpetra::EpetraIntVectorT::normInf
Teuchos::ScalarTraits< int >::magnitudeType normInf() const
Compute Inf-norm of this Vector.
Definition Xpetra_EpetraIntVector.hpp:107
Xpetra::EpetraIntVectorT::elementWiseMultiply
void elementWiseMultiply(int scalarAB, const Vector< int, int, GlobalOrdinal, Node > &A, const MultiVector< int, int, GlobalOrdinal, Node > &B, int scalarThis)
Element-wise multiply of a Vector A with a EpetraMultiVector B.
Definition Xpetra_EpetraIntVector.hpp:229
Xpetra::EpetraIntVectorT::update
void update(const int &alpha, const MultiVector< int, int, GlobalOrdinal, Node > &A, const int &beta, const MultiVector< int, int, GlobalOrdinal, Node > &B, const int &gamma)
Update multi-vector with scaled values of A and B, this = gamma*this + alpha*A + beta*B.
Definition Xpetra_EpetraIntVector.hpp:202
Xpetra::EpetraIntVectorT::isSameSize
bool isSameSize(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &vec) const
Checks to see if the local length, number of vectors and size of Scalar type match.
Definition Xpetra_EpetraIntVector.hpp:267
Xpetra::EpetraIntVectorT::replaceGlobalValue
void replaceGlobalValue(GlobalOrdinal globalRow, const Scalar &value)
Replace current value at the specified location with specified value.
Definition Xpetra_EpetraIntVector.hpp:122
Xpetra::EpetraIntVectorT::doExport
void doExport(const DistObject< int, LocalOrdinal, GlobalOrdinal, Node > &dest, const Export< int, GlobalOrdinal, Node > &exporter, CombineMode CM)
Export data into this object using an Export object ("forward mode").
Definition Xpetra_EpetraIntVector.hpp:304
Xpetra::EpetraIntVectorT::norm2
void norm2(const Teuchos::ArrayView< Teuchos::ScalarTraits< int >::magnitudeType > &norms) const
Compute 2-norm of each vector in multi-vector.
Definition Xpetra_EpetraIntVector.hpp:214
Xpetra::EpetraIntVectorT::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_EpetraIntVector.hpp:154
Xpetra::EpetraIntVectorT::update
void update(const int &alpha, const MultiVector< int, int, GlobalOrdinal, Node > &A, const int &beta)
Update multi-vector values with scaled values of A, this = beta*this + alpha*A.
Definition Xpetra_EpetraIntVector.hpp:194
Xpetra::EpetraIntVectorT::sumIntoLocalValue
void sumIntoLocalValue(LocalOrdinal myRow, size_t vectorIndex, const Scalar &value)
Add value to existing value, using local (row) index.
Definition Xpetra_EpetraIntVector.hpp:249
Xpetra::EpetraIntVectorT::getLocalLength
size_t getLocalLength() const
Returns the local vector length on the calling processor of vectors in the multi-vector.
Definition Xpetra_EpetraIntVector.hpp:261
Xpetra::EpetraIntVectorT::scale
void scale(Teuchos::ArrayView< const int > alpha)
Scale the current values of a multi-vector, this[j] = alpha[j]*this[j].
Definition Xpetra_EpetraIntVector.hpp:188
Xpetra::EpetraIntVectorT::setSeed
void setSeed(unsigned int seed)
Set seed for Random function.
Definition Xpetra_EpetraIntVector.hpp:145
Xpetra::EpetraIntVectorT::norm2
Teuchos::ScalarTraits< int >::magnitudeType norm2() const
Compute 2-norm of this Vector.
Definition Xpetra_EpetraIntVector.hpp:104
Xpetra::EpetraIntVectorT::maxValue
void maxValue(const Teuchos::ArrayView< int > &maxs) const
Compute max value of each vector in multi-vector.
Definition Xpetra_EpetraIntVector.hpp:223
Xpetra::EpetraIntVectorT::scale
void scale(const int &alpha)
Scale the current values of a multi-vector, this = alpha*this.
Definition Xpetra_EpetraIntVector.hpp:185
Xpetra::EpetraIntVectorT::dot
int dot(const Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &a) const
TODO missing comment.
Definition Xpetra_EpetraIntVector.hpp:96
Xpetra::EpetraIntVectorT::replaceLocalValue
void replaceLocalValue(LocalOrdinal myRow, const Scalar &value)
Replace current value at the specified location with specified values.
Definition Xpetra_EpetraIntVector.hpp:128
Xpetra::EpetraIntVectorT::meanValue
int meanValue() const
Compute mean (average) value of this Vector.
Definition Xpetra_EpetraIntVector.hpp:110
Xpetra::EpetraIntVectorT::describe
void describe(Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel=Teuchos::Describable::verbLevel_default) const
Print the object with some verbosity level to an FancyOStream object.
Definition Xpetra_EpetraIntVector.hpp:280
Xpetra::Exceptions::NotImplemented
Exception throws when you call an unimplemented method of Xpetra.
Definition Xpetra_Exceptions.hpp:95
Xpetra::MultiVector::getNumVectors
virtual size_t getNumVectors() const =0
Number of columns in the multivector.
TEUCHOS_TEST_FOR_EXCEPTION
#define TEUCHOS_TEST_FOR_EXCEPTION(throw_exception_test, Exception, msg)
Teuchos::as
TypeTo as(const TypeFrom &t)
TEUCHOS_UNREACHABLE_RETURN
#define TEUCHOS_UNREACHABLE_RETURN(dummyReturnVal)
Teuchos::FancyOStream
basic_FancyOStream< char > FancyOStream
Teuchos::EVerbosityLevel
EVerbosityLevel
Teuchos::VERB_MEDIUM
VERB_MEDIUM
Teuchos::VERB_HIGH
VERB_HIGH
Teuchos::VERB_EXTREME
VERB_EXTREME
Teuchos::VERB_NONE
VERB_NONE
Teuchos::VERB_DEFAULT
VERB_DEFAULT
Teuchos::VERB_LOW
VERB_LOW
Xpetra
Xpetra namespace
Definition Xpetra_BlockedCrsMatrix.hpp:89
Xpetra::global_size_t
size_t global_size_t
Global size_t object.
Definition Xpetra_ConfigDefs.hpp:174
Xpetra::toEpetra
const Epetra_CrsGraph & toEpetra(const RCP< const CrsGraph< int, GlobalOrdinal, Node > > &graph)
Definition Xpetra_EpetraCrsGraph.cpp:57
Xpetra::rcp
TEUCHOS_DEPRECATED RCP< T > rcp(T *p, Dealloc_T dealloc, bool owns_mem)
Xpetra::EpetraNode
Kokkos::Compat::KokkosSerialWrapperNode EpetraNode
Definition Xpetra_Map_decl.hpp:71
Xpetra::CombineMode
CombineMode
Xpetra::Combine Mode enumerable type.
Definition Xpetra_ConfigDefs.hpp:212