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 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174
|
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
<html xmlns="http://www.w3.org/1999/xhtml">
<head>
<meta http-equiv="Content-Type" content="text/xhtml;charset=UTF-8"/>
<title>Random123-1.09: Random123/conventional/gsl_cbrng.h Source File</title>
<link href="tabs.css" rel="stylesheet" type="text/css"/>
<link href="search/search.css" rel="stylesheet" type="text/css"/>
<script type="text/javaScript" src="search/search.js"></script>
<link href="doxygen.css" rel="stylesheet" type="text/css"/>
</head>
<body onload='searchBox.OnSelectItem(0);'>
<div class="tabs"><ul class="tablist"><li style="padding-left: 1.5em; font-weight: bold">Random123-1.09 Documentation</li></ul></div>
<!-- Generated by Doxygen 1.7.1 -->
<script type="text/javascript"><!--
var searchBox = new SearchBox("searchBox", "search",false,'Search');
--></script>
<div class="navigation" id="top">
<div class="tabs">
<ul class="tablist">
<li><a href="index.html"><span>Main Page</span></a></li>
<li><a href="pages.html"><span>Related Pages</span></a></li>
<li><a href="modules.html"><span>Modules</span></a></li>
<li><a href="namespaces.html"><span>Namespaces</span></a></li>
<li><a href="annotated.html"><span>Classes</span></a></li>
<li class="current"><a href="files.html"><span>Files</span></a></li>
<li id="searchli">
<div id="MSearchBox" class="MSearchBoxInactive">
<span class="left">
<img id="MSearchSelect" src="search/mag_sel.png"
onmouseover="return searchBox.OnSearchSelectShow()"
onmouseout="return searchBox.OnSearchSelectHide()"
alt=""/>
<input type="text" id="MSearchField" value="Search" accesskey="S"
onfocus="searchBox.OnSearchFieldFocus(true)"
onblur="searchBox.OnSearchFieldFocus(false)"
onkeyup="searchBox.OnSearchFieldChange(event)"/>
</span><span class="right">
<a id="MSearchClose" href="javascript:searchBox.CloseResultsWindow()"><img id="MSearchCloseImg" border="0" src="search/close.png" alt=""/></a>
</span>
</div>
</li>
</ul>
</div>
<div class="tabs2">
<ul class="tablist">
<li><a href="files.html"><span>File List</span></a></li>
<li><a href="globals.html"><span>File Members</span></a></li>
</ul>
</div>
<div class="header">
<div class="headertitle">
<h1>Random123/conventional/gsl_cbrng.h</h1> </div>
</div>
<div class="contents">
<a href="gsl__cbrng_8h.html">Go to the documentation of this file.</a><div class="fragment"><pre class="fragment"><a name="l00001"></a>00001 <span class="comment">/*</span>
<a name="l00002"></a>00002 <span class="comment">Copyright 2010-2011, D. E. Shaw Research.</span>
<a name="l00003"></a>00003 <span class="comment">All rights reserved.</span>
<a name="l00004"></a>00004 <span class="comment"></span>
<a name="l00005"></a>00005 <span class="comment">Redistribution and use in source and binary forms, with or without</span>
<a name="l00006"></a>00006 <span class="comment">modification, are permitted provided that the following conditions are</span>
<a name="l00007"></a>00007 <span class="comment">met:</span>
<a name="l00008"></a>00008 <span class="comment"></span>
<a name="l00009"></a>00009 <span class="comment">* Redistributions of source code must retain the above copyright</span>
<a name="l00010"></a>00010 <span class="comment"> notice, this list of conditions, and the following disclaimer.</span>
<a name="l00011"></a>00011 <span class="comment"></span>
<a name="l00012"></a>00012 <span class="comment">* Redistributions in binary form must reproduce the above copyright</span>
<a name="l00013"></a>00013 <span class="comment"> notice, this list of conditions, and the following disclaimer in the</span>
<a name="l00014"></a>00014 <span class="comment"> documentation and/or other materials provided with the distribution.</span>
<a name="l00015"></a>00015 <span class="comment"></span>
<a name="l00016"></a>00016 <span class="comment">* Neither the name of D. E. Shaw Research nor the names of its</span>
<a name="l00017"></a>00017 <span class="comment"> contributors may be used to endorse or promote products derived from</span>
<a name="l00018"></a>00018 <span class="comment"> this software without specific prior written permission.</span>
<a name="l00019"></a>00019 <span class="comment"></span>
<a name="l00020"></a>00020 <span class="comment">THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS</span>
<a name="l00021"></a>00021 <span class="comment">"AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT</span>
<a name="l00022"></a>00022 <span class="comment">LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR</span>
<a name="l00023"></a>00023 <span class="comment">A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT</span>
<a name="l00024"></a>00024 <span class="comment">OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,</span>
<a name="l00025"></a>00025 <span class="comment">SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT</span>
<a name="l00026"></a>00026 <span class="comment">LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,</span>
<a name="l00027"></a>00027 <span class="comment">DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY</span>
<a name="l00028"></a>00028 <span class="comment">THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT</span>
<a name="l00029"></a>00029 <span class="comment">(INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE</span>
<a name="l00030"></a>00030 <span class="comment">OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.</span>
<a name="l00031"></a>00031 <span class="comment">*/</span>
<a name="l00032"></a>00032 <span class="preprocessor">#ifndef __r123_compat_gslrng_dot_h__</span>
<a name="l00033"></a>00033 <span class="preprocessor"></span><span class="preprocessor">#define __r123_compat_gslrng_dot_h__</span>
<a name="l00034"></a>00034 <span class="preprocessor"></span>
<a name="l00035"></a>00035 <span class="preprocessor">#include <gsl/gsl_rng.h></span>
<a name="l00036"></a>00036 <span class="preprocessor">#include <string.h></span>
<a name="l00037"></a>00037
<a name="l00067"></a><a class="code" href="gsl__cbrng_8h.html#af561a004eef8e93cdfd6b255a8a1eb75">00067</a> <span class="preprocessor">#define GSL_CBRNG(NAME, CBRNGNAME) \</span>
<a name="l00068"></a>00068 <span class="preprocessor">const gsl_rng_type *gsl_rng_##NAME; \</span>
<a name="l00069"></a>00069 <span class="preprocessor"> \</span>
<a name="l00070"></a>00070 <span class="preprocessor">typedef struct{ \</span>
<a name="l00071"></a>00071 <span class="preprocessor"> CBRNGNAME##_ctr_t ctr; \</span>
<a name="l00072"></a>00072 <span class="preprocessor"> CBRNGNAME##_ctr_t r; \</span>
<a name="l00073"></a>00073 <span class="preprocessor"> CBRNGNAME##_key_t key; \</span>
<a name="l00074"></a>00074 <span class="preprocessor"> int elem; \</span>
<a name="l00075"></a>00075 <span class="preprocessor">} NAME##_state; \</span>
<a name="l00076"></a>00076 <span class="preprocessor"> \</span>
<a name="l00077"></a>00077 <span class="preprocessor">static unsigned long int NAME##_get(void *vstate){ \</span>
<a name="l00078"></a>00078 <span class="preprocessor"> NAME##_state *st = (NAME##_state *)vstate; \</span>
<a name="l00079"></a>00079 <span class="preprocessor"> const int N=sizeof(st->ctr.v)/sizeof(st->ctr.v[0]); \</span>
<a name="l00080"></a>00080 <span class="preprocessor"> if( st->elem == 0 ){ \</span>
<a name="l00081"></a>00081 <span class="preprocessor"> ++st->ctr.v[0]; \</span>
<a name="l00082"></a>00082 <span class="preprocessor"> if( N>1 && st->ctr.v[0] == 0 ) ++st->ctr.v[1]; \</span>
<a name="l00083"></a>00083 <span class="preprocessor"> if( N>2 && st->ctr.v[1] == 0 ) ++st->ctr.v[2]; \</span>
<a name="l00084"></a>00084 <span class="preprocessor"> if( N>3 && st->ctr.v[2] == 0 ) ++st->ctr.v[3]; \</span>
<a name="l00085"></a>00085 <span class="preprocessor"> st->r = CBRNGNAME(st->ctr, st->key); \</span>
<a name="l00086"></a>00086 <span class="preprocessor"> st->elem = N; \</span>
<a name="l00087"></a>00087 <span class="preprocessor"> } \</span>
<a name="l00088"></a>00088 <span class="preprocessor"> return 0xffffffffUL & st->r.v[--st->elem]; \</span>
<a name="l00089"></a>00089 <span class="preprocessor">} \</span>
<a name="l00090"></a>00090 <span class="preprocessor"> \</span>
<a name="l00091"></a>00091 <span class="preprocessor">static double \</span>
<a name="l00092"></a>00092 <span class="preprocessor">NAME##_get_double (void * vstate) \</span>
<a name="l00093"></a>00093 <span class="preprocessor">{ \</span>
<a name="l00094"></a>00094 <span class="preprocessor"> return NAME##_get (vstate)/4294967296.0; \</span>
<a name="l00095"></a>00095 <span class="preprocessor">} \</span>
<a name="l00096"></a>00096 <span class="preprocessor"> \</span>
<a name="l00097"></a>00097 <span class="preprocessor">static void NAME##_set(void *vstate, unsigned long int s){ \</span>
<a name="l00098"></a>00098 <span class="preprocessor"> NAME##_state *st = (NAME##_state *)vstate; \</span>
<a name="l00099"></a>00099 <span class="preprocessor"> st->elem = 0; \</span>
<a name="l00100"></a>00100 <span class="preprocessor"> </span><span class="comment">/* Assume that key and ctr have an array member, v, \</span>
<a name="l00101"></a>00101 <span class="comment"> as if they are r123arrayNxW. If not, this will fail \</span>
<a name="l00102"></a>00102 <span class="comment"> to compile. In particular, this macro fails to compile \</span>
<a name="l00103"></a>00103 <span class="comment"> when the underlying CBRNG requires use of keyinit */</span> \
<a name="l00104"></a>00104 memset(&st->ctr.v[0], 0, sizeof(st->ctr.v)); \
<a name="l00105"></a>00105 memset(&st->key.v[0], 0, sizeof(st->key.v)); \
<a name="l00106"></a>00106 <span class="comment">/* GSL 1.15 documentation says this about gsl_rng_set: \</span>
<a name="l00107"></a>00107 <span class="comment"> Note that the most generators only accept 32-bit seeds, with higher \</span>
<a name="l00108"></a>00108 <span class="comment"> values being reduced modulo 2^32. For generators with smaller \</span>
<a name="l00109"></a>00109 <span class="comment"> ranges the maximum seed value will typically be lower. \</span>
<a name="l00110"></a>00110 <span class="comment"> so we won't jump through any hoops here to deal with \</span>
<a name="l00111"></a>00111 <span class="comment"> high bits if sizeof(unsigned long) > sizeof(uint32_t). */</span> \
<a name="l00112"></a>00112 st->key.v[0] = s; \
<a name="l00113"></a>00113 } \
<a name="l00114"></a>00114 \
<a name="l00115"></a>00115 static const gsl_rng_type NAME##_type = { \
<a name="l00116"></a>00116 #NAME, \
<a name="l00117"></a>00117 0xffffffffUL, \
<a name="l00118"></a>00118 0, \
<a name="l00119"></a>00119 sizeof(NAME##_state), \
<a name="l00120"></a>00120 &NAME##_set, \
<a name="l00121"></a>00121 &NAME##_get, \
<a name="l00122"></a>00122 &NAME##_get_double \
<a name="l00123"></a>00123 }; \
<a name="l00124"></a>00124 \
<a name="l00125"></a>00125 const gsl_rng_type *gsl_rng_##NAME = &NAME##_type
<a name="l00126"></a>00126
<a name="l00127"></a>00127 <span class="preprocessor">#endif</span>
<a name="l00128"></a>00128 <span class="preprocessor"></span>
</pre></div></div>
</div>
<!--- window showing the filter options -->
<div id="MSearchSelectWindow"
onmouseover="return searchBox.OnSearchSelectShow()"
onmouseout="return searchBox.OnSearchSelectHide()"
onkeydown="return searchBox.OnSearchSelectKey(event)">
<a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(0)"><span class="SelectionMark"> </span>All</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(1)"><span class="SelectionMark"> </span>Classes</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(2)"><span class="SelectionMark"> </span>Namespaces</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(3)"><span class="SelectionMark"> </span>Files</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(4)"><span class="SelectionMark"> </span>Functions</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(5)"><span class="SelectionMark"> </span>Variables</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(6)"><span class="SelectionMark"> </span>Typedefs</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(7)"><span class="SelectionMark"> </span>Enumerations</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(8)"><span class="SelectionMark"> </span>Enumerator</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(9)"><span class="SelectionMark"> </span>Friends</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(10)"><span class="SelectionMark"> </span>Defines</a></div>
<!-- iframe showing the search results (closed by default) -->
<div id="MSearchResultsWindow">
<iframe src="" frameborder="0"
name="MSearchResults" id="MSearchResults">
</iframe>
</div>
<hr class="footer"/><address class="footer"><small>Generated on Mon Mar 7 2016 18:34:00 for Random123-1.09 by
<a href="http://www.doxygen.org/index.html">
<img class="footer" src="doxygen.png" alt="doxygen"/></a> 1.7.1 </small></address>
</body>
</html>
|