15#ifndef dealii_trilinos_vector_h
16#define dealii_trilinos_vector_h
21#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;
102 VectorReference(MPI::Vector &vector,
const size_type index);
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");
187 friend class ::TrilinosWrappers::MPI::Vector;
194# ifndef DEAL_II_WITH_64BIT_INDICES
199 gid(
const Epetra_BlockMap &map,
int i)
208 gid(
const Epetra_BlockMap &map,
int i)
495 template <
typename Number>
497 const ::Vector<Number> &v,
629 const std::shared_ptr<const Utilities::MPI::Partitioner> &partitioner,
719 template <
typename Number>
742 const ::TrilinosWrappers::SparseMatrix &matrix,
818 std::pair<size_type, size_type>
1021 std::vector<TrilinosScalar> &values)
const;
1057 template <
typename ForwardIterator,
typename OutputIterator>
1110 set(
const std::vector<size_type> &indices,
1111 const std::vector<TrilinosScalar> &values);
1118 set(
const std::vector<size_type> &indices,
1119 const ::Vector<TrilinosScalar> &values);
1136 add(
const std::vector<size_type> &indices,
1137 const std::vector<TrilinosScalar> &values);
1144 add(
const std::vector<size_type> &indices,
1145 const ::Vector<TrilinosScalar> &values);
1255 const Epetra_MultiVector &
1280 print(std::ostream &out,
1281 const unsigned int precision = 3,
1282 const bool scientific =
true,
1283 const bool across =
true)
const;
1324 <<
"An error with error number " <<
arg1
1325 <<
" occurred while calling a Trilinos function");
1336 <<
"You are trying to access element " <<
arg1
1337 <<
" of a distributed vector, but this element is not stored "
1338 <<
"on the current processor. Note: There are " <<
arg2
1339 <<
" elements stored "
1340 <<
"on the current processor from within the range [" <<
arg3 <<
','
1341 <<
arg4 <<
"] but Trilinos vectors need not store contiguous "
1342 <<
"ranges on each processor, and not every element in "
1343 <<
"this range may in fact be stored locally."
1345 <<
"A common source for this kind of problem is that you "
1346 <<
"are passing a 'fully distributed' vector into a function "
1347 <<
"that needs read access to vector elements that correspond "
1348 <<
"to degrees of freedom on ghost cells (or at least to "
1349 <<
"'locally active' degrees of freedom that are not also "
1350 <<
"'locally owned'). You need to pass a vector that has these "
1351 <<
"elements as ghost entries.");
1425 inline VectorReference::VectorReference(
MPI::Vector &vector,
1426 const size_type index)
1432 inline const VectorReference &
1433 VectorReference::operator=(
const VectorReference &r)
const
1446 inline VectorReference &
1447 VectorReference::operator=(
const VectorReference &r)
1456 inline const VectorReference &
1459 vector.set(1, &index, &value);
1465 inline const VectorReference &
1468 vector.add(1, &index, &value);
1474 inline const VectorReference &
1478 vector.add(1, &index, &new_value);
1484 inline const VectorReference &
1488 vector.set(1, &index, &new_value);
1494 inline const VectorReference &
1498 vector.set(1, &index, &new_value);
1510 return ((index >=
range.first) && (index <
range.second));
1520 "The locally owned elements have not been properly initialized!"
1521 " This happens for example if this object has been initialized"
1522 " with exactly one overlapping IndexSet."));
1542 inline internal::VectorReference
1550 inline internal::VectorReference
1568 std::vector<TrilinosScalar> &values)
const
1570 for (
size_type i = 0; i < indices.size(); ++i)
1571 values[i] =
operator()(indices[i]);
1581 for (
unsigned int i = 0; i < indices.
size(); ++i)
1584 elements[i] = (*this)[indices[i]];
1590 template <
typename ForwardIterator,
typename OutputIterator>
1593 const ForwardIterator indices_end,
1594 OutputIterator values_begin)
const
1639 Vector::set(
const std::vector<size_type> &indices,
1640 const std::vector<TrilinosScalar> &values)
1648 set(indices.size(), indices.data(),
values.data());
1654 Vector::set(
const std::vector<size_type> &indices,
1655 const ::Vector<TrilinosScalar> &values)
1663 set(indices.size(), indices.data(),
values.begin());
1670 const size_type *indices,
1686 for (
size_type i = 0; i < n_elements; ++i)
1695 const int ierr =
vector->ReplaceGlobalValues(1, &row, &values[i]);
1710 Vector::add(
const std::vector<size_type> &indices,
1711 const std::vector<TrilinosScalar> &values)
1718 add(indices.size(), indices.data(),
values.data());
1724 Vector::add(
const std::vector<size_type> &indices,
1725 const ::Vector<TrilinosScalar> &values)
1732 add(indices.size(), indices.data(),
values.begin());
1739 const size_type *indices,
1756 for (
size_type i = 0; i < n_elements; ++i)
1765 const int ierr =
vector->SumIntoGlobalValues(
1782 "Attempted to write into off-processor vector entry "
1783 "that has not be specified as being writable upon "
1796# ifndef DEAL_II_WITH_64BIT_INDICES
1797 return vector->Map().MaxAllGID() + 1 -
vector->Map().MinAllGID();
1799 return vector->Map().MaxAllGID64() + 1 -
vector->Map().MinAllGID64();
1813 inline std::pair<Vector::size_type, Vector::size_type>
1816# ifndef DEAL_II_WITH_64BIT_INDICES
1819 vector->Map().MaxMyGID() + 1;
1822 vector->Map().MinMyGID64();
1824 vector->Map().MaxMyGID64() + 1;
1830 "This function only makes sense if the elements that this "
1831 "vector stores on the current processor form a contiguous range. "
1832 "This does not appear to be the case for the current vector."));
2022 const int ierr =
vector->Update(1.0, *(v.vector), 1.0);
2037 const int ierr =
vector->Update(-1.0, *(v.vector), 1.0);
2070 const int ierr =
vector->Update(a, *(v.vector), 1.);
2091 const int ierr =
vector->Update(a, *(v.vector), b, *(
w.vector), 1.);
2113 Assert(this->
vector->Map().SameAs(v.vector->Map()) ==
true,
2115 const int ierr =
vector->Update(1., *(v.vector), s);
2144 Assert(this->
vector->Map().SameAs(v.vector->Map()) ==
true,
2146 const int ierr =
vector->Update(a, *(v.vector), s);
2154 this->
add(tmp,
true);
2183 if (
vector->Map().SameAs(v.vector->Map()) ==
false)
2185 this->
sadd(0., a, v);
2190 int ierr =
vector->Update(a, *v.vector, 0.0);
2199 inline const Epetra_MultiVector &
2202 return static_cast<const Epetra_MultiVector &
>(*vector);
2207 inline Epetra_FEVector &
2215 inline const Epetra_BlockMap &
2226 const Epetra_MpiComm *mpi_comm =
2227 dynamic_cast<const Epetra_MpiComm *
>(&
vector->Map().Comm());
2228 return mpi_comm->Comm();
2233 template <
typename number>
2235 const ::Vector<number> &v,
2272 namespace LinearOperatorImplementation
2285 template <
typename Matrix>
2289 bool omit_zeroing_entries)
2291 v.
reinit(matrix.locally_owned_range_indices(),
2292 matrix.get_mpi_communicator(),
2293 omit_zeroing_entries);
2296 template <
typename Matrix>
2300 bool omit_zeroing_entries)
2302 v.
reinit(matrix.locally_owned_domain_indices(),
2303 matrix.get_mpi_communicator(),
2304 omit_zeroing_entries);
size_type n_elements() const
real_type l1_norm() const
Vector & operator/=(const TrilinosScalar factor)
void compress(VectorOperation::values operation)
Vector(const IndexSet ¶llel_partitioning, const ::Vector< Number > &v, const MPI_Comm communicator=MPI_COMM_WORLD)
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)
void add(const std::vector< size_type > &indices, const std::vector< TrilinosScalar > &values)
std::unique_ptr< Epetra_MultiVector > nonlocal_vector
void add(const TrilinosScalar s)
void import_elements(const LinearAlgebra::ReadWriteVector< double > &rwv, const VectorOperation::values operation)
void update_ghost_values() const
MPI_Comm get_mpi_communicator() const
const Epetra_BlockMap & trilinos_partitioner() 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
size_type size() const override
bool in_local_range(const size_type index) 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
virtual void extract_subvector_to(const ArrayView< const size_type > &indices, ArrayView< TrilinosScalar > &elements) const override
Vector & operator-=(const Vector &V)
real_type l2_norm() const
Vector & operator+=(const Vector &V)
const Epetra_MultiVector & trilinos_vector() const
const internal::VectorReference const_reference
void reinit(const Vector &v, const bool omit_zeroing_entries=false, const bool allow_different_maps=false)
bool has_ghost_elements() const
std::pair< size_type, size_type > local_range() 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
Epetra_FEVector & trilinos_vector()
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)
bool operator==(const Vector &v) const
void sadd(const TrilinosScalar s, const TrilinosScalar a, const Vector &V)
Epetra_CombineMode last_action
reference operator[](const size_type index)
const value_type * const_iterator
TrilinosScalar value_type
size_type locally_owned_size() const
TrilinosScalar operator[](const size_type index) const
TrilinosScalar add_and_dot(const TrilinosScalar a, const Vector &V, const Vector &W)
Vector & operator=(const TrilinosScalar s)
Vector & operator*=(const TrilinosScalar factor)
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)
void extract_subvector_to(ForwardIterator indices_begin, const ForwardIterator indices_end, OutputIterator values_begin) const
std::size_t memory_consumption() const
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
Vector & operator=(const ::Vector< Number > &v)
const_iterator end() const
typename numbers::NumberTraits< Number >::real_type real_type
bool has_ghost_elements() const
const value_type * const_iterator
virtual size_type size() const override
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 & ExcGhostsPresent()
#define DeclException4(Exception4, type1, type2, type3, type4, outsequence)
static ::ExceptionBase & ExcDifferentParallelPartitioning()
static ::ExceptionBase & ExcAccessToNonLocalElement(size_type arg1, size_type arg2, size_type arg3, size_type arg4)
#define Assert(cond, exc)
#define AssertIsFinite(number)
#define AssertDimension(dim1, dim2)
#define AssertIndexRange(index, range)
static ::ExceptionBase & ExcTrilinosError(int arg1)
#define DeclException1(Exception1, type1, outsequence)
static ::ExceptionBase & ExcTrilinosError(int arg1)
static ::ExceptionBase & ExcMessage(std::string arg1)
#define AssertThrow(cond, exc)
types::global_dof_index size_type
double norm(const FEValuesBase< dim > &fe, const ArrayView< const std::vector< Tensor< 1, dim > > > &Du)
Tensor< 2, dim, Number > w(const Tensor< 2, dim, Number > &F, const Tensor< 2, dim, Number > &dF_dt)
SymmetricTensor< 2, dim, Number > d(const Tensor< 2, dim, Number > &F, const Tensor< 2, dim, Number > &dF_dt)
int gid(const Epetra_BlockMap &map, int i)
T sum(const T &t, const MPI_Comm mpi_communicator)
::VectorizedArray< Number, width > pow(const ::VectorizedArray< Number, width > &, const Number p)
unsigned int global_dof_index