/*
*
* Author: Dan Lynch
*
*/

#pragma hint switcher toggle
#pragma label switcher "Turbulence/Brownian"
#pragma label amp "Amplitude"
#pragma label lacunarity "Lacunarity"
#pragma label recursions "Recursions"
#pragma label freq "Frequency"
#pragma group Displace freq offset recursions lacunarity amp

displace
d_ndisp(
        int switcher = 0;
        vector freq = 1, offset = 0;
        int recursions = 5;
        float lacunarity = 2;
        float amp=1;
        )
{
    vector Pnz = wo_space(P+offset)*freq;
    vector Nn = normalize(N);
    float Nz = 0;
    float mlt = 1;
    float iF;
    int i;

    for(i=0;i<recursions;i++){
        iF = i*5.6789;

        if(switcher==0) {
                Nz += abs(noise(iF+Pnz*mlt)-.5)/mlt;
        } else {
                Nz += (noise(iF+Pnz*mlt)-.5)/mlt;
        }

        mlt *= lacunarity;
    }
    P += Nn * Nz * amp;
    N = computenormal(P);
}