Global Illumination Computations with Hierarchical Finite Element Methods