141 Assert(values.size() == this->n() || values.empty(),
143 Assert(grads.size() == this->n() || grads.empty(),
145 Assert(grad_grads.size() == this->n() || grad_grads.empty(),
147 Assert(third_derivatives.size() == this->n() || third_derivatives.empty(),
149 Assert(fourth_derivatives.size() == this->n() || fourth_derivatives.empty(),
152 std::vector<double> p_values;
153 std::vector<Tensor<1, dim>> p_grads;
154 std::vector<Tensor<2, dim>> p_grad_grads;
155 std::vector<Tensor<3, dim>> p_third_derivatives;
156 std::vector<Tensor<4, dim>> p_fourth_derivatives;
159 p_values.resize((values.empty()) ? 0 : n_sub);
160 p_grads.resize((grads.empty()) ? 0 : n_sub);
161 p_grad_grads.resize((grad_grads.empty()) ? 0 : n_sub);
162 p_third_derivatives.resize((third_derivatives.empty()) ? 0 : n_sub);
163 p_fourth_derivatives.resize((fourth_derivatives.empty()) ? 0 : n_sub);
165 for (
unsigned int d = 0; d < dim; ++d)
173 for (
unsigned int c = 0; c < dim; ++c)
174 p[c] = unit_point[(c + d) % dim];
181 p_fourth_derivatives);
183 for (
unsigned int i = 0; i < p_values.size(); ++i)
187 for (
unsigned int i = 0; i < p_grads.size(); ++i)
188 for (
unsigned int d1 = 0; d1 < dim; ++d1)
192 for (
unsigned int i = 0; i < p_grad_grads.size(); ++i)
193 for (
unsigned int d1 = 0; d1 < dim; ++d1)
194 for (
unsigned int d2 = 0; d2 < dim; ++d2)
196 [(d1 + d) % dim][(d2 + d) % dim] =
199 for (
unsigned int i = 0; i < p_third_derivatives.size(); ++i)
200 for (
unsigned int d1 = 0; d1 < dim; ++d1)
201 for (
unsigned int d2 = 0; d2 < dim; ++d2)
202 for (
unsigned int d3 = 0; d3 < dim; ++d3)
204 [(d1 + d) % dim][(d2 + d) % dim]
209 for (
unsigned int i = 0; i < p_fourth_derivatives.size(); ++i)
210 for (
unsigned int d1 = 0; d1 < dim; ++d1)
211 for (
unsigned int d2 = 0; d2 < dim; ++d2)
212 for (
unsigned int d3 = 0; d3 < dim; ++d3)
213 for (
unsigned int d4 = 0; d4 < dim; ++d4)
215 [d][(d1 + d) % dim][(d2 + d) % dim]
216 [(d3 + d) % dim][(d4 + d) % dim] =
260 const unsigned int n_dofs_face =
262 std::vector<unsigned int> lexicographic_numbering;
264 for (
unsigned int j = 0; j < n_dofs_face; ++j)
266 lexicographic_numbering.push_back(j);
268 for (
unsigned int i = n_dofs_face * 2 * dim;
271 lexicographic_numbering.push_back(i + j * (
normal_degree - 1));
272 lexicographic_numbering.push_back(n_dofs_face + j);
277 for (
unsigned int k = 0; k < layers; ++k)
280 for (
unsigned int j = n_dofs_face * 2;
283 lexicographic_numbering.push_back(j + k_add);
286 for (
unsigned int i = n_dofs_face * (2 * dim + (
normal_degree - 1));
293 for (
unsigned int j = n_dofs_face * 3;
296 lexicographic_numbering.push_back(j + k_add);
302 for (
unsigned int i = 4 * n_dofs_face; i < 5 * n_dofs_face; ++i)
303 lexicographic_numbering.push_back(i);
305 for (
unsigned int i =
309 lexicographic_numbering.push_back(i);
310 for (
unsigned int i = 5 * n_dofs_face; i < 6 * n_dofs_face; ++i)
311 lexicographic_numbering.push_back(i);
314 return lexicographic_numbering;
void evaluate(const Point< dim > &unit_point, std::vector< Tensor< 1, dim > > &values, std::vector< Tensor< 2, dim > > &grads, std::vector< Tensor< 3, dim > > &grad_grads, std::vector< Tensor< 4, dim > > &third_derivatives, std::vector< Tensor< 5, dim > > &fourth_derivatives) const override
SymmetricTensor< 2, dim, Number > d(const Tensor< 2, dim, Number > &F, const Tensor< 2, dim, Number > &dF_dt)