Не понимаю как работает алгоритм

typedef unsigned int uint;
#define CD static_cast<double>
#define CI static_cast<int>
double PerlinRand(uint M, uint A, uint C, uint &Z) {
    Z = (A * Z + C) % M;
    return CD(Z) / M;
}
void Draw(short x, short y, char c) {
    COORD p = { x, y };
    SetConsoleCursorPosition(GetStdHandle(STD_OUTPUT_HANDLE), p);
    cout << c;
}
int main() {
    default_random_engine rnd(time(NULL));
    uint M = (1 << 32) - 1;//4294967295
    uint A = 1'664'525;
    uint C = 1;
    uint Z = static_cast<uint>
        ((CD(rnd() % 1'000'000'000) / 1'000'000'000) * M);
    int w = 1024;
    int h = 256;
    
    int x = 0;
    int y = h / 2;
    int amp = 3;
    int wl = 5;
    double a = PerlinRand(M, A, C, Z);
    double b = PerlinRand(M, A, C, Z);
    
    while (x < w) {
        if (x % wl == 0) {
            a = b;
            b = PerlinRand(M, A, C, Z);
            y = CI(h / 4 + a * amp);
 
            Draw(x, y, '0');
        }
        x += 1;
    }
    _getch();
    return 0;
}

double PerlinRand(uint &Z) {
    uint64_t M = uint64_t(1) << 32;
    uint64_t A = 1'664'525;
    uint64_t C = 1;
    Z = static_cast<uint>((A * Z + C) % M);
    return static_cast<double>(Z) / M;
}