69       SUBROUTINE ab_sint_os(v, nmax, zetp, zetq, zetw, rho, rac2, omega, r_cutoff)
73       REAL(KIND=dp), INTENT(IN)                          :: zetp, zetq, zetw, rho, rac2, omega, &  local
129                                                             fcy, fcz, rho, zetp, zetq, zetw  local
164             zetw = 1.0_dp/(zeta(ipgf) + zetc(jpgf))
166             rho = zeta(ipgf)*zetc(jpgf)*zetw
170             CALL cps_operator2(v, nmax, zetp, zetq, zetw, rho, rac2, omega, r_cutoff)
179                rcw(:) = -zeta(ipgf)*zetw*rac(:)
247                f5 = 0.5_dp*zetw
249                raw(:) = zetc(jpgf)*zetw*rac(:)
440    SUBROUTINE cps_coulomb2(v, nmax, zetp, zetq, zetw, rho, rac2, omega, r_cutoff)  argument
443       REAL(KIND=dp), INTENT(IN)                          :: zetp, zetq, zetw, rho, rac2, omega, &  local
456       f0 = 2.0_dp*SQRT(pi**5*zetw)*zetp*zetq
485    SUBROUTINE cps_verf2(v, nmax, zetp, zetq, zetw, rho, rac2, omega, r_cutoff)  argument
488       REAL(KIND=dp), INTENT(IN)                          :: zetp, zetq, zetw, rho, rac2, omega, &  local
501       f0 = 2.0_dp*SQRT(pi**5*zetw*comega)*zetp*zetq
532    SUBROUTINE cps_verfc2(v, nmax, zetp, zetq, zetw, rho, rac2, omega, r_cutoff)  argument
535       REAL(KIND=dp), INTENT(IN)                          :: zetp, zetq, zetw, rho, rac2, omega, &  local
548       f0 = 2.0_dp*SQRT(pi**5*zetw)*zetp*zetq
581    SUBROUTINE cps_vgauss2(v, nmax, zetp, zetq, zetw, rho, rac2, omega, r_cutoff)  argument
584       REAL(KIND=dp), INTENT(IN)                          :: zetp, zetq, zetw, rho, rac2, omega, &  local
610       f0 = 2.0_dp*SQRT(pi**5*zetw**3)/(rho + omega)*expT
639    SUBROUTINE cps_gauss2(v, nmax, zetp, zetq, zetw, rho, rac2, omega, r_cutoff)  argument
642       REAL(KIND=dp), INTENT(IN)                          :: zetp, zetq, zetw, rho, rac2, omega, &  local
660       f0 = pi**3*SQRT(zetw**3/(rho + omega)**3)*expT
686    SUBROUTINE cps_truncated2(v, nmax, zetp, zetq, zetw, rho, rac2, omega, r_cutoff)  argument
689       REAL(KIND=dp), INTENT(IN)                          :: zetp, zetq, zetw, rho, rac2, omega, &  local
705       f0 = 2.0_dp*SQRT(pi**5*zetw)*zetp*zetq