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