In this paper, the preconditioned global least squares algorithm is applied to ap-
proximate the generalized inverse of nearly singular or rectangular matrices. This
preconditioner is based on the C-orthogonalization process, where C a symmetric pos-
itive definite matrix. Finally, some numerical experiments are given to illustrate the
efficiency of the new preconditioner.