https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/random-number-generator-vsrnggaussian-used-in-subroutine/m-p/874329#M8776 <BR /><SPAN style="font-size: small;">Your program creates and initializes random stream with the same basic random generator and seed 10 times in subroutine test_sub1. So, output of the Gaussian generator is not changed. Can you create/initialize random stream once in your main program andpass the stream in the subroutine test_sub1? Please, do not forget to remove initialization of the stream in the sub-routine.Also, please free resources allocated for random stream in the main programwhengeneration of random numbers is completed.Routine vsldeletestream can be used for this purpose. Thanks, Andrey</SPAN> Wed, 25 Mar 2009 06:49:09 GMT Andrey_N_Intel 2009-03-25T06:49:09Z https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/random-number-generator-vsrnggaussian-used-in-subroutine/m-p/874328#M8775 Hi, guys:<BR /> I want to create a random number with Gaussian distribution probability by function vsrnggaussian.<BR /><BR />The subroutine is:<BR />INCLUDE '/home/fengy/mkl/10.0.3.020/include/mkl_vsl.fi'<BR />INCLUDE '/home/fengy/mkl/10.0.3.020/examples/vslf/source/errcheck.inc'<BR /><BR /><BR /><BR />SUBROUTINE test_sub1<BR /> USE MKL_VSL_TYPE<BR /> USE MKL_VSL<BR /><BR /> IMPLICIT NONE<BR /> <BR /> REAL, DIMENSION(1) :: r<BR /><BR /> INTEGER(KIND=4) i, nn<BR /> INTEGER n<BR /> INTEGER(KIND=4) errcode<BR /><BR /> INTEGER brng, method, seed<BR /><BR /> TYPE(VSL_STREAM_STATE) :: stream<BR /> REAL a, sigma<BR /> <BR /><BR /><BR /> n=1<BR /> brng = VSL_BRNG_MCG31<BR /> method = 0<BR /> seed = 777<BR /> a = 0.0<BR /> sigma = 1.0<BR /> <BR /> errcode = vslnewstream(stream, brng, seed)<BR /> CALL CheckVslError(errcode)<BR /><BR /> <BR /> errcode = vsrnggaussian(method, stream, n, r, a, sigma)<BR /> CALL CheckVslError(errcode)<BR /><BR /> WRITE(*,*) 'now r is: ', r(1)<BR /><BR /> <BR /><BR /><BR /><BR /><BR />END SUBROUTINE test_sub1<BR /><BR /> The main program is:<BR />PROGRAM test_g_r<BR /><BR /><BR /> IMPLICIT NONE<BR /><BR /><BR /> INTEGER n<BR /><BR /> REAL, DIMENSION(1) :: r<BR /> <BR /><BR /><BR /><BR /> <BR /> DO n=1, 10<BR /> CALL test_sub1<BR /> ENDDO<BR /><BR /><BR /><BR />END PROGRAM test_g_r<BR /><BR />However, the output of r(1) does not change. It is a intrinsic problem for function vsrnggaussian? What can I do so that each time test_sub1 is called, it will give a random number which has Gaussian distribution density?<BR /><BR /> Thanks!<BR /> Arthur<BR /> Tue, 24 Mar 2009 17:34:01 GMT https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/random-number-generator-vsrnggaussian-used-in-subroutine/m-p/874328#M8775 fymaterials 2009-03-24T17:34:01Z https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/random-number-generator-vsrnggaussian-used-in-subroutine/m-p/874329#M8776 <BR /><SPAN style="font-size: small;">Your program creates and initializes random stream with the same basic random generator and seed 10 times in subroutine test_sub1. So, output of the Gaussian generator is not changed. Can you create/initialize random stream once in your main program andpass the stream in the subroutine test_sub1? Please, do not forget to remove initialization of the stream in the sub-routine.Also, please free resources allocated for random stream in the main programwhengeneration of random numbers is completed.Routine vsldeletestream can be used for this purpose. Thanks, Andrey</SPAN> Wed, 25 Mar 2009 06:49:09 GMT https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/random-number-generator-vsrnggaussian-used-in-subroutine/m-p/874329#M8776 Andrey_N_Intel 2009-03-25T06:49:09Z