1 /*
2  * Copyright (c) 1993 Martin Birgmeier
3  * All rights reserved.
4  *
5  * You may redistribute unmodified or modified versions of this source
6  * code provided that the above copyright notice and this and the
7  * following conditions are retained.
8  *
9  * This software is provided ``as is'', and comes with no warranties
10  * of any kind. I shall in no event be liable for anything that happens
11  * to anyone/anything when using this software.
12  */
13 
14 #include "rand48.h"
15 
16 double
_erand48_r(struct _rand48 * r,unsigned short xseed[3])17 _erand48_r (struct _rand48 *r,
18        unsigned short xseed[3])
19 {
20   __dorand48(r, xseed);
21   return ldexp((double) xseed[0], -48) +
22     ldexp((double) xseed[1], -32) +
23     ldexp((double) xseed[2], -16);
24 }
25 
26 #ifndef _REENT_ONLY
27 double
erand48(unsigned short xseed[3])28 erand48 (unsigned short xseed[3])
29 {
30   return _erand48_r (&_rand48, xseed);
31 }
32 #endif /* !_REENT_ONLY */
33