15#ifndef dealii_block_vector_h
16#define dealii_block_vector_h
36# ifdef DEAL_II_WITH_TRILINOS
69template <
typename Number>
128 template <typename OtherNumber>
131#ifdef DEAL_II_WITH_TRILINOS
160 template <
typename InputIterator>
162 const InputIterator
first,
163 const InputIterator
end);
217 template <
class Number2>
227#ifdef DEAL_II_WITH_TRILINOS
251 const bool omit_zeroing_entries =
false);
270 reinit(
const std::vector<size_type> &block_sizes,
271 const bool omit_zeroing_entries =
false);
284 const bool omit_zeroing_entries =
false);
299 template <
typename Number2>
302 const bool omit_zeroing_entries =
false);
308 template <
class BlockVector2>
331 const unsigned int precision = 3,
332 const bool scientific =
true,
333 const bool across =
true)
const;
376template <
typename Number>
377template <
typename InputIterator>
379 const InputIterator
first,
380 const InputIterator end)
386 reinit(block_sizes,
true);
387 InputIterator start =
first;
388 for (size_type b = 0; b < block_sizes.size(); ++b)
390 InputIterator end = start;
391 std::advance(end,
static_cast<signed int>(block_sizes[b]));
392 std::copy(start, end, this->block(b).begin());
395 Assert(start == end, ExcIteratorRangeDoesNotMatchVectorSize());
400template <
typename Number>
406 BaseClass::operator=(s);
412template <
typename Number>
417 BaseClass::operator=(v);
423template <
typename Number>
427 BaseClass::operator=(v);
433template <
typename Number>
434template <
typename Number2>
439 BaseClass::operator=(v);
443template <
typename Number>
447 for (size_type i = 0; i < this->
n_blocks(); ++i)
448 this->components[i].
compress(operation);
453template <
typename Number>
462template <
typename Number>
463template <
class BlockVector2>
480template <
typename Number>
490 namespace LinearOperatorImplementation
499 template <
typename number>
503 template <
typename Matrix>
507 bool omit_zeroing_entries)
509 v.
reinit(matrix.get_row_indices(), omit_zeroing_entries);
512 template <
typename Matrix>
516 bool omit_zeroing_entries)
518 v.
reinit(matrix.get_column_indices(), omit_zeroing_entries);
529template <
typename Number>
::internal::BlockVectorIterators::Iterator< BlockVectorBase, false > iterator
unsigned int n_blocks() const
const value_type * const_pointer
typename BlockType::const_reference const_reference
types::global_dof_index size_type
typename BlockType::value_type value_type
::internal::BlockVectorIterators::Iterator< BlockVectorBase, true > const_iterator
typename BlockType::reference reference
BlockIndices block_indices
typename BlockType::real_type real_type
Vector< Number > BlockType
BlockVector(const std::vector< size_type > &block_sizes)
typename BaseClass::const_reference const_reference
void reinit(const BlockIndices &block_indices, const bool omit_zeroing_entries=false)
void swap(BlockVector< Number > &v)
typename BaseClass::reference reference
void print(std::ostream &out, const unsigned int precision=3, const bool scientific=true, const bool across=true) const
void compress(VectorOperation::values operation=VectorOperation::unknown)
bool has_ghost_elements() const
BlockVector(const TrilinosWrappers::MPI::BlockVector &v)
BlockVector< Number > & operator=(const TrilinosWrappers::MPI::BlockVector &V)
BlockVector< Number > & operator=(const BlockVector< Number2 > &V)
BlockVector(const BlockIndices &block_indices)
BlockVector< Number > & operator=(BlockVector< Number > &&)=default
typename BaseClass::value_type value_type
typename BaseClass::const_pointer const_pointer
BlockVector(const BlockVector< Number > &V)
void block_write(std::ostream &out) const
typename BaseClass::pointer pointer
~BlockVector() override=default
typename BaseClass::BlockType BlockType
BlockVector(const unsigned int n_blocks=0, const size_type block_size=0)
void scale(const BlockVector2 &v)
BlockVector(BlockVector< Number > &&) noexcept=default
void reinit(const std::vector< size_type > &block_sizes, const bool omit_zeroing_entries=false)
typename BaseClass::iterator iterator
void block_read(std::istream &in)
BlockVector & operator=(const value_type s)
typename BaseClass::size_type size_type
BlockVector< Number > & operator=(const Vector< Number > &V)
BlockVector< Number > & operator=(const BlockVector< Number > &v)
BlockVector(const std::vector< size_type > &block_sizes, const InputIterator first, const InputIterator end)
void reinit(const unsigned int n_blocks, const size_type block_size=0, const bool omit_zeroing_entries=false)
void swap(BlockVector< Number > &u, BlockVector< Number > &v)
typename BaseClass::const_iterator const_iterator
void reinit(const BlockVector< Number2 > &V, const bool omit_zeroing_entries=false)
typename BaseClass::real_type real_type
static void reinit_range_vector(const Matrix &matrix, BlockVector< number > &v, bool omit_zeroing_entries)
static void reinit_domain_vector(const Matrix &matrix, BlockVector< number > &v, bool omit_zeroing_entries)
#define DEAL_II_NAMESPACE_OPEN
#define DEAL_II_NAMESPACE_CLOSE
#define DeclException0(Exception0)
#define Assert(cond, exc)
#define AssertIsFinite(number)
static ::ExceptionBase & ExcIteratorRangeDoesNotMatchVectorSize()
std::enable_if_t< IsBlockVector< VectorType >::value, unsigned int > n_blocks(const VectorType &vector)
std::string compress(const std::string &input)
void reinit(MatrixBlock< MatrixType > &v, const BlockSparsityPattern &p)