Abstract The method of paraxial complex geometrical optics (CGO) is presented, which describes Gaussian beam diffraction in arbitrary smoothly inhomogeneous media, including lens-like waveguides. By way of an example, the known analytical solution for Gaussian beam diffraction in free space is presented. Paraxial CGO reduces the problem of Gaussian beam diffraction in inhomogeneous media to the system of the first order ordinary differential equations, which can be readily solved numerically. As a result, CGO radically simplifies the description of Gaussian beam diffraction in inhomogeneous media as compared to the numerical methods of wave optics. For the paraxial on-axis Gaussian beam propagation in lens-like waveguide, we compare CGO solutions with numerical results for finite differences beam propagation method (FD-BPM). The CGO method is shown to provide 50-times higher rate of calculation then FD-BPM at comparable accuracy. Besides, paraxial eikonal-based complex geometrical optics is generalized for nonlinear Kerr type medium. This paper presents CGO analytical solutions for cylindrically symmetric Gaussian beam in Kerr type nonlinear medium and effective numerical solutions for the self-focusing effect of Gaussian beam with elliptic cross section. Both analytical and numerical solutions are shown to be in a good agreement with previous results, obtained by other methods.