Despite their numerous applications, efficiently rendering participating media remains a challenging task due to the intricacy of the radiative transport equation. As they provide a generic means of solving a wide variety of problems, numerical methods are most often used to solve the air-light integral even under simplifying assumptions. In this paper, we present a novel analytical approach to...

This paper surveys global illumination algorithms for environments including participating media and accounting for multiple scattering. The objective of this survey is the characterization of those methods: Identification of their base techniques, their assumptions, limitations and range of utilization. To this end, the algorithms are grouped into functional categories and each method is brief...

Several scalable many-light rendering methods have been proposed recently for the efficient computation of global illumination. However, gathering contributions of virtual lights in participating media remains an inefficient and time-consuming task. In this paper, we present a novel sparse sampling and reconstruction method to accelerate the gathering step of the many-light rendering for partic...

This paper presents a fast approximation for the solution of the radiative transfer equation in inhomogeneous participating media. This scheme traces rays from the point light creating spherical wavefronts. While we march on a particular ray, we assume that other rays find similar material properties, which allows the solution of the diffusion equation for this ray. The proposed approximation o...

For the rendering of multiple scattering effects in participating media, methods based on the diffusion approximation are an extremely efficient alternative to Monte Carlo path tracing. However, in sufficiently transparent regions, classical diffusion approximation suffers from non-physical radiative fluxes which leads to a poor match to correct light transport. In particular, this prevents the...

This paper presents a novel method that effectively combines both control variates and importance sampling in a sequential Monte Carlo context. The radiance estimates computed during the rendering process are cached in a 5D adaptive hierarchical structure that defines dynamic predicate functions for both variance reduction techniques and guarantees well-behaved PDFs, yielding continually increa...