deallog << "BDM" << poly.degree() << '<' << dim << '>' << '\t'
<< quadrature.point(k);
- poly.compute(quadrature.point(k), values, grads, grads2, thirds, fourths);
+ poly.evaluate(
+ quadrature.point(k), values, grads, grads2, thirds, fourths);
for (unsigned int i = 0; i < poly.n(); ++i)
for (unsigned int d = 0; d < dim; ++d)
deallog << "BDM" << poly.degree() << '<' << dim << '>' << points[k]
<< std::endl;
- poly.compute(points[k], values, grads, grads2, thirds, fourths);
+ poly.evaluate(points[k], values, grads, grads2, thirds, fourths);
for (unsigned int i = 0; i < poly.degree() + 1; ++i, start += dim)
for (unsigned int j = 0; j < dim; ++j)
deallog << ", " << points[i][d];
deallog << ")" << std::endl;
- poly.compute(points[i], values, grads, grads2, thirds, fourths);
+ poly.evaluate(points[i], values, grads, grads2, thirds, fourths);
// loop through the Q_1^d shape functions
for (unsigned int j = 0; j < dim * GeometryInfo<dim>::vertices_per_cell;
deallog << ", " << points[i][d];
deallog << ")" << std::endl;
- poly.compute(points[i], values, grads, grads2, thirds, fourths);
+ poly.evaluate(points[i], values, grads, grads2, thirds, fourths);
// loop through the bubble shape functions
for (unsigned int j = dim * GeometryInfo<dim>::vertices_per_cell;
deallog << "RT" << poly.degree() << '<' << dim << '>' << '\t'
<< quadrature.point(k);
- poly.compute(quadrature.point(k), values, grads, grads2, thirds, fourths);
+ poly.evaluate(
+ quadrature.point(k), values, grads, grads2, thirds, fourths);
for (unsigned int i = 0; i < poly.n(); ++i)
for (unsigned int d = 0; d < dim; ++d)
deallog << "RT_Bubbles" << poly.degree() << '<' << dim << '>' << '\t'
<< quadrature.point(k);
- poly.compute(quadrature.point(k), values, grads, grads2, thirds, fourths);
+ poly.evaluate(
+ quadrature.point(k), values, grads, grads2, thirds, fourths);
for (unsigned int i = 0; i < poly.n(); ++i)
for (unsigned int d = 0; d < dim; ++d)
std::vector<Tensor<4, dim>> thirds(0);
std::vector<Tensor<5, dim>> fourths(0);
- p.compute(x, values, gradients, seconds, thirds, fourths);
+ p.evaluate(x, values, gradients, seconds, thirds, fourths);
deallog << "Point " << x << std::endl;
for (unsigned int i = 0; i < n; ++i)