“…In the paper of Fedoseyev et al [8], the computational domain, Ω was extended slightly beyond the boundaries, ∂Ω k given by the line segments of the unit square: (0,1) (0,0), (0,1) (1,1) (0,1), The domain extended 0.020 beyond the unit square in the x 1 and x 2 directions. The parameters on the boundaries, see Wertz et al [9], were increased by a factor of 30. Consequently, numerical integration such as Gauss-Legendre (GL) needs to be performed, by first interpolating this integrand onto the appropriate GL zeros.…”