16 #ifndef dealii_trilinos_vector_h
17 #define dealii_trilinos_vector_h
22 #ifdef DEAL_II_WITH_TRILINOS
33 # include <Epetra_ConfigDefs.h>
34 # include <Epetra_FEVector.h>
35 # include <Epetra_LocalMap.h>
36 # include <Epetra_Map.h>
37 # include <Epetra_MpiComm.h>
50 template <
typename Number>
51 class ReadWriteVector;
108 VectorReference(
const VectorReference &) =
default;
121 const VectorReference &
122 operator=(
const VectorReference &r)
const;
128 operator=(
const VectorReference &r);
133 const VectorReference &
139 const VectorReference &
145 const VectorReference &
151 const VectorReference &
157 const VectorReference &
171 <<
"An error with error number " << arg1
172 <<
" occurred while calling a Trilinos function");
194 # ifndef DEAL_II_WITH_64BIT_INDICES
199 gid(
const Epetra_BlockMap &map,
int i)
208 gid(
const Epetra_BlockMap &map,
int i)
448 const MPI_Comm &communicator = MPI_COMM_WORLD);
463 const MPI_Comm &communicator = MPI_COMM_WORLD);
481 const MPI_Comm &communicator = MPI_COMM_WORLD);
495 template <
typename Number>
497 const ::Vector<Number> &v,
498 const MPI_Comm & communicator = MPI_COMM_WORLD);
547 const bool omit_zeroing_entries =
false,
548 const bool allow_different_maps =
false);
574 const MPI_Comm &communicator = MPI_COMM_WORLD,
575 const bool omit_zeroing_entries =
false);
605 const MPI_Comm &communicator = MPI_COMM_WORLD,
606 const bool vector_writable =
false);
614 const std::shared_ptr<const Utilities::MPI::Partitioner> &partitioner,
615 const bool vector_writable =
false);
679 template <
typename Number>
783 std::pair<size_type, size_type>
986 std::vector<TrilinosScalar> &
values)
const;
1015 template <
typename ForwardIterator,
typename OutputIterator>
1018 const ForwardIterator indices_end,
1019 OutputIterator values_begin)
const;
1070 set(
const std::vector<size_type> & indices,
1071 const std::vector<TrilinosScalar> &
values);
1078 set(
const std::vector<size_type> & indices,
1079 const ::Vector<TrilinosScalar> &
values);
1096 add(
const std::vector<size_type> & indices,
1097 const std::vector<TrilinosScalar> &
values);
1104 add(
const std::vector<size_type> & indices,
1105 const ::Vector<TrilinosScalar> &
values);
1161 add(
const Vector &
V,
const bool allow_different_maps =
false);
1215 const Epetra_MultiVector &
1229 const Epetra_BlockMap &
1240 print(std::ostream & out,
1241 const unsigned int precision = 3,
1242 const bool scientific =
true,
1243 const bool across =
true)
const;
1285 <<
"An error with error number " << arg1
1286 <<
" occurred while calling a Trilinos function");
1297 <<
"You are trying to access element " << arg1
1298 <<
" of a distributed vector, but this element is not stored "
1299 <<
"on the current processor. Note: There are " << arg2
1300 <<
" elements stored "
1301 <<
"on the current processor from within the range [" << arg3 <<
','
1302 << arg4 <<
"] but Trilinos vectors need not store contiguous "
1303 <<
"ranges on each processor, and not every element in "
1304 <<
"this range may in fact be stored locally."
1306 <<
"A common source for this kind of problem is that you "
1307 <<
"are passing a 'fully distributed' vector into a function "
1308 <<
"that needs read access to vector elements that correspond "
1309 <<
"to degrees of freedom on ghost cells (or at least to "
1310 <<
"'locally active' degrees of freedom that are not also "
1311 <<
"'locally owned'). You need to pass a vector that has these "
1312 <<
"elements as ghost entries.");
1386 inline VectorReference::VectorReference(
MPI::Vector & vector,
1393 inline const VectorReference &
1394 VectorReference::operator=(
const VectorReference &r)
const
1407 inline VectorReference &
1408 VectorReference::operator=(
const VectorReference &r)
1417 inline const VectorReference &
1420 vector.set(1, &index, &value);
1426 inline const VectorReference &
1429 vector.add(1, &index, &value);
1435 inline const VectorReference &
1439 vector.add(1, &index, &new_value);
1445 inline const VectorReference &
1449 vector.set(1, &index, &new_value);
1455 inline const VectorReference &
1459 vector.set(1, &index, &new_value);
1469 std::pair<size_type, size_type> range =
local_range();
1471 return ((index >= range.first) && (index < range.second));
1481 "The locally owned elements have not been properly initialized!"
1482 " This happens for example if this object has been initialized"
1483 " with exactly one overlapping IndexSet."));
1503 inline internal::VectorReference
1511 inline internal::VectorReference
1529 std::vector<TrilinosScalar> &
values)
const
1531 for (
size_type i = 0; i < indices.size(); ++i)
1532 values[i] =
operator()(indices[i]);
1537 template <
typename ForwardIterator,
typename OutputIterator>
1540 const ForwardIterator indices_end,
1541 OutputIterator values_begin)
const
1543 while (indices_begin != indices_end)
1586 Vector::set(
const std::vector<size_type> & indices,
1587 const std::vector<TrilinosScalar> &
values)
1595 set(indices.size(), indices.data(),
values.data());
1601 Vector::set(
const std::vector<size_type> & indices,
1602 const ::Vector<TrilinosScalar> &
values)
1610 set(indices.size(), indices.data(),
values.begin());
1626 const int ierr =
vector->GlobalAssemble(Add);
1633 for (
size_type i = 0; i < n_elements; ++i)
1638 if (local_row != -1)
1639 (*vector)[0][local_row] =
values[i];
1642 const int ierr =
vector->ReplaceGlobalValues(1, &row, &
values[i]);
1657 Vector::add(
const std::vector<size_type> & indices,
1658 const std::vector<TrilinosScalar> &
values)
1665 add(indices.size(), indices.data(),
values.data());
1671 Vector::add(
const std::vector<size_type> & indices,
1672 const ::Vector<TrilinosScalar> &
values)
1679 add(indices.size(), indices.data(),
values.begin());
1697 const int ierr =
vector->GlobalAssemble(Insert);
1703 for (
size_type i = 0; i < n_elements; ++i)
1708 if (local_row != -1)
1709 (*vector)[0][local_row] +=
values[i];
1712 const int ierr =
vector->SumIntoGlobalValues(
1729 "Attempted to write into off-processor vector entry "
1730 "that has not be specified as being writable upon "
1732 (*nonlocal_vector)[0][my_row] +=
values[i];
1743 # ifndef DEAL_II_WITH_64BIT_INDICES
1744 return vector->Map().MaxAllGID() + 1 -
vector->Map().MinAllGID();
1746 return vector->Map().MaxAllGID64() + 1 -
vector->Map().MinAllGID64();
1755 return vector->Map().NumMyElements();
1768 inline std::pair<Vector::size_type, Vector::size_type>
1771 # ifndef DEAL_II_WITH_64BIT_INDICES
1774 vector->Map().MaxMyGID() + 1;
1777 vector->Map().MinMyGID64();
1779 vector->Map().MaxMyGID64() + 1;
1785 "This function only makes sense if the elements that this "
1786 "vector stores on the current processor form a contiguous range. "
1787 "This does not appear to be the case for the current vector."));
1803 const int ierr =
vector->Dot(*(vec.vector), &result);
1838 const int ierr =
vector->MinValue(&min_value);
1850 const int ierr =
vector->MaxValue(&max_value);
1864 const int ierr =
vector->Norm1(&
d);
1878 const int ierr =
vector->Norm2(&
d);
1916 const int ierr =
vector->NormInf(&
d);
1945 const int ierr =
vector->Scale(a);
1962 const int ierr =
vector->Scale(factor);
1977 const int ierr =
vector->Update(1.0, *(v.vector), 1.0);
1992 const int ierr =
vector->Update(-1.0, *(v.vector), 1.0);
2025 const int ierr =
vector->Update(a, *(v.vector), 1.);
2046 const int ierr =
vector->Update(a, *(v.vector),
b, *(
w.vector), 1.);
2068 Assert(this->
vector->Map().SameAs(v.vector->Map()) ==
true,
2070 const int ierr =
vector->Update(1., *(v.vector), s);
2099 Assert(this->
vector->Map().SameAs(v.vector->Map()) ==
true,
2101 const int ierr =
vector->Update(a, *(v.vector), s);
2109 this->
add(tmp,
true);
2123 const int ierr =
vector->Multiply(1.0, *(factors.vector), *
vector, 0.0);
2138 if (
vector->Map().SameAs(v.vector->Map()) ==
false)
2140 this->
sadd(0., a, v);
2145 int ierr =
vector->Update(a, *v.vector, 0.0);
2154 inline const Epetra_MultiVector &
2157 return static_cast<const Epetra_MultiVector &
>(*vector);
2162 inline Epetra_FEVector &
2170 inline const Epetra_BlockMap &
2178 inline const MPI_Comm &
2181 static MPI_Comm
comm;
2183 const Epetra_MpiComm *mpi_comm =
2184 dynamic_cast<const Epetra_MpiComm *
>(&
vector->Map().Comm());
2185 comm = mpi_comm->Comm();
2190 template <
typename number>
2192 const ::Vector<number> &v,
2193 const MPI_Comm & communicator)
2207 int ierr =
vector->PutScalar(s);
2229 namespace LinearOperatorImplementation
2242 template <
typename Matrix>
2246 bool omit_zeroing_entries)
2249 matrix.get_mpi_communicator(),
2250 omit_zeroing_entries);
2253 template <
typename Matrix>
2257 bool omit_zeroing_entries)
2260 matrix.get_mpi_communicator(),
2261 omit_zeroing_entries);
size_type n_elements() const
Epetra_Map make_trilinos_map(const MPI_Comm &communicator=MPI_COMM_WORLD, const bool overlapping=false) const
real_type l1_norm() const
void compress(VectorOperation::values operation)
void sadd(const TrilinosScalar s, const Vector &V)
TrilinosScalar mean_value() const
void add(const size_type n_elements, const size_type *indices, const TrilinosScalar *values)
const Epetra_BlockMap & trilinos_partitioner() const
void add(const std::vector< size_type > &indices, const std::vector< TrilinosScalar > &values)
const MPI_Comm & get_mpi_communicator() const
std::unique_ptr< Epetra_MultiVector > nonlocal_vector
Epetra_FEVector & trilinos_vector()
void add(const TrilinosScalar s)
void update_ghost_values() const
reference operator()(const size_type index)
void import_nonlocal_data_for_fe(const ::TrilinosWrappers::SparseMatrix &matrix, const Vector &vector)
std::unique_ptr< Epetra_FEVector > vector
bool in_local_range(const size_type index) const
size_type local_size() const
void print(std::ostream &out, const unsigned int precision=3, const bool scientific=true, const bool across=true) const
friend class internal::VectorReference
real_type lp_norm(const TrilinosScalar p) const
IndexSet locally_owned_elements() const
real_type l2_norm() const
const internal::VectorReference const_reference
void reinit(const Vector &v, const bool omit_zeroing_entries=false, const bool allow_different_maps=false)
Vector & operator+=(const Vector &V)
bool has_ghost_elements() const
real_type norm_sqr() const
void extract_subvector_to(const std::vector< size_type > &indices, std::vector< TrilinosScalar > &values) const
bool operator!=(const Vector &v) const
~Vector() override=default
const_iterator begin() const
real_type linfty_norm() const
internal::VectorReference reference
TrilinosScalar min() const
void set(const std::vector< size_type > &indices, const ::Vector< TrilinosScalar > &values)
void set(const size_type n_elements, const size_type *indices, const TrilinosScalar *values)
void scale(const Vector &scaling_factors)
void equ(const TrilinosScalar a, const Vector &V)
const Epetra_MultiVector & trilinos_vector() const
bool operator==(const Vector &v) const
Vector & operator/=(const TrilinosScalar factor)
void sadd(const TrilinosScalar s, const TrilinosScalar a, const Vector &V)
Epetra_CombineMode last_action
reference operator[](const size_type index)
size_type locally_owned_size() const
std::pair< size_type, size_type > local_range() const
TrilinosScalar operator[](const size_type index) const
Vector & operator*=(const TrilinosScalar factor)
Vector & operator=(const TrilinosScalar s)
Vector & operator=(const ::Vector< Number > &v)
TrilinosScalar add_and_dot(const TrilinosScalar a, const Vector &V, const Vector &W)
bool is_non_negative() const
::types::global_dof_index size_type
TrilinosScalar operator*(const Vector &vec) const
void add(const std::vector< size_type > &indices, const ::Vector< TrilinosScalar > &values)
Vector & operator-=(const Vector &V)
void extract_subvector_to(ForwardIterator indices_begin, const ForwardIterator indices_end, OutputIterator values_begin) const
std::size_t memory_consumption() const
Vector(const IndexSet ¶llel_partitioning, const ::Vector< Number > &v, const MPI_Comm &communicator=MPI_COMM_WORLD)
void add(const TrilinosScalar a, const Vector &V, const TrilinosScalar b, const Vector &W)
void add(const TrilinosScalar a, const Vector &V)
void set(const std::vector< size_type > &indices, const std::vector< TrilinosScalar > &values)
void swap(Vector &u, Vector &v)
TrilinosScalar max() const
const_iterator end() const
typename numbers::NumberTraits< Number >::real_type real_type
bool has_ghost_elements() const
const value_type * const_iterator
types::global_dof_index size_type
size_type locally_owned_size() const
static void reinit_domain_vector(const Matrix &matrix, TrilinosWrappers::MPI::Vector &v, bool omit_zeroing_entries)
static void reinit_range_vector(const Matrix &matrix, TrilinosWrappers::MPI::Vector &v, bool omit_zeroing_entries)
#define DEAL_II_DEPRECATED
#define DEAL_II_NAMESPACE_OPEN
#define DEAL_II_NAMESPACE_CLOSE
#define DeclException0(Exception0)
static ::ExceptionBase & ExcTrilinosError(int arg1)
static ::ExceptionBase & ExcTrilinosError(int arg1)
#define DeclException4(Exception4, type1, type2, type3, type4, outsequence)
#define Assert(cond, exc)
#define AssertIsFinite(number)
#define AssertDimension(dim1, dim2)
static ::ExceptionBase & ExcDifferentParallelPartitioning()
#define DeclException1(Exception1, type1, outsequence)
static ::ExceptionBase & ExcGhostsPresent()
static ::ExceptionBase & ExcMessage(std::string arg1)
static ::ExceptionBase & ExcAccessToNonLocalElement(size_type arg1, size_type arg2, size_type arg3, size_type arg4)
#define AssertThrow(cond, exc)
Expression fabs(const Expression &x)
@ matrix
Contents is actually a matrix.
SparseMatrix< double > SparseMatrix
types::global_dof_index size_type
double norm(const FEValuesBase< dim > &fe, const ArrayView< const std::vector< Tensor< 1, dim >>> &Du)
SymmetricTensor< 2, dim, Number > d(const Tensor< 2, dim, Number > &F, const Tensor< 2, dim, Number > &dF_dt)
Tensor< 2, dim, Number > w(const Tensor< 2, dim, Number > &F, const Tensor< 2, dim, Number > &dF_dt)
SymmetricTensor< 2, dim, Number > b(const Tensor< 2, dim, Number > &F)
int gid(const Epetra_BlockMap &map, int i)
T sum(const T &t, const MPI_Comm &mpi_communicator)
unsigned int global_dof_index
::VectorizedArray< Number, width > pow(const ::VectorizedArray< Number, width > &, const Number p)