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
|
/* See LICENSE file for copyright and license details. */
#include "internals.h"
#include <fcntl.h>
#include <unistd.h>
#ifndef FAST_RANDOM_PATHNAME
# define FAST_RANDOM_PATHNAME "/dev/urandom"
#endif
#ifndef SECURE_RANDOM_PATHNAME
# define SECURE_RANDOM_PATHNAME "/dev/random"
#endif
static void
zrand_get_random_bits(z_t r, size_t bits, int fd)
{
size_t read_total = 0, n, chars = CEILING_BITS_TO_CHARS(bits);
ssize_t read_just;
zahl_char_t mask = 1;
char *buf;
ENSURE_SIZE(r, chars);
buf = (char *)(r->chars);
for (n = chars * sizeof(zahl_char_t); n;) {
read_just = read(fd, buf + read_total, n);
if (read_just < 0)
FAILURE(errno);
read_total += (size_t)read_just;
n -= (size_t)read_just;
}
bits = BITS_IN_LAST_CHAR(bits);
mask <<= bits;
mask -= 1;
r->chars[chars - 1] &= mask;
for (n = chars; n--;) {
if (r->chars[n]) {
r->used = n + 1;
SET_SIGNUM(r, 1);
return;
}
}
SET_SIGNUM(r, 0);
}
void
zrand(z_t r, enum zranddev dev, enum zranddist dist, z_t n)
{
const char *pathname = 0;
size_t bits;
int fd;
switch (dev) {
case FAST_RANDOM:
pathname = FAST_RANDOM_PATHNAME;
break;
case SECURE_RANDOM:
pathname = SECURE_RANDOM_PATHNAME;
break;
default:
abort();
}
if (EXPECT(zzero(n), 0)) {
SET_SIGNUM(r, 0);
return;
}
fd = open(pathname, O_RDONLY);
if (fd < 0)
FAILURE(errno);
switch (dist) {
case QUASIUNIFORM:
if (EXPECT(zsignum(n) < 0, 0))
FAILURE(EDOM); /* n must be non-negative. */
bits = zbits(n);
zrand_get_random_bits(r, bits, fd);
zadd(r, r, libzahl_const_1);
zmul(r, r, n);
zrsh(r, r, bits);
break;
case UNIFORM:
if (EXPECT(zsignum(n) < 0, 0))
FAILURE(EDOM); /* n must be non-negative. */
bits = zbits(n);
do
zrand_get_random_bits(r, bits, fd);
while (zcmpmag(r, n) > 0);
break;
default:
abort();
}
close(fd);
}
|
