Assert (dofs_per_quad==0, ExcInternalError());
const unsigned int dim=1;
- for (unsigned int i=0; i<(1<<dim); ++i)
+ for (unsigned int i=0; i<GeometryInfo<dim>::children_per_cell; ++i)
{
restriction[i].reinit (total_dofs, total_dofs);
prolongation[i].reinit (total_dofs, total_dofs);
dofs_per_quad)
{
const unsigned int dim=2;
- for (unsigned int i=0; i<(1<<dim); ++i)
+ for (unsigned int i=0; i<GeometryInfo<dim>::children_per_cell; ++i)
{
restriction[i].reinit (total_dofs, total_dofs);
prolongation[i].reinit (total_dofs, total_dofs);
template <int dim>
const dFMatrix &
FiniteElementBase<dim>::restrict (const unsigned int child) const {
- Assert (child<(1<<dim), ExcInvalidIndex(child));
+ Assert (child<GeometryInfo<dim>::children_per_cell, ExcInvalidIndex(child));
return restriction[child];
};
template <int dim>
const dFMatrix &
FiniteElementBase<dim>::prolongate (const unsigned int child) const {
- Assert (child<(1<<dim), ExcInvalidIndex(child));
+ Assert (child<GeometryInfo<dim>::children_per_cell, ExcInvalidIndex(child));
return prolongation[child];
};