GPU accelerated computation of Polarized Subsurface BRDF for Flat Particulate Layers

by   Charly Collin, et al.

BRDF of most real world materials has two components, the surface BRDF due to the light reflecting at the surface of the material and the subsurface BRDF due to the light entering and going through many scattering events inside the material. Each of these events modifies light's path, power, polarization state. Computing polarized subsurface BRDF of a material requires simulating the light transport inside the material. The transport of polarized light is modeled by the Vector Radiative Transfer Equation (VRTE), an integro-differential equation. Computing solution to that equation is expensive. The Discrete Ordinate Method (DOM) is a common approach to solving the VRTE. Such solvers are very time consuming for complex uses such as BRDF computation, where one must solve VRTE for surface radiance distribution due to light incident from every direction of the hemisphere above the surface. In this paper, we present a GPU based DOM solution of the VRTE to expedite the subsurface BRDF computation. As in other DOM based solutions, our solution is based on Fourier expansions of the phase function and the radiance function. This allows us to independently solve the VRTE for each order of expansion. We take advantage of those repetitions and of the repetitions in each of the sub-steps of the solution process. Our solver is implemented to run mainly on graphics hardware using the OpenCL library and runs up to seven times faster than its CPU equivalent, allowing the computation of subsurface BRDF in a matter of minutes. We compute and present the subsurface BRDF lobes due to powders and paints of a few materials. We also show the rendering of objects with the computed BRDF. The solver is available for public use through the authors' web site.


Haptic Rendering of Thin, Deformable Objects with Spatially Varying Stiffness

In the real world, we often come across soft objects having spatially va...

Learning Multiple-Scattering Solutions for Sphere-Tracing of Volumetric Subsurface Effects

Accurate subsurface scattering solutions require the integration of opti...

Accelerating Inverse Rendering By Using a GPU and Reuse of Light Paths

Inverse rendering seeks to estimate scene characteristics from a set of ...

An Extension of the Strain Transfer Principle for Fiber Reinforced Materials

Fiber optical strain sensors are used to measure the strain at a particu...

Position-free Multiple-bounce Computations for Smith Microfacet BSDFs

Bidirectional Scattering Distribution Functions (BSDFs) encode how a mat...

Frequency-Dependent Material Motion Benchmarks for Radiative Transfer

We present a general solution for the radiation intensity in front of a ...

Computing Light Transport Gradients using the Adjoint Method

This paper proposes a new equation from continuous adjoint theory to com...