02-May-2023 08:29:24 ziggurat_test(): MATLAB/Octave version 5.2.0. Test ziggurat(). ziggurat_test001(): shr3() returns pseudorandom uniformly distributed unsigned 32 bit integers. 0 123456789 1 -1579999413 2714967881 2838424670 2 -2056153898 2238813396 658813981 3 1250077441 1250077441 3488890837 4 -474866958 3820100336 775210481 5 -1117447608 3177519686 2702652726 6 -610828462 3684138832 2566691222 7 -1143879504 3151087790 2540259326 8 -632459186 3662508108 2518628602 9 -52590672 4242376622 3609917434 10 -920365316 3374601978 3322011304 0 987654321 1 248404469 248404469 1236058790 2 2078538413 2078538413 2326942882 3 -1837283784 2457683510 241254627 4 1841886731 1841886731 4602945 5 305946223 305946223 2147832954 6 -1479565191 2815402103 3121348326 7 736910199 736910199 3552312302 8 -1327526408 2967440886 3704351085 9 1833747846 1833747846 506221436 10 1249087608 1249087608 3082835454 0 123456789 1 -1579999413 2714967881 2838424670 2 -2056153898 2238813396 658813981 3 1250077441 1250077441 3488890837 4 -474866958 3820100336 775210481 5 -1117447608 3177519686 2702652726 6 -610828462 3684138832 2566691222 7 -1143879504 3151087790 2540259326 8 -632459186 3662508108 2518628602 9 -52590672 4242376622 3609917434 10 -920365316 3374601978 3322011304 ziggurat_test02(): r4_uni() returns pseudorandom uniformly distributed real numbers between 0 and 1. 0 123456789 1 -1579999413 2714967881 0.160872 2 -2056153898 2238813396 0.653392 3 1250077441 1250077441 0.312321 4 -474866958 3820100336 0.680493 5 -1117447608 3177519686 0.129260 6 -610828462 3684138832 0.097604 7 -1143879504 3151087790 0.091450 8 -632459186 3662508108 0.086414 9 -52590672 4242376622 0.340499 10 -920365316 3374601978 0.273466 0 987654321 1 248404469 248404469 0.787792 2 2078538413 2078538413 0.041784 3 -1837283784 2457683510 0.556171 4 1841886731 1841886731 0.501072 5 305946223 305946223 0.000081 6 -1479565191 2815402103 0.226745 7 736910199 736910199 0.327087 8 -1327526408 2967440886 0.362486 9 1833747846 1833747846 0.617864 10 1249087608 1249087608 0.217778 0 123456789 1 -1579999413 2714967881 0.160872 2 -2056153898 2238813396 0.653392 3 1250077441 1250077441 0.312321 4 -474866958 3820100336 0.680493 5 -1117447608 3177519686 0.129260 6 -610828462 3684138832 0.097604 7 -1143879504 3151087790 0.091450 8 -632459186 3662508108 0.086414 9 -52590672 4242376622 0.340499 10 -920365316 3374601978 0.273466 ziggurat_test03() r4_nor() returns pseudorandom normally distributed real numbers between 0 and 1. 0 123456789 1 -1579999413 2714967881 -1.348345 2 -2056153898 2238813396 0.321041 3 1250077441 1250077441 -0.689408 4 -474866958 3820100336 0.875903 5 -1117447608 3177519686 -1.036908 6 -610828462 3684138832 -0.749757 7 -1143879504 3151087790 -2.633581 8 -632459186 3662508108 -2.335211 9 -52590672 4242376622 -0.900580 10 -920365316 3374601978 -0.547212 0 987654321 1 248404469 248404469 0.678952 2 2078538413 2078538413 -1.028125 3 -1837283784 2457683510 0.233949 4 1841886731 1841886731 0.003315 5 736910199 736910199 -0.808394 6 -1327526408 2967440886 -0.635426 7 1833747846 1833747846 0.702161 8 1249087608 1249087608 -1.819247 9 2026661944 2026661944 -3.163224 10 1286162813 1286162813 -0.633432 0 123456789 1 -1579999413 2714967881 -1.348345 2 -2056153898 2238813396 0.321041 3 1250077441 1250077441 -0.689408 4 -474866958 3820100336 0.875903 5 -1117447608 3177519686 -1.036908 6 -610828462 3684138832 -0.749757 7 -1143879504 3151087790 -2.633581 8 -632459186 3662508108 -2.335211 9 -52590672 4242376622 -0.900580 10 -920365316 3374601978 -0.547212 ziggurat_test04() r4_exp() returns pseudorandom exponentially distributed real numbers between 0 and 1. 0 123456789 1 -1579999413 2714967881 0.863714 2 -2056153898 2238813396 0.171523 3 1250077441 1250077441 1.210186 4 -474866958 3820100336 1.645999 5 -1117447608 3177519686 0.624704 6 -610828462 3684138832 1.580322 7 -1143879504 3151087790 5.671516 8 -632459186 3662508108 4.687097 9 -52590672 4242376622 1.807592 10 -920365316 3374601978 1.014538 0 987654321 1 248404469 248404469 1.270067 2 2078538413 2078538413 1.963802 3 -1837283784 2457683510 0.419548 4 1841886731 1841886731 0.002059 5 736910199 736910199 1.498962 6 -1327526408 2967440886 0.396716 7 1833747846 1833747846 0.382060 8 1249087608 1249087608 3.917817 9 1286162813 1286162813 1.128685 10 -715281280 3579686014 1.563694 0 123456789 1 -1579999413 2714967881 0.863714 2 -2056153898 2238813396 0.171523 3 1250077441 1250077441 1.210186 4 -474866958 3820100336 1.645999 5 -1117447608 3177519686 0.624704 6 -610828462 3684138832 1.580322 7 -1143879504 3151087790 5.671516 8 -632459186 3662508108 4.687097 9 -52590672 4242376622 1.807592 10 -920365316 3374601978 1.014538 ziggurat_test05() Measure the time it takes shr3() to generate 100000 unsigned 32 bit integers. 8.511843 seconds. ziggurat_test06(): Measure the time it takes r4_uni() to generate 100000 uniform deviates. 10.251729 seconds. ziggurat_test07(): Measure the time it takes r4_nor() to generate 100000 normal deviates. 14.248554 seconds. ziggurat_test08(): Measure the time it takes r4_exp() to generate 100000 exponential deviates. 14.704635 seconds. ziggurat_test(): Normal end of execution. 02-May-2023 08:30:12