FORTRAN中截断正态分布的伪随机数生成

Question

我试图使用FORTRAN从左截断的正态分布中绘制伪随机数。我希望函数返回与输入具有相同维度的值：

功能（MU，SIGMA）; mu = N x 1且SIGMA = N X 1

有人可以帮忙吗？ 提前致谢

Answer 1

也许使用GNU Scientific Library和Fortran接口：http://www.lrz.de/services/software/mathematik/gsl/fortran/。 GSL有一个随机的正常发生器。如果数字偏错，请拒绝并重试。最简单的方法是将它全部放在一个循环中，并以串行方式使用底层随机数生成器。正确的多线程随机数是微妙的。