From a0d4e177a6b1e901fc9c18fa3bbd689bedc9d381 Mon Sep 17 00:00:00 2001 From: Martin Kronbichler Date: Sun, 19 Jul 2015 23:13:09 +0200 Subject: [PATCH] Fix the filtering of close-to-zero elements in interface_constraints --- source/fe/fe_q_base.cc | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/source/fe/fe_q_base.cc b/source/fe/fe_q_base.cc index 936e9ecf97..77839adcad 100644 --- a/source/fe/fe_q_base.cc +++ b/source/fe/fe_q_base.cc @@ -214,8 +214,8 @@ struct FE_Q_Base::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::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; } } } -- 2.39.5