Re: test randomizacyjny w C

Article View: pl.comp.lang.c

Article #308758

#308758

From: bartekltg
Date: Mon, 26 Sep 2016 21:17

149 lines
8700 bytes

On 26.09.2016 22:33, Borneq wrote:
> W dniu 26.09.2016 o 21:36, bartekltg pisze:
>> Nie, nie przechodzi testów; mt19937 przechodzi testy;p
>> dieharder -g 201 -f random10.dat  -a
>
> Ranlux w każdej wersji LUX nie przechodzi diehard_operm5, co więcej -
> fizyczny random nie przechodzi diehard_opso, podczas gdy mt19937
> przechodzi wszystkie ? Myślałem że Ranlux lepiej się spisuje niż mt19937.

Coś zepsułeś.
[Co? Któż to wie, znów nie napisałeś, _co_ dokąłdnie zrobiłeś.]

Sprawadźmy te ranluxy.

for i in {43..49}; do dieharder -g $i -d 1; echo;echo;  done
#=============================================================================#
#            dieharder version 3.31.1 Copyright 2003 Robert G. Brown 
      #
#=============================================================================#
    rng_name    |rands/second|   Seed   |
          ranlux|  1.78e+07  |2775139678|
#=============================================================================#
         test_name   |ntup| tsamples |psamples|  p-value |Assessment
#=============================================================================#
       diehard_operm5|   0|   1000000|     100|0.31240470|  PASSED


#=============================================================================#
#            dieharder version 3.31.1 Copyright 2003 Robert G. Brown 
      #
#=============================================================================#
    rng_name    |rands/second|   Seed   |
       ranlux389|  1.07e+07  |1235841629|
#=============================================================================#
         test_name   |ntup| tsamples |psamples|  p-value |Assessment
#=============================================================================#
       diehard_operm5|   0|   1000000|     100|0.85423943|  PASSED


#=============================================================================#
#            dieharder version 3.31.1 Copyright 2003 Robert G. Brown 
      #
#=============================================================================#
    rng_name    |rands/second|   Seed   |
         ranlxd1|  1.16e+07  |3090234114|
#=============================================================================#
         test_name   |ntup| tsamples |psamples|  p-value |Assessment
#=============================================================================#
       diehard_operm5|   0|   1000000|     100|0.99705872|   WEAK


#=============================================================================#
#            dieharder version 3.31.1 Copyright 2003 Robert G. Brown 
      #
#=============================================================================#
    rng_name    |rands/second|   Seed   |
         ranlxd2|  6.25e+06  |3221774477|
#=============================================================================#
         test_name   |ntup| tsamples |psamples|  p-value |Assessment
#=============================================================================#
       diehard_operm5|   0|   1000000|     100|0.92391596|  PASSED


#=============================================================================#
#            dieharder version 3.31.1 Copyright 2003 Robert G. Brown 
      #
#=============================================================================#
    rng_name    |rands/second|   Seed   |
         ranlxs0|  2.82e+07  |1766173486|
#=============================================================================#
         test_name   |ntup| tsamples |psamples|  p-value |Assessment
#=============================================================================#
       diehard_operm5|   0|   1000000|     100|0.60530086|  PASSED


#=============================================================================#
#            dieharder version 3.31.1 Copyright 2003 Robert G. Brown 
      #
#=============================================================================#
    rng_name    |rands/second|   Seed   |
         ranlxs1|  1.93e+07  |3269615808|
#=============================================================================#
         test_name   |ntup| tsamples |psamples|  p-value |Assessment
#=============================================================================#
       diehard_operm5|   0|   1000000|     100|0.03864195|  PASSED


#=============================================================================#
#            dieharder version 3.31.1 Copyright 2003 Robert G. Brown 
      #
#=============================================================================#
    rng_name    |rands/second|   Seed   |
         ranlxs2|  1.13e+07  |1751257048|
#=============================================================================#
         test_name   |ntup| tsamples |psamples|  p-value |Assessment
#=============================================================================#
       diehard_operm5|   0|   1000000|     100|0.41419599|  PASSED



ranlxd1 jest podejrzeny, ale miałem 0.6% sanszy na każym generratorze
uzyskać taki wynik przez przypadek, więć 1-(1-0.006)^7= 4% uzyskać
go rpzypadkiem.

Powtórzmy, ze trzy razy:


for i in {1..3}; do dieharder -g 45 -d 1; echo;echo;  done
#=============================================================================#
#            dieharder version 3.31.1 Copyright 2003 Robert G. Brown 
      #
#=============================================================================#
    rng_name    |rands/second|   Seed   |
         ranlxd1|  1.14e+07  |2458088494|
#=============================================================================#
         test_name   |ntup| tsamples |psamples|  p-value |Assessment
