Parallel Monte Carlo Approach for Integration of the Rendering Equation

This paper is addressed to the numerical solving of the rendering equation in realistic image creation. The rendering equation is integral equation describing the light propagation in a scene accordingly to a given illumination model. The used illumination model determines the kernel of the equation under consideration. Nowadays, widely used are the Monte Carlo methods for solving the rendering equation in order to create photorealistic images. In this work we consider the Monte Carlo solving of the rendering equation in the context of the parallel sampling scheme for hemisphere. Our aim is to apply this sampling scheme to stratified Monte Carlo integration method for parallel solving of the rendering equation. The domain for integration of the rendering equation is a hemisphere. We divide the hemispherical domain into a number of equal sub-domains of orthogonal spherical triangles. This domain partitioning allows to solve the rendering equation in parallel. It is known that the Neumann series represent the solution of the integral equation as a infinity sum of integrals. We approximate this sum with a desired truncation error (systematic error) receiving the fixed number of iteration. Then the rendering equation is solved iteratively using Monte Carlo approach. At each iteration we solve multi-dimensional integrals using uniform hemisphere partitioning scheme. An estimate of the rate of convergence is obtained using the stratified Monte Carlo method. This domain partitioning allows easy parallel realization and leads to convergence improvement of the Monte Carlo method. The high performance and Grid computing of the corresponding Monte Carlo scheme are discussed.

[1]  Anton Penzov.  Shading and illumination models in computer graphics : a literature survey , 1992 .

[2]  James T. Kajiya,et al.  The rendering equation , 1998 .

[3]  László Szirmay-Kalos,et al.  Monte-Carlo Global Illumination Methods State of the Art and New Developments , 2000 .

[4]  Ivan Dimov,et al.  Statistical modeling of pulse excited electron quantum kinetics in a one-band semiconductor , 1998 .

[5]  Alexander Keller.  Quasi-Monte Carlo Methods in Computer Graphics: The Global Illumination Problem , 1995 .

[6]  László Szirmay-Kalos,et al.  MONTE-CARLO METHODS IN GLOBAL ILLUMINATION , 2000 .

[7]  Åke Björck,et al.  Numerical Methods , 1995, Handbook of Marine Craft Hydrodynamics and Motion Control.

[8]  A. Penzov.  Shading and illumination models in computer graphics (Research report of the Computer Graphics Group CG 4/1992) , 1992 .

[9]  Leonidas J. Guibas,et al.  Robust Monte Carlo methods for light transport simulation , 1997 .

[10]  Stephen Gilmore,et al.  Evaluating the Performance of Skeleton-Based High Level Parallel Programs , 2004, International Conference on Computational Science.

[11]  Philip Dutré,et al.  Global Illumination Compendium , 2003 .

[12]  I. Dimov,et al.  Error analysis of an adaptive Monte Carlo method for numerical integration , 1998 .

[14]  Aneta Karaivanova,et al.  Parallel Importance Separation for Multiple Integrals and Integral Equations , 2004, International Conference on Computational Science.

[15]  P. Deuflhard,et al.  Large Scale Scientific Computing , 1987 .

[16]  Rayna Georgieva,et al.  Importance Separation for Solving Integral Equations , 2003, LSSC.

[17]  G. Marchuk,et al.  Numerical methods and applications , 1995 .

[18]  Ivan Tomov Dimov,et al.  A Monte Carlo Approach for the Cook-Torrance Model , 2004, NAA.

[19]  Rayna Georgieva,et al.  Parallel Importance Separation and Adaptive Monte Carlo Algorithms for Multiple Integrals , 2002, Numerical Methods and Application.

[20]  Ivan Tomov Dimov,et al.  Parallel Monte Carlo Sampling Scheme for Sphere and Hemisphere , 2006, Numerical Methods and Applications.