16 #ifndef dealii_la_parallel_block_vector_h
17 #define dealii_la_parallel_block_vector_h
38 # ifdef DEAL_II_WITH_PETSC
48 # ifdef DEAL_II_WITH_TRILINOS
84 template <
typename Number>
155 template <
typename OtherNumber>
169 const std::vector<IndexSet> &ghost_indices,
170 const MPI_Comm communicator);
176 const MPI_Comm communicator);
190 const std::vector<std::shared_ptr<const Utilities::MPI::Partitioner>>
192 const MPI_Comm &comm_sm = MPI_COMM_SELF);
222 template <
class Number2>
232 #ifdef DEAL_II_WITH_PETSC
244 #ifdef DEAL_II_WITH_TRILINOS
273 const bool omit_zeroing_entries =
false);
295 reinit(
const std::vector<size_type> &block_sizes,
296 const bool omit_zeroing_entries =
false);
312 template <
typename Number2>
315 const bool omit_zeroing_entries =
false);
335 reinit(
const std::vector<IndexSet> &local_ranges,
336 const std::vector<IndexSet> &ghost_indices,
337 const MPI_Comm communicator);
343 reinit(
const std::vector<IndexSet> &local_ranges,
344 const MPI_Comm communicator);
361 const std::vector<std::shared_ptr<const Utilities::MPI::Partitioner>>
363 const MPI_Comm &comm_sm = MPI_COMM_SELF);
373 const std::vector<std::shared_ptr<const Utilities::MPI::Partitioner>>
375 const bool make_ghosted,
376 const MPI_Comm &comm_sm = MPI_COMM_SELF);
452 template <
typename Number2>
460 template <
typename OtherNumber>
462 add(
const std::vector<size_type> &indices,
463 const ::Vector<OtherNumber> &
values);
555 const std::shared_ptr<const Utilities::MPI::CommunicationPatternBase>
556 &communication_pattern = {});
564 std::shared_ptr<const Utilities::MPI::CommunicationPatternBase>
565 communication_pattern = {})
592 template <
typename FullMatrixType>
614 template <
typename FullMatrixType>
629 template <
typename FullMatrixType>
632 const FullMatrixType &
matrix,
633 const Number s = Number(0.),
634 const Number
b = Number(1.))
const;
662 add(
const std::vector<size_type> &indices,
663 const std::vector<Number> &
values);
764 const unsigned int precision = 3,
765 const bool scientific =
true,
766 const bool across =
true)
const;
808 template <
typename Number>
821 template <
typename Number>
830 # include <deal.II/lac/la_parallel_block_vector.templates.h>
::internal::BlockVectorIterators::Iterator< BlockVectorBase, false > iterator
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
typename BlockType::real_type real_type
Vector< Number > BlockType
void swap(LinearAlgebra::distributed::BlockVector< Number > &u, LinearAlgebra::distributed::BlockVector< Number > &v)
BlockVector & operator=(const BlockVector &V)
BlockVector & operator=(const Vector< Number > &V)
typename BaseClass::value_type value_type
bool has_ghost_elements() const
BlockVector(const std::vector< std::shared_ptr< const Utilities::MPI::Partitioner >> &partitioners, const MPI_Comm &comm_sm=MPI_COMM_SELF)
void swap(BlockVector< Number > &v)
Number operator*(const BlockVector< Number > &V) const
BlockVector(const size_type num_blocks=0, const size_type block_size=0)
void add(const std::vector< size_type > &indices, const std::vector< Number > &values)
void reinit(const std::vector< std::shared_ptr< const Utilities::MPI::Partitioner >> &partitioners, const MPI_Comm &comm_sm=MPI_COMM_SELF)
void print(std::ostream &out, const unsigned int precision=3, const bool scientific=true, const bool across=true) const
BlockVector< Number > & operator=(const TrilinosWrappers::MPI::BlockVector &trilinos_vec)
void reinit(const std::vector< IndexSet > &local_ranges, const MPI_Comm communicator)
void sadd(const Number s, const Number a, const BlockVector< Number > &V)
static constexpr unsigned int communication_block_size
BlockVector & operator=(const value_type s)
void add(const Number a, const BlockVector< Number > &V, const Number b, const BlockVector< Number > &W)
real_type norm_sqr() const
real_type l2_norm() const
BlockVector & operator=(const BlockVector< Number2 > &V)
BlockVector(const std::vector< size_type > &block_sizes)
std::size_t memory_consumption() const
BlockVector(const BlockVector< OtherNumber > &v)
typename BaseClass::iterator iterator
BlockVector(const BlockVector< Number > &V)
void reinit(const size_type num_blocks, const size_type block_size=0, const bool omit_zeroing_entries=false)
void set_ghost_state(const bool ghosted) const
void add(const Number a, const BlockVector< Number > &V)
void add(const std::vector< size_type > &indices, const ::Vector< OtherNumber > &values)
BlockVector< Number > & operator+=(const BlockVector< Number > &V)
Number mean_value() const
BlockVector(const std::vector< IndexSet > &local_ranges, const MPI_Comm communicator)
void scale(const BlockVector< Number > &scaling_factors)
typename BaseClass::BlockType BlockType
void reinit(const std::vector< IndexSet > &local_ranges, const std::vector< IndexSet > &ghost_indices, const MPI_Comm communicator)
typename BaseClass::size_type size_type
Number multivector_inner_product_with_metric(const FullMatrixType &matrix, const BlockVector< Number > &V, const bool symmetric=false) const
void zero_out_ghost_values() const
BlockVector< Number > & operator-=(const BlockVector< Number > &V)
BlockVector< Number > & operator=(const PETScWrappers::MPI::BlockVector &petsc_vec)
void reinit(const std::vector< std::shared_ptr< const Utilities::MPI::Partitioner >> &partitioners, const bool make_ghosted, const MPI_Comm &comm_sm=MPI_COMM_SELF)
typename BaseClass::const_reference const_reference
real_type lp_norm(const real_type p) const
typename BaseClass::real_type real_type
void reinit(const BlockVector< Number2 > &V, const bool omit_zeroing_entries=false)
Number add_and_dot(const Number a, const BlockVector< Number > &V, const BlockVector< Number > &W)
void equ(const Number a, const BlockVector< Number > &V)
void sadd(const Number s, const BlockVector< Number > &V)
void copy_locally_owned_data_from(const BlockVector< Number2 > &src)
BlockVector(const std::vector< IndexSet > &local_ranges, const std::vector< IndexSet > &ghost_indices, const MPI_Comm communicator)
BlockVector< Number > & operator/=(const Number factor)
real_type linfty_norm() const
typename BaseClass::const_pointer const_pointer
typename BaseClass::reference reference
typename BaseClass::const_iterator const_iterator
void multivector_inner_product(FullMatrixType &matrix, const BlockVector< Number > &V, const bool symmetric=false) const
real_type l1_norm() const
void update_ghost_values() const
void mmult(BlockVector< Number > &V, const FullMatrixType &matrix, const Number s=Number(0.), const Number b=Number(1.)) const
BlockVector< Number > & operator*=(const Number factor)
void reinit(const std::vector< size_type > &block_sizes, const bool omit_zeroing_entries=false)
void compress(VectorOperation::values operation)
void import_elements(const LinearAlgebra::ReadWriteVector< Number > &V, const VectorOperation::values operation, const std::shared_ptr< const Utilities::MPI::CommunicationPatternBase > &communication_pattern={})
virtual size_type size() const override
typename BaseClass::pointer pointer
::IndexSet locally_owned_elements() const
#define DEAL_II_DEPRECATED
#define DEAL_II_NAMESPACE_OPEN
#define DEAL_II_NAMESPACE_CLOSE
static ::ExceptionBase & ExcIteratorRangeDoesNotMatchVectorSize()
#define DeclException0(Exception0)
static ::ExceptionBase & ExcVectorTypeNotCompatible()
@ matrix
Contents is actually a matrix.
@ symmetric
Matrix is symmetric.
BlockVector< double > BlockVector
SymmetricTensor< 2, dim, Number > b(const Tensor< 2, dim, Number > &F)