#=============================================================================#
       diehard_operm5|   0|   1000000|     100|0.21088550|  PASSED


#=============================================================================#
#            dieharder version 3.31.1 Copyright 2003 Robert G. Brown 
      #
#=============================================================================#
    rng_name    |rands/second|   Seed   |
         ranlxd1|  1.15e+07  |4264844718|
#=============================================================================#
         test_name   |ntup| tsamples |psamples|  p-value |Assessment
#=============================================================================#
       diehard_operm5|   0|   1000000|     100|0.58660050|  PASSED


#=============================================================================#
#            dieharder version 3.31.1 Copyright 2003 Robert G. Brown 
      #
#=============================================================================#
    rng_name    |rands/second|   Seed   |
         ranlxd1|  1.14e+07  |2473325856|
#=============================================================================#
         test_name   |ntup| tsamples |psamples|  p-value |Assessment
#=============================================================================#
       diehard_operm5|   0|   1000000|     100|0.31103629|  PASSED


Wszystkie generatory z rodziny ranlux przechodzą diehard_operm5 na
standardowych ustawieniach.

pzdr
bartekltg

Technical Details

Message-ID: <nsc38t$5t0$1@node1.news.atman.pl>

Path: mbox-mo1.pugleaf.net!news-archive.icm.edu.pl!news.icm.edu.pl!newsfeed2.atman.pl!newsfeed.atman.pl!.POSTED!not-for-mail

References:

<5777f44e$0$12548$65785112@news.neostrada.pl> <nlcdr5$rk4$1@portraits.wsisiz.edu.pl> <3fabf491-7526-4652-b460-a6c2b03318fe@googlegroups.com> <52620181-d1a6-4812-96a5-38eaed8f0860@googlegroups.com> <nlecfp$94i$1@portraits.wsisiz.edu.pl> <771ccedb-ab65-484c-95f6-fed3beab850d@googlegroups.com> <nlgdvu$7vo$1@node2.news.atman.pl> <nsbfoj$gk0$1@node1.news.atman.pl> <1121cd99-4f0c-42c9-af4e-bbb95846b221@googlegroups.com> <nsbi91$118$1@node2.news.atman.pl> <nsbja8$k91$1@node1.news.atman.pl> <nsbr2o$a64$1@node2.news.atman.pl> <nsbtc8$ch5$1@node2.news.atman.pl> <nsc0nj$fs7$1@node2.news.atman.pl>

Show Full Headers

Complete Article Headers

Path: mbox-mo1.pugleaf.net!news-archive.icm.edu.pl!news.icm.edu.pl!newsfeed2.atman.pl!newsfeed.atman.pl!.POSTED!not-for-mail
From: bartekltg <bartekltg@gmail.com>
Newsgroups: pl.comp.lang.c
Subject: Re: test randomizacyjny w C
Date: Mon, 26 Sep 2016 23:17:16 +0200
Organization: ATMAN - ATM S.A.
Lines: 153
Message-ID: <nsc38t$5t0$1@node1.news.atman.pl>
References: <5777f44e$0$12548$65785112@news.neostrada.pl> <nlcdr5$rk4$1@portraits.wsisiz.edu.pl> <3fabf491-7526-4652-b460-a6c2b03318fe@googlegroups.com> <52620181-d1a6-4812-96a5-38eaed8f0860@googlegroups.com> <nlecfp$94i$1@portraits.wsisiz.edu.pl> <771ccedb-ab65-484c-95f6-fed3beab850d@googlegroups.com> <nlgdvu$7vo$1@node2.news.atman.pl> <nsbfoj$gk0$1@node1.news.atman.pl> <1121cd99-4f0c-42c9-af4e-bbb95846b221@googlegroups.com> <nsbi91$118$1@node2.news.atman.pl> <nsbja8$k91$1@node1.news.atman.pl> <nsbr2o$a64$1@node2.news.atman.pl> <nsbtc8$ch5$1@node2.news.atman.pl> <nsc0nj$fs7$1@node2.news.atman.pl>
NNTP-Posting-Host: 89-70-119-159.dynamic.chello.pl
Mime-Version: 1.0
X-Trace: node1.news.atman.pl 1474924637 6048 89.70.119.159 (26 Sep 2016 21:17:17 GMT)
X-Complaints-To: usenet@atman.pl
NNTP-Posting-Date: Mon, 26 Sep 2016 21:17:17 +0000 (UTC)
User-Agent: Mozilla/5.0 (X11; Linux x86_64; rv:45.0) Gecko/20100101 Thunderbird/45.2.0
In-Reply-To: <nsc0nj$fs7$1@node2.news.atman.pl>
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

Reply to Parent Message