Fix two bugs in matrix-free code

dealii · Jan 5, 2020 · 2f3a0b9 · 2f3a0b9
1 parent d7e237b
commit 2f3a0b9
Show file tree

Hide file tree

Showing 2 changed files with 5 additions and 3 deletions.
diff --git a/include/deal.II/matrix_free/fe_evaluation.h b/include/deal.II/matrix_free/fe_evaluation.h
@@ -7417,11 +7417,11 @@ FEFaceEvaluation<dim,
                             .normal_vectors[offsets];
   this->jacobian = &this->matrix_info->get_mapping_info()
                       .face_data_by_cells[this->quad_no]
-                      .jacobians[!this->is_minus_face][offsets];
+                      .jacobians[!this->is_interior_face][offsets];
   this->normal_x_jacobian =
     &this->matrix_info->get_mapping_info()
        .face_data_by_cells[this->quad_no]
-       .normals_times_jacobians[!this->is_minus_face][offsets];
+       .normals_times_jacobians[!this->is_interior_face][offsets];
 
 #  ifdef DEBUG
   this->dof_values_initialized     = false;

diff --git a/include/deal.II/matrix_free/mapping_info.templates.h b/include/deal.II/matrix_free/mapping_info.templates.h
@@ -2049,7 +2049,9 @@ namespace internal
                   fe_val.reinit(cell_it, face);
 
                   const bool is_local =
-                    cell_it->is_locally_owned() &&
+                    (cell_it->active() ?
+                       cell_it->is_locally_owned() :
+                       cell_it->is_locally_owned_on_level()) &&
                     (!cell_it->at_boundary(face) ||
                      (cell_it->at_boundary(face) &&
                       cell_it->has_periodic_neighbor(face)));