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

Question

我想从使用FORTRAN的左截断正态分布中绘制伪随机数。我希望该函数返回与输入相同尺寸的值：

FUNCTION（MU，SIGMA）; mu = N x 1和SIGMA = N X 1

有人可以帮忙吗？ 在此先感谢

Answer 1

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