A fundamental imaging problem in microstructural analysis of metals is the reconstruction of local crystallographic orientations from X-ray diffraction measurements. This work deals with the computation of the 3D orientation distribution function for individual grains of the material in consideration. We present an iterative large-scale algorithm that uses preconditioned regularizing CGLS iterations with a stopping criterion based on the information available in the residual vectors.