-
Data: 2021-10-03 23:44:10
Temat: Sekwencje Weyl w generatorach liczb pseudolosowych
Od: "o...@g...com" <o...@g...com> szukaj wiadomości tego autora
[ pokaż wszystkie nagłówki ]Definicja Weyl sequence jest następująca:
https://en.wikipedia.org/wiki/Weyl_sequence
W artykule podano liczbę 362437 do utworzenia 32-bitowej sekwencji. Ale ogółem
wystarczy liczba względnie pierwsza z modulusem.
Pytanie dlaczego wzięto akurat taką liczbę? Czy można wziąć też liczbę 3? Również
jest względnie pierwsza z 2^32. Ile jest takich liczb, które nadadzą się do
utworzenia takiej sekwencji? Czy może to być dowolna liczba nieparzysta?
Nie do końca rozumiem też sens stosowania tej sekwencji. Wiem, że George Marsaglia
użył ich do udoskonalenia generatora Xorshift. Z jakichś powodów te sekwencje lepiej
nadają się mieszania z wynikami generatora niż zwykły licznik: 1,2,3,..., który też
moglibyśmy zastosować. Ale właściwie co czyni je użytecznymi? Bo same w sobie nie
produkują wystarczająco losowych wyników. Te sekwencje zostały też zastosowane do
ulepszenia generatora Middle Square:
https://arxiv.org/pdf/1704.00358.pdf
Wyniki sekwencji są po prostu dodawane do wyników klasyczne Middle Square. Ale nie do
końca rozumiem jaka teoria za tym stoi i znowu - dlaczego w publikacji przyjęto
akurat taką, a nie inną liczbę względnie pierwszą z modulusem? Zarzuty do Midlle
Square dobrze podsumowała Mellisa O'Neil:
https://www.pcg-random.org/posts/too-big-to-fail.htm
l
Generator ten tworzy random mapping, a w takim odwzorowaniu cykle są stosunkowo
krótkie i mamy bodaj niezerowe prawdopodobieństwo natknięcia się na jakiś fatalny,
krótki cykl (co dla generatora liczb pseudolosowych jest dyskwalifikujące). Widynski
zaproponował więc mieszanie wyników Middle Square z sekwencjami Weyl, co wydłuża cykl
całego generatora. Tylko jak to się dzieje, że zostaje on wydłużony.
Sam pracuję nad pewnym generatorem liczb pseudolosowych, który też działa jak random
mapping i 32-bitowa wersja ze względu oczekiwaną długość cyklu ~ 41067 () nie zdaje
testów. Ale dodawanie do wyników Weyl sequence (stała z artykułu 362437) poprawia
wyniki. Nie rozumiem jednak dlaczego.
Następne wpisy z tego wątku
- 03.10.21 23:55 o...@g...com
- 04.10.21 00:24 o...@g...com
Najnowsze wątki z tej grupy
- Bibl. Qt jest sztucznie ograniczona - jest nieprzydatna do celów komercyjnych
- Co sciaga kretynow
- AEiC 2024 - Ada-Europe conference - Deadlines Approaching
- Jakie są dobre zasady programowania programów opartych na wtyczkach?
- sprawdzanie słów kluczowych dot. zła
- Re: W czym sie teraz pisze programy??
- Re: (PDF) Surgical Pathology of Non-neoplastic Gastrointestinal Diseases by Lizhi Zhang
- CfC 28th Ada-Europe Int. Conf. Reliable Software Technologies
- Młodzi programiści i tajna policja
- Ada 2022 Language Reference Manual to be Published by Springer
- Press Release - AEiC 2023, Ada-Europe Reliable Softw. Technol.
- Ada-Europe - AEiC 2023 early registration deadline approaching
- Ada-Europe Int.Conf. Reliable Software Technologies, AEiC 2023
- Ile cykli zajmuje mnożenie liczb 64-bitowych?
- Ideologia Polskiego Programisty wer.3
Najnowsze wątki
- 2024-05-02 Silnik Stirlinga - energia mechaniczna wytwarzana z ciepła bez dostarczania paliwa!
- 2024-05-02 Kraków => Senior .Net Programmer <=
- 2024-05-02 Sprowadzenie pożaru zagrażającego mieniu w wielkich rozmiarch
- 2024-05-02 Zielona Góra => Spedytor międzynarodowy <=
- 2024-05-02 Katowice => Performance Analyst <=
- 2024-05-02 Gdańsk => Full Stack web developer (obszar .Net Core, Angular6+) <=
- 2024-05-02 Gdańsk => Business Development Manager - obszar bezpieczeństwa IT <=
- 2024-05-02 Warszawa => Starszy inżynier oprogramowania (Rust) <=
- 2024-05-02 Katowice => Java Developer <=
- 2024-05-02 'Studenci do nauki!'
- 2024-05-02 Dlaczego energia elektryczna jest droga i cały czas drożeje
- 2024-05-02 Dlaczego energia elektryczna jest droga i cały czas drożeje
- 2024-05-02 Dlaczego energia elektryczna jest droga i cały czas drożeje
- 2024-05-02 Wrocław => Lider Serwisu/Programista PLC z językiem niemieckim <=
- 2024-05-02 Łódź => Senior DevOps Engineer <=