Re: dieharder dla dużych wartości m, np. m=1000 - Grupy dyskusyjne w eGospodarka.pl

eGospodarka.pl › Grupy › pl.comp.programming › dieharder dla dużych wartości m, np. m=1000 › Re: dieharder dla dużych wartości m, np. m=1000

X-Received: by 10.31.180.138 with SMTP id d132mr834505vkf.9.1513187503986; Wed, 13
Dec 2017 09:51:43 -0800 (PST)
X-Received: by 10.31.180.138 with SMTP id d132mr834505vkf.9.1513187503986; Wed, 13
Dec 2017 09:51:43 -0800 (PST)
Path: news-archive.icm.edu.pl!news.icm.edu.pl!news.nask.pl!news.nask.org.pl!news.unit
0.net!peer01.am4!peer.am4.highwinds-media.com!peer03.iad!feed-me.highwinds-medi
a.com!news.highwinds-media.com!g35no425893qtk.1!news-out.google.com!v55ni255qtc
.0!nntp.google.com!m31no427256qtf.0!postnews.google.com!glegroupsg2000goo.googl
egroups.com!not-for-mail
Newsgroups: pl.comp.programming
Date: Wed, 13 Dec 2017 09:51:43 -0800 (PST)
In-Reply-To: <c...@g...com>
Complaints-To: g...@g...com
Injection-Info: glegroupsg2000goo.googlegroups.com; posting-host=31.0.33.98;
posting-account=xjvq9QoAAAATMPC2X3btlHd_LkaJo_rj
NNTP-Posting-Host: 31.0.33.98
References: <c...@g...com>
<c...@g...com>
User-Agent: G2/1.0
MIME-Version: 1.0
Message-ID: <b...@g...com>
Subject: Re: dieharder dla dużych wartości m, np. m=1000
From: "M.M." <m...@g...com>
Injection-Date: Wed, 13 Dec 2017 17:51:44 +0000
Content-Type: text/plain; charset="UTF-8"
Content-Transfer-Encoding: base64
X-Received-Body-CRC: 3215231526
X-Received-Bytes: 6257
Xref: news-archive.icm.edu.pl pl.comp.programming:211872
[ ukryj nagłówki ]
On Wednesday, December 13, 2017 at 6:38:16 PM UTC+1, bartekltg wrote:
> On Wednesday, December 13, 2017 at 3:03:54 PM UTC+1, M.M. wrote:
> > Dieharder chyba pada na dużych wartościach m.
> >
> > Testu dnia urodzin nie przechodzi generator ranlux48, ani mój generator
> > fibonacciego z poprzedniego wątku:
> >
> > FIBONACCI
> > time ./frnd | dieharder -d 0 -m 1000 -g 200
> > #===================================================
==========================#
> > # dieharder version 3.31.1 Copyright 2003 Robert G. Brown #
> > #===================================================
==========================#
> > rng_name |rands/second| Seed |
> > stdin_input_raw| 4.97e+06 |1997299816|
> > #===================================================
==========================#
> > test_name |ntup| tsamples |psamples| p-value |Assessment
> > #===================================================
==========================#
> > diehard_birthdays| 0| 100| 100000|0.00000014| FAILED
> >
> > real 34m8.515s
> > user 35m39.480s
> > sys 1m20.016s
> >
> >
> >
> > RANLUX48
> > time ./frnd | dieharder -d 0 -m 1000 -g 200
> > #===================================================
==========================#
> > # dieharder version 3.31.1 Copyright 2003 Robert G. Brown #
> > #===================================================
==========================#
> > rng_name |rands/second| Seed |
> > stdin_input_raw| 3.78e+05 |1785358803|
> > #===================================================
==========================#
> > test_name |ntup| tsamples |psamples| p-value |Assessment
> > #===================================================
==========================#
> > diehard_birthdays| 0| 100| 100000|0.00000013| FAILED
> >
> > real 93m39.541s
> > user 124m16.952s
> > sys 0m41.020s
> >
> > Nie wiem co jest grane, ale chyba brakuje precyzji dieharderowi, rozkłady
> > powinny być zaimplementowane chociaż na precyzji 1024 bitów.
> >
>
>
> Dodaj opcje -k 2
>
> Sprawdź w dokumentacji (tym pdfie z dieherder) czy ten test nie jest
> jednym z tych, które się wywalają dla zbyt duzych danych.
>
> pzdr
> bartekltg

Sprawdzałem, dodawałem. Faktycznie pomaga, niestety pomaga w bardzo
minimalnym stopniu.

Za to widzę, że w boost rozkład czi-kwadrat jest jako szablon zależny
od typu.

http://www.boost.org/doc/libs/1_38_0/libs/math/doc/s
f_and_dist/html/math_toolkit/dist/dist_ref/dists/chi
_squared_dist.html

Ciekawe czy można go zdefiniować na poniższym typie, a jak można, to czy to coś
pomaga:

http://www.boost.org/doc/libs/1_62_0/libs/multipreci
sion/doc/html/boost_multiprecision/tut/floats/cpp_de
c_float.html

Jest też gotowiec do odwrotnego chi-kwadrat

http://www.boost.org/doc/libs/1_52_0/libs/math/doc/s
f_and_dist/html/math_toolkit/dist/dist_ref/dists/inv
erse_chi_squared_dist.html

Jest też kilka przykładów użycia:
http://www.boost.org/doc/libs/1_38_0/libs/math/doc/s
f_and_dist/html/math_toolkit/dist/stat_tut/weg/cs_eg
.html

By trzeba się wczytać... jeśli to działa tak jak mi się wydaje, to napisanie
samemu (także dla dużej ilości punktów swobody) jest trywialne.

Ale np. to ( typ cpp_bin_float<N> ) u mnie się nawet nie kompiluje:

http://www.boost.org/doc/libs/1_62_0/libs/multipreci
sion/doc/html/boost_multiprecision/tut/floats/cpp_bi
n_float.html

Więc nie wiem czy jest tak kolorowo jak piszą :/

Pozdrawiam