Синтез гауссовского шума, Эффективный алгоритм синтеза псевдослучайной последовательности с норм Опции

[Pathfinder]

Кто-нибудь знает где взять эффективный алгоритм синтеза гауссовского шума в целочисленной арифметике? решение в лоб (2 рекуррентные последовательности максимальной длины+нелин. преобразования) требует вычисления функции sqrt(-ln(1-x)), а считается она очень медленно - для получения 16-битного корня требуется 16 операций деленья, и еще одна для вычисления логарифма. Может есть какие-то быстрые итеративные алгоритмы? Сразу уточню, период последовательности должен быть очень большим.

[SavSerG]

Вот алгоритм на основе суммы 20 равномерно распределенных чисел,
нелинейное преобразование для того, чтобы "хвосты" распределения
моделировались точнее. Уже при 12 числах нормальная гауссовская последовательность получается.
При сумме 12 чисел не надо вычислять корень, т.к. он равен 1.
long double gen()
{
long double x=0.0;
for(int i=0;i<20;i++){
x+=(double(rand())/double(RAND_MAX)-0.5);
};
x*=sqrtl(0.6); //sqrtl(12.0/N=20)
x+=0.01*x*(x*x-3.0);
return x;
}
Можно почитать книгу Numirical receipes on C там были алгоритмы генерирования
гауссовского шума