Конвейерная обработка генератора случайных чисел Зиггурата

Question

В настоящее время я использую версию генератора случайных чисел Зиггурата Джорджа Марсалья. Хотя это, предположительно, один из самых быстрых способов генерации нормально распределенных генераторов случайных чисел, он заполнен кодом управления контуром (т. Е. Операторы возврата в середине цикла, if-statements, branch и т. Д.), И это делает несколько вызовов стандартных функций C, таких как exp() и log(). Не говоря уже о бесконечном цикле.

Это делает код, который не может быть конвейеризован компилятором. В конечном счете, я чувствую, что базовый подход, например, использование центральной предельной теоремы напрямую, в конечном счете может быть быстрее, поскольку он может быть легко конвейерным. К сожалению, он не подходит для хвостов распределения Гаусса, и поэтому он неприемлем для моего приложения.

Есть ли у кого-нибудь идеи о том, как можно сократить код управления и вызовы функций. В настоящее время я использую Colin Green's implementation of the algorithm, который я портировал на C. Мой базовый унифицированный генератор - Tiny Mersenne Twister (так что, пожалуйста, не говорите мне, чтобы использовать MT, как я видел, другие люди, я уже здесь. для нормально распределенных RNG, неравномерных RNG).

Answer 1

Вы можете взглянуть на мою реализацию C here. Основная функция - это всего лишь 20 строк кода, поэтому должно быть легко развернуть цикл немного. Он также дает вам выбор использования целочисленного или поплавкового сравнения, в зависимости от того, что на вашем компьютере быстрее. Вы можете подключить любой базовый RNG.