ExcInternalError());
const Quadrature<dim> quadrature(
- GeometryInfo<dim>::project_to_unit_cell(cell_point.second));
+ cell_point.first->reference_cell().closest_point(cell_point.second));
FEValues<dim> fe_values(mapping, fe, quadrature, update_gradients);
fe_values.reinit(cell_point.first);
ExcInternalError());
const Quadrature<dim> quadrature(
- GeometryInfo<dim>::project_to_unit_cell(cell_point.second));
+ cell_point.first->reference_cell().closest_point(cell_point.second));
hp::FEValues<dim, spacedim> hp_fe_values(mapping,
fe,
hp::QCollection<dim>(quadrature),
ExcInternalError());
const Quadrature<dim> quadrature(
- GeometryInfo<dim>::project_to_unit_cell(cell_point.second));
+ cell_point.first->reference_cell().closest_point(cell_point.second));
FEValues<dim> fe_values(mapping, fe, quadrature, update_values);
fe_values.reinit(cell_point.first);
ExcInternalError());
const Quadrature<dim> quadrature(
- GeometryInfo<dim>::project_to_unit_cell(cell_point.second));
+ cell_point.first->reference_cell().closest_point(cell_point.second));
hp::FEValues<dim, spacedim> hp_fe_values(mapping,
fe,
hp::QCollection<dim>(quadrature),
ExcInternalError());
const Quadrature<dim> quadrature(
- GeometryInfo<dim>::project_to_unit_cell(cell_point.second));
+ cell_point.first->reference_cell().closest_point(cell_point.second));
FEValues<dim> fe_values(mapping, fe, quadrature, update_values);
fe_values.reinit(cell_point.first);
AssertThrow(cell_point.first.state() == IteratorState::valid,
ExcPointNotAvailableHere());
- Quadrature<dim> q(
- GeometryInfo<dim>::project_to_unit_cell(cell_point.second));
+ const Quadrature<dim> quadrature(
+ cell_point.first->reference_cell().closest_point(cell_point.second));
FEValues<dim, spacedim> fe_values(mapping,
dof_handler.get_fe(),
- q,
+ quadrature,
UpdateFlags(update_values));
fe_values.reinit(cell_point.first);
AssertThrow(cell_point.first.state() == IteratorState::valid,
ExcPointNotAvailableHere());
- Quadrature<dim> q(
- GeometryInfo<dim>::project_to_unit_cell(cell_point.second));
+ const Quadrature<dim> quadrature(
+ cell_point.first->reference_cell().closest_point(cell_point.second));
FEValues<dim> fe_values(mapping[cell_point.first->active_fe_index()],
cell_point.first->get_fe(),
- q,
+ quadrature,
UpdateFlags(update_values));
fe_values.reinit(cell_point.first);
AssertThrow(cell_point.first.state() == IteratorState::valid,
ExcPointNotAvailableHere());
- const Quadrature<dim> q(
- GeometryInfo<dim>::project_to_unit_cell(cell_point.second));
+ const Quadrature<dim> quadrature(
+ cell_point.first->reference_cell().closest_point(cell_point.second));
const FEValuesExtractors::Vector vec(0);
FEValues<dim, spacedim> fe_values(mapping,
dof_handler.get_fe(),
- q,
+ quadrature,
UpdateFlags(update_values));
fe_values.reinit(cell_point.first);
AssertThrow(cell_point.first.state() == IteratorState::valid,
ExcPointNotAvailableHere());
- Quadrature<dim> q(
- GeometryInfo<dim>::project_to_unit_cell(cell_point.second));
+ const Quadrature<dim> quadrature(
+ cell_point.first->reference_cell().closest_point(cell_point.second));
const FEValuesExtractors::Vector vec(0);
FEValues<dim> fe_values(mapping[cell_point.first->active_fe_index()],
cell_point.first->get_fe(),
- q,
+ quadrature,
UpdateFlags(update_values));
fe_values.reinit(cell_point.first);