Following some discussion with Brecht, submitting this diff here.
This is about adding support for anisotropy (phase) in the random walk subsurface, essentially implementing https://graphics.pixar.com/library/PathTracedSubsurface/.
The important things are:
1/ bssrdf->anisotropy and bssrdf->ior parameters (we have to raise the number of parameters accordingly in kernel_types.h);
2/ handling of the anisotropic walk and the corresponding albedo inversion (as published in our tech memo at Pixar) in integrator_subsurface.h.