In this take we explian how to use the conjugate gradient least squares (CGLS) algorithm in block forms for ill-condiioned linear system. We explore the applicaiton of CGLS so that multiple right-hand sides can be treated simultaneously, while maintaining orthogonality and minimization properties along iterations.  Numerical examples demonstrate effectiveness in handling rank deficiency and efficiency in convergence accelerations in these forms.