From e7be50b3af43ab9a57c6c70651ed37a102642735 Mon Sep 17 00:00:00 2001
From: d3v-null <dev.null@curtin.edu.au>
Date: Wed, 7 Aug 2024 08:50:56 +0000
Subject: [PATCH] bump up tolerance

---
 src/van_vleck.rs | 8 ++++----
 1 file changed, 4 insertions(+), 4 deletions(-)

diff --git a/src/van_vleck.rs b/src/van_vleck.rs
index f2db675..17f34eb 100644
--- a/src/van_vleck.rs
+++ b/src/van_vleck.rs
@@ -338,8 +338,8 @@ pub fn correct_van_vleck(
 
 /// Use Newton's method to solve the inverse of `sighat_vector`.
 fn van_vleck_auto(s: f64) -> Option<f64> {
-    let tol = 1e-10;
-    let niter = 10;
+    let tol = 1e-12;
+    let niter = 100;
     let mut guess = s;
     let mut delta = sighat(guess) - s;
     let mut count = 0;
@@ -728,8 +728,8 @@ fn van_vleck_cross_int(khat: f64, sigma_x: f64, sigma_y: f64) -> Option<f64> {
     let x_ = (0..7).map(|i| (i as f64 + 0.5) / sigma_x).collect_vec();
     let y_ = (0..7).map(|i| (i as f64 + 0.5) / sigma_y).collect_vec();
 
-    let tol = 1e-10; // it's 1e-10 for the autos
-    let niter = 10;
+    let tol = 1e-12; // it's 1e-10 for the autos
+    let niter = 100;
     let mut guess = khat / (sigma_x * sigma_y);
     debug_assert!(guess >= 0.0, "|ρ| must be >= 0: {guess:?}");
     debug_assert!(guess < 1.0, "|ρ| must be < 1: {guess:?}");