def METHOD_NAME(self):
    img_shape = [16, 16]
    mps_shape = [8, 16, 16]
    img = sp.randn(img_shape, dtype=complex)
    mps = sp.randn(mps_shape, dtype=complex)
    for coil_batch_size in [None, 1, 2, 3]:
        A = linop.Sense(mps, coil_batch_size=coil_batch_size)
        check_linop_adjoint(A, dtype=complex)
        npt.assert_allclose(sp.fft(img * mps, axes=[-1, -2]), A * img)