def test_shepp_logan_L1WaveletRecon(self): img, mps, ksp = self.shepp_logan_setup() lamda = 0 img_rec = app.L1WaveletRecon(ksp, mps, lamda, alg_name='GradientMethod').run() npt.assert_allclose(img, img_rec, atol=1e-3, rtol=1e-3) img_rec = app.L1WaveletRecon(ksp, mps, lamda, alg_name='PrimalDualHybridGradient', max_iter=1000).run() npt.assert_allclose(img, img_rec, atol=1e-3, rtol=1e-3)
def test_shepp_logan_L1WaveletRecon(self): img, mps, ksp = self.shepp_logan_setup() lamda = 0 for solver in ['GradientMethod', 'PrimalDualHybridGradient', 'ADMM']: with self.subTest(solver=solver): img_rec = app.L1WaveletRecon(ksp, mps, lamda, solver=solver, show_pbar=False).run() npt.assert_allclose(img, img_rec, atol=1e-2, rtol=1e-2)