Una respuesta reciente menciona el uso de los generadores de números aleatorios ( RNG ) Fortuna o Mersenne Twister para sembrar una simulación de Monte Carlo . No había oído hablar de Fortuna antes, así que lo busqué, parece que está destinado principalmente para uso criptográfico. Actualmente uso...
Supongamos que estoy usando un generador de números pseudoaleatorio congruencia lineal (PRNG). Dada una semilla , el factor de multiplicación (a), el factor de desplazamiento (c) y el factor de módulo (m), ¿cómo puedo determinar el período de mi PRNG? ¿Lo determino mediante algoritmos de...
Supongamos que supiera que un generador congruencial lineal genera una secuencia de números aleatorios. Es decir, Xn + 1= ( una Xnorte+ c ) mod mXnorte+1=(unaXnorte+C)modificaciónmetrox_{n+1}=(aX_n+c) \bmod m Si me dan todo el período (o al menos una gran subsecuencia contigua de la misma), ¿cómo...
Recientemente, me encontré con un comentario afirmando que casi todos los investigadores que hacen los métodos de Monte Carlo lo están haciendo mal. Continuó explicando que simplemente elegir diferentes semillas para diferentes instancias de un PRNG, como el Mersenne Twister, no es suficiente para...
Cada vez que uno usa un método cuasi-Monte Carlo para cubicación u optimización, parece que hay una gran variedad de secuencias de baja discrepancia para elegir, asociadas con los nombres de van der Corput, Halton, Hammersley, Faure, Niederreiter, Sobol ', y otros nombres que no recuerdo del todo....
Normalmente cuando siembro un generador secuencial de números aleatorios en C, uso la llamada srand(time(NULL)) luego usa rand() mod N para obtener un entero aleatorio entre 0 y N-1. Sin embargo, cuando hago esto en paralelo, las llamadas al tiempo (NULL) están tan cerca una de la otra que...
Quiero hacer que mi programa de Python sea rápido usando cython, ¡pero mi ciclo interno todavía está haciendo llamadas lentas de Python al generador de números aleatorios! Hace varios años, alguien con el apoyo de Sage planteó este mismo problema y parecía que no había una buena solución en ese...
Utilicé MT19937en un arnés de prueba para generar valores uniformes (sin signo) de 32 bits [0, - 1], basados en la implementación original mt19937.c de los Autores , para generar un suministro (esencialmente inagotable) de flujo estadísticamente aleatorio de bit-octetos. Esto es en lugar de un...