const unsigned int local_index,
dealii::internal::int2type<1>)
{
- return dof_handler.levels[obj_level]->lines.
+ return dof_handler.levels[obj_level]->dof_object.
get_dof_index (dof_handler,
obj_index,
fe_index,
dealii::internal::int2type<1>,
const unsigned int global_index)
{
- dof_handler.levels[obj_level]->lines.
+ dof_handler.levels[obj_level]->dof_object.
set_dof_index (dof_handler,
obj_index,
fe_index,
const unsigned int local_index,
dealii::internal::int2type<2>)
{
- return dof_handler.levels[obj_level]->quads.
+ return dof_handler.levels[obj_level]->dof_object.
get_dof_index (dof_handler,
obj_index,
fe_index,
dealii::internal::int2type<2>,
const unsigned int global_index)
{
- dof_handler.levels[obj_level]->quads.
+ dof_handler.levels[obj_level]->dof_object.
set_dof_index (dof_handler,
obj_index,
fe_index,
const unsigned int local_index,
dealii::internal::int2type<3>)
{
- return dof_handler.levels[obj_level]->hexes.
+ return dof_handler.levels[obj_level]->dof_object.
get_dof_index (dof_handler,
obj_index,
fe_index,
dealii::internal::int2type<3>,
const unsigned int global_index)
{
- dof_handler.levels[obj_level]->hexes.
+ dof_handler.levels[obj_level]->dof_object.
set_dof_index (dof_handler,
obj_index,
fe_index,
const unsigned int local_index,
const dealii::internal::int2type<1> &)
{
- return dof_handler.levels[obj_level]->lines.
+ return dof_handler.levels[obj_level]->dof_object.
get_dof_index (dof_handler,
obj_index,
fe_index,
const dealii::internal::int2type<1> &,
const unsigned int global_index)
{
- dof_handler.levels[obj_level]->lines.
+ dof_handler.levels[obj_level]->dof_object.
set_dof_index (dof_handler,
obj_index,
fe_index,
const unsigned int local_index,
const dealii::internal::int2type<2> &)
{
- return dof_handler.levels[obj_level]->quads.
+ return dof_handler.levels[obj_level]->dof_object.
get_dof_index (dof_handler,
obj_index,
fe_index,
const dealii::internal::int2type<2> &,
const unsigned int global_index)
{
- dof_handler.levels[obj_level]->quads.
+ dof_handler.levels[obj_level]->dof_object.
set_dof_index (dof_handler,
obj_index,
fe_index,
const unsigned int local_index,
const dealii::internal::int2type<3> &)
{
- return dof_handler.levels[obj_level]->hexes.
+ return dof_handler.levels[obj_level]->dof_object.
get_dof_index (dof_handler,
obj_index,
fe_index,
const dealii::internal::int2type<3> &,
const unsigned int global_index)
{
- dof_handler.levels[obj_level]->hexes.
+ dof_handler.levels[obj_level]->dof_object.
set_dof_index (dof_handler,
obj_index,
fe_index,
const unsigned int fe_index,
const dealii::internal::int2type<1> &)
{
- return dof_handler.levels[obj_level]->lines.fe_index_is_active(dof_handler,
+ return dof_handler.levels[obj_level]->dof_object.fe_index_is_active(dof_handler,
obj_index,
fe_index,
obj_level);
const unsigned int obj_index,
const dealii::internal::int2type<1> &)
{
- return dof_handler.levels[obj_level]->lines.n_active_fe_indices (dof_handler,
+ return dof_handler.levels[obj_level]->dof_object.n_active_fe_indices (dof_handler,
obj_index);
}
const unsigned int n,
const dealii::internal::int2type<1> &)
{
- return dof_handler.levels[obj_level]->lines.nth_active_fe_index (dof_handler,
+ return dof_handler.levels[obj_level]->dof_object.nth_active_fe_index (dof_handler,
obj_level,
obj_index,
n);
const unsigned int fe_index,
const dealii::internal::int2type<2> &)
{
- return dof_handler.levels[obj_level]->quads.fe_index_is_active(dof_handler,
+ return dof_handler.levels[obj_level]->dof_object.fe_index_is_active(dof_handler,
obj_index,
fe_index,
obj_level);
const unsigned int obj_index,
const dealii::internal::int2type<2> &)
{
- return dof_handler.levels[obj_level]->quads.n_active_fe_indices (dof_handler,
+ return dof_handler.levels[obj_level]->dof_object.n_active_fe_indices (dof_handler,
obj_index);
}
const unsigned int n,
const dealii::internal::int2type<2> &)
{
- return dof_handler.levels[obj_level]->quads.nth_active_fe_index (dof_handler,
+ return dof_handler.levels[obj_level]->dof_object.nth_active_fe_index (dof_handler,
obj_level,
obj_index,
n);
const unsigned int fe_index,
const dealii::internal::int2type<3> &)
{
- return dof_handler.levels[obj_level]->hexes.fe_index_is_active(dof_handler,
+ return dof_handler.levels[obj_level]->dof_object.fe_index_is_active(dof_handler,
obj_index,
fe_index,
obj_level);
const unsigned int obj_index,
const dealii::internal::int2type<3> &)
{
- return dof_handler.levels[obj_level]->hexes.n_active_fe_indices (dof_handler,
+ return dof_handler.levels[obj_level]->dof_object.n_active_fe_indices (dof_handler,
obj_index);
}
const unsigned int n,
const dealii::internal::int2type<3> &)
{
- return dof_handler.levels[obj_level]->hexes.nth_active_fe_index (dof_handler,
+ return dof_handler.levels[obj_level]->dof_object.nth_active_fe_index (dof_handler,
obj_level,
obj_index,
n);
* The object containing dof-indices
* and related access-functions
*/
- DoFObjects<1> lines;
+ DoFObjects<1> dof_object;
/**
* Determine an estimate for the
* The object containing dof-indices
* and related access-functions
*/
- internal::DoFHandler::DoFObjects<2> quads;
+ internal::DoFHandler::DoFObjects<2> dof_object;
/**
* Determine an estimate for the
* The object containing dof-indices
* and related access-functions
*/
- internal::DoFHandler::DoFObjects<3> hexes;
+ internal::DoFHandler::DoFObjects<3> dof_object;
/**
* Determine an estimate for the
const unsigned int version)
{
this->DoFLevel<0>::serialize (ar, version);
- ar & lines;
+ ar & dof_object;
}
const unsigned int version)
{
this->DoFLevel<0>::serialize (ar, version);
- ar & quads;
+ ar & dof_object;
}
const unsigned int version)
{
this->DoFLevel<0>::serialize (ar, version);
- ar & hexes;
+ ar & dof_object;
}
}
}
* store the dof-indices and related functions of
* lines
*/
- internal::hp::DoFObjects<1> lines;
+ internal::hp::DoFObjects<1> dof_object;
/**
* Determine an estimate for the
* store the dof-indices and related functions of
* quads
*/
- internal::hp::DoFObjects<2> quads;
+ internal::hp::DoFObjects<2> dof_object;
/**
* Determine an estimate for the
* store the dof-indices and related functions of
* hexes
*/
- internal::hp::DoFObjects<3> hexes;
+ internal::hp::DoFObjects<3> dof_object;
/**
* Determine an estimate for the
dof_handler.levels
.push_back (new internal::DoFHandler::DoFLevel<1>);
- dof_handler.levels.back()->lines.dofs
+ dof_handler.levels.back()->dof_object.dofs
.resize (dof_handler.tria->n_raw_lines(i) *
dof_handler.selected_fe->dofs_per_line,
DoFHandler<1,spacedim>::invalid_dof_index);
{
dof_handler.levels.push_back (new internal::DoFHandler::DoFLevel<2>);
- dof_handler.levels.back()->quads.dofs
+ dof_handler.levels.back()->dof_object.dofs
.resize (dof_handler.tria->n_raw_quads(i) *
dof_handler.selected_fe->dofs_per_quad,
DoFHandler<2,spacedim>::invalid_dof_index);
{
dof_handler.levels.push_back (new internal::DoFHandler::DoFLevel<3>);
- dof_handler.levels.back()->hexes.dofs
+ dof_handler.levels.back()->dof_object.dofs
.resize (dof_handler.tria->n_raw_hexs(i) *
dof_handler.selected_fe->dofs_per_hex,
DoFHandler<3,spacedim>::invalid_dof_index);
for (unsigned int level=0; level<dof_handler.levels.size(); ++level)
for (std::vector<unsigned int>::iterator
- i=dof_handler.levels[level]->lines.dofs.begin();
- i!=dof_handler.levels[level]->lines.dofs.end(); ++i)
+ i=dof_handler.levels[level]->dof_object.dofs.begin();
+ i!=dof_handler.levels[level]->dof_object.dofs.end(); ++i)
if (*i != DoFHandler<1,spacedim>::invalid_dof_index)
*i = new_numbers[*i];
for (unsigned int level=0; level<dof_handler.levels.size(); ++level)
{
for (std::vector<unsigned int>::iterator
- i=dof_handler.levels[level]->quads.dofs.begin();
- i!=dof_handler.levels[level]->quads.dofs.end(); ++i)
+ i=dof_handler.levels[level]->dof_object.dofs.begin();
+ i!=dof_handler.levels[level]->dof_object.dofs.end(); ++i)
if (*i != DoFHandler<2,spacedim>::invalid_dof_index)
*i = ((indices.n_elements() == 0) ?
new_numbers[*i] :
for (unsigned int level=0; level<dof_handler.levels.size(); ++level)
{
for (std::vector<unsigned int>::iterator
- i=dof_handler.levels[level]->hexes.dofs.begin();
- i!=dof_handler.levels[level]->hexes.dofs.end(); ++i)
+ i=dof_handler.levels[level]->dof_object.dofs.begin();
+ i!=dof_handler.levels[level]->dof_object.dofs.end(); ++i)
if (*i != DoFHandler<3,spacedim>::invalid_dof_index)
*i = ((indices.n_elements() == 0) ?
new_numbers[*i] :
DoFLevel<1>::memory_consumption () const
{
return (DoFLevel<0>::memory_consumption () +
- MemoryConsumption::memory_consumption (lines));
+ MemoryConsumption::memory_consumption (dof_object));
}
DoFLevel<2>::memory_consumption () const
{
return (DoFLevel<0>::memory_consumption () +
- MemoryConsumption::memory_consumption (quads));
+ MemoryConsumption::memory_consumption (dof_object));
}
DoFLevel<3>::memory_consumption () const
{
return (DoFLevel<0>::memory_consumption () +
- MemoryConsumption::memory_consumption (hexes));
+ MemoryConsumption::memory_consumption (dof_object));
}
}
// finite element is used for it
for (unsigned int level=0; level<dof_handler.tria->n_levels(); ++level)
{
- dof_handler.levels[level]->lines.dof_offsets
+ dof_handler.levels[level]->dof_object.dof_offsets
= std::vector<unsigned int> (dof_handler.tria->n_raw_lines(level),
DoFHandler<dim,spacedim>::invalid_dof_index);
cell!=dof_handler.end_active(level); ++cell)
if (!cell->has_children())
{
- dof_handler.levels[level]->lines.dof_offsets[cell->index()] = next_free_dof;
+ dof_handler.levels[level]->dof_object.dof_offsets[cell->index()] = next_free_dof;
next_free_dof += cell->get_fe().dofs_per_line;
}
- dof_handler.levels[level]->lines.dofs
+ dof_handler.levels[level]->dof_object.dofs
= std::vector<unsigned int> (next_free_dof,
DoFHandler<dim,spacedim>::invalid_dof_index);
}
if (!cell->has_children())
counter += cell->get_fe().dofs_per_line;
- Assert (dof_handler.levels[level]->lines.dofs.size() == counter,
+ Assert (dof_handler.levels[level]->dof_object.dofs.size() == counter,
ExcInternalError());
Assert (static_cast<unsigned int>
- (std::count (dof_handler.levels[level]->lines.dof_offsets.begin(),
- dof_handler.levels[level]->lines.dof_offsets.end(),
+ (std::count (dof_handler.levels[level]->dof_object.dof_offsets.begin(),
+ dof_handler.levels[level]->dof_object.dof_offsets.end(),
DoFHandler<dim,spacedim>::invalid_dof_index))
==
dof_handler.tria->n_raw_lines(level) - dof_handler.tria->n_active_lines(level),
// finite element is used for it
for (unsigned int level=0; level<dof_handler.tria->n_levels(); ++level)
{
- dof_handler.levels[level]->quads.dof_offsets
+ dof_handler.levels[level]->dof_object.dof_offsets
= std::vector<unsigned int> (dof_handler.tria->n_raw_quads(level),
DoFHandler<dim,spacedim>::invalid_dof_index);
cell!=dof_handler.end_active(level); ++cell)
if (!cell->has_children())
{
- dof_handler.levels[level]->quads.dof_offsets[cell->index()] = next_free_dof;
+ dof_handler.levels[level]->dof_object.dof_offsets[cell->index()] = next_free_dof;
next_free_dof += cell->get_fe().dofs_per_quad;
}
- dof_handler.levels[level]->quads.dofs
+ dof_handler.levels[level]->dof_object.dofs
= std::vector<unsigned int> (next_free_dof,
DoFHandler<dim,spacedim>::invalid_dof_index);
}
if (!cell->has_children())
counter += cell->get_fe().dofs_per_quad;
- Assert (dof_handler.levels[level]->quads.dofs.size() == counter,
+ Assert (dof_handler.levels[level]->dof_object.dofs.size() == counter,
ExcInternalError());
Assert (static_cast<unsigned int>
- (std::count (dof_handler.levels[level]->quads.dof_offsets.begin(),
- dof_handler.levels[level]->quads.dof_offsets.end(),
+ (std::count (dof_handler.levels[level]->dof_object.dof_offsets.begin(),
+ dof_handler.levels[level]->dof_object.dof_offsets.end(),
DoFHandler<dim,spacedim>::invalid_dof_index))
==
dof_handler.tria->n_raw_quads(level) - dof_handler.tria->n_active_quads(level),
// finite element is used for it
for (unsigned int level=0; level<dof_handler.tria->n_levels(); ++level)
{
- dof_handler.levels[level]->hexes.dof_offsets
+ dof_handler.levels[level]->dof_object.dof_offsets
= std::vector<unsigned int> (dof_handler.tria->n_raw_hexs(level),
DoFHandler<dim,spacedim>::invalid_dof_index);
cell!=dof_handler.end_active(level); ++cell)
if (!cell->has_children())
{
- dof_handler.levels[level]->hexes.dof_offsets[cell->index()] = next_free_dof;
+ dof_handler.levels[level]->dof_object.dof_offsets[cell->index()] = next_free_dof;
next_free_dof += cell->get_fe().dofs_per_hex;
}
- dof_handler.levels[level]->hexes.dofs
+ dof_handler.levels[level]->dof_object.dofs
= std::vector<unsigned int> (next_free_dof,
DoFHandler<dim,spacedim>::invalid_dof_index);
}
if (!cell->has_children())
counter += cell->get_fe().dofs_per_hex;
- Assert (dof_handler.levels[level]->hexes.dofs.size() == counter,
+ Assert (dof_handler.levels[level]->dof_object.dofs.size() == counter,
ExcInternalError());
Assert (static_cast<unsigned int>
- (std::count (dof_handler.levels[level]->hexes.dof_offsets.begin(),
- dof_handler.levels[level]->hexes.dof_offsets.end(),
+ (std::count (dof_handler.levels[level]->dof_object.dof_offsets.begin(),
+ dof_handler.levels[level]->dof_object.dof_offsets.end(),
DoFHandler<dim,spacedim>::invalid_dof_index))
==
dof_handler.tria->n_raw_hexs(level) - dof_handler.tria->n_active_hexs(level),
DoFLevel<1>::memory_consumption () const
{
return (DoFLevel<0>::memory_consumption() +
- MemoryConsumption::memory_consumption (lines));
+ MemoryConsumption::memory_consumption (dof_object));
}
DoFLevel<2>::memory_consumption () const
{
return (DoFLevel<0>::memory_consumption () +
- MemoryConsumption::memory_consumption (quads));
+ MemoryConsumption::memory_consumption (dof_object));
}
DoFLevel<3>::memory_consumption () const
{
return (DoFLevel<0>::memory_consumption () +
- MemoryConsumption::memory_consumption (hexes));
+ MemoryConsumption::memory_consumption (dof_object));
}
const unsigned int local_index,
const int2type<1>)
{
- return mg_level.lines.
+ return mg_level.dof_object.
get_dof_index (mg_dof_handler,
obj_index,
fe_index,
const unsigned int global_index,
const int2type<1>)
{
- mg_level.lines.
+ mg_level.dof_object.
set_dof_index (mg_dof_handler,
obj_index,
fe_index,
const unsigned int local_index,
const int2type<2>)
{
- return mg_level.quads.
+ return mg_level.dof_object.
get_dof_index (mg_dof_handler,
obj_index,
fe_index,
const unsigned int global_index,
const int2type<2>)
{
- mg_level.quads.
+ mg_level.dof_object.
set_dof_index (mg_dof_handler,
obj_index,
fe_index,
const unsigned int local_index,
const int2type<3>)
{
- return mg_level.hexes.
+ return mg_level.dof_object.
get_dof_index (mg_dof_handler,
obj_index,
fe_index,
const unsigned int global_index,
const int2type<3>)
{
- mg_level.hexes.
+ mg_level.dof_object.
set_dof_index (mg_dof_handler,
obj_index,
fe_index,
{
mg_dof_handler.mg_levels.push_back (new internal::DoFHandler::DoFLevel<1>);
- mg_dof_handler.mg_levels.back()->lines.dofs = std::vector<unsigned int>(mg_dof_handler.get_tria().n_raw_lines(i) *
+ mg_dof_handler.mg_levels.back()->dof_object.dofs = std::vector<unsigned int>(mg_dof_handler.get_tria().n_raw_lines(i) *
mg_dof_handler.get_fe().dofs_per_line,
DoFHandler<1>::invalid_dof_index);
}
{
mg_dof_handler.mg_levels.push_back (new internal::DoFHandler::DoFLevel<2>);
- mg_dof_handler.mg_levels.back()->quads.dofs = std::vector<unsigned int> (mg_dof_handler.get_tria().n_raw_quads(i) *
+ mg_dof_handler.mg_levels.back()->dof_object.dofs = std::vector<unsigned int> (mg_dof_handler.get_tria().n_raw_quads(i) *
mg_dof_handler.get_fe().dofs_per_quad,
DoFHandler<2>::invalid_dof_index);
}
{
mg_dof_handler.mg_levels.push_back (new internal::DoFHandler::DoFLevel<3>);
- mg_dof_handler.mg_levels.back()->hexes.dofs
+ mg_dof_handler.mg_levels.back()->dof_object.dofs
= std::vector<unsigned int> (mg_dof_handler.get_tria().n_raw_hexs(i) *
mg_dof_handler.get_fe().dofs_per_hex,
DoFHandler<3>::invalid_dof_index);
new_numbers[i->get_index (level, d,
this->get_fe().dofs_per_vertex)]);
- for (std::vector<unsigned int>::iterator i=mg_levels[level]->lines.dofs.begin();
- i!=mg_levels[level]->lines.dofs.end(); ++i)
+ for (std::vector<unsigned int>::iterator i=mg_levels[level]->dof_object.dofs.begin();
+ i!=mg_levels[level]->dof_object.dofs.end(); ++i)
{
if (*i != DoFHandler<1>::invalid_dof_index)
{
const_cast<Triangulation<2> &>(this->get_tria()).load_user_flags (user_flags);
}
- for (std::vector<unsigned int>::iterator i=mg_levels[level]->quads.dofs.begin();
- i!=mg_levels[level]->quads.dofs.end(); ++i)
+ for (std::vector<unsigned int>::iterator i=mg_levels[level]->dof_object.dofs.begin();
+ i!=mg_levels[level]->dof_object.dofs.end(); ++i)
{
if (*i != DoFHandler<2>::invalid_dof_index)
{
}
//HEX DoFs
- for (std::vector<unsigned int>::iterator i=mg_levels[level]->hexes.dofs.begin();
- i!=mg_levels[level]->hexes.dofs.end(); ++i)
+ for (std::vector<unsigned int>::iterator i=mg_levels[level]->dof_object.dofs.begin();
+ i!=mg_levels[level]->dof_object.dofs.end(); ++i)
{
if (*i != DoFHandler<3>::invalid_dof_index)
{