17 #include <deal.II/lac/block_sparse_matrix.h> 18 #include <deal.II/lac/la_parallel_block_vector.h> 19 #include <deal.II/lac/la_parallel_vector.h> 20 #include <deal.II/lac/la_vector.h> 21 #include <deal.II/lac/petsc_block_vector.h> 22 #include <deal.II/lac/petsc_vector.h> 23 #include <deal.II/lac/sparse_matrix.h> 24 #include <deal.II/lac/trilinos_parallel_block_vector.h> 25 #include <deal.II/lac/trilinos_vector.h> 26 #include <deal.II/lac/vector.h> 28 #include <deal.II/multigrid/mg_smoother.h> 29 #include <deal.II/multigrid/mg_transfer.h> 30 #include <deal.II/multigrid/mg_transfer_block.h> 31 #include <deal.II/multigrid/mg_transfer_block.templates.h> 32 #include <deal.II/multigrid/mg_transfer_component.h> 33 #include <deal.II/multigrid/mg_transfer_component.templates.h> 34 #include <deal.II/multigrid/multigrid.templates.h> 36 DEAL_II_NAMESPACE_OPEN
60 template <
typename number>
62 : memory(nullptr, typeid(*this).name())
66 template <
typename number>
74 template <
typename number>
84 template <
typename number>
100 for (
unsigned int b = 0; b < this->
mg_block.size(); ++b)
109 template <
typename number>
122 for (
unsigned int b = 0; b < this->
mg_block.size(); ++b)
132 *aux, src.
block(this->mg_block[b]));
151 std::size_t result =
sizeof(*this);
155 sizeof(mg_target_component);
158 sizeof(component_start);
160 sizeof(mg_component_start);
162 sizeof(prolongation_sparsities);
176 std::size_t result =
sizeof(*this);
177 result +=
sizeof(
unsigned int) *
sizes.size();
185 sizeof(prolongation_sparsities);
197 template <
typename number>
199 : selected_component(0)
200 , mg_selected_component(0)
204 template <
typename number>
213 template <
typename number>
230 template <
typename number>
242 .Tvmult_add(dst, src);
248 template <
typename number>
255 template <
typename number>
264 template <
typename number>
274 template <
typename number>
289 template <
typename number>
300 .Tvmult_add(dst, src);
307 #include "multigrid.inst" 317 DEAL_II_NAMESPACE_CLOSE
std::vector< std::shared_ptr< BlockSparseMatrix< double > > > prolongation_matrices
virtual VectorType * alloc()=0
virtual ~MGTransferBlock() override
std::size_t memory_consumption() const
static::ExceptionBase & ExcNotInitialized()
virtual void prolongate(const unsigned int to_level, Vector< number > &dst, const Vector< number > &src) const override
unsigned int n_blocks() const
std::vector< types::global_dof_index > block_start
static::ExceptionBase & ExcIndexRange(int arg1, int arg2, int arg3)
SmartPointer< const MGConstrainedDoFs, MGTransferBlockBase > mg_constrained_dofs
virtual void free(const VectorType *const)=0
std::vector< std::vector< types::global_dof_index > > sizes
std::vector< std::shared_ptr< BlockSparseMatrix< double > > > prolongation_matrices
#define Assert(cond, exc)
static::ExceptionBase & ExcDimensionMismatch(std::size_t arg1, std::size_t arg2)
virtual void restrict_and_add(const unsigned int from_level, Vector< number > &dst, const Vector< number > &src) const override
virtual void restrict_and_add(const unsigned int from_level, BlockVector< number > &dst, const BlockVector< number > &src) const override
virtual void restrict_and_add(const unsigned int from_level, Vector< number > &dst, const Vector< number > &src) const override
std::size_t memory_consumption() const
virtual void prolongate(const unsigned int to_level, Vector< number > &dst, const Vector< number > &src) const override
std::vector< unsigned int > mg_block
unsigned int mg_selected_component
std::vector< unsigned int > mg_target_component
std::vector< std::vector< types::global_dof_index > > mg_block_start
SmartPointer< const AffineConstraints< double > > constraints
virtual void reinit(const size_type N, const bool omit_zeroing_entries=false)
SmartPointer< VectorMemory< Vector< number > >, MGTransferBlock< number > > memory
std::vector< number > factors
BlockType & block(const unsigned int i)
virtual void prolongate(const unsigned int to_level, BlockVector< number > &dst, const BlockVector< number > &src) const override
void initialize(const std::vector< number > &factors, VectorMemory< Vector< number >> &memory)
std::vector< bool > selected
unsigned int selected_block
std::enable_if< std::is_fundamental< T >::value, std::size_t >::type memory_consumption(const T &t)
unsigned int selected_component