Section: Scientific Foundations
Keywords : Radiosity, Radiance, Finite Elements Modeling, Monte-Carlo, Physically-Based Rendering.
Numerical light simulation
Numerical simulation of light means solving for light intensity in the ``Rendering Equation'', an integral equation modeling energy transfers (or light intensity transfers). The Rendering Equation was first formalized by Kajiya  , and is given by:
Computing global illumination (i.e., solving for intensity in Equation 1 ) in general environments is a challenging task. Global illumination may be considered in terms of computing the interactions between the lighting signal and the geometric signal (i.e., the scene). These interactions occur at various scales . This issue belongs to the same class of problems encountered by Geometry Processing, described in the previous section. As a consequence, the signal processing family of approaches is again a well-suited formalism. As such, the multi-scale approach is a natural choice, which dramatically improves performances. Environments composed of a large number of primitives, such as highly tessellated models, show a high variability of these scales.
In addition, these methods are challenged with more and more complex materials that need to be taken into account in the simulation. The simple diffuse Lambert law has been replaced with much more complex reflection models. The goal is to create synthetic images that no longer have a synthetic aspect, in particular when human characters are considered.
One of the difficulties is finding efficient ways of evaluating the visibility term. This is typically a Computational Geometry problem, i.e., a matter of finding the right combinatorial data structure (the visibility complex ), studying its complexity and deriving algorithms to construct it. To deal with this issue, several teams (including VEGAS, ARTIS and REVES) study the visibility complex.
The other terms of the Rendering Equation cannot be solved analytically in general. Many different numerical resolution methods have been used. The main difficulties of the discipline is that each time a new physical effect should be simulated, the numerical resolution methods need to be adapted. In the worst case, it is even necessary to design a new ad-hoc numerical resolution method. For instance, in Monte-Carlo based solvers, several sampling maps are used, one for each effect (a map is used for the diffuse part of lighting, another map is used for caustics ...). As a consequence, the discipline becomes a collection of (sometimes mutually exclusive) techniques, where each of these technique can only simulate a specific lighting effect.
The other difficulty is to satisfy two somewhat antinomic objectives at the same time. On the one hand, we want to simulate complex physical phenomena (subsurface scattering, polarization, interferences, ...), responsible for subtle lighting effects. On the other hand, we want to visualize the result of the simulation in real-time.
We first experimented finite-element methods in parameter space, and developped the Virtual Mesh approach  and a parallel solution mechanism for the associated hierarchical finite element formulation  The initial method was dedicated to scenes composed of quadrics. We combined this method with our Geometry Processing methods to improve the visualization  . More recently, we developped a method to directly apply the finite element method on tesselated models  .
One of our goals is now to design new representations of lighting coupled with the geometric representation. These representations of lighting need to be general enough so as to be easily extended when multiple physical phenomena should be simulated. Moreover, we want to be able to use these representations of lighting in the frame of real-time visualization.