]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Fix the filtering of close-to-zero elements in interface_constraints
authorMartin Kronbichler <kronbichler@lnm.mw.tum.de>
Sun, 19 Jul 2015 21:13:09 +0000 (23:13 +0200)
committerMartin Kronbichler <kronbichler@lnm.mw.tum.de>
Mon, 20 Jul 2015 13:51:33 +0000 (15:51 +0200)
source/fe/fe_q_base.cc

index 936e9ecf970dff818dcf5f4b331efa023f32ac49..77839adcade303d3b50a92c5580dcdc593b7289c 100644 (file)
@@ -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;
           }
       }
   }

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.