16 #ifndef dealii_trilinos_vector_h
17 #define dealii_trilinos_vector_h
22 #ifdef DEAL_II_WITH_TRILINOS
34 # include <Epetra_ConfigDefs.h>
35 # include <Epetra_FEVector.h>
36 # include <Epetra_LocalMap.h>
37 # include <Epetra_Map.h>
38 # include <Epetra_MpiComm.h>
51 template <
typename Number>
52 class ReadWriteVector;
109 VectorReference(
const VectorReference &) =
default;
122 const VectorReference &
123 operator=(
const VectorReference &r)
const;
129 operator=(
const VectorReference &r);
134 const VectorReference &
140 const VectorReference &
146 const VectorReference &
152 const VectorReference &
158 const VectorReference &
172 <<
"An error with error number " << arg1
173 <<
" occurred while calling a Trilinos function");
195 # ifndef DEAL_II_WITH_64BIT_INDICES
200 gid(
const Epetra_BlockMap &map,
int i)
209 gid(
const Epetra_BlockMap &map,
int i)
449 const MPI_Comm communicator = MPI_COMM_WORLD);
464 const MPI_Comm communicator = MPI_COMM_WORLD);
482 const MPI_Comm communicator = MPI_COMM_WORLD);
496 template <
typename Number>
498 const ::Vector<Number> &v,
499 const MPI_Comm communicator = MPI_COMM_WORLD);
548 const bool omit_zeroing_entries =
false,
549 const bool allow_different_maps =
false);
575 const MPI_Comm communicator = MPI_COMM_WORLD,
576 const bool omit_zeroing_entries =
false);
614 const IndexSet &locally_relevant_or_ghost_entries,
615 const MPI_Comm communicator = MPI_COMM_WORLD,
616 const bool vector_writable =
false);
630 const std::shared_ptr<const Utilities::MPI::Partitioner> &partitioner,
631 const bool make_ghosted =
true,
632 const bool vector_writable =
false);
696 template <
typename Number>
795 std::pair<size_type, size_type>
998 std::vector<TrilinosScalar> &
values)
const;
1034 template <
typename ForwardIterator,
typename OutputIterator>
1037 const ForwardIterator indices_end,
1038 OutputIterator values_begin)
const;
1087 set(
const std::vector<size_type> &indices,
1088 const std::vector<TrilinosScalar> &
values);
1095 set(
const std::vector<size_type> &indices,
1096 const ::Vector<TrilinosScalar> &
values);
1113 add(
const std::vector<size_type> &indices,
1114 const std::vector<TrilinosScalar> &
values);
1121 add(
const std::vector<size_type> &indices,
1122 const ::Vector<TrilinosScalar> &
values);
1178 add(
const Vector &
V,
const bool allow_different_maps =
false);
1232 const Epetra_MultiVector &
1246 const Epetra_BlockMap &
1257 print(std::ostream &out,
1258 const unsigned int precision = 3,
1259 const bool scientific =
true,
1260 const bool across =
true)
const;
1301 <<
"An error with error number " << arg1
1302 <<
" occurred while calling a Trilinos function");
1313 <<
"You are trying to access element " << arg1
1314 <<
" of a distributed vector, but this element is not stored "
1315 <<
"on the current processor. Note: There are " << arg2
1316 <<
" elements stored "
1317 <<
"on the current processor from within the range [" << arg3 <<
','
1318 << arg4 <<
"] but Trilinos vectors need not store contiguous "
1319 <<
"ranges on each processor, and not every element in "
1320 <<
"this range may in fact be stored locally."
1322 <<
"A common source for this kind of problem is that you "
1323 <<
"are passing a 'fully distributed' vector into a function "
1324 <<
"that needs read access to vector elements that correspond "
1325 <<
"to degrees of freedom on ghost cells (or at least to "
1326 <<
"'locally active' degrees of freedom that are not also "
1327 <<
"'locally owned'). You need to pass a vector that has these "
1328 <<
"elements as ghost entries.");
1402 inline VectorReference::VectorReference(
MPI::Vector &vector,
1409 inline const VectorReference &
1410 VectorReference::operator=(
const VectorReference &r)
const
1423 inline VectorReference &
1424 VectorReference::operator=(
const VectorReference &r)
1433 inline const VectorReference &
1436 vector.set(1, &index, &value);
1442 inline const VectorReference &
1445 vector.add(1, &index, &value);
1451 inline const VectorReference &
1455 vector.add(1, &index, &new_value);
1461 inline const VectorReference &
1465 vector.set(1, &index, &new_value);
1471 inline const VectorReference &
1475 vector.set(1, &index, &new_value);
1485 std::pair<size_type, size_type> range =
local_range();
1487 return ((index >= range.first) && (index < range.second));
1497 "The locally owned elements have not been properly initialized!"
1498 " This happens for example if this object has been initialized"
1499 " with exactly one overlapping IndexSet."));
1519 inline internal::VectorReference
1527 inline internal::VectorReference
1545 std::vector<TrilinosScalar> &
values)
const
1547 for (
size_type i = 0; i < indices.size(); ++i)
1548 values[i] =
operator()(indices[i]);
1558 for (
unsigned int i = 0; i < indices.
size(); ++i)
1561 elements[i] = (*this)[indices[i]];
1567 template <
typename ForwardIterator,
typename OutputIterator>
1570 const ForwardIterator indices_end,
1571 OutputIterator values_begin)
const
1573 while (indices_begin != indices_end)
1616 Vector::set(
const std::vector<size_type> &indices,
1617 const std::vector<TrilinosScalar> &
values)
1625 set(indices.size(), indices.data(),
values.data());
1631 Vector::set(
const std::vector<size_type> &indices,
1632 const ::Vector<TrilinosScalar> &
values)
1640 set(indices.size(), indices.data(),
values.begin());
1656 const int ierr =
vector->GlobalAssemble(Add);
1663 for (
size_type i = 0; i < n_elements; ++i)
1668 if (local_row != -1)
1669 (*vector)[0][local_row] =
values[i];
1672 const int ierr =
vector->ReplaceGlobalValues(1, &row, &
values[i]);
1687 Vector::add(
const std::vector<size_type> &indices,
1688 const std::vector<TrilinosScalar> &
values)
1695 add(indices.size(), indices.data(),
values.data());
1701 Vector::add(
const std::vector<size_type> &indices,
1702 const ::Vector<TrilinosScalar> &
values)
1709 add(indices.size(), indices.data(),
values.begin());
1727 const int ierr =
vector->GlobalAssemble(Insert);
1733 for (
size_type i = 0; i < n_elements; ++i)
1738 if (local_row != -1)
1739 (*vector)[0][local_row] +=
values[i];
1742 const int ierr =
vector->SumIntoGlobalValues(
1759 "Attempted to write into off-processor vector entry "
1760 "that has not be specified as being writable upon "
1762 (*nonlocal_vector)[0][my_row] +=
values[i];
1773 # ifndef DEAL_II_WITH_64BIT_INDICES
1774 return vector->Map().MaxAllGID() + 1 -
vector->Map().MinAllGID();
1776 return vector->Map().MaxAllGID64() + 1 -
vector->Map().MinAllGID64();
1790 inline std::pair<Vector::size_type, Vector::size_type>
1793 # ifndef DEAL_II_WITH_64BIT_INDICES
1796 vector->Map().MaxMyGID() + 1;
1799 vector->Map().MinMyGID64();
1801 vector->Map().MaxMyGID64() + 1;
1807 "This function only makes sense if the elements that this "
1808 "vector stores on the current processor form a contiguous range. "
1809 "This does not appear to be the case for the current vector."));
1825 const int ierr =
vector->Dot(*(vec.vector), &result);
1860 const int ierr =
vector->MinValue(&min_value);
1872 const int ierr =
vector->MaxValue(&max_value);
1886 const int ierr =
vector->Norm1(&
d);
1900 const int ierr =
vector->Norm2(&
d);
1938 const int ierr =
vector->NormInf(&
d);
1967 const int ierr =
vector->Scale(a);
1984 const int ierr =
vector->Scale(factor);
1999 const int ierr =
vector->Update(1.0, *(v.vector), 1.0);
2014 const int ierr =
vector->Update(-1.0, *(v.vector), 1.0);
2047 const int ierr =
vector->Update(a, *(v.vector), 1.);
2068 const int ierr =
vector->Update(a, *(v.vector),
b, *(
w.vector), 1.);
2090 Assert(this->
vector->Map().SameAs(v.vector->Map()) ==
true,
2092 const int ierr =
vector->Update(1., *(v.vector), s);
2121 Assert(this->
vector->Map().SameAs(v.vector->Map()) ==
true,
2123 const int ierr =
vector->Update(a, *(v.vector), s);
2131 this->
add(tmp,
true);
2145 const int ierr =
vector->Multiply(1.0, *(factors.vector), *
vector, 0.0);
2160 if (
vector->Map().SameAs(v.vector->Map()) ==
false)
2162 this->
sadd(0., a, v);
2167 int ierr =
vector->Update(a, *v.vector, 0.0);
2176 inline const Epetra_MultiVector &
2179 return static_cast<const Epetra_MultiVector &
>(*vector);
2184 inline Epetra_FEVector &
2192 inline const Epetra_BlockMap &
2203 const Epetra_MpiComm *mpi_comm =
2204 dynamic_cast<const Epetra_MpiComm *
>(&
vector->Map().Comm());
2205 return mpi_comm->Comm();
2210 template <
typename number>
2212 const ::Vector<number> &v,
2213 const MPI_Comm communicator)
2227 int ierr =
vector->PutScalar(s);
2249 namespace LinearOperatorImplementation
2262 template <
typename Matrix>
2266 bool omit_zeroing_entries)
2269 matrix.get_mpi_communicator(),
2270 omit_zeroing_entries);
2273 template <
typename Matrix>
2277 bool omit_zeroing_entries)
2280 matrix.get_mpi_communicator(),
2281 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)
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)
const Epetra_BlockMap & trilinos_partitioner() const
void add(const std::vector< size_type > &indices, const std::vector< TrilinosScalar > &values)
std::unique_ptr< Epetra_MultiVector > nonlocal_vector
Epetra_FEVector & trilinos_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
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
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
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
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 & 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)
#define AssertIndexRange(index, range)
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)