quasirandom.h

Go to the documentation of this file.

 
 
#ifndef THEORETICA_QUASIRANDOM_H
#define THEORETICA_QUASIRANDOM_H
 
#include "../core/real_analysis.h"
#include "../algebra/algebra_types.h"
 
 
namespace theoretica {
 
 
    inline real qrand_weyl(unsigned int n, real alpha = INVPHI) {
        return fract(n * alpha);
    }
 
 
    inline real qrand_weyl_recurr(real prev = 0, real alpha = INVPHI) {
 
        if(prev == 0) {
            prev = qrand_weyl(1, alpha);
            return prev;
        }
 
        return fract(prev + alpha);
    }
 
 
    template<unsigned int N>
    inline vec<real, N> qrand_weyl_multi(unsigned int n, real alpha) {
        
        vec<real, N> r;
        for (unsigned int i = 0; i < N; ++i)
            r[i] = fract(n * pow(alpha, i + 1));
 
        return r;
    }
 
 
    inline vec2 qrand_weyl2(unsigned int n, real alpha = 0.7548776662466927) {
        return {fract(n * alpha), fract(n * square(alpha))};
    }
 
 
}
 
#endif