typename DoFAccessor<dim>::ExcInvalidObject());
Assert (&this->dof_handler->get_fe() != 0,
typename DoFAccessor<dim>::ExcInvalidObject());
- Assert (vertex<4, ExcIndexRange (i,0,4));
+ Assert (vertex<8, ExcIndexRange (vertex,0,8));
Assert (i<this->dof_handler->get_fe().dofs_per_vertex,
ExcIndexRange (i, 0, this->dof_handler->get_fe().dofs_per_vertex));
for (unsigned int line=0; line<12; ++line)
for (unsigned int d=0; d<dofs_per_line; ++d)
*next++ = this->line(line)->mg_dof_index(d);
- for (unsigned int quad=0; quad<12; ++quad)
+ for (unsigned int quad=0; quad<6; ++quad)
for (unsigned int d=0; d<dofs_per_quad; ++d)
*next++ = this->quad(quad)->mg_dof_index(d);
for (unsigned int d=0; d<dofs_per_hex; ++d)
// for the quads
if (this->selected_fe->dofs_per_quad > 0)
- for (unsigned int q=0; q<GeometryInfo<3>::lines_per_cell; ++q)
+ for (unsigned int q=0; q<GeometryInfo<3>::quads_per_cell; ++q)
{
quad_iterator quad = cell->quad(q);
// numbered (check only first dof)
if (quad->mg_dof_index(0) == DoFHandler<3>::invalid_dof_index)
// if not: distribute dofs
- for (unsigned int d=0; d<this->selected_fe->dofs_per_line; ++d)
+ for (unsigned int d=0; d<this->selected_fe->dofs_per_quad; ++d)
quad->set_mg_dof_index (d, next_free_dof++);
};