An iterative perturbation theory-based algorithm for calculation of diagonal force constants of the  modes in quasi-tetrahedral systems