Light field eikonal calculation for given scaling of illumination distribution

The problem of input illumination distribution scaling in non-paraxial case is considered. To obtain the light field eikonal function the numerical method is proposed. It is based on the representation of eikonal function as a system of bi-cubic splines. Spline parameters are derived numerically using the criterion of minimum scaling transformation error. The results of eikonal function computing for scaling of square and round areas with different scaling coefficients are presented.