The source code available here is provided without any guarantees. If you find bugs please send me a message. For a more detailed explanation of the method, see the corresponding paper on my publication page.
Error bounds for hierarchical radiosity.
Massively parallel progressive radiosity with parallel visibility computations and dynamic load balancing. Meshing for arbitrary convex polygons. Meshing with Delaunay triangulation and Voronoi diagrams. Adaptive discontinuity meshing. Hemispherical projection for the computation of exact and approximate local illumination.