49 const std::vector<unsigned int> &renumber)
51 Assert(renumber.size() == index_map.size(),
55 for (
unsigned int i = 0; i < index_map.size(); ++i)
56 index_map_inverse[index_map[i]] = i;
58 std::vector<unsigned int> renumber_base;
59 renumber_base.reserve(tensor_polys.n());
60 for (
unsigned int i = 0; i < tensor_polys.n(); ++i)
61 renumber_base.push_back(renumber[i]);
63 tensor_polys.set_numbering(renumber_base);
131 std::vector<double> &values,
137 Assert(values.size() == tensor_polys.n() + 1 || values.empty(),
139 Assert(grads.size() == tensor_polys.n() + 1 || grads.empty(),
141 Assert(grad_grads.size() == tensor_polys.n() + 1 || grad_grads.empty(),
143 Assert(third_derivatives.size() == tensor_polys.n() + 1 ||
144 third_derivatives.empty(),
146 tensor_polys.n() + 1,
148 Assert(fourth_derivatives.size() == tensor_polys.n() + 1 ||
149 fourth_derivatives.empty(),
151 tensor_polys.n() + 1,
156 bool do_values =
false, do_grads =
false, do_grad_grads =
false;
157 bool do_3rd_derivatives =
false, do_4th_derivatives =
false;
158 if (values.empty() ==
false)
163 if (grads.empty() ==
false)
168 if (grad_grads.empty() ==
false)
170 grad_grads.pop_back();
171 do_grad_grads =
true;
173 if (third_derivatives.empty() ==
false)
175 third_derivatives.resize(tensor_polys.n());
176 do_3rd_derivatives =
true;
178 if (fourth_derivatives.empty() ==
false)
180 fourth_derivatives.resize(tensor_polys.n());
181 do_4th_derivatives =
true;
184 tensor_polys.evaluate(
185 p, values, grads, grad_grads, third_derivatives, fourth_derivatives);
190 values.push_back(1.);
192 grads.emplace_back();
194 grad_grads.emplace_back();
195 if (do_3rd_derivatives)
196 third_derivatives.emplace_back();
197 if (do_4th_derivatives)
198 fourth_derivatives.emplace_back();
void evaluate(const Point< dim > &unit_point, std::vector< double > &values, std::vector< Tensor< 1, dim > > &grads, std::vector< Tensor< 2, dim > > &grad_grads, std::vector< Tensor< 3, dim > > &third_derivatives, std::vector< Tensor< 4, dim > > &fourth_derivatives) const override