X-ray computed tomography is a commonly used technique for noninvasive imaging at synchrotron facilities. Iterative tomographic reconstruction algorithms are often preferred for recovering high quality 3D volumetric images from 2D X-ray images, however, their use has been limited to small/medium datasets due to their computational requirements. In this paper, we propose a high-performance iterative reconstruction system for terabyte(s)-scale 3D volumes. Our design involves three novel optimizations: (1) optimization of (back)projection operators by extending the 2D memory-centric approach to 3D; (2) performing hierarchical communications by exploiting "fat-node" architecture with many GPUs; (3) utilization of mixed-precision types while preserving convergence rate and quality. We extensively evaluate the proposed optimizations and scaling on the Summit supercomputer. Our largest reconstruction is a mouse brain volume with 9Kx11Kx11K voxels, where the total reconstruction time is under three minutes using 24,576 GPUs, reaching 65 PFLOPS: 34% of Summit's peak performance.
@article{arxiv.2009.07226,
title = {Petascale XCT: 3D Image Reconstruction with Hierarchical Communications on Multi-GPU Nodes},
author = {Mert Hidayetoglu and Tekin Bicer and Simon Garcia de Gonzalo and Bin Ren and Vincent De Andrade and Doga Gursoy and Raj Kettimuthu and Ian T. Foster and Wen-mei W. Hwu},
journal= {arXiv preprint arXiv:2009.07226},
year = {2020}
}