From 7a812c9aad5c5707eaa6157a0327690782451c1e Mon Sep 17 00:00:00 2001
From: hartmann <hartmann@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Wed, 18 Jul 2001 09:38:49 +0000
Subject: [PATCH] Use relative error rather than absolute error for checking
 for an internal error.

git-svn-id: https://svn.dealii.org/trunk@4845 0785d39b-7218-0410-832d-ea1e28bc413d
---
 deal.II/deal.II/source/grid/tria_boundary_lib.cc | 8 +++++---
 1 file changed, 5 insertions(+), 3 deletions(-)

diff --git a/deal.II/deal.II/source/grid/tria_boundary_lib.cc b/deal.II/deal.II/source/grid/tria_boundary_lib.cc
index cadbfc8bad..5431ac6058 100644
--- a/deal.II/deal.II/source/grid/tria_boundary_lib.cc
+++ b/deal.II/deal.II/source/grid/tria_boundary_lib.cc
@@ -133,7 +133,7 @@ HyperBallBoundary<dim>::get_intermediate_points_between_points (
 		   v1=p1-center;
   const double length=std::sqrt((v1-v0).square());
   
-  double eps=1e-14;
+  double eps=1e-12;
   double r=0;
   if (compute_radius_automatically)
     {
@@ -143,8 +143,10 @@ HyperBallBoundary<dim>::get_intermediate_points_between_points (
   else
     r=radius;
 
-  Assert(std::fabs(v0.square()-r*r)<eps, ExcInternalError());
-  Assert(std::fabs(v1.square()-r*r)<eps, ExcInternalError());
+  
+  const double r2=r*r;
+  Assert(std::fabs(v0.square()-r2)<eps*r2, ExcInternalError());
+  Assert(std::fabs(v1.square()-r2)<eps*r2, ExcInternalError());
   
   const double alpha=std::acos((v0*v1)/std::sqrt(v0.square()*v1.square()));
   const double d_alpha=alpha/(n+1);
-- 
2.39.5