46 #ifndef XPETRA_EPETRAMULTIVECTOR_HPP
47 #define XPETRA_EPETRAMULTIVECTOR_HPP
51 #include <Kokkos_Core.hpp>
52 #include <Kokkos_DualView.hpp>
56 #include "Xpetra_MultiVector.hpp"
65 #include "Epetra_SerialComm.h"
67 #include <Epetra_MultiVector.h>
68 #include <Epetra_Vector.h>
73 template <
class GlobalOrdinal,
class Node>
74 const Epetra_MultiVector &
toEpetra(
const MultiVector<double, int, GlobalOrdinal, Node> &);
75 template <
class GlobalOrdinal,
class Node>
76 Epetra_MultiVector &
toEpetra(MultiVector<double, int, GlobalOrdinal, Node> &);
77 template <
class GlobalOrdinal,
class Node>
78 RCP<MultiVector<double, int, GlobalOrdinal, Node> >
toXpetra(RCP<Epetra_MultiVector> vec);
81 #ifndef DOXYGEN_SHOULD_SKIP_THIS
82 template <
class GlobalOrdinal,
class Node>
86 template <
class EpetraGlobalOrdinal,
class Node>
88 :
public virtual MultiVector<double, int, EpetraGlobalOrdinal, Node> {
100 "Xpetra::EpetraMultiVector only available for GO=int or GO=long long with EpetraNode (Serial or OpenMP depending on configuration)");
106 "Xpetra::EpetraMultiVector only available for GO=int or GO=long long with EpetraNode (Serial or OpenMP depending on configuration)");
112 "Xpetra::EpetraMultiVector only available for GO=int or GO=long long with EpetraNode (Serial or OpenMP depending on configuration)");
144 Teuchos::RCP<const Vector<double, int, GlobalOrdinal, Node> >
getVector(
size_t j)
const {
145 return Teuchos::null;
150 return Teuchos::null;
154 Teuchos::ArrayRCP<const Scalar>
getData(
size_t j)
const {
155 return ArrayRCP<const Scalar>();
160 return ArrayRCP<Scalar>();
181 void scale(Teuchos::ArrayView<const Scalar> alpha) {}
187 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) {}
190 void norm1(
const Teuchos::ArrayView<Teuchos::ScalarTraits<Scalar>::magnitudeType> &norms)
const {}
193 void norm2(
const Teuchos::ArrayView<Teuchos::ScalarTraits<Scalar>::magnitudeType> &norms)
const {}
196 void normInf(
const Teuchos::ArrayView<Teuchos::ScalarTraits<Scalar>::magnitudeType> &norms)
const {}
199 void meanValue(
const Teuchos::ArrayView<Scalar> &means)
const {}
202 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) {}
233 void describe(Teuchos::FancyOStream &out,
const Teuchos::EVerbosityLevel verbLevel = Teuchos::Describable::verbLevel_default)
const {}
238 void randomize(
bool bUseXpetraImplementation =
false) {}
247 Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> >
getMap()
const {
return Teuchos::null; }
272 "Xpetra::EpetraMultiVector only available for GO=int or GO=long long with EpetraNode (Serial or OpenMP depending on configuration)");
292 #ifndef XPETRA_EPETRA_NO_32BIT_GLOBAL_INDICES
295 :
public virtual MultiVector<double, int, int, EpetraNode> {
307 : vec_(Teuchos::rcp(new Epetra_MultiVector(
toEpetra<
GlobalOrdinal,
Node>(map), Teuchos::as<int>(NumVectors), zeroOut))) {}
311 if (copyOrView == Teuchos::Copy)
312 vec_ = Teuchos::rcp(
new Epetra_MultiVector(toEpetra<GlobalOrdinal, Node>(source)));
317 vec_ = Teuchos::rcp(
new Epetra_MultiVector(View, toEpetra<GlobalOrdinal, Node>(source), indices, source.
getNumVectors()));
326 const std::string tfecfFuncName(
"MultiVector(ArrayOfPtrs)");
327 TEUCHOS_TEST_FOR_EXCEPTION_CLASS_FUNC(NumVectors < 1 || NumVectors != Teuchos::as<size_t>(ArrayOfPtrs.size()), std::runtime_error,
328 ": ArrayOfPtrs.size() must be strictly positive and as large as ArrayOfPtrs.");
330 #ifdef HAVE_XPETRA_DEBUG
333 size_t localLength = map->getLocalNumElements();
334 for (
int j = 0; j < ArrayOfPtrs.size(); j++) {
335 TEUCHOS_TEST_FOR_EXCEPTION_CLASS_FUNC(Teuchos::as<size_t>(ArrayOfPtrs[j].size()) != localLength, std::runtime_error,
336 ": ArrayOfPtrs[" << j <<
"].size() (== " << ArrayOfPtrs[j].size() <<
") is not equal to getLocalLength() (== " << localLength);
342 Array<const double *> arrayOfRawPtrs(ArrayOfPtrs.size());
343 for (
int i = 0; i < ArrayOfPtrs.size(); i++) {
344 arrayOfRawPtrs[i] = ArrayOfPtrs[i].getRawPtr();
346 double **rawArrayOfRawPtrs =
const_cast<double **
>(arrayOfRawPtrs.getRawPtr());
348 vec_ = Teuchos::rcp(
new Epetra_MultiVector(Copy, toEpetra<GlobalOrdinal, Node>(map), rawArrayOfRawPtrs, static_cast<int>(NumVectors)));
362 vec_->ReplaceGlobalValue(globalRow, Teuchos::as<int>(vectorIndex), value);
368 vec_->SumIntoGlobalValue(globalRow, Teuchos::as<int>(vectorIndex), value);
374 vec_->ReplaceMyValue(myRow, Teuchos::as<int>(vectorIndex), value);
380 vec_->SumIntoMyValue(myRow, Teuchos::as<int>(vectorIndex), value);
386 vec_->PutScalar(value);
395 Teuchos::RCP<const Vector<double, int, int, EpetraNode> >
getVector(
size_t j)
const;
401 Teuchos::ArrayRCP<const Scalar>
getData(
size_t j)
const {
404 double **arrayOfPointers;
406 vec_->ExtractView(&arrayOfPointers);
408 double *data = arrayOfPointers[j];
409 int localLength = vec_->MyLength();
411 return ArrayRCP<double>(data, 0, localLength,
false);
418 double **arrayOfPointers;
420 vec_->ExtractView(&arrayOfPointers);
422 double *data = arrayOfPointers[j];
423 int localLength = vec_->MyLength();
425 return ArrayRCP<double>(data, 0, localLength,
false);
438 vec_->Dot(*eA.getEpetra_MultiVector(), dots.getRawPtr());
444 vec_->Abs(toEpetra<GlobalOrdinal, Node>(A));
450 vec_->Reciprocal(toEpetra<GlobalOrdinal, Node>(A));
460 void scale(Teuchos::ArrayView<const Scalar> alpha) {
465 for (
size_t j = 0; j < numVecs; ++j) {
466 Epetra_Vector *v = (*vec_)(j);
474 vec_->Update(alpha, toEpetra<GlobalOrdinal, Node>(A), beta);
478 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) {
480 vec_->Update(alpha, toEpetra<GlobalOrdinal, Node>(A), beta, toEpetra<GlobalOrdinal, Node>(B), gamma);
484 void norm1(
const Teuchos::ArrayView<Teuchos::ScalarTraits<Scalar>::magnitudeType> &norms)
const {
486 vec_->Norm1(norms.getRawPtr());
490 void norm2(
const Teuchos::ArrayView<Teuchos::ScalarTraits<Scalar>::magnitudeType> &norms)
const {
492 vec_->Norm2(norms.getRawPtr());
496 void normInf(
const Teuchos::ArrayView<Teuchos::ScalarTraits<Scalar>::magnitudeType> &norms)
const {
498 vec_->NormInf(norms.getRawPtr());
502 void meanValue(
const Teuchos::ArrayView<Scalar> &means)
const {
504 vec_->MeanValue(means.getRawPtr());
508 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) {
516 vec_->Multiply(scalarAB, toEpetra<GlobalOrdinal, Node>(A), toEpetra<GlobalOrdinal, Node>(B), scalarThis);
527 return vec_->NumVectors();
533 return vec_->MyLength();
539 return vec_->GlobalLength64();
545 auto vv = toEpetra<GlobalOrdinal, Node>(vec);
546 return ((vec_->MyLength() == vv.MyLength()) && (vec_->NumVectors() == vv.NumVectors()));
558 TEUCHOS_UNREACHABLE_RETURN(
"TODO");
562 void describe(Teuchos::FancyOStream &out,
const Teuchos::EVerbosityLevel = Teuchos::Describable::verbLevel_default)
const {
573 if (bUseXpetraImplementation)
583 if (bUseXpetraImplementation)
588 for (
size_t i = 0; i < numVectors; i++) {
592 for (
size_t j = 0; j < myLength; j++) {
593 datai[j] = 0.5 * (maxVal - minVal) * datai[j] + 0.5 * (maxVal + minVal);
603 Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> >
getMap()
const {
605 return toXpetra<GlobalOrdinal, Node>(vec_->Map());
615 RCP<Epetra_MultiVector> v = tSource.getEpetra_MultiVector();
617 TEUCHOS_TEST_FOR_EXCEPTION(err != 0, std::runtime_error,
"Catch error code returned by Epetra is " << err);
627 RCP<Epetra_MultiVector> v = tDest.getEpetra_MultiVector();
629 TEUCHOS_TEST_FOR_EXCEPTION(err != 0, std::runtime_error,
"Catch error code returned by Epetra.");
639 RCP<Epetra_MultiVector> v = tSource.getEpetra_MultiVector();
641 TEUCHOS_TEST_FOR_EXCEPTION(err != 0, std::runtime_error,
"Catch error code returned by Epetra.");
651 RCP<Epetra_MultiVector> v = tDest.getEpetra_MultiVector();
653 TEUCHOS_TEST_FOR_EXCEPTION(err != 0, std::runtime_error,
"Catch error code returned by Epetra.");
660 if (!map.is_null()) {
665 Epetra_SerialComm SComm;
669 TEUCHOS_TEST_FOR_EXCEPTION(err != 0, std::runtime_error,
"Catch error code returned by Epetra.");
688 Teuchos::ScalarTraits<Scalar>::seedrandom(seed);
703 typedef Kokkos::View<
typename dual_view_type::t_host::data_type,
705 typename dual_view_type::t_host::device_type,
706 Kokkos::MemoryUnmanaged>
712 vec_->ExtractView(&data, &myLDA);
713 int localLength = vec_->MyLength();
717 epetra_view_type test = epetra_view_type(data, localLength, numVectors);
718 typename dual_view_type::t_host_um ret = subview(test, Kokkos::ALL(), Kokkos::ALL());
733 const this_type *rhsPtr =
dynamic_cast<const this_type *
>(&rhs);
734 TEUCHOS_TEST_FOR_EXCEPTION(
735 rhsPtr == NULL, std::invalid_argument,
736 "Xpetra::MultiVector::operator=: "
737 "The left-hand side (LHS) of the assignment has a different type than "
738 "the right-hand side (RHS). The LHS has type Xpetra::EpetraMultiVectorT "
739 "(which means it wraps an Epetra_MultiVector), but the RHS has some "
740 "other type. This probably means that the RHS wraps a Tpetra::Multi"
741 "Vector. Xpetra::MultiVector does not currently implement assignment "
742 "from a Tpetra object to an Epetra object, though this could be added "
743 "with sufficient interest.");
745 RCP<const Epetra_MultiVector> rhsImpl = rhsPtr->getEpetra_MultiVector();
748 TEUCHOS_TEST_FOR_EXCEPTION(
749 rhsImpl.is_null(), std::logic_error,
750 "Xpetra::MultiVector::operator= "
751 "(in Xpetra::EpetraMultiVectorT::assign): *this (the right-hand side of "
752 "the assignment) has a null RCP<Epetra_MultiVector> inside. Please "
753 "report this bug to the Xpetra developers.");
754 TEUCHOS_TEST_FOR_EXCEPTION(
755 lhsImpl.is_null(), std::logic_error,
756 "Xpetra::MultiVector::operator= "
757 "(in Xpetra::EpetraMultiVectorT::assign): The left-hand side of the "
758 "assignment has a null RCP<Epetra_MultiVector> inside. Please report "
759 "this bug to the Xpetra developers.");
773 #ifndef XPETRA_EPETRA_NO_64BIT_GLOBAL_INDICES
776 :
public virtual MultiVector<double, int, long long, EpetraNode> {
788 : vec_(Teuchos::rcp(new Epetra_MultiVector(
toEpetra<
GlobalOrdinal,
Node>(map), Teuchos::as<int>(NumVectors), zeroOut))) {}
798 const std::string tfecfFuncName(
"MultiVector(ArrayOfPtrs)");
799 TEUCHOS_TEST_FOR_EXCEPTION_CLASS_FUNC(NumVectors < 1 || NumVectors != Teuchos::as<size_t>(ArrayOfPtrs.size()), std::runtime_error,
800 ": ArrayOfPtrs.size() must be strictly positive and as large as ArrayOfPtrs.");
802 #ifdef HAVE_XPETRA_DEBUG
805 size_t localLength = map->getLocalNumElements();
806 for (
int j = 0; j < ArrayOfPtrs.size(); j++) {
807 TEUCHOS_TEST_FOR_EXCEPTION_CLASS_FUNC(Teuchos::as<size_t>(ArrayOfPtrs[j].size()) != localLength, std::runtime_error,
808 ": ArrayOfPtrs[" << j <<
"].size() (== " << ArrayOfPtrs[j].size() <<
") is not equal to getLocalLength() (== " << localLength);
814 Array<const double *> arrayOfRawPtrs(ArrayOfPtrs.size());
815 for (
int i = 0; i < ArrayOfPtrs.size(); i++) {
816 arrayOfRawPtrs[i] = ArrayOfPtrs[i].getRawPtr();
818 double **rawArrayOfRawPtrs =
const_cast<double **
>(arrayOfRawPtrs.getRawPtr());
820 vec_ = Teuchos::rcp(
new Epetra_MultiVector(Copy, toEpetra<GlobalOrdinal, Node>(map), rawArrayOfRawPtrs, NumVectors));
834 vec_->ReplaceGlobalValue(globalRow, Teuchos::as<int>(vectorIndex), value);
840 vec_->SumIntoGlobalValue(globalRow, Teuchos::as<int>(vectorIndex), value);
846 vec_->ReplaceMyValue(myRow, Teuchos::as<int>(vectorIndex), value);
852 vec_->SumIntoMyValue(myRow, Teuchos::as<int>(vectorIndex), value);
858 vec_->PutScalar(value);
867 Teuchos::RCP<const Vector<double, int, long long, EpetraNode> >
getVector(
size_t j)
const;
870 Teuchos::RCP<Vector<double, int, long long, EpetraNode> >
getVectorNonConst(
size_t j);
873 Teuchos::ArrayRCP<const Scalar>
getData(
size_t j)
const {
876 double **arrayOfPointers;
878 vec_->ExtractView(&arrayOfPointers);
880 double *data = arrayOfPointers[j];
881 int localLength = vec_->MyLength();
883 return ArrayRCP<double>(data, 0, localLength,
false);
890 double **arrayOfPointers;
892 vec_->ExtractView(&arrayOfPointers);
894 double *data = arrayOfPointers[j];
895 int localLength = vec_->MyLength();
897 return ArrayRCP<double>(data, 0, localLength,
false);
910 vec_->Dot(*eA.getEpetra_MultiVector(), dots.getRawPtr());
916 vec_->Abs(toEpetra<GlobalOrdinal, Node>(A));
922 vec_->Reciprocal(toEpetra<GlobalOrdinal, Node>(A));
932 void scale(Teuchos::ArrayView<const Scalar> alpha) {
937 for (
size_t j = 0; j < numVecs; ++j) {
938 Epetra_Vector *v = (*vec_)(j);
946 vec_->Update(alpha, toEpetra<GlobalOrdinal, Node>(A), beta);
950 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) {
952 vec_->Update(alpha, toEpetra<GlobalOrdinal, Node>(A), beta, toEpetra<GlobalOrdinal, Node>(B), gamma);
956 void norm1(
const Teuchos::ArrayView<Teuchos::ScalarTraits<Scalar>::magnitudeType> &norms)
const {
958 vec_->Norm1(norms.getRawPtr());
962 void norm2(
const Teuchos::ArrayView<Teuchos::ScalarTraits<Scalar>::magnitudeType> &norms)
const {
964 vec_->Norm2(norms.getRawPtr());
968 void normInf(
const Teuchos::ArrayView<Teuchos::ScalarTraits<Scalar>::magnitudeType> &norms)
const {
970 vec_->NormInf(norms.getRawPtr());
974 void meanValue(
const Teuchos::ArrayView<Scalar> &means)
const {
976 vec_->MeanValue(means.getRawPtr());
980 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) {
988 vec_->Multiply(scalarAB, toEpetra<GlobalOrdinal, Node>(A), toEpetra<GlobalOrdinal, Node>(B), scalarThis);
999 return vec_->NumVectors();
1005 return vec_->MyLength();
1011 return vec_->GlobalLength64();
1017 auto vv = toEpetra<GlobalOrdinal, Node>(vec);
1018 return ((vec_->MyLength() == vv.MyLength()) && (vec_->NumVectors() == vv.NumVectors()));
1034 void describe(Teuchos::FancyOStream &out,
const Teuchos::EVerbosityLevel = Teuchos::Describable::verbLevel_default)
const {
1045 if (bUseXpetraImplementation)
1055 if (bUseXpetraImplementation)
1060 for (
size_t i = 0; i < numVectors; i++) {
1064 for (
size_t j = 0; j < myLength; j++) {
1065 datai[j] = 0.5 * (maxVal - minVal) * datai[j] + 0.5 * (maxVal + minVal);
1075 Teuchos::RCP<const Map<LocalOrdinal, GlobalOrdinal, Node> >
getMap()
const {
1077 return toXpetra<GlobalOrdinal, Node>(vec_->Map());
1087 RCP<Epetra_MultiVector> v = tSource.getEpetra_MultiVector();
1089 TEUCHOS_TEST_FOR_EXCEPTION(err != 0, std::runtime_error,
"Catch error code returned by Epetra is " << err);
1099 RCP<Epetra_MultiVector> v = tDest.getEpetra_MultiVector();
1101 TEUCHOS_TEST_FOR_EXCEPTION(err != 0, std::runtime_error,
"Catch error code returned by Epetra.");
1111 RCP<Epetra_MultiVector> v = tSource.getEpetra_MultiVector();
1113 TEUCHOS_TEST_FOR_EXCEPTION(err != 0, std::runtime_error,
"Catch error code returned by Epetra.");
1123 RCP<Epetra_MultiVector> v = tDest.getEpetra_MultiVector();
1125 TEUCHOS_TEST_FOR_EXCEPTION(err != 0, std::runtime_error,
"Catch error code returned by Epetra.");
1132 if (!map.is_null()) {
1137 Epetra_SerialComm SComm;
1141 TEUCHOS_TEST_FOR_EXCEPTION(err != 0, std::runtime_error,
"Catch error code returned by Epetra.");
1160 Teuchos::ScalarTraits<Scalar>::seedrandom(seed);
1161 vec_->SetSeed(seed);
1172 const this_type *rhsPtr =
dynamic_cast<const this_type *
>(&rhs);
1173 TEUCHOS_TEST_FOR_EXCEPTION(
1174 rhsPtr == NULL, std::invalid_argument,
1175 "Xpetra::MultiVector::operator=: "
1176 "The left-hand side (LHS) of the assignment has a different type than "
1177 "the right-hand side (RHS). The LHS has type Xpetra::EpetraMultiVectorT "
1178 "(which means it wraps an Epetra_MultiVector), but the RHS has some "
1179 "other type. This probably means that the RHS wraps a Tpetra::Multi"
1180 "Vector. Xpetra::MultiVector does not currently implement assignment "
1181 "from a Tpetra object to an Epetra object, though this could be added "
1182 "with sufficient interest.");
1184 RCP<const Epetra_MultiVector> rhsImpl = rhsPtr->getEpetra_MultiVector();
1187 TEUCHOS_TEST_FOR_EXCEPTION(
1188 rhsImpl.is_null(), std::logic_error,
1189 "Xpetra::MultiVector::operator= "
1190 "(in Xpetra::EpetraMultiVectorT::assign): *this (the right-hand side of "
1191 "the assignment) has a null RCP<Epetra_MultiVector> inside. Please "
1192 "report this bug to the Xpetra developers.");
1193 TEUCHOS_TEST_FOR_EXCEPTION(
1194 lhsImpl.is_null(), std::logic_error,
1195 "Xpetra::MultiVector::operator= "
1196 "(in Xpetra::EpetraMultiVectorT::assign): The left-hand side of the "
1197 "assignment has a null RCP<Epetra_MultiVector> inside. Please report "
1198 "this bug to the Xpetra developers.");
1201 *lhsImpl = *rhsImpl;
1215 #endif // XPETRA_EPETRAMULTIVECTOR_HPP
void dot(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const Teuchos::ArrayView< Scalar > &dots) const
Compute the dot product of each corresponding pair of vectors (columns) in A and B.
Teuchos::ArrayRCP< const Scalar > getData(size_t j) const
Const view of the local values in a particular vector of this multivector.
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: this = gamma*this + alpha*A + beta*B.
dual_view_type::t_dev_um getDeviceLocalView(Access::ReadWriteStruct) const override
EpetraMultiVectorT(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map, const Teuchos::ArrayView< const Teuchos::ArrayView< const Scalar > > &ArrayOfPtrs, size_t NumVectors)
Set multi-vector values from array of pointers using Teuchos memory management classes. (copy).
size_t getNumVectors() const
Number of columns in the multivector.
virtual void assign(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &rhs)
Implementation of the assignment operator (operator=); does a deep copy.
void update(const Scalar &alpha, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const Scalar &beta)
Update: this = beta*this + alpha*A.
dual_view_type::t_dev_um getDeviceLocalView(Access::OverwriteAllStruct) const override
std::string description() const
A simple one-line description of this object.
Teuchos::ArrayRCP< Scalar > getDataNonConst(size_t j)
View of the local values in a particular vector of this multivector.
void dot(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const Teuchos::ArrayView< Scalar > &dots) const
Compute the dot product of each corresponding pair of vectors (columns) in A and B.
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: this = gamma*this + alpha*A + beta*B.
EpetraMultiVectorT(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map, const Teuchos::ArrayView< const Teuchos::ArrayView< const Scalar > > &ArrayOfPtrs, size_t NumVectors)
Set multi-vector values from array of pointers using Teuchos memory management classes. (copy).
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: this = gamma*this + alpha*A + beta*B.
void scale(const Scalar &alpha)
Scale in place: this = alpha*this.
void scale(const Scalar &alpha)
Scale in place: this = alpha*this.
void doExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Export.
void update(const Scalar &alpha, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const Scalar &beta)
Update: this = beta*this + alpha*A.
void normInf(const Teuchos::ArrayView< Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const
Compute Inf-norm of each vector in multi-vector.
virtual ~EpetraMultiVectorT()
MultiVector destructor.
void putScalar(const Scalar &value)
Set all values in the multivector with the given value.
virtual void assign(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &rhs)
Implementation of the assignment operator (operator=); does a deep copy.
void putScalar(const Scalar &value)
Set all values in the multivector with the given value.
void randomize(bool bUseXpetraImplementation=false)
Set multi-vector values to random numbers.
void norm1(const Teuchos::ArrayView< Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const
Compute 1-norm of each vector in multi-vector.
void doImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Export< LocalOrdinal, GlobalOrdinal, Node > &exporter, CombineMode CM)
Import (using an Exporter).
size_t getLocalLength() const
Local number of rows on the calling process.
void doExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Export.
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...
virtual void assign(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &rhs)
Implementation of the assignment operator (operator=); does a deep copy.
EpetraMultiVectorT(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map, size_t NumVectors, bool zeroOut=true)
Basic MultiVector constuctor.
EpetraMultiVectorT(const RCP< Epetra_MultiVector > &vec)
EpetraMultiVectorT constructor to wrap a Epetra_MultiVector object.
Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > getMap() const
Implements DistObject interface.
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.
void doExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Export< LocalOrdinal, GlobalOrdinal, Node > &exporter, CombineMode CM)
Export (using an Importer).
Exception throws to report errors in the internal logical of the program.
virtual void Xpetra_randomize()
Set multi-vector values to random numbers. XPetra implementation.
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.
EpetraMultiVectorT(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map, size_t NumVectors, bool zeroOut=true)
Basic MultiVector constuctor.
Xpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node >::dual_view_type dual_view_type
Teuchos::ArrayRCP< const Scalar > getData(size_t j) const
Const view of the local values in a particular vector of this multivector.
void abs(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A)
Put element-wise absolute values of input Multi-vector in target: A = abs(this).
void describe(Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel=Teuchos::Describable::verbLevel_default) const
Print the object with the given verbosity level to a FancyOStream.
void scale(Teuchos::ArrayView< const Scalar > alpha)
Scale the current values of a multi-vector, this[j] = alpha[j]*this[j].
RCP< Epetra_MultiVector > getEpetra_MultiVector() const
Get the underlying Epetra multivector.
void randomize(bool bUseXpetraImplementation=false)
Set multi-vector values to random numbers.
void norm1(const Teuchos::ArrayView< Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const
Compute 1-norm of each vector in multi-vector.
void doImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Export< LocalOrdinal, GlobalOrdinal, Node > &exporter, CombineMode CM)
Import (using an Exporter).
void doExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Export.
void elementWiseMultiply(Scalar scalarAB, const Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &B, Scalar scalarThis)
Multiply a Vector A elementwise by a MultiVector B.
global_size_t getGlobalLength() const
Global number of rows in the multivector.
void norm1(const Teuchos::ArrayView< Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const
Compute 1-norm of each vector in multi-vector.
size_t getLocalLength() const
Local number of rows on the calling process.
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).
void replaceMap(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map)
Replace the underlying Map in place.
void sumIntoLocalValue(LocalOrdinal myRow, size_t vectorIndex, const Scalar &value)
Add value to existing value, using local (row) index.
EpetraMultiVectorT(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source)
MultiVector copy constructor.
RCP< Epetra_MultiVector > vec_
The Epetra_MultiVector which this class wraps.
void doImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Import.
dual_view_type::t_dev_const_um getDeviceLocalView(Access::ReadOnlyStruct) const override
Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > getMap() const
Implements DistObject interface.
void doExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Export< LocalOrdinal, GlobalOrdinal, Node > &exporter, CombineMode CM)
Export (using an Importer).
void elementWiseMultiply(Scalar scalarAB, const Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &B, Scalar scalarThis)
Multiply a Vector A elementwise by a MultiVector B.
EpetraGlobalOrdinal GlobalOrdinal
RCP< Epetra_MultiVector > getEpetra_MultiVector() const
Get the underlying Epetra multivector.
global_size_t getGlobalLength() const
Global number of rows in the multivector.
EpetraMultiVectorT(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Teuchos::DataAccess copyOrView=Teuchos::Copy)
MultiVector copy constructor.
void norm2(const Teuchos::ArrayView< Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const
RCP< Epetra_MultiVector > vec_
The Epetra_MultiVector which this class wraps.
void randomize(const Scalar &minVal, const Scalar &maxVal, bool bUseXpetraImplementation=false)
Set multi-vector values to random numbers.
std::string description() const
A simple one-line description of this object.
Teuchos::ArrayRCP< Scalar > getDataNonConst(size_t j)
View of the local values in a particular vector of this multivector.
Teuchos::ArrayRCP< const Scalar > getData(size_t j) const
Const view of the local values in a particular vector of this multivector.
size_t getNumVectors() const
Number of columns in the multivector.
void setSeed(unsigned int seed)
Set seed for Random function.
void abs(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A)
Put element-wise absolute values of input Multi-vector in target: A = abs(this).
void replaceGlobalValue(GlobalOrdinal globalRow, size_t vectorIndex, const Scalar &value)
Replace value, using global (row) index.
void randomize(bool bUseXpetraImplementation=false)
Set multi-vector values to random numbers.
void sumIntoLocalValue(LocalOrdinal myRow, size_t vectorIndex, const Scalar &value)
Add value to existing value, using local (row) index.
EpetraMultiVectorT(const RCP< Epetra_MultiVector > &vec)
EpetraMultiVectorT constructor to wrap a Epetra_MultiVector object.
virtual ~EpetraMultiVectorT()
MultiVector destructor.
dual_view_type::t_host_um getHostLocalView(Access::OverwriteAllStruct) const override
void replaceGlobalValue(GlobalOrdinal globalRow, size_t vectorIndex, const Scalar &value)
Replace value, using global (row) index.
EpetraMultiVectorT(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map, const Teuchos::ArrayView< const Teuchos::ArrayView< const Scalar > > &ArrayOfPtrs, size_t NumVectors)
Set multi-vector values from array of pointers using Teuchos memory management classes. (copy).
void replaceLocalValue(LocalOrdinal myRow, size_t vectorIndex, const Scalar &value)
Replace value, using local (row) index.
EpetraMultiVectorT(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map, size_t NumVectors, bool zeroOut=true)
Basic MultiVector constuctor.
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...
void doImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Export< LocalOrdinal, GlobalOrdinal, Node > &exporter, CombineMode CM)
Import (using an Exporter).
size_t getNumVectors() const
Number of columns in the multivector.
void norm2(const Teuchos::ArrayView< Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const
void replaceMap(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map)
Replace the underlying Map in place.
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).
global_size_t getGlobalLength() const
Global number of rows in the multivector.
Exception throws when you call an unimplemented method of Xpetra.
#define XPETRA_DYNAMIC_CAST(type, obj, newObj, exceptionMsg)
size_t global_size_t
Global size_t object.
void doExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Export< LocalOrdinal, GlobalOrdinal, Node > &exporter, CombineMode CM)
Export (using an Importer).
dual_view_type::t_host_um getHostLocalView(Access::ReadWriteStruct) const override
void update(const Scalar &alpha, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const Scalar &beta)
Update: this = beta*this + alpha*A.
void replaceLocalValue(LocalOrdinal myRow, size_t vectorIndex, const Scalar &value)
Replace value, using local (row) index.
void sumIntoGlobalValue(GlobalOrdinal globalRow, size_t vectorIndex, const Scalar &value)
Add value to existing value, using global (row) index.
size_t getLocalLength() const
Local number of rows on the calling process.
void sumIntoLocalValue(LocalOrdinal myRow, size_t vectorIndex, const Scalar &value)
Add value to existing value, using local (row) index.
Kokkos::DualView< impl_scalar_type **, Kokkos::LayoutStride, typename node_type::device_type, Kokkos::MemoryUnmanaged > dual_view_type
void abs(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A)
Put element-wise absolute values of input Multi-vector in target: A = abs(this).
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.
void sumIntoGlobalValue(GlobalOrdinal globalRow, size_t vectorIndex, const Scalar &value)
Add value to existing value, using global (row) index.
void setSeed(unsigned int seed)
Set seed for Random function.
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).
RCP< const CrsGraph< int, GlobalOrdinal, Node > > toXpetra(const Epetra_CrsGraph &g)
Tpetra::KokkosCompat::KokkosSerialWrapperNode EpetraNode
Teuchos::RCP< const Vector< double, int, GlobalOrdinal, Node > > getVector(size_t j) const
Return a Vector which is a const view of column j.
void normInf(const Teuchos::ArrayView< Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const
Compute Inf-norm of each vector in multi-vector.
void scale(Teuchos::ArrayView< const Scalar > alpha)
Scale the current values of a multi-vector, this[j] = alpha[j]*this[j].
void replaceGlobalValue(GlobalOrdinal globalRow, size_t vectorIndex, const Scalar &value)
Replace value, using global (row) index.
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...
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).
EpetraMultiVectorT(const RCP< Epetra_MultiVector > &vec)
EpetraMultiVectorT constructor to wrap a Epetra_MultiVector object.
void putScalar(const Scalar &value)
Set all values in the multivector with the given value.
void scale(Teuchos::ArrayView< const Scalar > alpha)
Scale the current values of a multi-vector, this[j] = alpha[j]*this[j].
EpetraMultiVectorT(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Teuchos::DataAccess copyOrView=Teuchos::Copy)
MultiVector copy constructor.
void setSeed(unsigned int seed)
Set seed for Random function.
void replaceLocalValue(LocalOrdinal myRow, size_t vectorIndex, const Scalar &value)
Replace value, using local (row) index.
void scale(const Scalar &alpha)
Scale in place: this = alpha*this.
Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > getMap() const
Implements DistObject interface.
virtual ~EpetraMultiVectorT()
MultiVector destructor.
void doImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Import.
bool isSameSize(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &vec) const
void describe(Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel=Teuchos::Describable::verbLevel_default) const
Print the object with the given verbosity level to a FancyOStream.
void replaceMap(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map)
Replace the underlying Map in place.
RCP< Epetra_MultiVector > getEpetra_MultiVector() const
Get the underlying Epetra multivector.
CombineMode
Xpetra::Combine Mode enumerable type.
Teuchos::ArrayRCP< Scalar > getDataNonConst(size_t j)
View of the local values in a particular vector of this multivector.
constexpr struct ReadWriteStruct ReadWrite
#define XPETRA_MONITOR(funcName)
void normInf(const Teuchos::ArrayView< Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const
Compute Inf-norm of each vector in multi-vector.
virtual dual_view_type::t_host_const_um getHostLocalView(Access::ReadOnlyStruct) const
void elementWiseMultiply(Scalar scalarAB, const Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &B, Scalar scalarThis)
Multiply a Vector A elementwise by a MultiVector B.
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).
void doImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Import.
virtual size_t getNumVectors() const =0
Number of columns in the multivector.
void randomize(const Scalar &minVal, const Scalar &maxVal, bool bUseXpetraImplementation=false)
Set multi-vector values to random numbers.
void sumIntoGlobalValue(GlobalOrdinal globalRow, size_t vectorIndex, const Scalar &value)
Add value to existing value, using global (row) index.
void randomize(const Scalar &minVal, const Scalar &maxVal, bool bUseXpetraImplementation=false)
Set multi-vector values to random numbers.
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).
std::string description() const
A simple one-line description of this object.
const Epetra_CrsGraph & toEpetra(const RCP< const CrsGraph< int, GlobalOrdinal, Node > > &graph)
Teuchos::RCP< Vector< double, int, GlobalOrdinal, Node > > getVectorNonConst(size_t j)
Return a Vector which is a nonconst view of column j.
void dot(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A, const Teuchos::ArrayView< Scalar > &dots) const
Compute the dot product of each corresponding pair of vectors (columns) in A and B.
void norm2(const Teuchos::ArrayView< Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const
dual_view_type::t_host_const_um getHostLocalView(Access::ReadOnlyStruct) const override