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)