In summary, we have implemented the combination of finite difference method and conjugate gradient method in the positron annihilation lifetime calculations. From the perfect and defect calculations, the analysis shows that the positron wave functions could be optimized by conjugate gradient method efficiently without loss of accuracy. From the tests on mesh points and supercell sizes, the numerical results show that our calculations are extremely stable for different mesh points and supercells.