Papers › Decomposed Diffusion Sampler for Accelerating Large-Scale Inverse Problems

Decomposed Diffusion Sampler for Accelerating Large-Scale Inverse Problems

10 Mar 2023arXiv:2303.05754archive 2025-07-28

Hyungjin Chung, Suhyeon Lee, Jong Chul Ye

Krylov subspace, which is generated by multiplying a given vector by the matrix of a linear transformation and its successive powers, has been extensively studied in classical optimization literature to design algorithms that converge quickly for large linear inverse problems. For example, the conjugate gradient method (CG), one of the most popular Krylov subspace methods, is based on the idea of minimizing the residual error in the Krylov subspace. However, with the recent advancement of high-performance diffusion solvers for inverse problems, it is not clear how classical wisdom can be synergistically combined with modern diffusion models. In this study, we propose a novel and efficient diffusion sampling strategy that synergistically combines the diffusion sampling and Krylov subspace methods. Specifically, we prove that if the tangent space at a denoised sample by Tweedie's formula forms a Krylov subspace, then the CG initialized with the denoised data ensures the data consistency update to remain in the tangent space. This negates the need to compute the manifold-constrained gradient (MCG), leading to a more efficient diffusion sampling method. Our method is applicable regardless of the parametrization and setting (i.e., VE, VP). Notably, we achieve state-of-the-art reconstruction quality on challenging real-world medical inverse imaging problems, including multi-coil MRI reconstruction and 3D CT reconstruction. Moreover, our proposed method achieves more than 80 times faster inference time than the previous state-of-the-art method. Code is available at https://github.com/HJ-harry/DDS

PaperPDFCodeCode Syntology ran

In Syntology Open this paper in Syntology's Atlas, the map of the papers in Syntology's graph and their citations.

For agents, Syntology's MCP tool lists every function and class Syntology harvested from this paper and whether it ran (how to connect): get_harvested_code_for_paper(arxiv_id="2303.05754")

Code

Syntology Ran 16 of 27 code samples harvested from 1 repository linked to this paper; 11 have no recorded run. Of those that ran: 3 ran · honoured contract; 1 ran · violated contract; 6 ran · our draft was wrong; 3 ran · fixture could not drive it; 3 ran with no contract checked.

By repository: official repository: 26 samples from 1 repository, 15 ran; 1 identical to code first harvested elsewhere. The run record, sample by sample. “Ran” means executed on a synthesized input, not that the code is correct or reproduces the paper.

hj-harry/dds officialmentioned in papermentioned on GitHubpytorch report
hj-harry/ddip3d mentioned on GitHubpytorch report

Repository list and official/mentioned flags are the archive's, frozen 2025-07-28. Reachability, where shown, is from one Syntology probe window (2026-09-16 to 2026-09-18); repositories not probed show nothing. GitHub stars are not tracked.

Code Syntology ran Syntology

27 samples harvested; 16 ran; 3 honoured the contract we drafted; 11 have no recorded run. Read from Syntology's graph 2026-09-24; that is when this build read the record, not when the samples ran.

3ran · honoured contract
1ran · violated contract
6ran · our draft was wrong
3ran · fixture could not drive it
3ran
11unverified

Licence: 27 of the 27 samples are pointer only, meaning Syntology does not serve that copy's text. This page shows no code text for any sample; each one links to its file in the repository.

Harvested from hj-harry/dds. Some samples are identical code Syntology first harvested from another repository; for those, this paper's copy is not located and its licence is not recorded. “Ran” means the sample executed on a synthesized input. It does not mean the output is correct, and nothing here reproduces the paper's results. “Honoured” and “violated” refer to a contract Syntology drafted from the code itself; “our draft was wrong” and “fixture could not drive it” are failures of Syntology's instrument, not of the code.

Each sample ends with its code_sha256, Syntology's identity for that exact code. An agent fetches the stored sample with Syntology's MCP tool get_code(code_sha256="…") (how to connect); click an identity to copy that call.

Repository labels, per sample. official repository: The archive marks this repository official for the paper. named in the paper: The archive records that the paper mentions this repository; it is not marked official. community (archive-listed): In the archive's code links for this paper, not marked official and not recorded as mentioned in the paper. found in paper text by Syntology: Syntology found this repository in the paper's own text; whether it is the authors' implementation is not asserted. community: Not in the archive's code links for this paper; a community repository Syntology harvested. Samples from a repository marked official are listed first. Licence labels name the repository's licence as recorded at harvest. “Pointer only” means Syntology does not serve that copy's text, for one of four reasons: no licence file was found; the licence was not identified; the licence is recorded as permissive but that copy's record is not marked cleared; or the licence is outside the permissive list Syntology serves text under (MIT, Apache-2.0, BSD and similar). Some licences outside that list permit redistribution, such as WTFPL, and GPL-3.0 under its conditions; they are simply not on the list. Hover a licence label for the reason. File links open the file on GitHub at the default branch, which may have changed since the harvest.

