libcrbn/sampler/vanDerCorput.hpp

#ifndef _VAN_DER_CORPUT_HPP_
#define _VAN_DER_CORPUT_HPP_

// Generalized Van Der Corput function
float vdC(int i, int base, int** permutation) {
        float u, f, factor;
        int   digit;

        u = 0.0;
        f = factor = 1.0 / (float)base;
        
        while(i > 0) {
                digit = i % base;

                u += (float)(permutation[base-1][digit]) * factor;
                i /= base;
                factor *= f;
        }

        return u;
}

#endif

---