16 #ifndef dealii_read_write_vector_h
17 #define dealii_read_write_vector_h
37 #ifdef DEAL_II_WITH_TRILINOS
42 # include <Epetra_MultiVector.h>
59 template <
typename,
typename>
64 # ifdef DEAL_II_WITH_PETSC
74 # ifdef DEAL_II_WITH_TRILINOS
84 # ifdef DEAL_II_WITH_CUDA
135 template <
typename Number>
204 template <
typename Number2>
207 const bool omit_zeroing_entries =
false);
220 const bool omit_zeroing_entries =
false);
223 #ifdef DEAL_II_WITH_TRILINOS
252 template <
typename Functor>
280 template <
typename Number2>
302 const ::Vector<Number> &vec,
304 const std::shared_ptr<const Utilities::MPI::CommunicationPatternBase>
305 &communication_pattern = {});
312 import(const ::Vector<Number> &
V,
314 const std::shared_ptr<const Utilities::MPI::CommunicationPatternBase>
315 &communication_pattern = {})
328 template <
typename MemorySpace>
333 const std::shared_ptr<const Utilities::MPI::CommunicationPatternBase>
334 &communication_pattern = {});
339 template <
typename MemorySpace>
343 const std::shared_ptr<const Utilities::MPI::CommunicationPatternBase>
344 &communication_pattern = {})
350 #ifdef DEAL_II_WITH_PETSC
363 const std::shared_ptr<const Utilities::MPI::CommunicationPatternBase>
364 &communication_pattern = {});
373 const std::shared_ptr<const Utilities::MPI::CommunicationPatternBase>
374 &communication_pattern = {})
380 #ifdef DEAL_II_WITH_TRILINOS
395 const std::shared_ptr<const Utilities::MPI::CommunicationPatternBase>
396 &communication_pattern = {});
405 const std::shared_ptr<const Utilities::MPI::CommunicationPatternBase>
406 &communication_pattern = {})
411 # ifdef DEAL_II_TRILINOS_WITH_TPETRA
420 template <
typename Dummy = Number>
421 std::enable_if_t<std::is_same_v<Dummy, Number> &&
426 const std::shared_ptr<const Utilities::MPI::CommunicationPatternBase>
427 &communication_pattern = {});
432 template <
typename Dummy = Number>
437 const std::shared_ptr<const Utilities::MPI::CommunicationPatternBase>
438 &communication_pattern = {})
456 const std::shared_ptr<const Utilities::MPI::CommunicationPatternBase>
457 &communication_pattern = {});
466 const std::shared_ptr<const Utilities::MPI::CommunicationPatternBase>
467 &communication_pattern = {})
473 #ifdef DEAL_II_WITH_CUDA
484 const std::shared_ptr<const Utilities::MPI::CommunicationPatternBase>
485 &communication_pattern = {});
495 const std::shared_ptr<const Utilities::MPI::CommunicationPatternBase>
496 &communication_pattern = {})
613 template <
typename Number2>
616 std::vector<Number2> &
values)
const;
653 template <
typename ForwardIterator,
typename OutputIterator>
656 const ForwardIterator indices_end,
657 OutputIterator values_begin)
const;
696 template <
typename Number2>
698 add(
const std::vector<size_type> &indices,
699 const std::vector<Number2> &
values);
705 template <
typename Number2>
707 add(
const std::vector<size_type> &indices,
715 template <
typename Number2>
726 const unsigned int precision = 3,
727 const bool scientific =
true)
const;
737 #ifdef DEAL_II_WITH_TRILINOS
738 # ifdef DEAL_II_TRILINOS_WITH_TPETRA
744 template <
typename Dummy = Number>
745 std::enable_if_t<std::is_same_v<Dummy, Number> &&
748 const Tpetra::Vector<Number, int, types::signed_global_dof_index>
750 const IndexSet &locally_owned_elements,
752 const MPI_Comm mpi_comm,
753 const std::shared_ptr<const Utilities::MPI::CommunicationPatternBase>
754 &communication_pattern);
764 const Epetra_MultiVector &multivector,
765 const IndexSet &locally_owned_elements,
767 const MPI_Comm mpi_comm,
768 const std::shared_ptr<const Utilities::MPI::CommunicationPatternBase>
769 &communication_pattern);
784 #ifdef DEAL_II_WITH_TRILINOS
785 # ifdef DEAL_II_TRILINOS_WITH_TPETRA
792 const MPI_Comm mpi_comm);
801 const MPI_Comm mpi_comm);
818 std::shared_ptr<Utilities::MPI::CommunicationPatternBase>
comm_pattern;
823 std::unique_ptr<Number[], void (*)(
void *)>
830 mutable std::shared_ptr<::parallel::internal::TBBPartitioner>
834 template <
typename Number2>
843 template <
typename Functor>
879 template <
typename Number>
892 template <
typename Number>
894 const ReadWriteVector<Number> &v)
903 template <
typename Number>
916 template <
typename Number>
918 const IndexSet &locally_stored_indices)
930 template <
typename Number>
939 template <
typename Number>
948 template <
typename Number>
957 template <
typename Number>
966 template <
typename Number>
975 template <
typename Number>
984 template <
typename Number>
993 template <
typename Number>
1002 template <
typename Number>
1011 template <
typename Number>
1020 template <
typename Number>
1029 template <
typename Number>
1030 template <
typename Number2>
1033 const std::vector<size_type> &indices,
1034 std::vector<Number2> &extracted_values)
const
1036 for (
size_type i = 0; i < indices.size(); ++i)
1037 extracted_values[i] =
operator()(indices[i]);
1042 template <
typename Number>
1049 for (
unsigned int i = 0; i < indices.
size(); ++i)
1050 entries[i] = (*
this)[indices[i]];
1055 template <
typename Number>
1056 template <
typename ForwardIterator,
typename OutputIterator>
1059 ForwardIterator indices_begin,
1060 const ForwardIterator indices_end,
1061 OutputIterator values_begin)
const
1063 while (indices_begin != indices_end)
1073 template <
typename Number>
1079 return values[local_index];
1084 template <
typename Number>
1090 return values[local_index];
1095 template <
typename Number>
1096 template <
typename Number2>
1099 const std::vector<Number2> &
values)
1102 add(indices.size(), indices.data(),
values.data());
1107 template <
typename Number>
1108 template <
typename Number2>
1111 const ReadWriteVector<Number2> &
values)
1119 "The given value is not finite but either infinite or Not A Number (NaN)"));
1126 template <
typename Number>
1127 template <
typename Number2>
1131 const Number2 *values_to_add)
1133 for (
size_type i = 0; i < n_indices; ++i)
1138 "The given value is not finite but either infinite or Not A Number (NaN)"));
1139 this->
operator()(indices[i]) += values_to_add[i];
1145 template <
typename Number>
1152 return static_cast<unsigned int>(
1158 template <
typename Number>
1159 template <
typename Functor>
1161 ReadWriteVector<Number> &parent,
1162 const Functor &functor)
1169 template <
typename Number>
1170 template <
typename Functor>
1177 functor(parent.values[i]);
1193 template <
typename Number>
size_type index_within_set(const size_type global_index) const
size_type n_elements() const
FunctorTemplate(ReadWriteVector< Number > &parent, const Functor &functor)
virtual void operator()(const size_type begin, const size_type end)
void import_elements(const ::Vector< Number > &vec, VectorOperation::values operation, const std::shared_ptr< const Utilities::MPI::CommunicationPatternBase > &communication_pattern={})
ReadWriteVector< Number > & operator=(const ReadWriteVector< Number2 > &in_vector)
void reinit(const ReadWriteVector< Number2 > &in_vector, const bool omit_zeroing_entries=false)
ReadWriteVector(const ReadWriteVector< Number > &in_vector)
ReadWriteVector(const size_type size)
void import_elements(const TrilinosWrappers::MPI::Vector &trilinos_vec, VectorOperation::values operation, const std::shared_ptr< const Utilities::MPI::CommunicationPatternBase > &communication_pattern={})
friend class ReadWriteVector
size_type size() const override
void add(const std::vector< size_type > &indices, const std::vector< Number2 > &values)
Number operator()(const size_type global_index) const
void apply(const Functor &func)
std::size_t memory_consumption() const
void import_elements(const CUDAWrappers::Vector< Number > &cuda_vec, VectorOperation::values operation, const std::shared_ptr< const Utilities::MPI::CommunicationPatternBase > &communication_pattern={})
void print(std::ostream &out, const unsigned int precision=3, const bool scientific=true) const
Number & operator()(const size_type global_index)
void add(const std::vector< size_type > &indices, const ReadWriteVector< Number2 > &values)
void import_elements(const EpetraWrappers::Vector &epetra_vec, VectorOperation::values operation, const std::shared_ptr< const Utilities::MPI::CommunicationPatternBase > &communication_pattern={})
ReadWriteVector< Number > & operator=(const ReadWriteVector< Number > &in_vector)
ReadWriteVector< Number > & operator=(const Number s)
unsigned int global_to_local(const types::global_dof_index global_index) const
void resize_val(const size_type new_allocated_size)
Number operator[](const size_type global_index) const
Number local_element(const size_type local_index) const
void add(const size_type n_elements, const size_type *indices, const Number2 *values)
virtual void extract_subvector_to(const ArrayView< const types::global_dof_index > &indices, ArrayView< Number > &entries) const override
void import_elements(const PETScWrappers::MPI::Vector &petsc_vec, VectorOperation::values operation, const std::shared_ptr< const Utilities::MPI::CommunicationPatternBase > &communication_pattern={})
std::enable_if_t< std::is_same_v< Dummy, Number > &&::is_tpetra_type< Number >::value > import_elements(const TpetraWrappers::Vector< Number > &tpetra_vec, VectorOperation::values operation, const std::shared_ptr< const Utilities::MPI::CommunicationPatternBase > &communication_pattern={})
TpetraWrappers::CommunicationPattern create_tpetra_comm_pattern(const IndexSet &source_index_set, const MPI_Comm mpi_comm)
Number & local_element(const size_type local_index)
void reinit(const TrilinosWrappers::MPI::Vector &trilinos_vec)
virtual void reinit(const IndexSet &locally_stored_indices, const bool omit_zeroing_entries=false)
const value_type * const_pointer
void extract_subvector_to(ForwardIterator indices_begin, const ForwardIterator indices_end, OutputIterator values_begin) const
virtual void reinit(const size_type size, const bool omit_zeroing_entries=false)
const IndexSet & get_stored_elements() const
size_type locally_owned_size() const
const_iterator end() const
IndexSet source_stored_elements
void extract_subvector_to(const std::vector< size_type > &indices, std::vector< Number2 > &values) const
const value_type * const_iterator
const_iterator begin() const
const value_type & const_reference
~ReadWriteVector() override=default
ReadWriteVector(const IndexSet &locally_stored_indices)
std::enable_if_t< std::is_same_v< Dummy, Number > &&::is_tpetra_type< Number >::value > import_elements(const Tpetra::Vector< Number, int, types::signed_global_dof_index > &tpetra_vector, const IndexSet &locally_owned_elements, VectorOperation::values operation, const MPI_Comm mpi_comm, const std::shared_ptr< const Utilities::MPI::CommunicationPatternBase > &communication_pattern)
EpetraWrappers::CommunicationPattern create_epetra_comm_pattern(const IndexSet &source_index_set, const MPI_Comm mpi_comm)
std::unique_ptr< Number[], void(*)(void *)> values
types::global_dof_index size_type
void import_elements(const distributed::Vector< Number, MemorySpace > &vec, VectorOperation::values operation, const std::shared_ptr< const Utilities::MPI::CommunicationPatternBase > &communication_pattern={})
void swap(ReadWriteVector< Number > &v)
std::shared_ptr< Utilities::MPI::CommunicationPatternBase > comm_pattern
std::shared_ptr<::parallel::internal::TBBPartitioner > thread_loop_partitioner
typename numbers::NumberTraits< Number >::real_type real_type
void import_elements(const Epetra_MultiVector &multivector, const IndexSet &locally_owned_elements, VectorOperation::values operation, const MPI_Comm mpi_comm, const std::shared_ptr< const Utilities::MPI::CommunicationPatternBase > &communication_pattern)
Number & operator[](const size_type global_index)
types::global_dof_index size_type
void swap(LinearAlgebra::ReadWriteVector< Number > &u, LinearAlgebra::ReadWriteVector< Number > &v)
#define DEAL_II_DEPRECATED
#define DEAL_II_NAMESPACE_OPEN
#define DEAL_II_NAMESPACE_CLOSE
#define Assert(cond, exc)
#define AssertDimension(dim1, dim2)
#define AssertIndexRange(index, range)
static ::ExceptionBase & ExcMessage(std::string arg1)
TrilinosWrappers::types::int_type global_index(const Epetra_BlockMap &map, const ::types::global_dof_index i)
bool is_finite(const double x)
unsigned int global_dof_index