17 #include <deal.II/lac/vector.h> 18 #include <deal.II/lac/la_vector.h> 19 #include <deal.II/lac/la_parallel_vector.h> 20 #include <deal.II/lac/la_parallel_block_vector.h> 21 #include <deal.II/lac/petsc_vector.h> 22 #include <deal.II/lac/petsc_block_vector.h> 23 #include <deal.II/lac/trilinos_vector.h> 24 #include <deal.II/lac/trilinos_block_vector.h> 25 #include <deal.II/lac/sparse_matrix.h> 26 #include <deal.II/lac/block_sparse_matrix.h> 27 #include <deal.II/multigrid/mg_transfer.h> 28 #include <deal.II/multigrid/mg_transfer_block.h> 29 #include <deal.II/multigrid/mg_transfer_component.h> 30 #include <deal.II/multigrid/mg_smoother.h> 31 #include <deal.II/multigrid/mg_transfer_block.templates.h> 32 #include <deal.II/multigrid/mg_transfer_component.templates.h> 33 #include <deal.II/multigrid/multigrid.templates.h> 35 DEAL_II_NAMESPACE_OPEN
61 template <
typename number>
64 memory(0, typeid(*this).name())
68 template <
typename number>
75 template <
typename number>
85 template <
typename number>
87 const unsigned int to_level,
101 for (
unsigned int b=0; b<this->
mg_block.size(); ++b)
110 template <
typename number>
112 const unsigned int from_level,
123 for (
unsigned int b=0; b<this->
mg_block.size(); ++b)
133 *aux, src.
block(this->mg_block[b]));
152 std::size_t result =
sizeof(*this);
156 -
sizeof(mg_target_component);
160 -
sizeof(component_start);
162 -
sizeof(mg_component_start);
164 -
sizeof(prolongation_sparsities);
178 std::size_t result =
sizeof(*this);
179 result +=
sizeof(
unsigned int) *
sizes.size();
189 -
sizeof(prolongation_sparsities);
201 template<
typename number>
204 selected_component (0),
205 mg_selected_component (0)
209 template<
typename number>
217 template <
typename number>
222 template <
typename number>
224 const unsigned int to_level,
238 template <
typename number>
240 const unsigned int from_level,
249 .Tvmult_add (dst, src);
255 template <
typename number>
263 template <
typename number>
272 template <
typename number>
282 template <
typename number>
288 template <
typename number>
303 template <
typename number>
314 .Tvmult_add (dst, src);
321 #include "multigrid.inst" 331 DEAL_II_NAMESPACE_CLOSE
virtual VectorType * alloc()=0
std::size_t memory_consumption() const
std::vector< std_cxx11::shared_ptr< BlockSparseMatrix< double > > > prolongation_matrices
static::ExceptionBase & ExcNotInitialized()
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_cxx11::shared_ptr< BlockSparseMatrix< double > > > prolongation_matrices
#define Assert(cond, exc)
virtual void restrict_and_add(const unsigned int from_level, Vector< number > &dst, const Vector< number > &src) const
virtual void prolongate(const unsigned int to_level, Vector< number > &dst, const Vector< number > &src) const
static::ExceptionBase & ExcDimensionMismatch(std::size_t arg1, std::size_t arg2)
virtual ~MGTransferBlock()
std::size_t memory_consumption() const
void initialize(const std::vector< number > &factors, VectorMemory< Vector< number > > &memory)
virtual void prolongate(const unsigned int to_level, Vector< number > &dst, const Vector< number > &src) const
virtual void prolongate(const unsigned int to_level, BlockVector< number > &dst, const BlockVector< number > &src) const
std::vector< unsigned int > mg_block
unsigned int mg_selected_component
std_cxx11::enable_if< std_cxx11::is_fundamental< T >::value, std::size_t >::type memory_consumption(const T &t)
std::vector< unsigned int > mg_target_component
std::vector< std::vector< types::global_dof_index > > mg_block_start
SmartPointer< VectorMemory< Vector< number > >, MGTransferBlock< number > > memory
virtual void restrict_and_add(const unsigned int from_level, BlockVector< number > &dst, const BlockVector< number > &src) const
SmartPointer< const ConstraintMatrix > constraints
virtual void reinit(const size_type N, const bool omit_zeroing_entries=false)
virtual ~MGTransferSelect()
std::vector< number > factors
virtual ~MGTransferBlockSelect()
virtual void restrict_and_add(const unsigned int from_level, Vector< number > &dst, const Vector< number > &src) const
BlockType & block(const unsigned int i)
std::vector< bool > selected
unsigned int selected_block
unsigned int selected_component