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
|
<!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/gsl_microrng.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/gsl_microrng.h</h1> </div>
</div>
<div class="contents">
<a href="gsl__microrng_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_gslmicrorng_dot_h__</span>
<a name="l00033"></a>00033 <span class="preprocessor"></span><span class="preprocessor">#define __r123_gslmicrorng_dot_h__</span>
<a name="l00034"></a>00034 <span class="preprocessor"></span>
<a name="l00035"></a>00035
<a name="l00036"></a>00036 <span class="preprocessor">#include <gsl/gsl_rng.h></span>
<a name="l00037"></a>00037 <span class="preprocessor">#include <string.h></span>
<a name="l00038"></a>00038
<a name="l00079"></a><a class="code" href="gsl__microrng_8h.html#a21c7bb64a536a1704c6dc96856b78297">00079</a> <span class="preprocessor">#define GSL_MICRORNG(NAME, CBRNGNAME) \</span>
<a name="l00080"></a>00080 <span class="preprocessor">const gsl_rng_type *gsl_rng_##NAME; \</span>
<a name="l00081"></a>00081 <span class="preprocessor"> \</span>
<a name="l00082"></a>00082 <span class="preprocessor">typedef struct{ \</span>
<a name="l00083"></a>00083 <span class="preprocessor"> CBRNGNAME##_ctr_t ctr; \</span>
<a name="l00084"></a>00084 <span class="preprocessor"> CBRNGNAME##_ctr_t r; \</span>
<a name="l00085"></a>00085 <span class="preprocessor"> CBRNGNAME##_key_t key; \</span>
<a name="l00086"></a>00086 <span class="preprocessor"> R123_ULONG_LONG n; \</span>
<a name="l00087"></a>00087 <span class="preprocessor"> int elem; \</span>
<a name="l00088"></a>00088 <span class="preprocessor">} NAME##_state; \</span>
<a name="l00089"></a>00089 <span class="preprocessor"> \</span>
<a name="l00090"></a>00090 <span class="preprocessor">static unsigned long int NAME##_get(void *vstate){ \</span>
<a name="l00091"></a>00091 <span class="preprocessor"> NAME##_state *st = (NAME##_state *)vstate; \</span>
<a name="l00092"></a>00092 <span class="preprocessor"> const int N=sizeof(st->ctr.v)/sizeof(st->ctr.v[0]); \</span>
<a name="l00093"></a>00093 <span class="preprocessor"> if( st->elem == 0 ){ \</span>
<a name="l00094"></a>00094 <span class="preprocessor"> CBRNGNAME##_ctr_t c = st->ctr; \</span>
<a name="l00095"></a>00095 <span class="preprocessor"> c.v[N-1] |= st->n<<(R123_W(CBRNGNAME##_ctr_t)-32); \</span>
<a name="l00096"></a>00096 <span class="preprocessor"> st->n++; \</span>
<a name="l00097"></a>00097 <span class="preprocessor"> st->r = CBRNGNAME(c, st->key); \</span>
<a name="l00098"></a>00098 <span class="preprocessor"> st->elem = N; \</span>
<a name="l00099"></a>00099 <span class="preprocessor"> } \</span>
<a name="l00100"></a>00100 <span class="preprocessor"> return 0xffffffff & st->r.v[--st->elem]; \</span>
<a name="l00101"></a>00101 <span class="preprocessor">} \</span>
<a name="l00102"></a>00102 <span class="preprocessor"> \</span>
<a name="l00103"></a>00103 <span class="preprocessor">static double \</span>
<a name="l00104"></a>00104 <span class="preprocessor">NAME##_get_double (void * vstate) \</span>
<a name="l00105"></a>00105 <span class="preprocessor">{ \</span>
<a name="l00106"></a>00106 <span class="preprocessor"> return NAME##_get (vstate)/4294967296.; \</span>
<a name="l00107"></a>00107 <span class="preprocessor">} \</span>
<a name="l00108"></a>00108 <span class="preprocessor"> \</span>
<a name="l00109"></a>00109 <span class="preprocessor">static void NAME##_set(void *vstate, unsigned long int s){ \</span>
<a name="l00110"></a>00110 <span class="preprocessor"> NAME##_state *st = (NAME##_state *)vstate; \</span>
<a name="l00111"></a>00111 <span class="preprocessor"> (void)s; </span><span class="comment">/* ignored */</span> \
<a name="l00112"></a>00112 st->elem = 0; \
<a name="l00113"></a>00113 st->n = ~0; <span class="comment">/* will abort if _reset is not called */</span> \
<a name="l00114"></a>00114 } \
<a name="l00115"></a>00115 \
<a name="l00116"></a>00116 static const gsl_rng_type NAME##_type = { \
<a name="l00117"></a>00117 #NAME, \
<a name="l00118"></a>00118 0xffffffffUL, \
<a name="l00119"></a>00119 0, \
<a name="l00120"></a>00120 sizeof(NAME##_state), \
<a name="l00121"></a>00121 &NAME##_set, \
<a name="l00122"></a>00122 &NAME##_get, \
<a name="l00123"></a>00123 &NAME##_get_double \
<a name="l00124"></a>00124 }; \
<a name="l00125"></a>00125 \
<a name="l00126"></a>00126 R123_STATIC_INLINE void NAME##_reset(const gsl_rng* gr, CBRNGNAME##_ctr_t c, CBRNGNAME##_key_t k) { \
<a name="l00127"></a>00127 NAME##_state* state = (NAME##_state *)gr->state; \
<a name="l00128"></a>00128 state->ctr = c; \
<a name="l00129"></a>00129 state->key = k; \
<a name="l00130"></a>00130 state->n = 0; \
<a name="l00131"></a>00131 state->elem = 0; \
<a name="l00132"></a>00132 } \
<a name="l00133"></a>00133 \
<a name="l00134"></a>00134 const gsl_rng_type *gsl_rng_##NAME = &NAME##_type
<a name="l00135"></a>00135
<a name="l00136"></a>00136 <span class="preprocessor">#endif</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>
|