template <>
void FEQuadraticSub<2>::get_support_points (const typename DoFHandler<2>::cell_iterator &cell,
- const Boundary<2>&,
- vector<Point<2> > &support_points) const {
+ const Boundary<2>&,
+ vector<Point<2> > &support_points) const {
Assert (support_points.size() == total_dofs,
ExcWrongFieldDimension (support_points.size(), total_dofs));
// to the mean coordinates of the vertices
for (unsigned int line=0; line<4; ++line)
support_points[4+line] = (cell->line(line)->vertex(0) +
- cell->line(line)->vertex(1)) / 2;
+ cell->line(line)->vertex(1)) / 2;
// same for the center of the square:
// since all four linear basis functions
// take on the value 1/4 at the center,
// the center is mapped to the mean
// coordinates of the four vertices
support_points[8] = (support_points[0] +
- support_points[1] +
- support_points[2] +
- support_points[3]) / 4;
+ support_points[1] +
+ support_points[2] +
+ support_points[3]) / 4;
};
for (unsigned int face=0; face<GeometryInfo<dim>::faces_per_cell; ++face)
QProjector<dim>::project_to_face (quadrature, face, unit_quadrature_points[face]);
- for (unsigned int face=0; face<GeometryInfo<dim>::faces_per_cell; ++face)
- for (unsigned int i=0; i<fe.total_dofs; ++i)
- for (unsigned int j=0; j<n_quadrature_points; ++j)
+ for (unsigned int i=0; i<fe.total_dofs; ++i)
+ for (unsigned int j=0; j<n_quadrature_points; ++j)
+ for (unsigned int face=0; face<GeometryInfo<dim>::faces_per_cell; ++face)
{
shape_values[face](i,j)
= fe.shape_value(i, unit_quadrature_points[face][j]);
= fe.shape_grad(i, unit_quadrature_points[face][j]);
};
- for (unsigned int face=0; face<GeometryInfo<dim>::faces_per_cell; ++face)
- for (unsigned int i=0; i<n_transform_functions; ++i)
- for (unsigned int j=0; j<n_quadrature_points; ++j)
+ for (unsigned int i=0; i<n_transform_functions; ++i)
+ for (unsigned int j=0; j<n_quadrature_points; ++j)
+ for (unsigned int face=0; face<GeometryInfo<dim>::faces_per_cell; ++face)
{
shape_values_transform[face] (i,j)
= fe.shape_value_transform (i, unit_quadrature_points[face][j]);
unit_shape_gradients_transform[face][i][j]
= fe.shape_grad_transform(i, unit_quadrature_points[face][j]);
- };
+ };
};
face, subface,
unit_quadrature_points[face*(1<<(dim-1))+subface]);
- for (unsigned int face=0; face<GeometryInfo<dim>::faces_per_cell; ++face)
- for (unsigned int subface=0; subface<GeometryInfo<dim>::subfaces_per_face; ++subface)
- for (unsigned int i=0; i<fe.total_dofs; ++i)
- for (unsigned int j=0; j<n_quadrature_points; ++j)
+ for (unsigned int i=0; i<fe.total_dofs; ++i)
+ for (unsigned int j=0; j<n_quadrature_points; ++j)
+ for (unsigned int face=0; face<GeometryInfo<dim>::faces_per_cell; ++face)
+ for (unsigned int subface=0; subface<GeometryInfo<dim>::subfaces_per_face; ++subface)
{
shape_values[face*GeometryInfo<dim>::subfaces_per_face+subface](i,j)
- = fe.shape_value(i, unit_quadrature_points[face*GeometryInfo<dim>::subfaces_per_face+subface][j]);
+ = fe.shape_value(i, unit_quadrature_points[face *
+ GeometryInfo<dim>::
+ subfaces_per_face+subface][j]);
unit_shape_gradients[face*GeometryInfo<dim>::subfaces_per_face+subface][i][j]
- = fe.shape_grad(i, unit_quadrature_points[face*GeometryInfo<dim>::subfaces_per_face+subface][j]);
+ = fe.shape_grad(i, unit_quadrature_points[face *
+ GeometryInfo<dim>::
+ subfaces_per_face+subface][j]);
};
- for (unsigned int face=0; face<GeometryInfo<dim>::faces_per_cell; ++face)
- for (unsigned int subface=0; subface<GeometryInfo<dim>::subfaces_per_face; ++subface)
- for (unsigned int i=0; i<n_transform_functions; ++i)
- for (unsigned int j=0; j<n_quadrature_points; ++j)
+ for (unsigned int i=0; i<n_transform_functions; ++i)
+ for (unsigned int j=0; j<n_quadrature_points; ++j)
+ for (unsigned int face=0; face<GeometryInfo<dim>::faces_per_cell; ++face)
+ for (unsigned int subface=0; subface<GeometryInfo<dim>::subfaces_per_face; ++subface)
{
shape_values_transform[face*GeometryInfo<dim>::subfaces_per_face+subface] (i,j)
- = fe.shape_value_transform (i, unit_quadrature_points[face*GeometryInfo<dim>::subfaces_per_face+subface][j]);
- unit_shape_gradients_transform[face*GeometryInfo<dim>::subfaces_per_face+subface][i][j]
- = fe.shape_grad_transform(i, unit_quadrature_points[face*GeometryInfo<dim>::subfaces_per_face+subface][j]);
+ = fe.shape_value_transform (i, unit_quadrature_points[face *
+ GeometryInfo<dim>::
+ subfaces_per_face +
+ subface][j]);
+ unit_shape_gradients_transform[face *
+ GeometryInfo<dim>::subfaces_per_face +
+ subface][i][j]
+ = fe.shape_grad_transform(i, unit_quadrature_points[face *
+ GeometryInfo<dim>::
+ subfaces_per_face +
+ subface][j]);
};
};