def METHOD_NAME(self):
    target, sens = self.problem_2d(8)
    # makes dim*dim*2 trajectory
    traj = sp.mri.radial(
        (sens.shape[1], sens.shape[1], 2),
        target.shape,
        golden=True,
        dtype=np.float32,
    )
    # reshape to be Nt*2 trajectory
    traj = np.reshape(traj, [traj.shape[0] * traj.shape[1], 2])
    A = linop.Sense(sens, coord=traj, weights=None, ishape=target.shape).H
    pulses = rf.stspa(
        target,
        sens,
        traj,
        dt=4e-6,
        alpha=1,
        b0=None,
        st=None,
        explicit=False,
        max_iter=100,
        tol=1e-4,
    )
    npt.assert_array_almost_equal(A * pulses, target, 1e-3)