MulticoilMRI hj-harry/dds/solver_2d.py official repository ran no licence file found · pointer only · 0b34c7d318e99dda · report
QKVAttention hj-harry/dds/solver_2d.py official repository ran no licence file found · pointer only · 4cea67dfff70dfe2 · report
SinglecoilMRI_comp hj-harry/dds/solver_2d.py official repository ran no licence file found · pointer only · 5258121713055e95 · report
clear hj-harry/dds/solver_2d.py official repository ran · honoured contract fingerprinted no licence file found · pointer only · 2c52eaf728ba40f9 · report
fft2c_new hj-harry/dds/solver_2d.py official repository ran · fixture could not drive it no licence file found · pointer only · 4285b3d174aef9e0 · report
fftshift hj-harry/dds/solver_2d.py official repository ran · our draft was wrong fingerprinted no licence file found · pointer only · b41fc4b36259eb57 · report
get_mask hj-harry/dds/solver_2d.py official repository ran · our draft was wrong no licence file found · pointer only · b4db25e6c5279fcc · report
ifft2_m hj-harry/dds/solver_2d.py official repository ran · our draft was wrong no licence file found · pointer only · 55164bfcc9459982 · report
ifft2c_new hj-harry/dds/solver_2d.py official repository ran · fixture could not drive it no licence file found · pointer only · fe5581b42c91e639 · report
ifftshift hj-harry/dds/solver_2d.py official repository ran · our draft was wrong fingerprinted no licence file found · pointer only · 9869a096979dafe0 · report
nchw_comp_to_real hj-harry/dds/solver_2d.py official repository ran · fixture could not drive it no licence file found · pointer only · f4750dccd85a13f3 · report
normalize_np hj-harry/dds/solver_2d.py official repository ran · honoured contract fingerprinted no licence file found · pointer only · 135303c090dee252 · report
real_to_nchw_comp hj-harry/dds/solver_2d.py official repository ran · violated contract fingerprinted no licence file found · pointer only · 8d87d11aa5f9cde2 · report
roll hj-harry/dds/solver_2d.py official repository ran · our draft was wrong no licence file found · pointer only · 96ca075a8ebe94a7 · report
roll_one_dim hj-harry/dds/solver_2d.py official repository ran · our draft was wrong no licence file found · pointer only · d61fe2184194de9f · report
AttentionBlock hj-harry/dds/solver_2d.py official repository unverified no licence file found · pointer only · 15b10af89289e00c · report
CG hj-harry/dds/solver_2d.py official repository unverified no licence file found · pointer only · abc9b70e3146ed20 · report
CT hj-harry/dds/solver_2d.py official repository unverified no licence file found · pointer only · 0fe729cb5016e5b1 · report
Diffusion hj-harry/dds/solver_2d.py official repository unverified no licence file found · pointer only · 14da7a4d709b89ba · report
Downsample hj-harry/dds/solver_2d.py official repository unverified no licence file found · pointer only · 72f641bded2a58d5 · report
QKVAttentionLegacy hj-harry/dds/solver_2d.py official repository unverified no licence file found · pointer only · 6503094226d94604 · report
ResBlock hj-harry/dds/solver_2d.py official repository unverified no licence file found · pointer only · f49a281fab7d9cff · report
UNetModel hj-harry/dds/solver_2d.py official repository unverified no licence file found · pointer only · 4fa4885d061af603 · report
Upsample hj-harry/dds/solver_2d.py official repository unverified no licence file found · pointer only · 6e246e792b44ae83 · report
create_model hj-harry/dds/solver_2d.py official repository unverified no licence file found · pointer only · 22a4f88b2aa37adf · report
fft2_m hj-harry/dds/solver_2d.py official repository unverified no licence file found · pointer only · 9c0580953d54c536 · report
get_beta_schedule identical code first harvested elsewhere ran · honoured contract licence of this copy not recorded · ca0cea1ca1eeec89 · report

Tasks

CT ReconstructionMRI Reconstructionsubspace methods

Results from the paper archive 2025-07-28

No leaderboard rows for this paper in the archive.

Methods

Diffusion

Report a problem or propose a change · a person checks every report against the paper or source before anything changes; decisions are listed on /corrections