From: Martin Kronbichler <kronbichler@lnm.mw.tum.de>
Date: Sun, 5 Jan 2020 16:37:55 +0000 (+0100)
Subject: Fix two bugs in matrix-free code
X-Git-Tag: v9.2.0-rc1~741^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=2f3a0b99aba4e5bc6dfbceac14da8b2ffc24a564;p=dealii.git

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 c6e2604582..b2c1182e24 100644
--- 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 e7357c76bb..9d8d568229 100644
--- 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)));