Uma resposta recente mencionou o uso de Fortuna ou Mersenne Twister Random Number Generators ( RNGs ) para semear uma simulação de Monte Carlo . Eu nunca tinha ouvido falar de Fortuna antes, então procurei - parece que ele se destina principalmente ao uso criptográfico. Atualmente, uso um Mersenne...
Suponha que eu esteja usando um gerador de números pseudo-aleatórios lineares congruentes (PRNG). Dada uma semente x0x0x_0 , o fator multiplicador (a), o fator de deslocamento (c) e o fator de módulo (m), como posso determinar o período do meu PRNG? Eu o determino por algoritmos de detecção de...
Suponha que eu soubesse que uma sequência numérica aleatória foi gerada por um gerador congruencial linear. Isso é, xn+1=(aXn+c)modmxn+1=(aXn+c)modmx_{n+1}=(aX_n+c) \bmod m Se eu sou dado o período inteiro (ou, pelo menos, uma grande subsequência contígua do mesmo), como pode reconstruir os...
Recentemente, me deparei com um comentário afirmando que quase todos os pesquisadores que usam os métodos de Monte Carlo estão fazendo errado. Em seguida, foi elaborado que a simples escolha de sementes diferentes para diferentes instâncias de um PRNG, como o Mersenne Twister, não é suficiente para...
Sempre que se usa um método quase-Monte Carlo para cubagem ou otimização, parece que há uma grande variedade de sequências de baixa discrepância para escolher, associadas aos nomes de van der Corput, Halton, Hammersley, Faure, Niederreiter, Sobol ', e outros nomes dos quais não me lembro. Existem...
Normalmente, quando eu proponho um gerador seqüencial de números aleatórios em C, uso a chamada srand(time(NULL)) então use rand() mod N para obter um número inteiro aleatório entre 0 e N-1. No entanto, quando faço isso em paralelo, as chamadas ao tempo (NULL) são tão próximas uma da outra...
Quero acelerar meu programa python usando cython, mas meu loop interno ainda está fazendo chamadas lentas de python para o gerador de números aleatórios! Vários anos atrás, esse mesmo problema foi levantado por alguém que apoiava os sábios e parecia não haver uma boa solução naquele momento. Não é...
Eu usei MT19937em um equipamento de teste para gerar valores [0, - 1] de 32 bits uniformemente (não assinados) , com base na implementação mt19937.c do Authors original , para gerar um suprimento (essencialmente inesgotável) de fluxo estatisticamente aleatório de octetos de bits. Isso substitui um...