Ahoj, potřeboval bych doporučit vhodný algoritmus generování PRND.
Řeším fyzikální úlohu (proudění) kde se jistý člen řešených rovnic z určitého důvodu trochu pošťouchne náhodnou fluktuací.
Jde o to, že v každé buňce výpočetní sítě je potřeba 50 PRND, výpočetní síť má 1e6 až 1e7 buněk a provádí se 1e4 až 1e5 iterací, tzn. náhodné číslo je potřeba získat 5e11-krát až 5e13-krát.
Výpočetní čas se přitom nesmí zvýšit o více než 20-30% oproti modelu bez náhodných flutuací (dny až týdny).
RAND_MAX stačí 1e3 až 1a4, posloupnost PRND by se (asi) neměla opakovat dříve než po 1e8 realizacích.
Nesmí vadit omp parallel.
Díky