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
114
115
116
117
118
119
|
/*
urandom.c:
Copyright (C) 2010 by John ffitch
This file is part of Csound.
The Csound Library is free software; you can redistribute it
and/or modify it under the terms of the GNU Lesser General Public
License as published by the Free Software Foundation; either
version 2.1 of the License, or (at your option) any later version.
Csound is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
GNU Lesser General Public License for more details.
You should have received a copy of the GNU Lesser General Public
License along with Csound; if not, write to the Free Software
Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA
02110-1301 USA
*/
#include "csdl.h"
//#include <ieee754.h>
#ifdef __HAIKU__
#include <fcntl.h>
#endif
#ifdef MACOSX
#include <unistd.h>
#endif
typedef struct {
OPDS h;
MYFLT *ar; /* Output */
MYFLT *imin;
MYFLT *imax;
int32_t ur;
MYFLT mul;
MYFLT add;
} URANDOM;
static int32_t urand_deinit(CSOUND *csound, URANDOM *p)
{
IGN(csound);
close(p->ur);
return OK;
}
static int32_t urand_init(CSOUND *csound, URANDOM *p)
{
int32_t ur = open("/dev/urandom", O_RDONLY);
if (UNLIKELY(ur<0)) return NOTOK;
p->ur = ur;
csound->RegisterDeinitCallback(csound, p,
(int32_t (*)(CSOUND *, void *)) urand_deinit);
p->mul = FL(0.5)*(*p->imax - *p->imin);
p->add = FL(0.5)*(*p->imax + *p->imin);
return OK;
}
static int32_t urand_run(CSOUND *csound, URANDOM *p)
{
IGN(csound);
int32_t ur = p->ur;
/* union ieee754_double x; */
int64_t x;
if (UNLIKELY(read(ur, &x, sizeof(int64_t))!=sizeof(int64_t))) return NOTOK;
/* x.ieee.exponent = x.ieee.exponent& 0x377; */
/* printf("Debug: %s(%d): %g %d %03x %05x %08x\n", __FILE__, __LINE__, x.d, */
/* x.ieee.negative, x.ieee.exponent, x.ieee.mantissa0, x.ieee.mantissa1); */
*p->ar = p->mul *((MYFLT)x/(MYFLT)INT64_MAX) + p->add;
return OK;
}
static int32_t urand_irate(CSOUND *csound, URANDOM *p)
{
if (LIKELY(urand_init(csound,p)==OK)) return urand_run(csound,p);
else return NOTOK;
}
static int32_t urand_arun(CSOUND *csound, URANDOM *p)
{
IGN(csound);
int32_t ur = p->ur;
/* union ieee754_double x; */
int64_t x;
MYFLT *ar = p->ar;
uint32_t offset = p->h.insdshead->ksmps_offset;
uint32_t early = p->h.insdshead->ksmps_no_end;
uint32_t n, nsmps = CS_KSMPS;
if (UNLIKELY(offset)) memset(ar, '\0', offset*sizeof(MYFLT));
if (UNLIKELY(early)) {
nsmps -= early;
memset(&ar[nsmps], '\0', early*sizeof(MYFLT));
}
for (n=offset; n<nsmps; n++) {
if (UNLIKELY(read(ur, &x, sizeof(int64_t))!= sizeof(int64_t))) return NOTOK;
ar[n] = p->mul *((MYFLT)x/(MYFLT)INT64_MAX) + p->add;
}
return OK;
}
#define S(x) sizeof(x)
static OENTRY urandom_localops[] = {
{ "urandom.i", S(URANDOM), 0, 1, "i", "jp", (SUBR) urand_irate },
{ "urandom.k", S(URANDOM), 0, 3, "k", "jp", (SUBR) urand_init, (SUBR) urand_run},
{ "urandom.a", S(URANDOM), 0, 3, "a", "jp",
(SUBR) urand_init, (SUBR) urand_arun}
};
LINKAGE_BUILTIN(urandom_localops)
|
