mirror of
https://github.com/jart/cosmopolitan.git
synced 2025-02-01 03:53:33 +00:00
f8b9bd2b47
Things are a little better. The LLD that comes with Linux seems to work. Old versions like LLVM 8 haven't been supported since Cosmopolitan v0.2. Running Clang on Windows with --target=x86_64-pc-linux-gnu doesn't seem to work. It has something to do with the recently added .zip section in the linker script. But even if that's removed, LLD on Windows thinks it is building an EFI application for some reason. Linker scripts are such a brittle house of cards, even for just ld.bfd alone.. We should just find a way to run our one true musl-cross-make linux gcc toolchain under Blinkenlights on non-Linux because GCC and Clang are so nondeterministic, inconsistent, and unreproducible when built for other operating systems. We need an actually portable compiler/linker that'll always behave the same way no matter what. See #180
242 lines
8.1 KiB
C
242 lines
8.1 KiB
C
/*-*- mode:c;indent-tabs-mode:nil;c-basic-offset:2;tab-width:8;coding:utf-8 -*-│
|
|
│vi: set net ft=c ts=2 sts=2 sw=2 fenc=utf-8 :vi│
|
|
╞══════════════════════════════════════════════════════════════════════════════╡
|
|
│ Copyright 2021 Justine Alexandra Roberts Tunney │
|
|
│ │
|
|
│ Permission to use, copy, modify, and/or distribute this software for │
|
|
│ any purpose with or without fee is hereby granted, provided that the │
|
|
│ above copyright notice and this permission notice appear in all copies. │
|
|
│ │
|
|
│ THE SOFTWARE IS PROVIDED "AS IS" AND THE AUTHOR DISCLAIMS ALL │
|
|
│ WARRANTIES WITH REGARD TO THIS SOFTWARE INCLUDING ALL IMPLIED │
|
|
│ WARRANTIES OF MERCHANTABILITY AND FITNESS. IN NO EVENT SHALL THE │
|
|
│ AUTHOR BE LIABLE FOR ANY SPECIAL, DIRECT, INDIRECT, OR CONSEQUENTIAL │
|
|
│ DAMAGES OR ANY DAMAGES WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR │
|
|
│ PROFITS, WHETHER IN AN ACTION OF CONTRACT, NEGLIGENCE OR OTHER │
|
|
│ TORTIOUS ACTION, ARISING OUT OF OR IN CONNECTION WITH THE USE OR │
|
|
│ PERFORMANCE OF THIS SOFTWARE. │
|
|
╚─────────────────────────────────────────────────────────────────────────────*/
|
|
#include "libc/log/check.h"
|
|
#include "libc/math.h"
|
|
#include "libc/nexgen32e/x86feature.h"
|
|
#include "libc/rand/lcg.internal.h"
|
|
#include "libc/rand/rand.h"
|
|
#include "libc/stdio/stdio.h"
|
|
#include "libc/sysv/consts/grnd.h"
|
|
#include "libc/testlib/testlib.h"
|
|
|
|
/* JustReturnZero */
|
|
/* entropy: 0 */
|
|
/* chi-square: 2.55e+07 */
|
|
/* chi-square percent: 0 */
|
|
/* mean: 0 */
|
|
/* monte-carlo-pi: 27.324 */
|
|
/* serial-correlation: -100000 */
|
|
|
|
/* JustIncrement */
|
|
/* entropy: 2.63951 */
|
|
/* chi-square: 1.443e+07 */
|
|
/* chi-square percent: 0 */
|
|
/* mean: 18.8803 */
|
|
/* monte-carlo-pi: 27.324 */
|
|
/* serial-correlation: 0.0092003 */
|
|
|
|
/* UNIX Sixth Edition */
|
|
/* entropy: 8 */
|
|
/* chi-square: 0.1536 */
|
|
/* chi-square percent: 1 */
|
|
/* mean: 127.502 */
|
|
/* monte-carlo-pi: 3.4192 */
|
|
/* serial-correlation: -0.470645 */
|
|
|
|
/* UNIX Seventh Edition */
|
|
/* entropy: 7.99818 */
|
|
/* chi-square: 251.843 */
|
|
/* chi-square percent: 0.544128 */
|
|
/* mean: 127.955 */
|
|
/* monte-carlo-pi: 0.675703 */
|
|
/* serial-correlation: -0.00207669 */
|
|
|
|
/* KnuthLcg */
|
|
/* entropy: 7.99835 */
|
|
/* chi-square: 228.383 */
|
|
/* chi-square percent: 0.883476 */
|
|
/* mean: 127.1 */
|
|
/* monte-carlo-pi: 0.561935 */
|
|
/* serial-correlation: -0.0038954 */
|
|
|
|
/* rand64 */
|
|
/* entropy: 7.99832 */
|
|
/* chi-square: 233.267 */
|
|
/* chi-square percent: 0.831821 */
|
|
/* mean: 127.427 */
|
|
/* monte-carlo-pi: 0.0271532 */
|
|
/* serial-correlation: -0.00255319 */
|
|
|
|
/* Rand64LowByte */
|
|
/* entropy: 7.99798 */
|
|
/* chi-square: 278.344 */
|
|
/* chi-square percent: 0.150796 */
|
|
/* mean: 127.88 */
|
|
/* monte-carlo-pi: 0.00340573 */
|
|
/* serial-correlation: 0.00162231 */
|
|
|
|
/* GetRandomNoSystem */
|
|
/* entropy: 7.99819 */
|
|
/* chi-square: 249.743 */
|
|
/* chi-square percent: 0.58114 */
|
|
/* mean: 127.124 */
|
|
/* monte-carlo-pi: 0.293716 */
|
|
/* serial-correlation: 0.00198516 */
|
|
|
|
/* GetRandomNoRdrrnd */
|
|
/* entropy: 7.99816 */
|
|
/* chi-square: 254.797 */
|
|
/* chi-square percent: 0.491811 */
|
|
/* mean: 127.308 */
|
|
/* monte-carlo-pi: 0.0118738 */
|
|
/* serial-correlation: 0.000197669 */
|
|
|
|
/* GetRandom */
|
|
/* entropy: 7.99808 */
|
|
/* chi-square: 266.737 */
|
|
/* chi-square percent: 0.294131 */
|
|
/* mean: 127.178 */
|
|
/* monte-carlo-pi: 0.0577122 */
|
|
/* serial-correlation: 0.00598793 */
|
|
|
|
typedef uint64_t (*random_f)(void);
|
|
|
|
static uint32_t randx = 1;
|
|
|
|
uint64_t JustReturnZero(void) {
|
|
return 0;
|
|
}
|
|
|
|
uint64_t JustIncrement(void) {
|
|
static uint64_t x;
|
|
return x++;
|
|
}
|
|
|
|
uint16_t SixthEditionRand(void) {
|
|
static int16_t gorp;
|
|
gorp = (gorp + 625) & 077777;
|
|
return gorp;
|
|
}
|
|
|
|
uint64_t SixthEditionLowByte(void) {
|
|
unsigned i;
|
|
uint64_t x;
|
|
for (x = i = 0; i < 8; ++i) {
|
|
x <<= 8;
|
|
x |= SixthEditionRand() & 255;
|
|
}
|
|
return x;
|
|
}
|
|
|
|
uint32_t SeventhEditionRand(void) {
|
|
return ((randx = randx * 1103515245 + 12345) >> 16) & 077777;
|
|
}
|
|
|
|
uint64_t SeventhEditionLowByte(void) {
|
|
unsigned i;
|
|
uint64_t x;
|
|
for (x = i = 0; i < 8; ++i) {
|
|
x <<= 8;
|
|
x |= SeventhEditionRand() & 255;
|
|
}
|
|
return x;
|
|
}
|
|
|
|
uint64_t KnuthLcg(void) {
|
|
unsigned i;
|
|
uint64_t x;
|
|
for (x = i = 0; i < 8; ++i) {
|
|
x <<= 8;
|
|
x |= rand() & 255;
|
|
}
|
|
return x;
|
|
}
|
|
|
|
uint64_t Rand64LowByte(void) {
|
|
unsigned i;
|
|
uint64_t x;
|
|
for (x = i = 0; i < 8; ++i) {
|
|
x <<= 8;
|
|
x |= rand64() & 255;
|
|
}
|
|
return x;
|
|
}
|
|
|
|
uint64_t GetRandomNoRdrrnd(void) {
|
|
uint64_t x;
|
|
ASSERT_EQ(8, getrandom(&x, 8, GRND_NORDRND));
|
|
return x;
|
|
}
|
|
|
|
uint64_t GetRandomNoSystem(void) {
|
|
uint64_t x;
|
|
if (X86_HAVE(RDRND) || X86_HAVE(RDSEED)) {
|
|
ASSERT_EQ(8, getrandom(&x, 8, GRND_NOSYSTEM));
|
|
} else {
|
|
ASSERT_EQ(8, getrandom(&x, 8, 0));
|
|
}
|
|
return x;
|
|
}
|
|
|
|
uint64_t GetRandom(void) {
|
|
uint64_t x;
|
|
ASSERT_EQ(8, getrandom(&x, 8, 0));
|
|
return x;
|
|
}
|
|
|
|
static const struct RandomFunction {
|
|
const char *s;
|
|
random_f f;
|
|
bool r;
|
|
} kRandomFunctions[] = {
|
|
{"JustReturnZero", JustReturnZero, false}, //
|
|
{"JustIncrement", JustIncrement, false}, //
|
|
{"SixthEditionLowByte", SixthEditionLowByte, false}, //
|
|
{"SeventhEditionLowByte", SeventhEditionLowByte, false}, //
|
|
{"KnuthLcg", KnuthLcg, false}, //
|
|
{"rand64", rand64, true}, //
|
|
{"Rand64LowByte", Rand64LowByte, true}, //
|
|
{"GetRandomNoRdrrnd", GetRandomNoRdrrnd, true}, //
|
|
{"GetRandomNoSystem", GetRandomNoSystem, true}, //
|
|
{"GetRandom", GetRandom, true}, //
|
|
};
|
|
|
|
TEST(getrandom, sanityTest) {
|
|
uint64_t q;
|
|
size_t i, j, k;
|
|
double montepi, chip, scc, ent, mean, chisq;
|
|
for (k = 0; k < 1; ++k) {
|
|
for (j = 0; j < ARRAYLEN(kRandomFunctions); ++j) {
|
|
rt_init(0);
|
|
for (i = 0; i + 8 <= 100000; i += 8) {
|
|
q = kRandomFunctions[j].f();
|
|
rt_add(&q, 8);
|
|
}
|
|
rt_end(&ent, &chisq, &mean, &montepi, &scc);
|
|
chip = pochisq(chisq, 255);
|
|
#if 0
|
|
fprintf(stderr, "\n");
|
|
fprintf(stderr, "/* %-32s */\n", kRandomFunctions[j].s);
|
|
fprintf(stderr, "/* entropy: %-12g */\n", ent);
|
|
fprintf(stderr, "/* chi-square: %-12g */\n", chisq);
|
|
fprintf(stderr, "/* chi-square percent: %-12g */\n", chip);
|
|
fprintf(stderr, "/* mean: %-12g */\n", mean);
|
|
fprintf(stderr, "/* monte-carlo-pi: %-12g */\n",
|
|
100 * fabs(M_PI - montepi) / M_PI);
|
|
fprintf(stderr, "/* serial-correlation: %-12g */\n", scc);
|
|
#endif
|
|
if (kRandomFunctions[j].r) {
|
|
CHECK_GE(chisq, 180, "%s", kRandomFunctions[j].s);
|
|
CHECK_GE(ent * 10, 78, "%s", kRandomFunctions[j].s);
|
|
CHECK_LT(fabs(scc) * 100, 5, "%s", kRandomFunctions[j].s);
|
|
CHECK_LT(fabs(128 - mean), 3, "%s", kRandomFunctions[j].s);
|
|
}
|
|
}
|
|
}
|
|
}
|