[cig-commits] commit 2338 by gassmoeller to /var/svn/dealii/aspect

[dealii.demon at gmail.com]

Revision 2338

Removed debug code and loosened Assert condition in gplates plugin to hold true for coarse meshes.

U   trunk/aspect/source/velocity_boundary_conditions/gplates.cc


http://www.dealii.org/websvn/revision.php?repname=Aspect+Repository&path=%2F&rev=2338&peg=2338

Diff:
Modified: trunk/aspect/source/velocity_boundary_conditions/gplates.cc
===================================================================
--- trunk/aspect/source/velocity_boundary_conditions/gplates.cc	2014-03-19 14:30:43 UTC (rev 2337)
+++ trunk/aspect/source/velocity_boundary_conditions/gplates.cc	2014-03-19 14:52:37 UTC (rev 2338)
@@ -263,13 +263,7 @@
                             // visited all possible interpolation points. Return current surf_vel.
                             const double residual_normal_velocity = (surf_vel * position) / (surf_vel.norm() * position.norm());
 
-                            // TODO: Debug hack, can be removed
-                            if (!(residual_normal_velocity < 1e-12))
-                              {
-                                std::cout << residual_normal_velocity;
-                              }
-
-                            AssertThrow(residual_normal_velocity < 1e-12,
+                            AssertThrow(residual_normal_velocity < 1e-11,
                                         ExcMessage("Error in velocity boundary module interpolation. "
                                                    "Radial component of velocity should be zero, but is not."));
 
@@ -297,13 +291,7 @@
 
         const double residual_normal_velocity = (surf_vel * position) / (surf_vel.norm() * position.norm());
 
-        //TODO: Debug hack. Can be removed
-        if (!(residual_normal_velocity < 1e-12))
-          {
-            std::cout << residual_normal_velocity;
-          }
-
-        AssertThrow( residual_normal_velocity < 1e-12,
+        AssertThrow( residual_normal_velocity < 1e-11,
                      ExcMessage("Error in velocity boundary module interpolation. "
                                 "Radial component of velocity should be zero, but is not."));