Method init

openpgl/directional/vmm/VMMPhaseFunctions.cpp:11–2455

int VMMSingleLobeHenyeyGreensteinOracle::numRepresentations;

void VMMSingleLobeHenyeyGreensteinOracle::init()
{
    /* VMM representations of the single lobe HG phase function.
       The 3-lobe VMMs are fitted using a least squares optimization framework
       as described by S. Herholz:
                - "Path Guiding in Production"
                - - "Volumetric Zero-variance Based Path Guiding" (Chapter 11.5)
    */

    minMeanCosine = 0.0f;
    maxMeanCosine = 0.99f;
    numRepresentations = 128;

    representations[0].K = 3;
    representations[0].g = 0.0f;
    representations[0].weights[0] = 0.01688512612073299f;
    representations[0].meanCosines[0] = 0.0f;
    representations[0].kappas[0] = 0.0f;
    representations[0].normalizations[0] = 0.07957747154594767f;
    representations[0].weights[1] = 0.34940739552081895f;
    representations[0].meanCosines[1] = 0.0f;
    representations[0].kappas[1] = 0.0f;
    representations[0].normalizations[1] = 0.07957747154594767f;
    representations[0].weights[2] = 0.6337074783584481f;
    representations[0].meanCosines[2] = 0.0f;
    representations[0].kappas[2] = 0.0f;
    representations[0].normalizations[2] = 0.07957747154594767f;
    representations[0].weights[3] = 0.0f;
    representations[0].meanCosines[3] = 0.0f;
    representations[0].kappas[3] = 0.0f;
    representations[0].normalizations[3] = 0.0f;

    representations[1].K = 3;
    representations[1].g = 0.007795275590551181f;
    representations[1].weights[0] = 0.01688512612254555f;
    representations[1].meanCosines[0] = 0.03588398182493312f;
    representations[1].kappas[0] = 0.10774447737138962f;
    representations[1].normalizations[0] = 0.08845920147738665f;
    representations[1].weights[1] = 0.34940739552017397f;
    representations[1].meanCosines[1] = 0.014362686373439479f;
    representations[1].kappas[1] = 0.043093986007026464f;
    representations[1].normalizations[1] = 0.08305603677895068f;
    representations[1].weights[2] = 0.6337074783572806f;
    representations[1].meanCosines[2] = 0.003425584458792068f;
    representations[1].kappas[2] = 0.010276833773244303f;
    representations[1].normalizations[2] = 0.08039807745351672f;
    representations[1].weights[3] = 0.0f;
    representations[1].meanCosines[3] = 0.0f;
    representations[1].kappas[3] = 0.0f;
    representations[1].normalizations[3] = 0.0f;

    representations[2].K = 3;
    representations[2].g = 0.015590551181102362f;
    representations[2].weights[0] = 0.01688437903886321f;
    representations[2].meanCosines[0] = 0.07165186574299764f;
    representations[2].kappas[0] = 0.21569511380458395f;
    representations[2].normalizations[0] = 0.09797223006527885f;
    representations[2].weights[1] = 0.34940741896923394f;