mg_additional_data);
MatrixFree<dim, LevelNumberType> mg_level;
- mg_level.reinit(dof_handler,
+ mg_level.reinit(MappingQ1<dim>{},
+ dof_handler,
level_constraints,
quadrature_formula,
mg_additional_data);
mg_additional_data);
MatrixFree<dim, LevelNumberType> mg_level;
- mg_level.reinit(dof_handler,
+ mg_level.reinit(MappingQ1<dim>{},
+ dof_handler,
level_constraints,
quadrature_formula,
mg_additional_data);
const QGauss<1> quad(fe_degree + 1);
typename MatrixFree<dim, number>::AdditionalData data;
data.tasks_parallel_scheme = MatrixFree<dim, number>::AdditionalData::none;
- mf_data.reinit(std::vector<const DoFHandler<dim> *>({&dof2, &dof}),
+ mf_data.reinit(MappingQ1<dim>{},
+ std::vector<const DoFHandler<dim> *>({&dof2, &dof}),
std::vector<const AffineConstraints<double> *>(
{&constraints, &constraints}),
std::vector<Quadrature<1>>({quad}),
std::shared_ptr<MatrixFree<dim, double>> stokes_mf_storage(
new MatrixFree<dim, double>());
- stokes_mf_storage->reinit(stokes_dofs,
+ stokes_mf_storage->reinit(MappingQ1<dim>{},
+ stokes_dofs,
stokes_constraints,
QGauss<1>(degree_u + 1),
additional_data_stokes);
additional_data.mapping_update_flags =
update_values | update_JxW_values | update_quadrature_points;
std::shared_ptr<MatrixFree<dim, double>> data(new MatrixFree<dim, double>());
- data->reinit(dof_handlers_mf,
+ data->reinit(MappingQ1<dim>{},
+ dof_handlers_mf,
constraints_mf,
quadrature_formula_1d,
additional_data);