The required time and memory for gravity data inversion can be decreased by priorconditioned lanczos bidiagonalization. The orthonormal basis vectors of the discrete cosine transform can be applied for solution space priorconditioning. We performed a number of inversions of synthetic and real cases, the proposed method works for fast inversion of the data.