Fix the filtering of close-to-zero elements in interface_constraints

diff --git a/source/fe/fe_q_base.cc b/source/fe/fe_q_base.cc
index 936e9ecf970dff818dcf5f4b331efa023f32ac49..77839adcade303d3b50a92c5580dcdc593b7289c 100644 (file)
--- a/source/fe/fe_q_base.cc
+++ b/source/fe/fe_q_base.cc
@@ -214,8 +214,8 @@ struct FE_Q_Base<POLY,xdim,xspacedim>::Implementation
           // to avoid unwanted fill-in of the constraint matrices (which would
           // then increase the number of other DoFs a constrained DoF would
           // couple to)
-          if (std::fabs(fe.interface_constraints(i,j)) < 1e-14)
-            fe.interface_constraints(i,j) = 0;
+          if (std::fabs(fe.interface_constraints(i,face_index_map[j])) < 1e-13)
+            fe.interface_constraints(i,face_index_map[j]) = 0;
         }
   }
 
@@ -390,8 +390,8 @@ struct FE_Q_Base<POLY,xdim,xspacedim>::Implementation
             // zero to avoid unwanted fill-in of the constraint matrices
             // (which would then increase the number of other DoFs a
             // constrained DoF would couple to)
-            if (std::fabs(fe.interface_constraints(i,j)) < 1e-14)
-              fe.interface_constraints(i,j) = 0;
+            if (std::fabs(fe.interface_constraints(i,face_index_map[j])) < 1e-13)
+              fe.interface_constraints(i,face_index_map[j]) = 0;
           }
       }
   }