Compute a hierarchical radiosity solution with some initial tolerance on the error in the links using something like CreateLinks. Your error term should probably ignore energy. Solve the system using something like the GatherSweep routine.
