46 #ifndef XPETRA_EPETRAMULTIVECTOR_HPP 47 #define XPETRA_EPETRAMULTIVECTOR_HPP 51 #ifdef HAVE_XPETRA_KOKKOS_REFACTOR 52 #include <Kokkos_Core.hpp> 53 #include <Kokkos_DualView.hpp> 68 #include "Epetra_SerialComm.h" 70 #include <Epetra_MultiVector.h> 75 template<
class GlobalOrdinal,
class Node>
76 const Epetra_MultiVector &
toEpetra(
const MultiVector<double,int,GlobalOrdinal,Node> &);
77 template<
class GlobalOrdinal,
class Node>
78 Epetra_MultiVector &
toEpetra(MultiVector<double, int,GlobalOrdinal,Node> &);
79 template<
class GlobalOrdinal,
class Node>
80 RCP<MultiVector<double, int, GlobalOrdinal, Node> >
toXpetra(RCP<Epetra_MultiVector> vec);
83 #ifndef DOXYGEN_SHOULD_SKIP_THIS 84 template<
class GlobalOrdinal,
class Node>
class EpetraVectorT;
87 template<
class EpetraGlobalOrdinal,
class Node>
89 :
public virtual MultiVector<double, int, EpetraGlobalOrdinal, Node>
103 "Xpetra::EpetraMultiVector only available for GO=int or GO=long long with EpetraNode (Serial or OpenMP depending on configuration)");
109 "Xpetra::EpetraMultiVector only available for GO=int or GO=long long with EpetraNode (Serial or OpenMP depending on configuration)");
115 "Xpetra::EpetraMultiVector only available for GO=int or GO=long long with EpetraNode (Serial or OpenMP depending on configuration)");
147 Teuchos::RCP< const Vector< double, int, GlobalOrdinal, Node > >
getVector(
size_t j)
const {
148 return Teuchos::null;
153 return Teuchos::null;
157 Teuchos::ArrayRCP< const Scalar >
getData(
size_t j)
const {
158 return ArrayRCP<const Scalar>();
163 return ArrayRCP<Scalar>();
184 void scale (Teuchos::ArrayView< const Scalar > alpha) { }
190 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) { }
193 void norm1(
const Teuchos::ArrayView< Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms)
const { }
196 void norm2(
const Teuchos::ArrayView< Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms)
const { }
199 void normInf(
const Teuchos::ArrayView< Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms)
const { }
202 void meanValue(
const Teuchos::ArrayView< Scalar > &means)
const { }
205 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) { }
244 Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > >
getMap()
const {
return Teuchos::null; }
269 "Xpetra::EpetraMultiVector only available for GO=int or GO=long long with EpetraNode (Serial or OpenMP depending on configuration)");
278 #ifdef HAVE_XPETRA_KOKKOS_REFACTOR 292 template<
class TargetDeviceType>
293 typename Kokkos::Impl::if_c<
294 Kokkos::Impl::is_same<
295 typename dual_view_type::t_dev_um::execution_space::memory_space,
296 typename TargetDeviceType::memory_space>::value,
297 typename dual_view_type::t_dev_um,
298 typename dual_view_type::t_host_um>::type
299 getLocalView ()
const {
300 typename Kokkos::Impl::if_c<
301 Kokkos::Impl::is_same<
302 typename dual_view_type::t_dev_um::execution_space::memory_space,
303 typename TargetDeviceType::memory_space>::value,
304 typename dual_view_type::t_dev_um,
305 typename dual_view_type::t_host_um>::type dummy;
309 typename dual_view_type::t_host_um getHostLocalView ()
const {
310 return typename dual_view_type::t_host_um();
313 typename dual_view_type::t_dev_um getDeviceLocalView()
const {
314 throw std::runtime_error(
"Epetra does not support device views!");
315 typename dual_view_type::t_dev_um ret;
332 #ifndef XPETRA_EPETRA_NO_32BIT_GLOBAL_INDICES 335 :
public virtual MultiVector<double, int, int, EpetraNode>
349 : vec_(Teuchos::rcp(new Epetra_MultiVector(
toEpetra<
GlobalOrdinal,
Node>(map), Teuchos::as<int>(NumVectors), zeroOut))) { }
359 const std::string tfecfFuncName(
"MultiVector(ArrayOfPtrs)");
360 TEUCHOS_TEST_FOR_EXCEPTION_CLASS_FUNC(NumVectors < 1 || NumVectors != Teuchos::as<size_t>(ArrayOfPtrs.size()), std::runtime_error,
361 ": ArrayOfPtrs.size() must be strictly positive and as large as ArrayOfPtrs.");
363 #ifdef HAVE_XPETRA_DEBUG 366 size_t localLength = map->getNodeNumElements();
367 for(
int j=0; j<ArrayOfPtrs.size(); j++) {
368 TEUCHOS_TEST_FOR_EXCEPTION_CLASS_FUNC(Teuchos::as<size_t>(ArrayOfPtrs[j].size()) != localLength, std::runtime_error,
369 ": ArrayOfPtrs[" << j <<
"].size() (== " << ArrayOfPtrs[j].size() <<
370 ") is not equal to getLocalLength() (== " << localLength);
377 Array<const double*> arrayOfRawPtrs(ArrayOfPtrs.size());
378 for(
int i=0; i<ArrayOfPtrs.size(); i++) {
379 arrayOfRawPtrs[i] = ArrayOfPtrs[i].getRawPtr();
381 double** rawArrayOfRawPtrs =
const_cast<double**
>(arrayOfRawPtrs.getRawPtr());
383 vec_ = Teuchos::rcp(
new Epetra_MultiVector(Copy, toEpetra<GlobalOrdinal,Node>(map), rawArrayOfRawPtrs, NumVectors));
415 Teuchos::RCP< const Vector< double, int, int, EpetraNode > >
getVector(
size_t j)
const;
418 Teuchos::RCP< Vector< double, int, int, EpetraNode > >
getVectorNonConst(
size_t j);
421 Teuchos::ArrayRCP< const Scalar >
getData(
size_t j)
const {
424 double ** arrayOfPointers;
426 vec_->ExtractView(&arrayOfPointers);
428 double * data = arrayOfPointers[j];
429 int localLength = vec_->MyLength();
431 return ArrayRCP<double>(data, 0, localLength,
false);
438 double ** arrayOfPointers;
440 vec_->ExtractView(&arrayOfPointers);
442 double * data = arrayOfPointers[j];
443 int localLength = vec_->MyLength();
445 return ArrayRCP<double>(data, 0, localLength,
false);
458 vec_->Dot(*eA.getEpetra_MultiVector(), dots.getRawPtr());
471 void scale (Teuchos::ArrayView< const Scalar > alpha) {
476 for (
size_t j = 0; j < numVecs; ++j) {
477 vec_->Scale (alpha[j]);
485 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) {
XPETRA_MONITOR(
"EpetraMultiVectorT::update"); vec_->Update(alpha, toEpetra<GlobalOrdinal,Node>(A), beta, toEpetra<GlobalOrdinal,Node>(B), gamma); }
488 void norm1(
const Teuchos::ArrayView< Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms)
const {
XPETRA_MONITOR(
"EpetraMultiVectorT::norm1"); vec_->Norm1(norms.getRawPtr()); }
491 void norm2(
const Teuchos::ArrayView< Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms)
const {
XPETRA_MONITOR(
"EpetraMultiVectorT::norm2"); vec_->Norm2(norms.getRawPtr()); }
494 void normInf(
const Teuchos::ArrayView< Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms)
const {
XPETRA_MONITOR(
"EpetraMultiVectorT::normInf"); vec_->NormInf(norms.getRawPtr()); }
497 void meanValue(
const Teuchos::ArrayView< Scalar > &means)
const {
XPETRA_MONITOR(
"EpetraMultiVectorT::meanValue"); vec_->MeanValue(means.getRawPtr()); }
500 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) {
XPETRA_MONITOR(
"EpetraMultiVectorT::multiply"); vec_->Multiply(
toEpetra(transA),
toEpetra(transB), alpha,
toEpetra(A),
toEpetra(B), beta); }
503 void elementWiseMultiply(
Scalar scalarAB,
const Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A,
const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &B,
Scalar scalarThis) {
XPETRA_MONITOR(
"EpetraMultiVectorT::elementWiseMultiply"); vec_->Multiply(scalarAB, toEpetra<GlobalOrdinal,Node>(A), toEpetra<GlobalOrdinal,Node>(B), scalarThis); }
532 void describe(Teuchos::FancyOStream &out,
const Teuchos::EVerbosityLevel verbLevel=Teuchos::Describable::verbLevel_default)
const {
543 if (bUseXpetraImplementation)
553 Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > >
getMap()
const {
XPETRA_MONITOR(
"EpetraMultiVectorT::getMap");
return toXpetra<GlobalOrdinal,Node>(vec_->Map()); }
562 RCP<Epetra_MultiVector> v = tSource.getEpetra_MultiVector();
564 TEUCHOS_TEST_FOR_EXCEPTION(err != 0, std::runtime_error,
"Catch error code returned by Epetra is " << err);
574 RCP<Epetra_MultiVector> v = tDest.getEpetra_MultiVector();
576 TEUCHOS_TEST_FOR_EXCEPTION(err != 0, std::runtime_error,
"Catch error code returned by Epetra.");
586 RCP<Epetra_MultiVector> v = tSource.getEpetra_MultiVector();
588 TEUCHOS_TEST_FOR_EXCEPTION(err != 0, std::runtime_error,
"Catch error code returned by Epetra.");
598 RCP<Epetra_MultiVector> v = tDest.getEpetra_MultiVector();
600 TEUCHOS_TEST_FOR_EXCEPTION(err != 0, std::runtime_error,
"Catch error code returned by Epetra.");
607 if (!map.is_null()) {
612 Epetra_SerialComm SComm;
616 TEUCHOS_TEST_FOR_EXCEPTION(err != 0, std::runtime_error,
"Catch error code returned by Epetra.");
634 Teuchos::ScalarTraits< Scalar >::seedrandom(seed);
638 #ifdef HAVE_XPETRA_KOKKOS_REFACTOR 652 template<
class TargetDeviceType>
653 typename Kokkos::Impl::if_c<
654 Kokkos::Impl::is_same<
655 typename dual_view_type::t_dev_um::execution_space::memory_space,
656 typename TargetDeviceType::memory_space>::value,
657 typename dual_view_type::t_dev_um,
658 typename dual_view_type::t_host_um>::type
659 getLocalView ()
const {
663 typename dual_view_type::t_host_um getHostLocalView ()
const {
664 typedef Kokkos::View<
typename dual_view_type::t_host::data_type ,
666 typename dual_view_type::t_host::device_type ,
667 Kokkos::MemoryUnmanaged> epetra_view_type;
672 vec_->ExtractView(&data, &myLDA);
673 int localLength = vec_->MyLength();
677 epetra_view_type test = epetra_view_type(data, localLength, numVectors);
678 typename dual_view_type::t_host_um ret = subview(test, Kokkos::ALL(), Kokkos::ALL());
683 typename dual_view_type::t_dev_um getDeviceLocalView()
const {
684 throw std::runtime_error(
"Epetra does not support device views!");
685 typename dual_view_type::t_dev_um ret;
699 const this_type* rhsPtr =
dynamic_cast<const this_type*
> (&rhs);
700 TEUCHOS_TEST_FOR_EXCEPTION(
701 rhsPtr == NULL, std::invalid_argument,
"Xpetra::MultiVector::operator=: " 702 "The left-hand side (LHS) of the assignment has a different type than " 703 "the right-hand side (RHS). The LHS has type Xpetra::EpetraMultiVectorT " 704 "(which means it wraps an Epetra_MultiVector), but the RHS has some " 705 "other type. This probably means that the RHS wraps a Tpetra::Multi" 706 "Vector. Xpetra::MultiVector does not currently implement assignment " 707 "from a Tpetra object to an Epetra object, though this could be added " 708 "with sufficient interest.");
710 RCP<const Epetra_MultiVector> rhsImpl = rhsPtr->getEpetra_MultiVector ();
713 TEUCHOS_TEST_FOR_EXCEPTION(
714 rhsImpl.is_null (), std::logic_error,
"Xpetra::MultiVector::operator= " 715 "(in Xpetra::EpetraMultiVectorT::assign): *this (the right-hand side of " 716 "the assignment) has a null RCP<Epetra_MultiVector> inside. Please " 717 "report this bug to the Xpetra developers.");
718 TEUCHOS_TEST_FOR_EXCEPTION(
719 lhsImpl.is_null (), std::logic_error,
"Xpetra::MultiVector::operator= " 720 "(in Xpetra::EpetraMultiVectorT::assign): The left-hand side of the " 721 "assignment has a null RCP<Epetra_MultiVector> inside. Please report " 722 "this bug to the Xpetra developers.");
730 RCP< Epetra_MultiVector >
vec_;
736 #ifndef XPETRA_EPETRA_NO_64BIT_GLOBAL_INDICES 739 :
public virtual MultiVector<double, int, long long, EpetraNode>
753 : vec_(Teuchos::rcp(new Epetra_MultiVector(
toEpetra<
GlobalOrdinal,
Node>(map), Teuchos::as<int>(NumVectors), zeroOut))) { }
763 const std::string tfecfFuncName(
"MultiVector(ArrayOfPtrs)");
764 TEUCHOS_TEST_FOR_EXCEPTION_CLASS_FUNC(NumVectors < 1 || NumVectors != Teuchos::as<size_t>(ArrayOfPtrs.size()), std::runtime_error,
765 ": ArrayOfPtrs.size() must be strictly positive and as large as ArrayOfPtrs.");
767 #ifdef HAVE_XPETRA_DEBUG 770 size_t localLength = map->getNodeNumElements();
771 for(
int j=0; j<ArrayOfPtrs.size(); j++) {
772 TEUCHOS_TEST_FOR_EXCEPTION_CLASS_FUNC(Teuchos::as<size_t>(ArrayOfPtrs[j].size()) != localLength, std::runtime_error,
773 ": ArrayOfPtrs[" << j <<
"].size() (== " << ArrayOfPtrs[j].size() <<
774 ") is not equal to getLocalLength() (== " << localLength);
781 Array<const double*> arrayOfRawPtrs(ArrayOfPtrs.size());
782 for(
int i=0; i<ArrayOfPtrs.size(); i++) {
783 arrayOfRawPtrs[i] = ArrayOfPtrs[i].getRawPtr();
785 double** rawArrayOfRawPtrs =
const_cast<double**
>(arrayOfRawPtrs.getRawPtr());
787 vec_ = Teuchos::rcp(
new Epetra_MultiVector(Copy, toEpetra<GlobalOrdinal,Node>(map), rawArrayOfRawPtrs, NumVectors));
819 Teuchos::RCP< const Vector< double, int, long long, EpetraNode > >
getVector(
size_t j)
const;
822 Teuchos::RCP< Vector< double, int, long long, EpetraNode > >
getVectorNonConst(
size_t j);
825 Teuchos::ArrayRCP< const Scalar >
getData(
size_t j)
const {
828 double ** arrayOfPointers;
830 vec_->ExtractView(&arrayOfPointers);
832 double * data = arrayOfPointers[j];
833 int localLength = vec_->MyLength();
835 return ArrayRCP<double>(data, 0, localLength,
false);
842 double ** arrayOfPointers;
844 vec_->ExtractView(&arrayOfPointers);
846 double * data = arrayOfPointers[j];
847 int localLength = vec_->MyLength();
849 return ArrayRCP<double>(data, 0, localLength,
false);
862 vec_->Dot(*eA.getEpetra_MultiVector(), dots.getRawPtr());
875 void scale (Teuchos::ArrayView< const Scalar > alpha) {
880 for (
size_t j = 0; j < numVecs; ++j) {
881 vec_->Scale (alpha[j]);
889 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) {
XPETRA_MONITOR(
"EpetraMultiVectorT::update"); vec_->Update(alpha, toEpetra<GlobalOrdinal,Node>(A), beta, toEpetra<GlobalOrdinal,Node>(B), gamma); }
892 void norm1(
const Teuchos::ArrayView< Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms)
const {
XPETRA_MONITOR(
"EpetraMultiVectorT::norm1"); vec_->Norm1(norms.getRawPtr()); }
895 void norm2(
const Teuchos::ArrayView< Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms)
const {
XPETRA_MONITOR(
"EpetraMultiVectorT::norm2"); vec_->Norm2(norms.getRawPtr()); }
898 void normInf(
const Teuchos::ArrayView< Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms)
const {
XPETRA_MONITOR(
"EpetraMultiVectorT::normInf"); vec_->NormInf(norms.getRawPtr()); }
901 void meanValue(
const Teuchos::ArrayView< Scalar > &means)
const {
XPETRA_MONITOR(
"EpetraMultiVectorT::meanValue"); vec_->MeanValue(means.getRawPtr()); }
904 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) {
XPETRA_MONITOR(
"EpetraMultiVectorT::multiply"); vec_->Multiply(
toEpetra(transA),
toEpetra(transB), alpha,
toEpetra(A),
toEpetra(B), beta); }
907 void elementWiseMultiply(
Scalar scalarAB,
const Vector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A,
const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &B,
Scalar scalarThis) {
XPETRA_MONITOR(
"EpetraMultiVectorT::elementWiseMultiply"); vec_->Multiply(scalarAB, toEpetra<GlobalOrdinal,Node>(A), toEpetra<GlobalOrdinal,Node>(B), scalarThis); }
936 void describe(Teuchos::FancyOStream &out,
const Teuchos::EVerbosityLevel verbLevel=Teuchos::Describable::verbLevel_default)
const {
947 if (bUseXpetraImplementation)
957 Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > >
getMap()
const {
XPETRA_MONITOR(
"EpetraMultiVectorT::getMap");
return toXpetra<GlobalOrdinal,Node>(vec_->Map()); }
966 RCP<Epetra_MultiVector> v = tSource.getEpetra_MultiVector();
968 TEUCHOS_TEST_FOR_EXCEPTION(err != 0, std::runtime_error,
"Catch error code returned by Epetra is " << err);
978 RCP<Epetra_MultiVector> v = tDest.getEpetra_MultiVector();
980 TEUCHOS_TEST_FOR_EXCEPTION(err != 0, std::runtime_error,
"Catch error code returned by Epetra.");
990 RCP<Epetra_MultiVector> v = tSource.getEpetra_MultiVector();
992 TEUCHOS_TEST_FOR_EXCEPTION(err != 0, std::runtime_error,
"Catch error code returned by Epetra.");
1002 RCP<Epetra_MultiVector> v = tDest.getEpetra_MultiVector();
1004 TEUCHOS_TEST_FOR_EXCEPTION(err != 0, std::runtime_error,
"Catch error code returned by Epetra.");
1011 if (!map.is_null()) {
1016 Epetra_SerialComm SComm;
1020 TEUCHOS_TEST_FOR_EXCEPTION(err != 0, std::runtime_error,
"Catch error code returned by Epetra.");
1038 Teuchos::ScalarTraits< Scalar >::seedrandom(seed);
1039 vec_->SetSeed(seed);
1042 #ifdef HAVE_XPETRA_KOKKOS_REFACTOR 1056 template<
class TargetDeviceType>
1057 typename Kokkos::Impl::if_c<
1058 Kokkos::Impl::is_same<
1059 typename dual_view_type::t_dev_um::execution_space::memory_space,
1060 typename TargetDeviceType::memory_space>::value,
1061 typename dual_view_type::t_dev_um,
1062 typename dual_view_type::t_host_um>::type
1063 getLocalView ()
const {
1067 typename dual_view_type::t_host_um getHostLocalView ()
const {
1068 typedef Kokkos::View<
typename dual_view_type::t_host::data_type ,
1070 typename dual_view_type::t_host::device_type ,
1071 Kokkos::MemoryUnmanaged> epetra_view_type;
1074 double* data = NULL;
1076 vec_->ExtractView(&data, &myLDA);
1077 int localLength = vec_->MyLength();
1081 epetra_view_type test = epetra_view_type(data, localLength, numVectors);
1082 typename dual_view_type::t_host_um ret = subview(test, Kokkos::ALL(), Kokkos::ALL());
1087 typename dual_view_type::t_dev_um getDeviceLocalView()
const {
1088 throw std::runtime_error(
"Epetra does not support device views!");
1089 typename dual_view_type::t_dev_um ret;
1103 const this_type* rhsPtr =
dynamic_cast<const this_type*
> (&rhs);
1104 TEUCHOS_TEST_FOR_EXCEPTION(
1105 rhsPtr == NULL, std::invalid_argument,
"Xpetra::MultiVector::operator=: " 1106 "The left-hand side (LHS) of the assignment has a different type than " 1107 "the right-hand side (RHS). The LHS has type Xpetra::EpetraMultiVectorT " 1108 "(which means it wraps an Epetra_MultiVector), but the RHS has some " 1109 "other type. This probably means that the RHS wraps a Tpetra::Multi" 1110 "Vector. Xpetra::MultiVector does not currently implement assignment " 1111 "from a Tpetra object to an Epetra object, though this could be added " 1112 "with sufficient interest.");
1114 RCP<const Epetra_MultiVector> rhsImpl = rhsPtr->getEpetra_MultiVector ();
1117 TEUCHOS_TEST_FOR_EXCEPTION(
1118 rhsImpl.is_null (), std::logic_error,
"Xpetra::MultiVector::operator= " 1119 "(in Xpetra::EpetraMultiVectorT::assign): *this (the right-hand side of " 1120 "the assignment) has a null RCP<Epetra_MultiVector> inside. Please " 1121 "report this bug to the Xpetra developers.");
1122 TEUCHOS_TEST_FOR_EXCEPTION(
1123 lhsImpl.is_null (), std::logic_error,
"Xpetra::MultiVector::operator= " 1124 "(in Xpetra::EpetraMultiVectorT::assign): The left-hand side of the " 1125 "assignment has a null RCP<Epetra_MultiVector> inside. Please report " 1126 "this bug to the Xpetra developers.");
1129 *lhsImpl = *rhsImpl;
1143 #endif // XPETRA_EPETRAMULTIVECTOR_HPP 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).
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.
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.
RCP< Epetra_MultiVector > getEpetra_MultiVector() const
Get the underlying Epetra 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.
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.
size_t getNumVectors() const
Number of columns in the multivector.
global_size_t getGlobalLength() const
Global number of rows in the multivector.
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 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 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)
Update: this = beta*this + alpha*A.
global_size_t getGlobalLength() const
Global number of rows in the multivector.
void norm2(const Teuchos::ArrayView< Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const
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.
RCP< Epetra_MultiVector > vec_
The Epetra_MultiVector which this class wraps.
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.
virtual void Xpetra_randomize()
Set multi-vector values to random numbers. XPetra implementation.
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 describe(Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel=Teuchos::Describable::verbLevel_default) const
Print the object with the given verbosity level to a FancyOStream.
size_t getLocalLength() const
Local number of rows on the calling process.
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.
void doExport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &dest, const Export< LocalOrdinal, GlobalOrdinal, Node > &exporter, CombineMode CM)
Export (using an Importer).
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 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.
Exception throws to report errors in the internal logical of the program.
void norm1(const Teuchos::ArrayView< Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const
Compute 1-norm of each vector in multi-vector.
const Epetra_CrsGraph & toEpetra(const RCP< const CrsGraph< int, GlobalOrdinal, Node > > &graph)
EpetraMultiVectorT(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map, size_t NumVectors, bool zeroOut=true)
Basic MultiVector constuctor.
Teuchos::ArrayRCP< Scalar > getDataNonConst(size_t j)
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 scale(Teuchos::ArrayView< const Scalar > alpha)
Scale the current values of a multi-vector, this[j] = alpha[j]*this[j].
void randomize(bool bUseXpetraImplementation=false)
Set multi-vector values to random numbers.
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.
RCP< Epetra_MultiVector > getEpetra_MultiVector() const
Get the underlying Epetra multivector.
Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > getMap() const
Implements DistObject interface.
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 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).
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 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.
void norm1(const Teuchos::ArrayView< Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const
Compute 1-norm of each vector in multi-vector.
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 doImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Import.
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
std::string description() const
A simple one-line description of this object.
RCP< Epetra_MultiVector > getEpetra_MultiVector() const
Get the underlying Epetra multivector.
void setSeed(unsigned int seed)
Set seed for Random function.
Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > getMap() const
Implements DistObject interface.
size_t getLocalLength() const
Local number of rows on the calling process.
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 abs(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A)
Put element-wise absolute values of input Multi-vector in target: A = abs(this).
std::string description() const
A simple one-line description of this object.
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.
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 doImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Export< LocalOrdinal, GlobalOrdinal, Node > &exporter, CombineMode CM)
Import (using an Exporter).
void norm2(const Teuchos::ArrayView< Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const
void normInf(const Teuchos::ArrayView< Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const
Compute Inf-norm of each vector in multi-vector.
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).
void normInf(const Teuchos::ArrayView< Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const
Compute Inf-norm of each vector in multi-vector.
Exception throws when you call an unimplemented method of Xpetra.
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.
#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).
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.
void sumIntoLocalValue(LocalOrdinal myRow, size_t vectorIndex, const Scalar &value)
Add value to existing value, using local (row) index.
size_t getNumVectors() const
Number of columns in the multivector.
void abs(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &A)
Put element-wise absolute values of input Multi-vector in target: A = abs(this).
Teuchos::ArrayRCP< const Scalar > getData(size_t j) const
Const view of the local values in a particular vector of this multivector.
void sumIntoGlobalValue(GlobalOrdinal globalRow, size_t vectorIndex, const Scalar &value)
Add value to existing value, using global (row) index.
Teuchos::ArrayRCP< Scalar > getDataNonConst(size_t j)
View of the local values in a particular vector of this multivector.
global_size_t getGlobalLength() const
Global number of rows in the multivector.
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)
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 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].
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.
void normInf(const Teuchos::ArrayView< Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const
Compute Inf-norm of each vector in multi-vector.
RCP< Epetra_MultiVector > vec_
The Epetra_MultiVector which this class wraps.
virtual ~EpetraMultiVectorT()
MultiVector destructor.
size_t getLocalLength() const
Local number of rows on the calling process.
void doImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Import.
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 norm1(const Teuchos::ArrayView< Teuchos::ScalarTraits< Scalar >::magnitudeType > &norms) const
Compute 1-norm of each vector in multi-vector.
void replaceMap(const Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > &map)
Replace the underlying Map in place.
EpetraMultiVectorT(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source)
MultiVector copy constructor.
CombineMode
Xpetra::Combine Mode enumerable type.
#define XPETRA_MONITOR(funcName)
Teuchos::RCP< const Map< LocalOrdinal, GlobalOrdinal, Node > > getMap() const
Implements DistObject interface.
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).
std::string description() const
A simple one-line description of this object.
void doImport(const DistObject< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source, const Import< LocalOrdinal, GlobalOrdinal, Node > &importer, CombineMode CM)
Import.
EpetraMultiVectorT(const MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &source)
MultiVector copy constructor.
void sumIntoGlobalValue(GlobalOrdinal globalRow, size_t vectorIndex, const Scalar &value)
Add value to existing value, using global (row) index.
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).
Teuchos::RCP< Vector< double, int, GlobalOrdinal, Node > > getVectorNonConst(size_t j)
Return a Vector which is a nonconst view of column j.
Teuchos::ArrayRCP< const Scalar > getData(size_t j) const
Const view of the local values in a particular vector of this multivector.