Fix two bugs in matrix-free code

diff --git a/include/deal.II/matrix_free/fe_evaluation.h b/include/deal.II/matrix_free/fe_evaluation.h
index c6e2604582ed9e985e4f8338aaf0f16e5c25d7ba..b2c1182e2427e41ed5603f60f2c0634f19acda1f 100644 (file)
--- 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
index e7357c76bbb10bb9758f977b7ec3682e4a174d57..9d8d56822906c0501c2fb84acb7ac4422f747adb 100644 (file)
--- 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)));