xref: /dflybsd-src/contrib/gcc-8.0/libstdc++-v3/src/c++11/random.cc (revision 38fd149817dfbff97799f62fcb70be98c4e32523)
1*38fd1498Szrj // random -*- C++ -*-
2*38fd1498Szrj 
3*38fd1498Szrj // Copyright (C) 2012-2018 Free Software Foundation, Inc.
4*38fd1498Szrj //
5*38fd1498Szrj // This file is part of the GNU ISO C++ Library.  This library is free
6*38fd1498Szrj // software; you can redistribute it and/or modify it under the
7*38fd1498Szrj // terms of the GNU General Public License as published by the
8*38fd1498Szrj // Free Software Foundation; either version 3, or (at your option)
9*38fd1498Szrj // any later version.
10*38fd1498Szrj 
11*38fd1498Szrj // This library is distributed in the hope that it will be useful,
12*38fd1498Szrj // but WITHOUT ANY WARRANTY; without even the implied warranty of
13*38fd1498Szrj // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
14*38fd1498Szrj // GNU General Public License for more details.
15*38fd1498Szrj 
16*38fd1498Szrj // Under Section 7 of GPL version 3, you are granted additional
17*38fd1498Szrj // permissions described in the GCC Runtime Library Exception, version
18*38fd1498Szrj // 3.1, as published by the Free Software Foundation.
19*38fd1498Szrj 
20*38fd1498Szrj // You should have received a copy of the GNU General Public License and
21*38fd1498Szrj // a copy of the GCC Runtime Library Exception along with this program;
22*38fd1498Szrj // see the files COPYING3 and COPYING.RUNTIME respectively.  If not, see
23*38fd1498Szrj // <http://www.gnu.org/licenses/>.
24*38fd1498Szrj 
25*38fd1498Szrj #define _GLIBCXX_USE_CXX11_ABI 1
26*38fd1498Szrj #include <random>
27*38fd1498Szrj 
28*38fd1498Szrj #ifdef  _GLIBCXX_USE_C99_STDINT_TR1
29*38fd1498Szrj 
30*38fd1498Szrj #if defined __i386__ || defined __x86_64__
31*38fd1498Szrj # include <cpuid.h>
32*38fd1498Szrj #endif
33*38fd1498Szrj 
34*38fd1498Szrj #include <cerrno>
35*38fd1498Szrj #include <cstdio>
36*38fd1498Szrj 
37*38fd1498Szrj #ifdef _GLIBCXX_HAVE_UNISTD_H
38*38fd1498Szrj # include <unistd.h>
39*38fd1498Szrj #endif
40*38fd1498Szrj 
41*38fd1498Szrj #ifdef _GLIBCXX_HAVE_SYS_IOCTL_H
42*38fd1498Szrj # include <sys/ioctl.h>
43*38fd1498Szrj #endif
44*38fd1498Szrj 
45*38fd1498Szrj #ifdef _GLIBCXX_HAVE_LINUX_TYPES_H
46*38fd1498Szrj # include <linux/types.h>
47*38fd1498Szrj #endif
48*38fd1498Szrj 
49*38fd1498Szrj #ifdef _GLIBCXX_HAVE_LINUX_RANDOM_H
50*38fd1498Szrj # include <linux/random.h>
51*38fd1498Szrj #endif
52*38fd1498Szrj 
53*38fd1498Szrj namespace std _GLIBCXX_VISIBILITY(default)
54*38fd1498Szrj {
55*38fd1498Szrj   namespace
56*38fd1498Szrj   {
57*38fd1498Szrj     static unsigned long
_M_strtoul(const std::string & __str)58*38fd1498Szrj     _M_strtoul(const std::string& __str)
59*38fd1498Szrj     {
60*38fd1498Szrj       unsigned long __ret = 5489UL;
61*38fd1498Szrj       if (__str != "mt19937")
62*38fd1498Szrj 	{
63*38fd1498Szrj 	  const char* __nptr = __str.c_str();
64*38fd1498Szrj 	  char* __endptr;
65*38fd1498Szrj 	  __ret = std::strtoul(__nptr, &__endptr, 0);
66*38fd1498Szrj 	  if (*__nptr == '\0' || *__endptr != '\0')
67*38fd1498Szrj 	    std::__throw_runtime_error(__N("random_device::_M_strtoul"
68*38fd1498Szrj 					   "(const std::string&)"));
69*38fd1498Szrj 	}
70*38fd1498Szrj       return __ret;
71*38fd1498Szrj     }
72*38fd1498Szrj 
73*38fd1498Szrj #if (defined __i386__ || defined __x86_64__) && defined _GLIBCXX_X86_RDRAND
74*38fd1498Szrj     unsigned int
75*38fd1498Szrj     __attribute__ ((target("rdrnd")))
__x86_rdrand(void)76*38fd1498Szrj     __x86_rdrand(void)
77*38fd1498Szrj     {
78*38fd1498Szrj       unsigned int retries = 100;
79*38fd1498Szrj       unsigned int val;
80*38fd1498Szrj 
81*38fd1498Szrj       while (__builtin_ia32_rdrand32_step(&val) == 0)
82*38fd1498Szrj 	if (--retries == 0)
83*38fd1498Szrj 	  std::__throw_runtime_error(__N("random_device::__x86_rdrand(void)"));
84*38fd1498Szrj 
85*38fd1498Szrj       return val;
86*38fd1498Szrj     }
87*38fd1498Szrj #endif
88*38fd1498Szrj   }
89*38fd1498Szrj 
90*38fd1498Szrj   void
_M_init(const std::string & token)91*38fd1498Szrj   random_device::_M_init(const std::string& token)
92*38fd1498Szrj   {
93*38fd1498Szrj     const char *fname = token.c_str();
94*38fd1498Szrj 
95*38fd1498Szrj     if (token == "default")
96*38fd1498Szrj       {
97*38fd1498Szrj #if (defined __i386__ || defined __x86_64__) && defined _GLIBCXX_X86_RDRAND
98*38fd1498Szrj 	unsigned int eax, ebx, ecx, edx;
99*38fd1498Szrj 	// Check availability of cpuid and, for now at least, also the
100*38fd1498Szrj 	// CPU signature for Intel's
101*38fd1498Szrj 	if (__get_cpuid_max(0, &ebx) > 0 && ebx == signature_INTEL_ebx)
102*38fd1498Szrj 	  {
103*38fd1498Szrj 	    __cpuid(1, eax, ebx, ecx, edx);
104*38fd1498Szrj 	    if (ecx & bit_RDRND)
105*38fd1498Szrj 	      {
106*38fd1498Szrj 		_M_file = nullptr;
107*38fd1498Szrj 		return;
108*38fd1498Szrj 	      }
109*38fd1498Szrj 	  }
110*38fd1498Szrj #endif
111*38fd1498Szrj 
112*38fd1498Szrj 	fname = "/dev/urandom";
113*38fd1498Szrj       }
114*38fd1498Szrj     else if (token != "/dev/urandom" && token != "/dev/random")
115*38fd1498Szrj     fail:
116*38fd1498Szrj       std::__throw_runtime_error(__N("random_device::"
117*38fd1498Szrj 				     "random_device(const std::string&)"));
118*38fd1498Szrj 
119*38fd1498Szrj     _M_file = static_cast<void*>(std::fopen(fname, "rb"));
120*38fd1498Szrj     if (!_M_file)
121*38fd1498Szrj       goto fail;
122*38fd1498Szrj   }
123*38fd1498Szrj 
124*38fd1498Szrj   void
_M_init_pretr1(const std::string & token)125*38fd1498Szrj   random_device::_M_init_pretr1(const std::string& token)
126*38fd1498Szrj   {
127*38fd1498Szrj     _M_mt.seed(_M_strtoul(token));
128*38fd1498Szrj   }
129*38fd1498Szrj 
130*38fd1498Szrj   void
_M_fini()131*38fd1498Szrj   random_device::_M_fini()
132*38fd1498Szrj   {
133*38fd1498Szrj     if (_M_file)
134*38fd1498Szrj       std::fclose(static_cast<FILE*>(_M_file));
135*38fd1498Szrj   }
136*38fd1498Szrj 
137*38fd1498Szrj   random_device::result_type
_M_getval()138*38fd1498Szrj   random_device::_M_getval()
139*38fd1498Szrj   {
140*38fd1498Szrj #if (defined __i386__ || defined __x86_64__) && defined _GLIBCXX_X86_RDRAND
141*38fd1498Szrj     if (!_M_file)
142*38fd1498Szrj       return __x86_rdrand();
143*38fd1498Szrj #endif
144*38fd1498Szrj 
145*38fd1498Szrj     result_type __ret;
146*38fd1498Szrj     void* p = &__ret;
147*38fd1498Szrj     size_t n = sizeof(result_type);
148*38fd1498Szrj #ifdef _GLIBCXX_HAVE_UNISTD_H
149*38fd1498Szrj     do
150*38fd1498Szrj       {
151*38fd1498Szrj 	const int e = read(fileno(static_cast<FILE*>(_M_file)), p, n);
152*38fd1498Szrj 	if (e > 0)
153*38fd1498Szrj 	  {
154*38fd1498Szrj 	    n -= e;
155*38fd1498Szrj 	    p = static_cast<char*>(p) + e;
156*38fd1498Szrj 	  }
157*38fd1498Szrj 	else if (e != -1 || errno != EINTR)
158*38fd1498Szrj 	  __throw_runtime_error(__N("random_device could not be read"));
159*38fd1498Szrj       }
160*38fd1498Szrj     while (n > 0);
161*38fd1498Szrj #else
162*38fd1498Szrj     const size_t e = std::fread(p, n, 1, static_cast<FILE*>(_M_file));
163*38fd1498Szrj     if (e != 1)
164*38fd1498Szrj       __throw_runtime_error(__N("random_device could not be read"));
165*38fd1498Szrj #endif
166*38fd1498Szrj 
167*38fd1498Szrj     return __ret;
168*38fd1498Szrj   }
169*38fd1498Szrj 
170*38fd1498Szrj   random_device::result_type
_M_getval_pretr1()171*38fd1498Szrj   random_device::_M_getval_pretr1()
172*38fd1498Szrj   {
173*38fd1498Szrj     return _M_mt();
174*38fd1498Szrj   }
175*38fd1498Szrj 
176*38fd1498Szrj   double
_M_getentropy() const177*38fd1498Szrj   random_device::_M_getentropy() const noexcept
178*38fd1498Szrj   {
179*38fd1498Szrj #if defined _GLIBCXX_HAVE_SYS_IOCTL_H && defined RNDGETENTCNT
180*38fd1498Szrj     if (!_M_file)
181*38fd1498Szrj       return 0.0;
182*38fd1498Szrj 
183*38fd1498Szrj     const int fd = fileno(static_cast<FILE*>(_M_file));
184*38fd1498Szrj     if (fd < 0)
185*38fd1498Szrj       return 0.0;
186*38fd1498Szrj 
187*38fd1498Szrj     int ent;
188*38fd1498Szrj     if (ioctl(fd, RNDGETENTCNT, &ent) < 0)
189*38fd1498Szrj       return 0.0;
190*38fd1498Szrj 
191*38fd1498Szrj     if (ent < 0)
192*38fd1498Szrj       return 0.0;
193*38fd1498Szrj 
194*38fd1498Szrj     const int max = sizeof(result_type) * __CHAR_BIT__;
195*38fd1498Szrj     if (ent > max)
196*38fd1498Szrj       ent = max;
197*38fd1498Szrj 
198*38fd1498Szrj     return static_cast<double>(ent);
199*38fd1498Szrj #else
200*38fd1498Szrj     return 0.0;
201*38fd1498Szrj #endif
202*38fd1498Szrj   }
203*38fd1498Szrj 
204*38fd1498Szrj   template class mersenne_twister_engine<
205*38fd1498Szrj     uint_fast32_t,
206*38fd1498Szrj     32, 624, 397, 31,
207*38fd1498Szrj     0x9908b0dfUL, 11,
208*38fd1498Szrj     0xffffffffUL, 7,
209*38fd1498Szrj     0x9d2c5680UL, 15,
210*38fd1498Szrj     0xefc60000UL, 18, 1812433253UL>;
211*38fd1498Szrj }
212*38fd1498Szrj #endif
213