49 const unsigned int start_index2d[6] = {0, 1, 4, 10, 20, 35};
50 const double points2d[35][2] = {
51 {0, 0}, {0, 0}, {1, 0}, {0, 1}, {0, 0},
52 {1, 0}, {0, 1}, {1, 1}, {0.5, 0}, {0, 0.5},
53 {0, 0}, {1, 0}, {0, 1}, {1, 1}, {1. / 3., 0},
54 {2. / 3., 0}, {0, 1. / 3.}, {0, 2. / 3.}, {0.5, 1}, {1, 0.5},
55 {0, 0}, {1, 0}, {0, 1}, {1, 1}, {0.25, 0},
56 {0.5, 0}, {0.75, 0}, {0, 0.25}, {0, 0.5}, {0, 0.75},
57 {1. / 3., 1}, {2. / 3., 1}, {1, 1. / 3.}, {1, 2. / 3.}, {0.5, 0.5}};
65 const unsigned int start_index3d[6] = {0, 1, 5, 15 };
66 const double points3d[35][3] = {{0, 0, 0},
85 generate_unit_points(
const unsigned int, std::vector<
Point<dim>> &);
89 generate_unit_points(
const unsigned int k, std::vector<
Point<1>> &p)
92 const double h = 1. / k;
93 for (
unsigned int i = 0; i < p.size(); ++i)
99 generate_unit_points(
const unsigned int k, std::vector<
Point<2>> &p)
102 Assert(p.size() == start_index2d[k + 1] - start_index2d[k],
104 for (
unsigned int i = 0; i < p.size(); ++i)
106 p[i](0) = points2d[start_index2d[k] + i][0];
107 p[i](1) = points2d[start_index2d[k] + i][1];
113 generate_unit_points(
const unsigned int k, std::vector<
Point<3>> &p)
116 Assert(p.size() == start_index3d[k + 1] - start_index3d[k],
118 for (
unsigned int i = 0; i < p.size(); ++i)
120 p[i](0) = points3d[start_index3d[k] + i][0];
121 p[i](1) = points3d[start_index3d[k] + i][1];
122 p[i](2) = points3d[start_index3d[k] + i][2];
203 if (source_dgp_monomial)
208 const unsigned int m = interpolation_matrix.
m();
209 const unsigned int n = interpolation_matrix.
n();
212 Assert(m == this->n_dofs_per_cell(),
217 const unsigned int min_mn =
218 interpolation_matrix.
m() < interpolation_matrix.
n() ?
219 interpolation_matrix.
m() :
220 interpolation_matrix.
n();
222 for (
unsigned int i = 0; i < min_mn; ++i)
223 interpolation_matrix(i, i) = 1.;
227 std::vector<Point<dim>> unit_points(this->n_dofs_per_cell());
228 internal::FE_DGPMonomial::generate_unit_points(this->degree, unit_points);
233 for (
unsigned int k = 0; k < unit_points.size(); ++k)
234 source_fe_matrix(k, j) = source_fe.
shape_value(j, unit_points[k]);
237 this->n_dofs_per_cell());
238 for (
unsigned int j = 0; j < this->n_dofs_per_cell(); ++j)
239 for (
unsigned int k = 0; k < unit_points.size(); ++k)
241 this->poly_space->compute_value(j, unit_points[k]);
245 this_matrix.
mmult(interpolation_matrix, source_fe_matrix);