The authors present a stabilized approximate inverse algorithm for arbitrary symmetric positive definite matrices to be used in preconditioned conjugate gradient methods. They also investigate another approach to prevent breakdowns that is based on the technique of diagonally compensated reduction of positive off-diagonal entries. Numerical results for problems arising from finite element discretizations of elasticity and diffusion problems illustrate the performance of the algorithms.