diff --git a/test/test_operators/test_adjoint.py b/test/test_operators/test_adjoint.py
index 2e8bcf271940081e2ec8b033d6558049b9243197..c6fa4d2eb5dd38eec2488b7f9efea33fbd13554b 100644
--- a/test/test_operators/test_adjoint.py
+++ b/test/test_operators/test_adjoint.py
@@ -35,6 +35,15 @@ _pow_spaces = [ift.PowerSpace(ift.RGSpace((17, 38), harmonic=True))]
 
 
 class Consistency_Tests(unittest.TestCase):
+    @expand(product(_p_RG_spaces, [np.float64, np.complex128]))
+    def testLOSResponse(self, sp, dtype):
+        starts = np.random.randn(len(sp.shape), 10)
+        ends = np.random.randn(len(sp.shape), 10)
+        sigma_low = 1e-4*np.random.randn(10)
+        sigma_ups = 1e-5*np.random.randn(10)
+        op = ift.LOSResponse(sp, starts, ends, sigma_low, sigma_ups)
+        ift.extra.consistency_check(op, dtype, dtype)
+
     @expand(product(_h_spaces + _p_spaces + _pow_spaces,
                     [np.float64, np.complex128]))
     def testOperatorCombinations(self, sp, dtype):