Random number generation using ranlib

danielpollmann — Mon, 09 Jul 2012 12:32:36 GMT

Hi,

for an application, I am generating random variates from a standard exponential distribution using the function sexpo() included in the popular ranlib.f (e.g., http://people.sc.fsu.edu/~jburkardt/f77_src/ranlib/ranlib.html).

It works perfectly fine using the compiler version 11.1.072. However, if I use version 12.3.174 (with MKL included), the function returns very large integer values instead which are definitely not drawn from a standard exponential distribution. Unfortunately, I need to use the newer version because it has better support for allocatable arrays.

I was wondering whether anyone had come across either this or a similar problem, or has an idea what might go wrong in the newer version given its new features or any particularities of MKL.

Thanks!

Random number generation using ranlib

Steven_L_Intel1 — Mon, 09 Jul 2012 14:52:03 GMT

I wonder if you have declared the procedure properly in the caller. It should be returning reals, not integers. Since you have the source code, you should be able to step through it with the debugger and see what is going wrong.

topic Random number generation using ranlib in Intel® Fortran Compiler

Random number generation using ranlib

Random number generation using ranlib