Merge pull request #216 from seismic-anisotropy/fix-lsoda-tolerances

Improve default LSODA tolerances
seismic-anisotropy · Sep 16, 2024 · 56940fa · 56940fa
2 parents d4fb04b + 4b62dbf
commit 56940fa
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/pydrex/minerals.py b/src/pydrex/minerals.py
@@ -520,7 +520,7 @@ def perform_step(solver):
             time_start,
             y_start,
             time_end,
-            atol=kwargs.pop("atol", np.abs(y_start * 1e-6) + 1e-12),
+            atol=kwargs.pop("atol", np.abs(y_start * 1e-6) + 1e-4),
             rtol=kwargs.pop("rtol", 1e-6),
             first_step=kwargs.pop("first_step", np.abs(time_end - time_start) * 1e-1),
             # max_step=kwargs.pop("max_step", np.abs(time_end - time_start)),