Cryo-electron microscopy single particle analysis (SPA) rep-resents a vital tool for structure determination of macro-molecules. Discrete inverse problems arising in this field are extremely large and seriously contaminated by noise.
The model matrix is highly structured and can not be stored ex-plicitly. Iterative regularization methods are used here only rarely, since it is believed that they are computationally ex-pensive and require complicated stopping criteria.
In this paper, we overcome these difficulties and demonstrate that Hybrid Krylov subspace methods can be used to solve SPA inverse problems efficiently. We propose an application-driven regularization parameter selection approach and present a matrix-free implementation of the hybrid solver for GPU com-putations in single precision arithmetic.
In comparison to Fourier-based techniques, the hybrid approach allows to com-pute reconstructions from full (nonreduced) SPA models even for highly noisy data sets.(c) 2022 Elsevier Inc. All rights reserved.