The problem of generating equally weighted test particles from the one way flux of drifting Maxwellian is tackled. This paper extends previous work on the subject by presenting a simple and efficient rejection sampling algorithm together with C++ source files. The properties of the underlying probability distribution function, having the form of a normal distribution times x with positive support, are also disseminated. The method presented in this paper has been successfully used to combine fluid and kinetic models for trace impurity problems in plasma physics.
- Monte Carlo