From: Martin Kronbichler Date: Sun, 19 Jul 2015 21:13:09 +0000 (+0200) Subject: Fix the filtering of close-to-zero elements in interface_constraints X-Git-Tag: v8.4.0-rc2~741^2~2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=a0d4e177a6b1e901fc9c18fa3bbd689bedc9d381;p=dealii.git 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 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; } } }