SART algorithm stuck in "set_w" when the input image is large
- Dominant language
- MATLAB
- Stars
- 803
- Forks
- 262
- Avg merge
- 9d 13h
- Merged PRs (30d)
- 8
Description
## Behavior
When I ran iterative methods like SART on a CT scan (512×512×463, can be downloaded [here](https://drive.google.com/file/d/1BGGxrHR984JRAF0AYlLAWksPivCU9hKw/view?usp=drive_link)), it was stuck in "set_w" (w=Ax(...)) when initializing IterativeReconAlg class ([link to the code](https://github.com/CERN/TIGRE/blob/4598aaa32a9aa3ae31aef45f5fe2f5daa01854ea/Python/tigre/algorithms/iterative_recon_alg.py#L222)). Do you know the possible reason and how to fix it?
Smaller images (like 256×256×256) work well. Running FDK on the 512×512×463 CT also works well.
## Code to reproduce the problem (If applicable)
```
image_pred = algs.sart(projections, geo, angles, niter=4)
```
where "projections, geo, angles" can be obtained from the CT pickle file I attached above.
## Specifications
- MATLAB/python version: python 3.9
- OS: Linux
- CUDA version: 11.3
(Data reference: Zha et. al, [NAF](https://github.com/Ruyi-Zha/naf_cbct))
Thanks!
Contributor guide
Assessment
This issue has not been assessed yet.