Testar a qualidade dos PRNGs

Estou brincando com PRNGs (como Mersenne Twister erand() stdlib) e eu gostaria de um bom teste que me ajudasse a verificar a qualidade dos dados aleatórios produzidos pelos PRNGs. Calculei o valor de Pi usando números aleatórios gerados pelos ...

Encontrando dígitos PI usando Monte Carlo

Tentei muitos algoritmos para encontrar π usando Monte Carlo. Uma das soluções (em Python) é esta: def calc_PI(): n_points = 1000000 hits = 0 for i in range(1, n_points): x, y = uniform(0.0, 1.0), uniform(0.0, 1.0) if (x**2 + y**2) <= 1.0: hits ...

Calculando π usando as limitações de uma simulação de Monte Carlo

Fiz uma pergunta muito semelhante a essa, então mencionarei as soluções anteriores no final, tenho umalocal na rede Internet [https://monte-carlo.sjorsvanholst.nl]que calcula π com a CPU do cliente enquanto a armazena em um servidor, até agora eu ...

Cálculo de Monte Carlo de Pi em Scala

Suponha que eu gostaria de calcular Pi com a simulação de Monte Carlo como um exercício. Estou escrevendo uma função, que escolhe um ponto em um quadrado(0, 1), (1, 0) aleatoriamente e testa se o ponto está dentro do círculo. import ...

Simulação de preço de ações R code - Slow - Monte Carlo

Distribuição aleatória uniforme (Monte-Carlo) na esfera unitária

Preciso de um esclarecimento com o algoritmo que gera valores aleatórios para o meu traçador de raios de estimação. Emito raios de um ponto. E eu tenho o problema com a distribuição desses raios: preciso que a distribuição seja uniforme, mas não ...

Acelerar a amostragem da estimativa do kernel

Aqui está um

Distribuição uniforme de pontos Python em esfera de 4 dimensões

Eu preciso de uma distribuição uniforme de pontos em uma esfera de 4 dimensões. Eu sei que isso não é tão trivial quanto escolher 3 ângulos e usar coordenada...

Pode o Generator de números aleatórios do Fortran 90 ser confiável para a integração de Monte Carlo?

Escrevi um algoritmo de integração curto monte carlo para calcular uma integral no Fortran 90. Certa vez, comparei o resultado obtido resolvendo a integral com respeito a alguns parâmetros usando o gerador de números aleatórios intrínseco com o ...

Modelo Dinâmico Simples em PyMC3

Estou tentando montar um modelo de sistema dinâmico em PyMC3, para inferir dois parâmetros. O modelo é o SIR básico, comumente usado em epidemiologia: dS / dt = - r0 * g * S * I dI / dt = g * I (r * S - 1) onde r0 e g são parâmetros a serem ...