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 83 template <
typename Number>
100 static constexpr
unsigned int communication_block_size = 20;
161 template <
typename OtherNumber>
168 BlockVector(
const std::vector<size_type> &block_sizes);
174 BlockVector(
const std::vector<IndexSet> &local_ranges,
175 const std::vector<IndexSet> &ghost_indices,
181 BlockVector(
const std::vector<IndexSet> &local_ranges,
212 template <
class Number2>
222 #ifdef DEAL_II_WITH_PETSC 234 #ifdef DEAL_II_WITH_TRILINOS 263 const bool omit_zeroing_entries =
false);
285 reinit(
const std::vector<size_type> &
N,
286 const bool omit_zeroing_entries =
false);
302 template <
typename Number2>
305 const bool omit_zeroing_entries =
false);
342 update_ghost_values()
const;
353 zero_out_ghosts()
const;
359 has_ghost_elements()
const;
365 template <
typename OtherNumber>
367 add(
const std::vector<size_type> & indices,
368 const ::Vector<OtherNumber> &
values);
383 all_zero()
const override;
389 mean_value()
const override;
430 const bool omit_zeroing_entries =
false)
override;
436 operator*=(
const Number factor)
override;
442 operator/=(
const Number factor)
override;
468 std::shared_ptr<const CommunicationPatternBase> communication_pattern =
469 std::shared_ptr<const CommunicationPatternBase>())
override;
493 template <
typename FullMatrixType>
495 multivector_inner_product(FullMatrixType &
matrix,
515 template <
typename FullMatrixType>
517 multivector_inner_product_with_metric(
const FullMatrixType &
matrix,
530 template <
typename FullMatrixType>
533 const FullMatrixType &
matrix,
534 const Number s = Number(0.),
535 const Number
b = Number(1.))
const;
541 add(
const Number a)
override;
563 add(
const std::vector<size_type> &indices,
564 const std::vector<Number> &
values);
646 size()
const override;
660 locally_owned_elements()
const override;
666 print(std::ostream & out,
667 const unsigned int precision = 3,
668 const bool scientific =
true,
669 const bool across =
true)
const override;
711 template <
typename Number>
724 template <
typename Number>
733 # include <deal.II/lac/la_parallel_block_vector.templates.h>
void reinit(MatrixBlock< MatrixType > &v, const BlockSparsityPattern &p)
typename BaseClass::BlockType BlockType
Contents is actually a matrix.
types::global_dof_index size_type
std::vector< value_type > l2_norm(const typename ::Triangulation< dim, spacedim >::cell_iterator &parent, const value_type parent_value)
typename BaseClass::reference reference
typename BaseClass::const_reference const_reference
__global__ void add_and_dot(Number *res, Number *v1, const Number *v2, const Number *v3, const Number a, const size_type N)
typename BaseClass::const_iterator const_iterator
typename BaseClass::value_type value_type
void swap(BlockVector &u, BlockVector &v)
std::string compress(const std::string &input)
BlockVector< double > BlockVector
Number linfty_norm(const Tensor< 2, dim, Number > &t)
#define DeclException0(Exception0)
#define DEAL_II_NAMESPACE_CLOSE
void swap(BlockVector< Number > &v)
__global__ void equ(Number *val, const Number a, const Number *V_val, const size_type N)
__global__ void sadd(const Number s, Number *val, const Number a, const Number *V_val, const size_type N)
SymmetricTensor< 2, dim, Number > b(const Tensor< 2, dim, Number > &F)
Vector< Number > BlockType
typename BaseClass::real_type real_type
typename BaseClass::pointer pointer
typename BlockType::value_type value_type
#define DEAL_II_NAMESPACE_OPEN
typename BlockType::real_type real_type
Number l1_norm(const Tensor< 2, dim, Number > &t)
typename BaseClass::size_type size_type
typename BlockType::const_reference const_reference
typename BaseClass::iterator iterator
std::enable_if< std::is_floating_point< T >::value &&std::is_floating_point< U >::value, typename ProductType< std::complex< T >, std::complex< U > >::type >::type operator*(const std::complex< T > &left, const std::complex< U > &right)
typename BaseClass::const_pointer const_pointer
typename BlockType::reference reference
std::enable_if< std::is_fundamental< T >::value, std::size_t >::type memory_consumption(const T &t)
const value_type * const_pointer