diff --git a/src/qmckl_distance.org b/src/qmckl_distance.org
index b34b939..5bdbb53 100644
--- a/src/qmckl_distance.org
+++ b/src/qmckl_distance.org
@@ -217,7 +217,7 @@ integer(c_int32_t) function test_qmckl_distance_sq(context) bind(C)
 	x =  (A(i,1)-B(j,1))**2 + &
 	     (A(i,2)-B(j,2))**2 + &
 	     (A(i,3)-B(j,3))**2
-	if ( dabs(1.d0 - C(i,j)/x) > 1.d-12 ) return
+	if ( dabs(1.d0 - C(i,j)/x) > 1.d-14 ) return
      end do
   end do
   test_qmckl_distance_sq = 0