Improvement of the rejection technique for sampling

in Monte Carlo simulation

 

Amélioration de la technique de rejet pour l’échantillonnage

dans une simulation par Monte Carlo

 

Z. Faik Ouahab, A. Jehouani*, J. Ghassoun,  A. Oubaiddou

Laboratoire de Physique et Techniques Nucléaires, Département de Physique, Faculté des Sciences Semlalia,

B.P.2390, Université Cadi Ayyad, 40000 Marrakech, Maroc.

* Corresponding author. E-mail:  jehouani@ucam.ac.ma

Received: 02 April 2006; revised version accepted:22 May 2008

 

Abstract

     Basically, the calculations by the Monte Carlo simulation use several techniques to sample the density functions. The rejection technique is used when the density function is impossible to invert but it requires huge random numbers to converge. In this paper we propose an improvement of the rejection technique in order to reduce the amount of these random numbers and make this technique more efficient. Many numerical tests are done for fission spectrum and Gauss distribution.

 

Keywords: Monte Carlo ; Rejet ; Echantillonnage; Fission.

 

Résumé

     La simulation des phénomènes physiques par la méthode Monte Carlo est basée sur l’échantillonnage des différents paramètres mis en jeu. La technique de rejet permet de faire l’échantillonnage en imposant des critères sur l’acceptation de la valeur tirée aléatoirement. La masse des nombres aléatoires rejetés est souvent énorme.

     Pour remédier à ce problème nous proposons une amélioration de la technique de rejet pour réduire la perte des nombres aléatoires et du temps de calculs. Cette amélioration est testée avec succès pour des cas réels du spectre de fission et la distribution de Gauss.

 

Mots clés: Monte Carlo ; Rejet ; Echantillonnage; Fission.

© 2015