16 #ifndef dealii_vector_h
17 #define dealii_vector_h
32 #include <boost/serialization/split_member.hpp>
35 #include <initializer_list>
44 # ifdef DEAL_II_WITH_PETSC
51 # ifdef DEAL_II_WITH_TRILINOS
61 template <
typename number>
108 template <
typename Number>
122 "The Vector class only supports basic numeric types. In particular, it "
123 "does not support automatically differentiated numbers.");
183 template <
typename OtherNumber>
201 template <
typename OtherNumber>
202 explicit Vector(
const std::initializer_list<OtherNumber> &v);
204 #ifdef DEAL_II_WITH_PETSC
219 #ifdef DEAL_II_WITH_TRILINOS
252 template <
typename InputIterator>
340 template <
typename Number2>
396 template <
typename Number2>
406 #ifdef DEAL_II_WITH_PETSC
423 #ifdef DEAL_II_WITH_TRILINOS
448 template <
typename Number2>
457 template <
typename Number2>
482 template <
typename Number2>
661 template <
typename OtherNumber>
664 std::vector<OtherNumber> &
values)
const;
700 template <
typename ForwardIterator,
typename OutputIterator>
703 const ForwardIterator indices_end,
704 OutputIterator values_begin)
const;
733 template <
typename OtherNumber>
735 add(
const std::vector<size_type> &indices,
736 const std::vector<OtherNumber> &
values);
742 template <
typename OtherNumber>
751 template <
typename OtherNumber>
755 const OtherNumber *
values);
832 template <
typename Number2>
847 template <
typename Number2>
872 const unsigned int precision = 3,
873 const bool scientific =
true,
874 const bool across =
true)
const;
903 template <
class Archive>
905 save(Archive &ar,
const unsigned int version)
const;
912 template <
class Archive>
914 load(Archive &ar,
const unsigned int version);
922 template <
class Archive>
924 serialize(Archive &archive,
const unsigned int version);
928 BOOST_SERIALIZATION_SPLIT_MEMBER()
1045 const bool omit_zeroing_entries,
1046 const bool reset_partitioner);
1052 mutable std::shared_ptr<parallel::internal::TBBPartitioner>
1056 template <
typename Number2>
1075 template <
typename Number>
1086 template <
typename Number>
1087 template <
typename OtherNumber>
1094 template <
typename Number>
1095 template <
typename InputIterator>
1106 template <
typename Number>
1117 template <
typename Number>
1126 template <
typename Number>
1135 template <
typename Number>
1144 template <
typename Number>
1153 template <
typename Number>
1162 template <
typename Number>
1171 template <
typename Number>
1180 template <
typename Number>
1189 template <
typename Number>
1198 template <
typename Number>
1208 template <
typename Number>
1218 template <
typename Number>
1227 template <
typename Number>
1236 template <
typename Number>
1237 template <
typename OtherNumber>
1240 std::vector<OtherNumber> &
values)
const
1242 for (
size_type i = 0; i < indices.size(); ++i)
1243 values[i] =
operator()(indices[i]);
1248 template <
typename Number>
1249 template <
typename ForwardIterator,
typename OutputIterator>
1252 const ForwardIterator indices_end,
1253 OutputIterator values_begin)
const
1255 while (indices_begin != indices_end)
1265 template <
typename Number>
1278 template <
typename Number>
1279 template <
typename OtherNumber>
1282 const std::vector<OtherNumber> &
values)
1286 add(indices.size(), indices.data(),
values.data());
1291 template <
typename Number>
1292 template <
typename OtherNumber>
1299 add(indices.size(), indices.data(),
values.values.begin());
1304 template <
typename Number>
1305 template <
typename OtherNumber>
1309 const OtherNumber *
values)
1311 for (
size_type i = 0; i < n_indices; ++i)
1317 "The given value is not finite but either infinite or Not A Number (NaN)"));
1319 this->values[indices[i]] +=
values[i];
1325 template <
typename Number>
1326 template <
typename Number2>
1330 return !(*
this == v);
1335 template <
typename Number>
1342 template <
typename Number>
1351 template <
typename Number>
1358 template <
typename Number>
1365 template <
typename Number>
1366 template <
typename Number2>
1369 const bool omit_zeroing_entries)
1373 if (!omit_zeroing_entries ||
size() != v.
size())
1385 template <
typename Number>
1395 template <
typename Number>
1396 template <
class Archive>
1407 template <
typename Number>
1408 template <
class Archive>
1434 template <
typename Number>
1449 template <
typename number>
1450 inline std::ostream &
1458 out << v(v.
size() - 1);
1473 template <
typename Number>
OutputOperator< VectorType > & operator<<(OutputOperator< VectorType > &out, unsigned int step)
types::global_dof_index size_type
typename numbers::NumberTraits< Number >::real_type real_type
void extract_subvector_to(ForwardIterator indices_begin, const ForwardIterator indices_end, OutputIterator values_begin) const
bool operator!=(const Vector< Number2 > &v) const
const value_type * const_pointer
void scale(const Vector< Number2 > &scaling_factors)
Vector< Number > & operator+=(const Vector< Number > &V)
bool has_ghost_elements() const
void add(const std::vector< size_type > &indices, const std::vector< OtherNumber > &values)
Number & operator()(const size_type i)
void add(const Number a, const Vector< Number > &V, const Number b, const Vector< Number > &W)
Vector(const TrilinosWrappers::MPI::Vector &v)
bool in_local_range(const size_type global_index) const
Number operator*(const Vector< Number2 > &V) const
const_pointer data() const
void do_reinit(const size_type new_size, const bool omit_zeroing_entries, const bool reset_partitioner)
void sadd(const Number s, const Number a, const Vector< Number > &V)
Vector(const size_type n)
const value_type * const_iterator
void block_write(std::ostream &out) const
Number mean_value() const
Number add_and_dot(const Number a, const Vector< Number > &V, const Vector< Number > &W)
void load(Archive &ar, const unsigned int version)
const value_type & const_reference
Vector< Number > & operator/=(const Number factor)
void reinit(const Vector< Number2 > &V, const bool omit_zeroing_entries=false)
Vector< Number > & operator=(const TrilinosWrappers::MPI::Vector &v)
void print(std::ostream &out, const unsigned int precision=3, const bool scientific=true, const bool across=true) const
Number & operator[](const size_type i)
Vector< Number > & operator*=(const Number factor)
Vector(const Vector< Number > &v)
std::shared_ptr< parallel::internal::TBBPartitioner > thread_loop_partitioner
void block_read(std::istream &in)
types::global_dof_index size_type
Vector< Number > & operator=(const BlockVector< Number > &v)
Number operator[](const size_type i) const
Vector< Number > & operator=(const Vector< Number > &v)
virtual size_type size() const override
void equ(const Number a, const Vector< Number > &u)
void zero_out_ghost_values() const
Vector< Number > & operator-=(const Vector< Number > &V)
void serialize(Archive &archive, const unsigned int version)
Vector(const PETScWrappers::VectorBase &v)
real_type lp_norm(const real_type p) const
const_iterator begin() const
void sadd(const Number s, const Vector< Number > &V)
Vector< Number > & operator=(const PETScWrappers::VectorBase &v)
bool operator==(const Vector< Number2 > &v) const
real_type l2_norm() const
virtual void swap(Vector< Number > &v)
void grow_or_shrink(const size_type N)
virtual ~Vector() override=default
Vector(const Vector< OtherNumber > &v)
void save(Archive &ar, const unsigned int version) const
Vector< Number > & operator=(const Number s)
real_type linfty_norm() const
Vector< Number > & operator=(const Vector< Number2 > &v)
void scale(const Vector< Number > &scaling_factors)
virtual void extract_subvector_to(const ArrayView< const types::global_dof_index > &indices, ArrayView< Number > &elements) const override
void maybe_reset_thread_partitioner()
IndexSet locally_owned_elements() const
real_type norm_sqr() const
void add(const Number a, const Vector< Number > &V)
AlignedVector< Number > values
void equ(const Number a, const Vector< Number2 > &u)
void apply_givens_rotation(const std::array< Number, 3 > &csr, const size_type i, const size_type k)
const_iterator end() const
Vector(Vector< Number > &&v) noexcept=default
virtual void reinit(const size_type N, const bool omit_zeroing_entries=false)
void add(const size_type n_elements, const size_type *indices, const OtherNumber *values)
size_type locally_owned_size() const
Number operator()(const size_type i) const
void compress(VectorOperation::values operation=VectorOperation::unknown) const
bool is_non_negative() const
void extract_subvector_to(const std::vector< size_type > &indices, std::vector< OtherNumber > &values) const
void update_ghost_values() const
Vector< Number > & operator=(Vector< Number > &&v) noexcept=default
std::size_t memory_consumption() const
real_type l1_norm() const
void add(const std::vector< size_type > &indices, const Vector< OtherNumber > &values)
Vector(const InputIterator first, const InputIterator last)
Vector(const std::initializer_list< OtherNumber > &v)
#define DEAL_II_NAMESPACE_OPEN
#define DEAL_II_NAMESPACE_CLOSE
static ::ExceptionBase & ExcZero()
static ::ExceptionBase & ExcDimensionMismatch(std::size_t arg1, std::size_t arg2)
#define Assert(cond, exc)
#define AssertIsFinite(number)
static ::ExceptionBase & ExcEmptyObject()
#define AssertIndexRange(index, range)
static ::ExceptionBase & ExcIO()
static ::ExceptionBase & ExcMessage(std::string arg1)
#define AssertThrow(cond, exc)
void swap(Vector< Number > &u, Vector< Number > &v)
types::global_dof_index size_type
SymmetricTensor< 2, dim, Number > b(const Tensor< 2, dim, Number > &F)
VectorType::value_type * begin(VectorType &V)
VectorType::value_type * end(VectorType &V)
TrilinosWrappers::types::int_type global_index(const Epetra_BlockMap &map, const ::types::global_dof_index i)
void copy(const T *begin, const T *end, U *dest)
void reinit(MatrixBlock< MatrixType > &v, const BlockSparsityPattern &p)
bool is_finite(const double x)
unsigned int global_dof_index