for (unsigned int d=0;d<dim;++d)
{
- // guard the following statement
- Assert (dim <= 3, ExcNotImplemented());
- const QAnisotropic<dim>
- quadrature (dim == 1 ?
- QAnisotropic<dim>(high) :
- (dim == 2 ?
- QAnisotropic<dim>(((d==0) ? low : high),
- ((d==1) ? low : high)) :
- QAnisotropic<dim>(((d==0) ? low : high),
- ((d==1) ? low : high),
- ((d==2) ? low : high))));
+ QAnisotropic<dim>* quadrature;
+ if (dim == 1) quadrature = new QAnisotropic<dim>(high);
+ if (dim == 2) quadrature = new QAnisotropic<dim>(((d==0) ? low : high),
+ ((d==1) ? low : high));
+ if (dim == 3) quadrature = new QAnisotropic<dim>(((d==0) ? low : high),
+ ((d==1) ? low : high),
+ ((d==2) ? low : high));
+ Assert(dim<=3, ExcNotImplemented());
- for (unsigned int k=0;k<quadrature.n_quadrature_points;++k)
- this->generalized_support_points[current++] = quadrature.point(k);
+ for (unsigned int k=0;k<quadrature->n_quadrature_points;++k)
+ this->generalized_support_points[current++] = quadrature->point(k);
+ delete quadrature;
}
Assert (current == this->dofs_per_cell, ExcInternalError());
}