Reference documentation for deal.II version Git 1f3f4df3cd 2020-02-17 16:54:31 -0500
\(\newcommand{\vcentcolon}{\mathrel{\mathop{:}}}\) \(\newcommand{\dealcoloneq}{\vcentcolon\mathrel{\mkern-1.2mu}=}\) \(\newcommand{\jump}[1]{\left[\!\left[ #1 \right]\!\right]}\) \(\newcommand{\average}[1]{\left\{\!\left\{ #1 \right\}\!\right\}}\)
sparse_matrix_collection.h
1 // ---------------------------------------------------------------------
2 //
3 // Copyright (C) 2003 - 2018 by the deal.II authors
4 //
5 // This file is part of the deal.II library.
6 //
7 // The deal.II library is free software; you can use it, redistribute
8 // it, and/or modify it under the terms of the GNU Lesser General
9 // Public License as published by the Free Software Foundation; either
10 // version 2.1 of the License, or (at your option) any later version.
11 // The full text of the license can be found in the file LICENSE.md at
12 // the top level directory of deal.II.
13 //
14 // ---------------------------------------------------------------------
15 
16 #ifndef dealii_mg_sparse_matrix_collection_h
17 #define dealii_mg_sparse_matrix_collection_h
18 
19 #include <deal.II/base/config.h>
20 
21 #include <deal.II/base/mg_level_object.h>
22 
23 #include <deal.II/lac/dynamic_sparsity_pattern.h>
24 #include <deal.II/lac/sparse_matrix.h>
25 #include <deal.II/lac/vector.h>
26 
27 #include <deal.II/multigrid/mg_base.h>
28 #include <deal.II/multigrid/mg_tools.h>
29 
30 #include <memory>
31 
32 DEAL_II_NAMESPACE_OPEN
33 
34 namespace mg
35 {
43  template <typename number>
45  {
46  public:
47  void
48  resize(const unsigned int minlevel, const unsigned int maxlevel);
49 
50  template <typename DoFHandlerType>
51  void
52  reinit(const DoFHandlerType &dof_handler);
53 
54  void
55  set_zero();
56 
58  MGLevelObject<SparsityPattern> sparsity_edge;
59 
65  };
66 
67 
68  template <typename number>
69  void
70  SparseMatrixCollection<number>::resize(const unsigned int minlevel,
71  const unsigned int maxlevel)
72  {
73  matrix.resize(minlevel, maxlevel);
74  matrix.clear_elements();
75  matrix_up.resize(minlevel + 1, maxlevel);
76  matrix_up.clear_elements();
77  matrix_down.resize(minlevel + 1, maxlevel);
78  matrix_down.clear_elements();
79  matrix_in.resize(minlevel, maxlevel);
80  matrix_in.clear_elements();
81  matrix_out.resize(minlevel, maxlevel);
82  matrix_out.clear_elements();
83  sparsity.resize(minlevel, maxlevel);
84  sparsity_edge.resize(minlevel, maxlevel);
85  }
86 
87 
88  template <typename number>
89  template <typename DoFHandlerType>
90  void
91  SparseMatrixCollection<number>::reinit(const DoFHandlerType &dof_handler)
92  {
93  AssertIndexRange(sparsity.max_level(),
94  dof_handler.get_triangulation().n_levels());
95 
96  for (unsigned int level = sparsity.min_level();
97  level <= sparsity.max_level();
98  ++level)
99  {
100  DynamicSparsityPattern dsp(dof_handler.n_dofs(level));
101  MGTools::make_flux_sparsity_pattern(dof_handler, dsp, level);
102  sparsity[level].copy_from(dsp);
103  matrix[level].reinit(sparsity[level]);
104  matrix_in[level].reinit(sparsity[level]);
105  matrix_out[level].reinit(sparsity[level]);
106  if (level > 0)
107  {
108  DynamicSparsityPattern ci_sparsity;
109  ci_sparsity.reinit(dof_handler.n_dofs(level - 1),
110  dof_handler.n_dofs(level));
112  ci_sparsity,
113  level);
114  sparsity_edge[level].copy_from(ci_sparsity);
115  matrix_up[level].reinit(sparsity_edge[level]);
116  matrix_down[level].reinit(sparsity_edge[level]);
117  }
118  }
119  }
120 
121  template <typename number>
122  void
124  {
125  matrix = 0.;
126  matrix_in = 0.;
127  matrix_out = 0.;
128  matrix_up = 0.;
129  matrix_down = 0.;
130  }
131 
132 } // namespace mg
133 
134 DEAL_II_NAMESPACE_CLOSE
135 
136 #endif
unsigned int max_level() const
Definition: mg.h:81
void make_flux_sparsity_pattern(const DoFHandler< dim, spacedim > &dof_handler, SparsityPatternType &sparsity, const unsigned int level)
Definition: mg_tools.cc:599
#define AssertIndexRange(index, range)
Definition: exceptions.h:1641
unsigned int min_level() const
void reinit(const size_type m, const size_type n, const IndexSet &rowset=IndexSet())
void resize(const unsigned int new_minlevel, const unsigned int new_maxlevel)
void make_flux_sparsity_pattern_edge(const DoFHandler< dim, spacedim > &dof_handler, SparsityPatternType &sparsity, const unsigned int level)
Definition: mg_tools.cc:674