diff --git a/src/filter/find_shocks.jl b/src/filter/find_shocks.jl
index a3232a44..05ee398b 100644
--- a/src/filter/find_shocks.jl
+++ b/src/filter/find_shocks.jl
@@ -347,7 +347,7 @@ function find_shocks(::Val{:LagrangeNewton},
     # println("Norm: $(ℒ.norm(x̂) / max(norm1,norm2))")
     # println(ℒ.norm(Δxλ))
     # println(ℒ.norm(x̂) / max(norm1,norm2) < tol && ℒ.norm(Δxλ) / ℒ.norm(xλ) < tol)
-    return x, ℒ.norm(x̂) / max(norm1,norm2) < 1e-14 && ℒ.norm(Δxλ) / ℒ.norm(xλ) < tol
+    return x, ℒ.norm(x̂) / max(norm1,norm2) < tol && ℒ.norm(Δxλ) / ℒ.norm(xλ) < sqrt(tol)
 end
 
 
diff --git a/test/test_3rd_order_estimation.jl b/test/test_3rd_order_estimation.jl
index 8f8beaaf..ab3e79fa 100644
--- a/test/test_3rd_order_estimation.jl
+++ b/test/test_3rd_order_estimation.jl
@@ -83,6 +83,8 @@ if isfinite(LLH)
 
     pt = Pigeons.pigeons(target = Caldara_lp, n_rounds = 0, n_chains = 1)
 else
+    pt = Pigeons.pigeons(target = Caldara_lp, n_rounds = 0, n_chains = 1)
+
     replica = pt.replicas[end]
     XMAX = deepcopy(replica.state)
     LPmax = Caldara_lp(XMAX)