The latter author, together with collaborators, proposed a numerical scheme
to calculate the price of barrier options. The scheme is based on a
symmetrization of diffusion process. The present paper aims to give a
mathematical credit to the use of the numerical scheme for Heston or SABR type
stochastic volatility models. This will be done by showing a fairly general
result on the symmetrization (in multi-dimension/multi-reflections). Further
applications (to time-inhomogeneous diffusions/ to time dependent boundaries/to
curved boundaries) are also discussed