From 0f249d98b84ceceaa88e797c0fa80c775a2985d4 Mon Sep 17 00:00:00 2001 From: Magdalena Schreter Date: Wed, 5 Jul 2023 11:53:55 +0200 Subject: [PATCH] fix ComponentMask --- source/navier_stokes.cc | 2 +- source/navier_stokes_matrix.cc | 2 +- source/navier_stokes_preconditioner.cc | 6 +++--- 3 files changed, 5 insertions(+), 5 deletions(-) diff --git a/source/navier_stokes.cc b/source/navier_stokes.cc index 7234090..a33e898 100644 --- a/source/navier_stokes.cc +++ b/source/navier_stokes.cc @@ -1415,7 +1415,7 @@ NavierStokes::interpolate_pressure_field( // set DG0 components to zero std::vector> constant_modes; DoFTools::extract_constant_modes(dof_handler_p, - std::vector(1, true), + ComponentMask(std::vector(1, true)), constant_modes); AssertDimension(constant_modes.size(), 2); for (unsigned int i = 0; i < pressure_vector.locally_owned_size(); ++i) diff --git a/source/navier_stokes_matrix.cc b/source/navier_stokes_matrix.cc index 82b0c99..698ddc1 100644 --- a/source/navier_stokes_matrix.cc +++ b/source/navier_stokes_matrix.cc @@ -132,7 +132,7 @@ NavierStokesMatrix::initialize(const MatrixFree &matrix_free_in, std::vector> constant_modes; DoFTools::extract_constant_modes(matrix_free->get_dof_handler(dof_index_p), - std::vector(1, true), + ComponentMask(std::vector(1, true)), constant_modes); for (unsigned int mode = 0; mode < 2; ++mode) diff --git a/source/navier_stokes_preconditioner.cc b/source/navier_stokes_preconditioner.cc index 5a3e035..168e24e 100644 --- a/source/navier_stokes_preconditioner.cc +++ b/source/navier_stokes_preconditioner.cc @@ -1152,7 +1152,7 @@ NavierStokesPreconditioner::initialize_matrices( { std::vector> constant_modes; DoFTools::extract_constant_modes(dof_handler_p, - std::vector(1, true), + ComponentMask(std::vector(1, true)), constant_modes); AssertDimension(constant_modes.size(), 2); int min_local_index = std::numeric_limits::min(); @@ -1387,12 +1387,12 @@ NavierStokesPreconditioner::initialize_matrices( if (parameters.precondition_velocity == FlowParameters::u_amg || parameters.precondition_velocity == FlowParameters::u_amg_linear) DoFTools::extract_constant_modes(dof_handler_u, - std::vector(dim, true), + ComponentMask(std::vector(dim, true)), constant_modes_u); if (parameters.density > 0 || parameters.augmented_taylor_hood) DoFTools::extract_constant_modes(dof_handler_p, - std::vector(1, true), + ComponentMask(std::vector(1, true)), constant_modes_p);