Chadia, the first method that comes to me is rejection. This means that you keep generating regular Gaussion variables, using the GSL generator, with the prescribed mean and std dev until the outcome is in the interval where you want it to be. That is, you accept permissible values and you toss the other values away! Is this what you have used before? Abbas Alhakim ________________________________
From: [email protected] on behalf of chadia kanaan Sent: Mon 12/22/2008 1:45 PM To: [email protected] Subject: [Help-gsl] Truncated normal distribution using gsl_rng_gaussian Hello, Could someone, please, guide me how to generate truncated Gaussian probability distribution numbers between [min, max] using the already existing gsl_gaussian generators . lets say between 0 and infinity for example. In earlier gsl archive, this point was discussed and refered to other method like rejection. I have used this method before, but gsl_gaussian generator is at least 10 times faster and reliable than generators based on rejection/acceptance method. so I need a solution based on gsl_generator , maybe through simple transformation .. Many thanks for your help, regards, C. K _______________________________________________ Help-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/help-gsl _______________________________________________ Help-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/help-gsl
