fix rv determination that would sometimes crash

AWehrhahn · Nov 26, 2021 · 986be74 · 986be74
1 parent 1c8d178
commit 986be74
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/src/pysme/continuum_and_radial_velocity.py b/src/pysme/continuum_and_radial_velocity.py
@@ -948,7 +948,9 @@ def func(rv):
             return resid
 
         interpolator = lambda x: np.interp(x, x_syn, y_syn)
-        res = least_squares(func, x0=rvel, loss="soft_l1", bounds=rv_bounds)
+        res = least_squares(
+            func, x0=rvel, loss="soft_l1", bounds=rv_bounds, jac="3-point"
+        )
         rvel = res.x[0]
 
     return rvel