diff --git a/lisainstrument/shift_inversion_dask.py b/lisainstrument/shift_inversion_dask.py
index afac2ad0afe786353623fc80dec18003ad39ceb6..f41f22686358109e301e0ea50241c921e55fccf0 100644
--- a/lisainstrument/shift_inversion_dask.py
+++ b/lisainstrument/shift_inversion_dask.py
@@ -17,6 +17,7 @@ from lisainstrument.regular_interpolators import (
     RegularInterpolator,
     make_regular_interpolator_lagrange,
 )
+from lisainstrument.dynamic_delay_dsp import make_regular_interpolator_dsp
 from lisainstrument.shift_inversion_numpy import fixed_point_iter
 
 
@@ -145,3 +146,23 @@ def make_shift_inverse_lagrange_dask(
     """
     interp = make_regular_interpolator_lagrange(order)
     return ShiftInverseDask(max_abs_shift, interp, max_iter, tolerance)
+
+def make_shift_inverse_dsp_dask(
+    order: int,
+    max_abs_shift: float,
+    max_iter: int,
+    tolerance: float,
+) -> ShiftInverseDask:
+    """Set up ShiftInverseDask instance using dsp.timeshift as interpolation method.
+
+    Arguments:
+        order: Order of the Lagrange polynomials
+        max_abs_shift: Upper limit for absolute difference between coordinate frames w.r.t index space
+        max_iter: Maximum iterations before fail
+        tolerance: Maximum absolute error of result
+
+    Returns:
+        Inversion function of type ShiftInverseNumpy
+    """
+    interp = make_regular_interpolator_dsp(order)
+    return ShiftInverseDask(max_abs_shift, interp, max_iter, tolerance)