stats-claw 0.1.0

Data science on the hot path: in-process, zero-dependency statistical computing for Rust (distributions, hypothesis tests, resampling) validated against scipy
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113

// @generated ziggurat layer tables for the standard Normal (Doornik 2005 form).
//
// Equal-area 128-layer tables. `ZIG_K[i]` is the 31-bit acceptance threshold,
// `ZIG_W[i]` scales the magnitude to the layer's x, `ZIG_F[i]` is the density at
// the layer edge. Regenerate via the equal-area layer construction (Doornik 2005);
// verified by `ziggurat::tests::table_construction_invariant_holds`.

/// 31-bit acceptance thresholds per layer.
#[allow(clippy::unreadable_literal)]
static ZIG_K: [u32; 128] = [
    1991057938, 0, 1611602771, 1826899878, 1918584482, 1969227037, 2001281515, 2023368125,
    2039498179, 2051788381, 2061460127, 2069267110, 2075699398, 2081089314, 2085670119,
    2089610331, 2093034710, 2096037586, 2098691595, 2101053571, 2103168620, 2105072996,
    2106796166, 2108362327, 2109791536, 2111100552, 2112303493, 2113412330, 2114437283,
    2115387130, 2116269447, 2117090813, 2117856962, 2118572919, 2119243101, 2119871411,
    2120461303, 2121015852, 2121537798, 2122029592, 2122493434, 2122931299, 2123344971,
    2123736059, 2124106020, 2124456175, 2124787725, 2125101763, 2125399283, 2125681194,
    2125948325, 2126201433, 2126441213, 2126668298, 2126883268, 2127086657, 2127278949,
    2127460589, 2127631985, 2127793506, 2127945490, 2128088244, 2128222044, 2128347141,
    2128463758, 2128572095, 2128672327, 2128764606, 2128849065, 2128925811, 2128994934,
    2129056501, 2129110560, 2129157136, 2129196237, 2129227847, 2129251929, 2129268426,
    2129277255, 2129278312, 2129271467, 2129256561, 2129233410, 2129201800, 2129161480,
    2129112170, 2129053545, 2128985244, 2128906855, 2128817916, 2128717911, 2128606255,
    2128482298, 2128345305, 2128194452, 2128028813, 2127847342, 2127648860, 2127432031,
    2127195339, 2126937058, 2126655214, 2126347546, 2126011445, 2125643893, 2125241376,
    2124799783, 2124314271, 2123779094, 2123187386, 2122530867, 2121799464, 2120980787,
    2120059418, 2119015917, 2117825402, 2116455471, 2114863093, 2112989789, 2110753906,
    2108037662, 2104664315, 2100355223, 2094642347, 2086670106, 2074676188, 2054300022,
    2010539237,
];
/// Magnitude-to-x scale per layer (edge = `ZIG_W[i] * 2^31`).
#[allow(clippy::unreadable_literal)]
static ZIG_W: [f64; 128] = [
    1.729040521542798e-09, 1.2680928447002762e-10, 1.689751777318455e-10,
    1.9862688442479051e-10, 2.2232431792499955e-10, 2.424493612544893e-10,
    2.6016131900632064e-10, 2.7611988711703956e-10, 2.907396281771598e-10,
    3.0429970414376596e-10, 3.1699795213954273e-10, 3.2898020527113064e-10,
    3.4035738121834064e-10, 3.512160221366471e-10, 3.616250995056517e-10,
    3.7164057634959785e-10, 3.813085643110598e-10, 3.906675680994882e-10,
    3.997501186997691e-10, 4.0858398615984403e-10, 4.1719309640160654e-10,
    4.2559823534592626e-10, 4.3381759739255105e-10, 4.418672181252886e-10,
    4.497613196266582e-10, 4.5751258894588287e-10, 4.65132404814001e-10,
    4.726310238481176e-10, 4.800177347232567e-10, 4.873009867798748e-10,
    4.944884980538973e-10, 5.015873466119616e-10, 5.08604048242456e-10,
    5.15544622919539e-10, 5.224146519706316e-10, 5.292193275006305e-10,
    5.35963495331289e-10, 5.426516924820619e-10, 5.492881800346021e-10,
    5.558769720760773e-10, 5.624218612983588e-10, 5.68926441734655e-10,
    5.753941290375603e-10, 5.818281786390898e-10, 5.88231702081217e-10,
    5.946076817624996e-10, 6.009589843108302e-10, 6.072883727627885e-10,
    6.135985177054135e-10, 6.198920075155922e-10, 6.261713578149429e-10,
    6.324390202435402e-10, 6.386973906435736e-10, 6.449488167337383e-10,
    6.511956053464698e-10, 6.574400292928599e-10, 6.636843339139875e-10,
    6.699307433723302e-10, 6.761814667327444e-10, 6.824387038791137e-10,
    6.887046513100733e-10, 6.949815078551667e-10, 7.012714803513155e-10,
    7.07576789318556e-10, 7.138996746735849e-10, 7.202424015197486e-10,
    7.266072660527047e-10, 7.329966016220864e-10, 7.394127849911228e-10,
    7.458582428383539e-10, 7.523354585483488e-10, 7.588469793417652e-10,
    7.653954237992263e-10, 7.7198348983844e-10, 7.786139632098381e-10,
    7.852897265828997e-10, 7.920137693034098e-10, 7.987891979113536e-10,
    8.05619247520217e-10, 8.125072941713968e-10, 8.194568682925745e-10,
    8.264716694066625e-10, 8.335555822587845e-10, 8.407126945532991e-10,
    8.479473165218372e-10, 8.552640025776094e-10, 8.626675753519363e-10,
    8.701631524574424e-10, 8.777561763803284e-10, 8.854524479737278e-10,
    8.932581641080369e-10, 9.011799601356605e-10, 9.092249579511381e-10,
    9.174008205786005e-10, 9.257158144040126e-10, 9.341788803988472e-10,
    9.427997159666314e-10, 9.515888693998883e-10, 9.605578493831253e-10,
    9.697192525453944e-10, 9.7908691279089e-10, 9.886760770687724e-10,
    9.985036134535425e-10, 1.0085882589914473e-09, 1.0189509168621382e-09,
    1.0296150152006668e-09, 1.0406069436999874e-09, 1.0519565892728039e-09,
    1.0636979991930871e-09, 1.0758702101645819e-09, 1.0885182960607283e-09,
    1.1016947078135044e-09, 1.1154610095597163e-09, 1.1298901613493216e-09,
    1.1450695700067237e-09, 1.1611052426022348e-09, 1.178127560945613e-09,
    1.1962995053850756e-09, 1.2158286983295564e-09, 1.2369856290804966e-09,
    1.2601323300608525e-09, 1.2857696844205153e-09, 1.3146201849677183e-09,
    1.3477839562210855e-09, 1.3870635315067043e-09, 1.435740319181638e-09,
    1.5008659030222993e-09, 1.6030947938091123e-09,
];
/// Normal density at each layer edge.
#[allow(clippy::unreadable_literal)]
static ZIG_F: [f64; 128] = [
    1.0, 0.9635996931270862, 0.9362826816850596, 0.9130436479717402, 0.8922816507840261,
    0.8732430489100695, 0.8555006078694506, 0.8387836052959896, 0.822907211381409,
    0.8077382946829605, 0.7931770117713051, 0.7791460859296877, 0.7655841738977045,
    0.7524415591746114, 0.7396772436726473, 0.7272569183441848, 0.7151515074104986,
    0.7033360990161581, 0.6917891434366751, 0.6804918409973341, 0.6694276673488904,
    0.658582000050088, 0.6479418211102225, 0.6374954773350423, 0.6272324852499273,
    0.6171433708188809, 0.6072195366251203, 0.5974531509445167, 0.5878370544347066,
    0.5783646811197631, 0.5690299910679509, 0.5598274127040869, 0.5507517931146045,
    0.5417983550254255, 0.5329626593838361, 0.5242405726729841, 0.5156282382440018,
    0.507122051075569, 0.4987186354709795, 0.4904148252838441, 0.4822076463294852,
    0.47409430069301695, 0.4660721526894561, 0.45813871626787206, 0.4502916436820392,
    0.44252871527546844, 0.4348478302499909, 0.4272469983049961, 0.4197243320495744,
    0.412278040102661, 0.40490642080722294, 0.3976078564938733, 0.3903808082373146,
    0.3832238110559012, 0.3761354695105626, 0.3691144536644722, 0.3621594953693176,
    0.3552693848479171, 0.3484429675463266, 0.3416791412315504, 0.3349768533135892,
    0.3283350983728503, 0.3217529158759849, 0.3152293880650109, 0.3087636380061811,
    0.30235482778648354, 0.296002156846933, 0.28970486044295984, 0.283462208223233,
    0.2772735029191881, 0.2711380791383846, 0.2650553022555892, 0.25902456739620483,
    0.25304529850732577, 0.2471169475123214, 0.24123899354543982, 0.23541094226347908,
    0.22963232523211613, 0.22390269938500842, 0.2182216465543054, 0.2125887730717303,
    0.20700370943992652, 0.20146611007431367, 0.19597565311627774, 0.19053204031913715,
    0.1851349970089922, 0.17978427212329545, 0.1744796383307895, 0.169220892237365,
    0.16400785468342038, 0.1588403711394793, 0.15371831220818166, 0.14864157424234226,
    0.14361008009062776, 0.1386237799845946, 0.13368265258343937, 0.1287867061959432,
    0.12393598020286782, 0.11913054670765083, 0.11437051244886601, 0.10965602101484027,
    0.10498725540942132, 0.10036444102865587, 0.09578784912173144, 0.09125780082683026,
    0.08677467189478018, 0.08233889824223566, 0.0779509825139734, 0.0736115018841134,
    0.06932111739357791, 0.06508058521306807, 0.060890770348040406, 0.05675266348104985,
    0.05266740190305101, 0.048636295859867805, 0.044660862200491425, 0.040742868074444175,
    0.0368843887866562, 0.03308788614622575, 0.02935631744000685, 0.02569329193593427,
    0.022103304615927098, 0.018592102737011288, 0.015167298010546568, 0.011839478657884862,
    0.008624484412859885, 0.005548995220771345, 0.002669629083880923,
];