AssertIndexRange(q_index, q_collections.size());
- // set the triple of indices
- // that we want to work with
- present_fe_values_index = TableIndices<3>(fe_index, mapping_index, q_index);
-
- // first check whether we
- // already have an object for
- // this particular combination
- // of indices
+ // set the triple of indices that we want to work with
+ present_fe_values_index = {fe_index, mapping_index, q_index};
+
+ // first check whether we already have an object for this particular
+ // combination of indices
if (fe_values_table(present_fe_values_index).get() == nullptr)
fe_values_table(present_fe_values_index) =
std::make_unique<FEValuesType>((*mapping_collection)[mapping_index],
task_group +=
Threads::new_task([&, fe_index, mapping_index, q_index]() {
- fe_values_table(TableIndices<3>(fe_index, mapping_index, q_index)) =
+ fe_values_table[fe_index][mapping_index][q_index] =
std::make_unique<FEValuesType>(
(*mapping_collection)[mapping_index],
(*fe_collection)[fe_index],
const unsigned int mapping_index,
const unsigned int fe_index)
{
- // determine which indices we
- // should actually use
+ // determine which indices we should actually use
unsigned int real_q_index = q_index, real_mapping_index = mapping_index,
real_fe_index = fe_index;
AssertIndexRange(real_mapping_index, this->mapping_collection->size());
AssertIndexRange(real_fe_index, this->fe_collection->size());
- // now finally actually get the
- // corresponding object and
- // initialize it
+ // now finally actually get the corresponding object and initialize it
this->select_fe_values(real_fe_index, real_mapping_index, real_q_index)
.reinit(cell);
}
const unsigned int mapping_index,
const unsigned int fe_index)
{
- // determine which indices we
- // should actually use
+ // determine which indices we should actually use
unsigned int real_q_index = q_index, real_mapping_index = mapping_index,
real_fe_index = fe_index;
AssertIndexRange(real_mapping_index, this->mapping_collection->size());
AssertIndexRange(real_fe_index, this->fe_collection->size());
- // now finally actually get the
- // corresponding object and
- // initialize it
+ // now finally actually get the corresponding object and initialize it
this->select_fe_values(real_fe_index, real_mapping_index, real_q_index)
.reinit(cell);
}
AssertIndexRange(real_mapping_index, this->mapping_collection->size());
AssertIndexRange(real_fe_index, this->fe_collection->size());
- // now finally actually get the
- // corresponding object and
- // initialize it
+ // now finally actually get the corresponding object and initialize it
this->select_fe_values(real_fe_index, real_mapping_index, real_q_index)
.reinit(cell, face_no);
}
AssertIndexRange(real_mapping_index, this->mapping_collection->size());
AssertIndexRange(real_fe_index, this->fe_collection->size());
- // now finally actually get the
- // corresponding object and
- // initialize it
+ // now finally actually get the corresponding object and initialize it
this->select_fe_values(real_fe_index, real_mapping_index, real_q_index)
.reinit(cell, face_no);
}
const unsigned int mapping_index,
const unsigned int fe_index)
{
- // determine which indices we
- // should actually use
+ // determine which indices we should actually use
unsigned int real_q_index = q_index, real_mapping_index = mapping_index,
real_fe_index = fe_index;
AssertIndexRange(real_mapping_index, this->mapping_collection->size());
AssertIndexRange(real_fe_index, this->fe_collection->size());
- // now finally actually get the
- // corresponding object and
- // initialize it
+ // now finally actually get the corresponding object and initialize it
this->select_fe_values(real_fe_index, real_mapping_index, real_q_index)
.reinit(cell, face_no, subface_no);
}
const unsigned int mapping_index,
const unsigned int fe_index)
{
- // determine which indices we
- // should actually use
+ // determine which indices we should actually use
unsigned int real_q_index = q_index, real_mapping_index = mapping_index,
real_fe_index = fe_index;
AssertIndexRange(real_mapping_index, this->mapping_collection->size());
AssertIndexRange(real_fe_index, this->fe_collection->size());
- // now finally actually get the
- // corresponding object and
- // initialize it
+ // now finally actually get the corresponding object and initialize it
this->select_fe_values(real_fe_index, real_mapping_index, real_q_index)
.reinit(cell, face_no, subface_no);
}