diff --git a/src/pysme/continuum_and_radial_velocity.py b/src/pysme/continuum_and_radial_velocity.py
index eeae6c3c..c9a2dc15 100644
--- 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