Set minimum allowed skew to 1e-12

diff --git a/reference.c b/reference.c
index 19cb8adbd7bab82744a2e9ca647e055ee4a2a7ea..01b53ef599dd18e8886c87af156c04dde2a55c72 100644 (file)
--- a/reference.c
+++ b/reference.c
@@ -541,10 +541,10 @@ REF_SetReference(int stratum,
   assert(initialised);
 
   /* Avoid getting NaNs */
-  if (skew == 0.0)
-    skew = 1e-10;
-  if (our_skew == 0.0)
-    our_skew = 1e-10;
+  if (skew < 1e-12)
+    skew = 1e-12;
+  if (our_skew < 1e-12)
+    our_skew = 1e-12;
 
   /* If we get a serious rounding error in the source stats regression
      processing, there is a remote chance that the skew argument is a