DataOut determines whether something is a DoF or cell vector by the
size of the vector. This will fail if someone is using DG0 elements,
since there the number of elements of both types of vectors is the
- same, but the ordering will usually differ.
+ same, but the ordering will usually differ. Errors cannot be plotted
+ anyway, since they are float vectors.
+
+
+Move GeometryInfo to base
vector<Tensor<2,dim> >(quadrature.n_quadrature_points)),
unit_shape_gradients_transform(fe.n_transform_functions(),
vector<Tensor<1,dim> >(quadrature.n_quadrature_points)),
- unit_quadrature_points(quadrature.get_quad_points())
+ unit_quadrature_points(quadrature.get_points())
{
Assert ((update_flags & update_normal_vectors) == false,
ExcInvalidUpdateFlag());
update_flags,
fe)
{
- unit_face_quadrature_points = quadrature.get_quad_points();
+ unit_face_quadrature_points = quadrature.get_points();
weights = quadrature.get_weights ();
// set up an array of the unit points
Assert ((update_flags & update_support_points) == false,
ExcInvalidUpdateFlag());
- unit_face_quadrature_points = quadrature.get_quad_points();
+ unit_face_quadrature_points = quadrature.get_points();
weights = quadrature.get_weights ();
// set up an array of the unit points
for (unsigned int n=2; n<9; ++n)
{
Quadrature<1> *quadrature=gauss_quadratures[n];
- const vector<Point<1> > &points=quadrature->get_quad_points();
+ const vector<Point<1> > &points=quadrature->get_points();
const vector<double> &weights=quadrature->get_weights();
deallog << "Gauss" << n;