Just the indirect specular pass by importance sampling.
With all layers.
Manually traced by 3D Hammersley sequence with 8 samples. Link to the paper.
Actually this should be done in a pure Path Tracer with 3 random number from multiple dimension rng, that keeps numerical correction.
