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

[Pathfinder]

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

[Pathfinder]

Спасибо за ответы, но хотел бы сделать несколько уточнений:

У каждого понятие о "хорошем" распределении свое.

Есть объективные численные критерии - период последовательности, хи квадрат, критерий Андерсена-Дарлинга, Колмогорова-Смирнова и др.
С суммированием и центральной предельной теоремой все понятно, только тут куча недостатков есть:
1. Нужен многопоточный генератор равномерного распределения - следовательно еще придется алгоритм вычисления начальных условий для него реализовать
2. Период формируемой последовательности уменьшится пропорционально числу суммируемых потоков, чего очень не хочется - желательно уложиться в 32 бита и получить период порядка 10^9-10^10
3. Число суммируемых потоков - боюсь, 20 все же с трудом хватит...
4. Суммируемые числа надо сгенерить, потом их просуммировать и отнормировать (тут еще небольшая потеря точности). При числе потоков порядка 10 на сигнальном проце все удовольствие займет порядка 40 инструкций на одну выборку - многовато

2 gab, спасибо за ссылку, еще нашел несколько статей этих же авторов. Похоже, предложенный ими метод наиболее эффективен.