*
* Computed once.
*/
- std::vector<std::vector<double> > cell_manifold_quadratures_weights;
+ std::vector<std::vector<double> > cell_manifold_quadrature_weights;
// /**
// * Values of shape function derivatives. Access by function @p derivative.
void
MappingManifold<dim,spacedim>::InternalData::compute_manifold_quadrature_weights (const Quadrature<dim> &quad)
{
- static FE_Q<dim> fe_q(1);
- cell_manifold_quadratures_weights.resize(quad.size(), std::vector<double>(GeometryInfo<dim>::vertices_per_cell));
+ FE_Q<dim> fe_q(1);
+ cell_manifold_quadrature_weights.resize(quad.size(), std::vector<double>(GeometryInfo<dim>::vertices_per_cell));
for (unsigned int q=0; q<quad.size(); ++q)
{
for (unsigned int i=0; i<GeometryInfo<dim>::vertices_per_cell; ++i)
{
- cell_manifold_quadratures_weights[q][i] = fe_q.shape_value(i, quad.point(q));
+ cell_manifold_quadrature_weights[q][i] = fe_q.shape_value(i, quad.point(q));
}
}
}
{
const UpdateFlags update_flags = data.update_each;
- AssertDimension(quadrature_points.size(),
- data.cell_manifold_quadratures_weights.size());
-
- std::vector<Point<spacedim> > vertices;
- for (unsigned int v=0; v<GeometryInfo<dim>::vertices_per_cell; ++v)
- vertices[v] = cell->vertex(v);
if (update_flags & update_quadrature_points)
{
+
+ AssertDimension(quadrature_points.size(),
+ data.cell_manifold_quadrature_weights.size());
+
+ std::vector<Point<spacedim> > vertices;
+ for (unsigned int v=0; v<GeometryInfo<dim>::vertices_per_cell; ++v)
+ vertices.push_back(cell->vertex(v));
+
for (unsigned int point=0; point<quadrature_points.size(); ++point)
{
quadrature_points[point] = cell->get_manifold().
- get_new_point(Quadrature<spacedim>(vertices, data.cell_manifold_quadratures_weights[point]));
+ get_new_point(Quadrature<spacedim>(vertices,
+ data.cell_manifold_quadrature_weights[point]));
}
}
}
cell = triangulation.begin_active();
+ std::ofstream two_be_plotted("plot.txt");
for (; cell!=endc; ++cell)
{
fe_values.reinit(cell);
deallog << "Expected: " << pq << ", got: "
<< q_points_from_fe_values[q] << std::endl;
}
+ two_be_plotted << pq << " " << q_points_from_fe_values[q] << std::endl;
}
}
deallog << "OK" << std::endl;