“…One attractive feature of the method is that the computational complexity per iteration does not depend on the data size n, and thus it is directly scalable to large data volume, which is especially attractive in the era of big data. SGD type methods have found applications in several inverse problems, e.g., randomized Kaczmarz method [12,32] in computed tomography, ordered subset expectation maximization [13,21] for positron emission tomography, and more recently also some nonlinear inverse problems, e.g., optical tomography [4] and phonon transmission coefficient recovery [8]. However, the relevant mathematical theory for inverse problems in the lens of regularization theory [7,14,20] is still not fully understood.…”