diff --git a/test/test.jl b/test/test.jl
index 1d598d1c..5ce51c31 100644
--- a/test/test.jl
+++ b/test/test.jl
@@ -462,7 +462,7 @@ end
     random = Metida.VarEffect(Metida.@covstr(1 + r2 * r1|subject), Metida.DIAG; coding=Dict(:r1 => DummyCoding(), :r2 => DummyCoding()))
     )
     Metida.fit!(lmm)
-    @test Metida.theta(lmm)  ≈ [2.796694409004289, 2.900485570555582, 3.354913215348968, 2.0436114769223237, 1.8477830405766895, 2.0436115732330955, 1.0131934233937254] atol=1E-6 # atol=1E-8 !
+    @test Metida.theta(lmm)  ≈ [2.796694409004289, 2.900485570555582, 3.354913215348968, 2.0436114769223237, 1.8477830405766895, 2.0436115732330955, 1.0131934233937254] atol=1E-5 # atol=1E-8 !
     @test Metida.m2logreml(lmm)  ≈ 713.0655862252027 atol=1E-8
 end
 @testset "  Model: &, DIAG/SI                                        " begin
@@ -470,7 +470,7 @@ end
     random = Metida.VarEffect(Metida.@covstr(r1&r2|subject), Metida.DIAG),
     )
     Metida.fit!(lmm)
-    @test Metida.theta(lmm)  ≈ [3.0325005960015985, 3.343826588448401, 1.8477830405766895, 1.8477830405766895, 1.8477830405766895, 4.462942536844632, 1.0082345219318216] atol=1E-6 # atol=1E-8 !
+    @test Metida.theta(lmm)  ≈ [3.0325005960015985, 3.343826588448401, 1.8477830405766895, 1.8477830405766895, 1.8477830405766895, 4.462942536844632, 1.0082345219318216] atol=1E-5 # atol=1E-8 !
     @test Metida.m2logreml(lmm)  ≈ 719.9413776641368 atol=1E-8
 end
 @testset "  Model: INT, +,  TOEPHP(3)/SI                             " begin