Modified: websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/internal/source32/AbstractWell.html ============================================================================== --- websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/internal/source32/AbstractWell.html (original) +++ websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/internal/source32/AbstractWell.html Fri Sep 2 16:37:56 2016 @@ -89,127 +89,124 @@ <a class="jxr_linenumber" name="L81" href="#L81">81</a> <em class="jxr_javadoccomment">/**</em> <a class="jxr_linenumber" name="L82" href="#L82">82</a> <em class="jxr_javadoccomment"> * Initializes the generator with the given {@code seed}.</em> <a class="jxr_linenumber" name="L83" href="#L83">83</a> <em class="jxr_javadoccomment"> *</em> -<a class="jxr_linenumber" name="L84" href="#L84">84</a> <em class="jxr_javadoccomment"> * <p>The state of the generator is exactly the same as a new generator built</em> -<a class="jxr_linenumber" name="L85" href="#L85">85</a> <em class="jxr_javadoccomment"> * with the same seed.</p></em> -<a class="jxr_linenumber" name="L86" href="#L86">86</a> <em class="jxr_javadoccomment"> *</em> -<a class="jxr_linenumber" name="L87" href="#L87">87</a> <em class="jxr_javadoccomment"> * @param seed Seed. Cannot be null.</em> -<a class="jxr_linenumber" name="L88" href="#L88">88</a> <em class="jxr_javadoccomment"> */</em> -<a class="jxr_linenumber" name="L89" href="#L89">89</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">void</strong> setSeedInternal(<strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong>[] seed) { -<a class="jxr_linenumber" name="L90" href="#L90">90</a> System.arraycopy(seed, 0, v, 0, Math.min(seed.length, v.length)); -<a class="jxr_linenumber" name="L91" href="#L91">91</a> -<a class="jxr_linenumber" name="L92" href="#L92">92</a> <strong class="jxr_keyword">if</strong> (seed.length < v.length) { -<a class="jxr_linenumber" name="L93" href="#L93">93</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> i = seed.length; i < v.length; ++i) { -<a class="jxr_linenumber" name="L94" href="#L94">94</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">long</strong> current = v[i - seed.length]; -<a class="jxr_linenumber" name="L95" href="#L95">95</a> v[i] = (<strong class="jxr_keyword">int</strong>) ((1812433253L * (current ^ (current >> 30)) + i) & 0xffffffffL); -<a class="jxr_linenumber" name="L96" href="#L96">96</a> } -<a class="jxr_linenumber" name="L97" href="#L97">97</a> } +<a class="jxr_linenumber" name="L84" href="#L84">84</a> <em class="jxr_javadoccomment"> * @param seed Seed. Cannot be null.</em> +<a class="jxr_linenumber" name="L85" href="#L85">85</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L86" href="#L86">86</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">void</strong> setSeedInternal(<strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong>[] seed) { +<a class="jxr_linenumber" name="L87" href="#L87">87</a> System.arraycopy(seed, 0, v, 0, Math.min(seed.length, v.length)); +<a class="jxr_linenumber" name="L88" href="#L88">88</a> +<a class="jxr_linenumber" name="L89" href="#L89">89</a> <strong class="jxr_keyword">if</strong> (seed.length < v.length) { +<a class="jxr_linenumber" name="L90" href="#L90">90</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> i = seed.length; i < v.length; ++i) { +<a class="jxr_linenumber" name="L91" href="#L91">91</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">long</strong> current = v[i - seed.length]; +<a class="jxr_linenumber" name="L92" href="#L92">92</a> v[i] = (<strong class="jxr_keyword">int</strong>) ((1812433253L * (current ^ (current >> 30)) + i) & 0xffffffffL); +<a class="jxr_linenumber" name="L93" href="#L93">93</a> } +<a class="jxr_linenumber" name="L94" href="#L94">94</a> } +<a class="jxr_linenumber" name="L95" href="#L95">95</a> +<a class="jxr_linenumber" name="L96" href="#L96">96</a> index = 0; +<a class="jxr_linenumber" name="L97" href="#L97">97</a> } <a class="jxr_linenumber" name="L98" href="#L98">98</a> -<a class="jxr_linenumber" name="L99" href="#L99">99</a> index = 0; -<a class="jxr_linenumber" name="L100" href="#L100">100</a> } -<a class="jxr_linenumber" name="L101" href="#L101">101</a> -<a class="jxr_linenumber" name="L102" href="#L102">102</a> <em class="jxr_javadoccomment">/**</em> -<a class="jxr_linenumber" name="L103" href="#L103">103</a> <em class="jxr_javadoccomment"> * Calculate the number of 32-bits blocks.</em> -<a class="jxr_linenumber" name="L104" href="#L104">104</a> <em class="jxr_javadoccomment"> *</em> -<a class="jxr_linenumber" name="L105" href="#L105">105</a> <em class="jxr_javadoccomment"> * @param k Number of bits in the pool (not necessarily a multiple of 32).</em> -<a class="jxr_linenumber" name="L106" href="#L106">106</a> <em class="jxr_javadoccomment"> * @return the number of 32-bits blocks.</em> -<a class="jxr_linenumber" name="L107" href="#L107">107</a> <em class="jxr_javadoccomment"> */</em> -<a class="jxr_linenumber" name="L108" href="#L108">108</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">static</strong> <strong class="jxr_keyword">int</strong> calculateBlockCount(<strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> k) { -<a class="jxr_linenumber" name="L109" href="#L109">109</a> <em class="jxr_comment">// the bits pool contains k bits, k = r w - p where r is the number</em> -<a class="jxr_linenumber" name="L110" href="#L110">110</a> <em class="jxr_comment">// of w bits blocks, w is the block size (always 32 in the original paper)</em> -<a class="jxr_linenumber" name="L111" href="#L111">111</a> <em class="jxr_comment">// and p is the number of unused bits in the last block</em> -<a class="jxr_linenumber" name="L112" href="#L112">112</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> w = 32; -<a class="jxr_linenumber" name="L113" href="#L113">113</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> r = (k + w - 1) / w; -<a class="jxr_linenumber" name="L114" href="#L114">114</a> <strong class="jxr_keyword">return</strong> r; -<a class="jxr_linenumber" name="L115" href="#L115">115</a> } -<a class="jxr_linenumber" name="L116" href="#L116">116</a> -<a class="jxr_linenumber" name="L117" href="#L117">117</a> <em class="jxr_javadoccomment">/**</em> -<a class="jxr_linenumber" name="L118" href="#L118">118</a> <em class="jxr_javadoccomment"> * Inner class used to store the indirection index table which is fixed for a given</em> -<a class="jxr_linenumber" name="L119" href="#L119">119</a> <em class="jxr_javadoccomment"> * type of WELL class of pseudo-random number generator.</em> -<a class="jxr_linenumber" name="L120" href="#L120">120</a> <em class="jxr_javadoccomment"> */</em> -<a class="jxr_linenumber" name="L121" href="#L121">121</a> <strong class="jxr_keyword">protected</strong> <strong class="jxr_keyword">static</strong> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">class</strong> <a href="../../../../../../org/apache/commons/rng/internal/source32/AbstractWell.html">IndexTable</a> { -<a class="jxr_linenumber" name="L122" href="#L122">122</a> <em class="jxr_javadoccomment">/** Index indirection table giving for each index its predecessor taking table size into account. */</em> -<a class="jxr_linenumber" name="L123" href="#L123">123</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong>[] iRm1; -<a class="jxr_linenumber" name="L124" href="#L124">124</a> <em class="jxr_javadoccomment">/** Index indirection table giving for each index its second predecessor taking table size into account. */</em> -<a class="jxr_linenumber" name="L125" href="#L125">125</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong>[] iRm2; -<a class="jxr_linenumber" name="L126" href="#L126">126</a> <em class="jxr_javadoccomment">/** Index indirection table giving for each index the value index + m1 taking table size into account. */</em> -<a class="jxr_linenumber" name="L127" href="#L127">127</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong>[] i1; -<a class="jxr_linenumber" name="L128" href="#L128">128</a> <em class="jxr_javadoccomment">/** Index indirection table giving for each index the value index + m2 taking table size into account. */</em> -<a class="jxr_linenumber" name="L129" href="#L129">129</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong>[] i2; -<a class="jxr_linenumber" name="L130" href="#L130">130</a> <em class="jxr_javadoccomment">/** Index indirection table giving for each index the value index + m3 taking table size into account. */</em> -<a class="jxr_linenumber" name="L131" href="#L131">131</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong>[] i3; -<a class="jxr_linenumber" name="L132" href="#L132">132</a> -<a class="jxr_linenumber" name="L133" href="#L133">133</a> <em class="jxr_javadoccomment">/** Creates a new pre-calculated indirection index table.</em> -<a class="jxr_linenumber" name="L134" href="#L134">134</a> <em class="jxr_javadoccomment"> * @param k number of bits in the pool (not necessarily a multiple of 32)</em> -<a class="jxr_linenumber" name="L135" href="#L135">135</a> <em class="jxr_javadoccomment"> * @param m1 first parameter of the algorithm</em> -<a class="jxr_linenumber" name="L136" href="#L136">136</a> <em class="jxr_javadoccomment"> * @param m2 second parameter of the algorithm</em> -<a class="jxr_linenumber" name="L137" href="#L137">137</a> <em class="jxr_javadoccomment"> * @param m3 third parameter of the algorithm</em> -<a class="jxr_linenumber" name="L138" href="#L138">138</a> <em class="jxr_javadoccomment"> */</em> -<a class="jxr_linenumber" name="L139" href="#L139">139</a> <strong class="jxr_keyword">public</strong> <a href="../../../../../../org/apache/commons/rng/internal/source32/AbstractWell.html">IndexTable</a>(<strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> k, <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> m1, <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> m2, <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> m3) { -<a class="jxr_linenumber" name="L140" href="#L140">140</a> -<a class="jxr_linenumber" name="L141" href="#L141">141</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> r = calculateBlockCount(k); -<a class="jxr_linenumber" name="L142" href="#L142">142</a> -<a class="jxr_linenumber" name="L143" href="#L143">143</a> <em class="jxr_comment">// precompute indirection index tables. These tables are used for optimizing access</em> -<a class="jxr_linenumber" name="L144" href="#L144">144</a> <em class="jxr_comment">// they allow saving computations like "(j + r - 2) % r" with costly modulo operations</em> -<a class="jxr_linenumber" name="L145" href="#L145">145</a> iRm1 = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">int</strong>[r]; -<a class="jxr_linenumber" name="L146" href="#L146">146</a> iRm2 = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">int</strong>[r]; -<a class="jxr_linenumber" name="L147" href="#L147">147</a> i1 = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">int</strong>[r]; -<a class="jxr_linenumber" name="L148" href="#L148">148</a> i2 = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">int</strong>[r]; -<a class="jxr_linenumber" name="L149" href="#L149">149</a> i3 = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">int</strong>[r]; -<a class="jxr_linenumber" name="L150" href="#L150">150</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> j = 0; j < r; ++j) { -<a class="jxr_linenumber" name="L151" href="#L151">151</a> iRm1[j] = (j + r - 1) % r; -<a class="jxr_linenumber" name="L152" href="#L152">152</a> iRm2[j] = (j + r - 2) % r; -<a class="jxr_linenumber" name="L153" href="#L153">153</a> i1[j] = (j + m1) % r; -<a class="jxr_linenumber" name="L154" href="#L154">154</a> i2[j] = (j + m2) % r; -<a class="jxr_linenumber" name="L155" href="#L155">155</a> i3[j] = (j + m3) % r; -<a class="jxr_linenumber" name="L156" href="#L156">156</a> } -<a class="jxr_linenumber" name="L157" href="#L157">157</a> } -<a class="jxr_linenumber" name="L158" href="#L158">158</a> -<a class="jxr_linenumber" name="L159" href="#L159">159</a> <em class="jxr_javadoccomment">/**</em> -<a class="jxr_linenumber" name="L160" href="#L160">160</a> <em class="jxr_javadoccomment"> * Returns the predecessor of the given index modulo the table size.</em> -<a class="jxr_linenumber" name="L161" href="#L161">161</a> <em class="jxr_javadoccomment"> * @param index the index to look at</em> -<a class="jxr_linenumber" name="L162" href="#L162">162</a> <em class="jxr_javadoccomment"> * @return (index - 1) % table size</em> -<a class="jxr_linenumber" name="L163" href="#L163">163</a> <em class="jxr_javadoccomment"> */</em> -<a class="jxr_linenumber" name="L164" href="#L164">164</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">int</strong> getIndexPred(<strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> index) { -<a class="jxr_linenumber" name="L165" href="#L165">165</a> <strong class="jxr_keyword">return</strong> iRm1[index]; -<a class="jxr_linenumber" name="L166" href="#L166">166</a> } -<a class="jxr_linenumber" name="L167" href="#L167">167</a> -<a class="jxr_linenumber" name="L168" href="#L168">168</a> <em class="jxr_javadoccomment">/**</em> -<a class="jxr_linenumber" name="L169" href="#L169">169</a> <em class="jxr_javadoccomment"> * Returns the second predecessor of the given index modulo the table size.</em> -<a class="jxr_linenumber" name="L170" href="#L170">170</a> <em class="jxr_javadoccomment"> * @param index the index to look at</em> -<a class="jxr_linenumber" name="L171" href="#L171">171</a> <em class="jxr_javadoccomment"> * @return (index - 2) % table size</em> -<a class="jxr_linenumber" name="L172" href="#L172">172</a> <em class="jxr_javadoccomment"> */</em> -<a class="jxr_linenumber" name="L173" href="#L173">173</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">int</strong> getIndexPred2(<strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> index) { -<a class="jxr_linenumber" name="L174" href="#L174">174</a> <strong class="jxr_keyword">return</strong> iRm2[index]; -<a class="jxr_linenumber" name="L175" href="#L175">175</a> } -<a class="jxr_linenumber" name="L176" href="#L176">176</a> -<a class="jxr_linenumber" name="L177" href="#L177">177</a> <em class="jxr_javadoccomment">/**</em> -<a class="jxr_linenumber" name="L178" href="#L178">178</a> <em class="jxr_javadoccomment"> * Returns index + M1 modulo the table size.</em> -<a class="jxr_linenumber" name="L179" href="#L179">179</a> <em class="jxr_javadoccomment"> * @param index the index to look at</em> -<a class="jxr_linenumber" name="L180" href="#L180">180</a> <em class="jxr_javadoccomment"> * @return (index + M1) % table size</em> -<a class="jxr_linenumber" name="L181" href="#L181">181</a> <em class="jxr_javadoccomment"> */</em> -<a class="jxr_linenumber" name="L182" href="#L182">182</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">int</strong> getIndexM1(<strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> index) { -<a class="jxr_linenumber" name="L183" href="#L183">183</a> <strong class="jxr_keyword">return</strong> i1[index]; -<a class="jxr_linenumber" name="L184" href="#L184">184</a> } -<a class="jxr_linenumber" name="L185" href="#L185">185</a> -<a class="jxr_linenumber" name="L186" href="#L186">186</a> <em class="jxr_javadoccomment">/**</em> -<a class="jxr_linenumber" name="L187" href="#L187">187</a> <em class="jxr_javadoccomment"> * Returns index + M2 modulo the table size.</em> -<a class="jxr_linenumber" name="L188" href="#L188">188</a> <em class="jxr_javadoccomment"> * @param index the index to look at</em> -<a class="jxr_linenumber" name="L189" href="#L189">189</a> <em class="jxr_javadoccomment"> * @return (index + M2) % table size</em> -<a class="jxr_linenumber" name="L190" href="#L190">190</a> <em class="jxr_javadoccomment"> */</em> -<a class="jxr_linenumber" name="L191" href="#L191">191</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">int</strong> getIndexM2(<strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> index) { -<a class="jxr_linenumber" name="L192" href="#L192">192</a> <strong class="jxr_keyword">return</strong> i2[index]; -<a class="jxr_linenumber" name="L193" href="#L193">193</a> } -<a class="jxr_linenumber" name="L194" href="#L194">194</a> -<a class="jxr_linenumber" name="L195" href="#L195">195</a> <em class="jxr_javadoccomment">/**</em> -<a class="jxr_linenumber" name="L196" href="#L196">196</a> <em class="jxr_javadoccomment"> * Returns index + M3 modulo the table size.</em> -<a class="jxr_linenumber" name="L197" href="#L197">197</a> <em class="jxr_javadoccomment"> * @param index the index to look at</em> -<a class="jxr_linenumber" name="L198" href="#L198">198</a> <em class="jxr_javadoccomment"> * @return (index + M3) % table size</em> -<a class="jxr_linenumber" name="L199" href="#L199">199</a> <em class="jxr_javadoccomment"> */</em> -<a class="jxr_linenumber" name="L200" href="#L200">200</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">int</strong> getIndexM3(<strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> index) { -<a class="jxr_linenumber" name="L201" href="#L201">201</a> <strong class="jxr_keyword">return</strong> i3[index]; -<a class="jxr_linenumber" name="L202" href="#L202">202</a> } -<a class="jxr_linenumber" name="L203" href="#L203">203</a> } -<a class="jxr_linenumber" name="L204" href="#L204">204</a> } +<a class="jxr_linenumber" name="L99" href="#L99">99</a> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L100" href="#L100">100</a> <em class="jxr_javadoccomment"> * Calculate the number of 32-bits blocks.</em> +<a class="jxr_linenumber" name="L101" href="#L101">101</a> <em class="jxr_javadoccomment"> *</em> +<a class="jxr_linenumber" name="L102" href="#L102">102</a> <em class="jxr_javadoccomment"> * @param k Number of bits in the pool (not necessarily a multiple of 32).</em> +<a class="jxr_linenumber" name="L103" href="#L103">103</a> <em class="jxr_javadoccomment"> * @return the number of 32-bits blocks.</em> +<a class="jxr_linenumber" name="L104" href="#L104">104</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L105" href="#L105">105</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">static</strong> <strong class="jxr_keyword">int</strong> calculateBlockCount(<strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> k) { +<a class="jxr_linenumber" name="L106" href="#L106">106</a> <em class="jxr_comment">// the bits pool contains k bits, k = r w - p where r is the number</em> +<a class="jxr_linenumber" name="L107" href="#L107">107</a> <em class="jxr_comment">// of w bits blocks, w is the block size (always 32 in the original paper)</em> +<a class="jxr_linenumber" name="L108" href="#L108">108</a> <em class="jxr_comment">// and p is the number of unused bits in the last block</em> +<a class="jxr_linenumber" name="L109" href="#L109">109</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> w = 32; +<a class="jxr_linenumber" name="L110" href="#L110">110</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> r = (k + w - 1) / w; +<a class="jxr_linenumber" name="L111" href="#L111">111</a> <strong class="jxr_keyword">return</strong> r; +<a class="jxr_linenumber" name="L112" href="#L112">112</a> } +<a class="jxr_linenumber" name="L113" href="#L113">113</a> +<a class="jxr_linenumber" name="L114" href="#L114">114</a> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L115" href="#L115">115</a> <em class="jxr_javadoccomment"> * Inner class used to store the indirection index table which is fixed for a given</em> +<a class="jxr_linenumber" name="L116" href="#L116">116</a> <em class="jxr_javadoccomment"> * type of WELL class of pseudo-random number generator.</em> +<a class="jxr_linenumber" name="L117" href="#L117">117</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L118" href="#L118">118</a> <strong class="jxr_keyword">protected</strong> <strong class="jxr_keyword">static</strong> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">class</strong> <a href="../../../../../../org/apache/commons/rng/internal/source32/AbstractWell.html">IndexTable</a> { +<a class="jxr_linenumber" name="L119" href="#L119">119</a> <em class="jxr_javadoccomment">/** Index indirection table giving for each index its predecessor taking table size into account. */</em> +<a class="jxr_linenumber" name="L120" href="#L120">120</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong>[] iRm1; +<a class="jxr_linenumber" name="L121" href="#L121">121</a> <em class="jxr_javadoccomment">/** Index indirection table giving for each index its second predecessor taking table size into account. */</em> +<a class="jxr_linenumber" name="L122" href="#L122">122</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong>[] iRm2; +<a class="jxr_linenumber" name="L123" href="#L123">123</a> <em class="jxr_javadoccomment">/** Index indirection table giving for each index the value index + m1 taking table size into account. */</em> +<a class="jxr_linenumber" name="L124" href="#L124">124</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong>[] i1; +<a class="jxr_linenumber" name="L125" href="#L125">125</a> <em class="jxr_javadoccomment">/** Index indirection table giving for each index the value index + m2 taking table size into account. */</em> +<a class="jxr_linenumber" name="L126" href="#L126">126</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong>[] i2; +<a class="jxr_linenumber" name="L127" href="#L127">127</a> <em class="jxr_javadoccomment">/** Index indirection table giving for each index the value index + m3 taking table size into account. */</em> +<a class="jxr_linenumber" name="L128" href="#L128">128</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong>[] i3; +<a class="jxr_linenumber" name="L129" href="#L129">129</a> +<a class="jxr_linenumber" name="L130" href="#L130">130</a> <em class="jxr_javadoccomment">/** Creates a new pre-calculated indirection index table.</em> +<a class="jxr_linenumber" name="L131" href="#L131">131</a> <em class="jxr_javadoccomment"> * @param k number of bits in the pool (not necessarily a multiple of 32)</em> +<a class="jxr_linenumber" name="L132" href="#L132">132</a> <em class="jxr_javadoccomment"> * @param m1 first parameter of the algorithm</em> +<a class="jxr_linenumber" name="L133" href="#L133">133</a> <em class="jxr_javadoccomment"> * @param m2 second parameter of the algorithm</em> +<a class="jxr_linenumber" name="L134" href="#L134">134</a> <em class="jxr_javadoccomment"> * @param m3 third parameter of the algorithm</em> +<a class="jxr_linenumber" name="L135" href="#L135">135</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L136" href="#L136">136</a> <strong class="jxr_keyword">public</strong> <a href="../../../../../../org/apache/commons/rng/internal/source32/AbstractWell.html">IndexTable</a>(<strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> k, <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> m1, <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> m2, <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> m3) { +<a class="jxr_linenumber" name="L137" href="#L137">137</a> +<a class="jxr_linenumber" name="L138" href="#L138">138</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> r = calculateBlockCount(k); +<a class="jxr_linenumber" name="L139" href="#L139">139</a> +<a class="jxr_linenumber" name="L140" href="#L140">140</a> <em class="jxr_comment">// precompute indirection index tables. These tables are used for optimizing access</em> +<a class="jxr_linenumber" name="L141" href="#L141">141</a> <em class="jxr_comment">// they allow saving computations like "(j + r - 2) % r" with costly modulo operations</em> +<a class="jxr_linenumber" name="L142" href="#L142">142</a> iRm1 = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">int</strong>[r]; +<a class="jxr_linenumber" name="L143" href="#L143">143</a> iRm2 = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">int</strong>[r]; +<a class="jxr_linenumber" name="L144" href="#L144">144</a> i1 = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">int</strong>[r]; +<a class="jxr_linenumber" name="L145" href="#L145">145</a> i2 = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">int</strong>[r]; +<a class="jxr_linenumber" name="L146" href="#L146">146</a> i3 = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">int</strong>[r]; +<a class="jxr_linenumber" name="L147" href="#L147">147</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> j = 0; j < r; ++j) { +<a class="jxr_linenumber" name="L148" href="#L148">148</a> iRm1[j] = (j + r - 1) % r; +<a class="jxr_linenumber" name="L149" href="#L149">149</a> iRm2[j] = (j + r - 2) % r; +<a class="jxr_linenumber" name="L150" href="#L150">150</a> i1[j] = (j + m1) % r; +<a class="jxr_linenumber" name="L151" href="#L151">151</a> i2[j] = (j + m2) % r; +<a class="jxr_linenumber" name="L152" href="#L152">152</a> i3[j] = (j + m3) % r; +<a class="jxr_linenumber" name="L153" href="#L153">153</a> } +<a class="jxr_linenumber" name="L154" href="#L154">154</a> } +<a class="jxr_linenumber" name="L155" href="#L155">155</a> +<a class="jxr_linenumber" name="L156" href="#L156">156</a> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L157" href="#L157">157</a> <em class="jxr_javadoccomment"> * Returns the predecessor of the given index modulo the table size.</em> +<a class="jxr_linenumber" name="L158" href="#L158">158</a> <em class="jxr_javadoccomment"> * @param index the index to look at</em> +<a class="jxr_linenumber" name="L159" href="#L159">159</a> <em class="jxr_javadoccomment"> * @return (index - 1) % table size</em> +<a class="jxr_linenumber" name="L160" href="#L160">160</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L161" href="#L161">161</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">int</strong> getIndexPred(<strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> index) { +<a class="jxr_linenumber" name="L162" href="#L162">162</a> <strong class="jxr_keyword">return</strong> iRm1[index]; +<a class="jxr_linenumber" name="L163" href="#L163">163</a> } +<a class="jxr_linenumber" name="L164" href="#L164">164</a> +<a class="jxr_linenumber" name="L165" href="#L165">165</a> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L166" href="#L166">166</a> <em class="jxr_javadoccomment"> * Returns the second predecessor of the given index modulo the table size.</em> +<a class="jxr_linenumber" name="L167" href="#L167">167</a> <em class="jxr_javadoccomment"> * @param index the index to look at</em> +<a class="jxr_linenumber" name="L168" href="#L168">168</a> <em class="jxr_javadoccomment"> * @return (index - 2) % table size</em> +<a class="jxr_linenumber" name="L169" href="#L169">169</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L170" href="#L170">170</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">int</strong> getIndexPred2(<strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> index) { +<a class="jxr_linenumber" name="L171" href="#L171">171</a> <strong class="jxr_keyword">return</strong> iRm2[index]; +<a class="jxr_linenumber" name="L172" href="#L172">172</a> } +<a class="jxr_linenumber" name="L173" href="#L173">173</a> +<a class="jxr_linenumber" name="L174" href="#L174">174</a> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L175" href="#L175">175</a> <em class="jxr_javadoccomment"> * Returns index + M1 modulo the table size.</em> +<a class="jxr_linenumber" name="L176" href="#L176">176</a> <em class="jxr_javadoccomment"> * @param index the index to look at</em> +<a class="jxr_linenumber" name="L177" href="#L177">177</a> <em class="jxr_javadoccomment"> * @return (index + M1) % table size</em> +<a class="jxr_linenumber" name="L178" href="#L178">178</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L179" href="#L179">179</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">int</strong> getIndexM1(<strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> index) { +<a class="jxr_linenumber" name="L180" href="#L180">180</a> <strong class="jxr_keyword">return</strong> i1[index]; +<a class="jxr_linenumber" name="L181" href="#L181">181</a> } +<a class="jxr_linenumber" name="L182" href="#L182">182</a> +<a class="jxr_linenumber" name="L183" href="#L183">183</a> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L184" href="#L184">184</a> <em class="jxr_javadoccomment"> * Returns index + M2 modulo the table size.</em> +<a class="jxr_linenumber" name="L185" href="#L185">185</a> <em class="jxr_javadoccomment"> * @param index the index to look at</em> +<a class="jxr_linenumber" name="L186" href="#L186">186</a> <em class="jxr_javadoccomment"> * @return (index + M2) % table size</em> +<a class="jxr_linenumber" name="L187" href="#L187">187</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L188" href="#L188">188</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">int</strong> getIndexM2(<strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> index) { +<a class="jxr_linenumber" name="L189" href="#L189">189</a> <strong class="jxr_keyword">return</strong> i2[index]; +<a class="jxr_linenumber" name="L190" href="#L190">190</a> } +<a class="jxr_linenumber" name="L191" href="#L191">191</a> +<a class="jxr_linenumber" name="L192" href="#L192">192</a> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L193" href="#L193">193</a> <em class="jxr_javadoccomment"> * Returns index + M3 modulo the table size.</em> +<a class="jxr_linenumber" name="L194" href="#L194">194</a> <em class="jxr_javadoccomment"> * @param index the index to look at</em> +<a class="jxr_linenumber" name="L195" href="#L195">195</a> <em class="jxr_javadoccomment"> * @return (index + M3) % table size</em> +<a class="jxr_linenumber" name="L196" href="#L196">196</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L197" href="#L197">197</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">int</strong> getIndexM3(<strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> index) { +<a class="jxr_linenumber" name="L198" href="#L198">198</a> <strong class="jxr_keyword">return</strong> i3[index]; +<a class="jxr_linenumber" name="L199" href="#L199">199</a> } +<a class="jxr_linenumber" name="L200" href="#L200">200</a> } +<a class="jxr_linenumber" name="L201" href="#L201">201</a> } </pre> <hr/> <div id="footer">Copyright © 2016 <a href="https://www.apache.org/">The Apache Software Foundation</a>. All rights reserved.</div>
Added: websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/internal/source32/KISSRandom.html ============================================================================== --- websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/internal/source32/KISSRandom.html (added) +++ websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/internal/source32/KISSRandom.html Fri Sep 2 16:37:56 2016 @@ -0,0 +1,144 @@ +<!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" xml:lang="en" lang="en"> +<head><meta http-equiv="content-type" content="text/html; charset=UTF-8" /> +<title>KISSRandom xref</title> +<link type="text/css" rel="stylesheet" href="../../../../../../stylesheet.css" /> +</head> +<body> +<div id="overview"><a href="../../../../../../../apidocs/org/apache/commons/rng/internal/source32/KISSRandom.html">View Javadoc</a></div><pre> +<a class="jxr_linenumber" name="L1" href="#L1">1</a> <em class="jxr_comment">/*</em> +<a class="jxr_linenumber" name="L2" href="#L2">2</a> <em class="jxr_comment"> * Licensed to the Apache Software Foundation (ASF) under one or more</em> +<a class="jxr_linenumber" name="L3" href="#L3">3</a> <em class="jxr_comment"> * contributor license agreements. See the NOTICE file distributed with</em> +<a class="jxr_linenumber" name="L4" href="#L4">4</a> <em class="jxr_comment"> * this work for additional information regarding copyright ownership.</em> +<a class="jxr_linenumber" name="L5" href="#L5">5</a> <em class="jxr_comment"> * The ASF licenses this file to You under the Apache License, Version 2.0</em> +<a class="jxr_linenumber" name="L6" href="#L6">6</a> <em class="jxr_comment"> * (the "License"); you may not use this file except in compliance with</em> +<a class="jxr_linenumber" name="L7" href="#L7">7</a> <em class="jxr_comment"> * the License. You may obtain a copy of the License at</em> +<a class="jxr_linenumber" name="L8" href="#L8">8</a> <em class="jxr_comment"> *</em> +<a class="jxr_linenumber" name="L9" href="#L9">9</a> <em class="jxr_comment"> * <a href="http://www.apache.org/licenses/LICENSE-2." target="alexandria_uri">http://www.apache.org/licenses/LICENSE-2.</a>0</em> +<a class="jxr_linenumber" name="L10" href="#L10">10</a> <em class="jxr_comment"> *</em> +<a class="jxr_linenumber" name="L11" href="#L11">11</a> <em class="jxr_comment"> * Unless required by applicable law or agreed to in writing, software</em> +<a class="jxr_linenumber" name="L12" href="#L12">12</a> <em class="jxr_comment"> * distributed under the License is distributed on an "AS IS" BASIS,</em> +<a class="jxr_linenumber" name="L13" href="#L13">13</a> <em class="jxr_comment"> * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.</em> +<a class="jxr_linenumber" name="L14" href="#L14">14</a> <em class="jxr_comment"> * See the License for the specific language governing permissions and</em> +<a class="jxr_linenumber" name="L15" href="#L15">15</a> <em class="jxr_comment"> * limitations under the License.</em> +<a class="jxr_linenumber" name="L16" href="#L16">16</a> <em class="jxr_comment"> */</em> +<a class="jxr_linenumber" name="L17" href="#L17">17</a> <strong class="jxr_keyword">package</strong> org.apache.commons.rng.internal.source32; +<a class="jxr_linenumber" name="L18" href="#L18">18</a> +<a class="jxr_linenumber" name="L19" href="#L19">19</a> <strong class="jxr_keyword">import</strong> org.apache.commons.rng.internal.util.NumberFactory; +<a class="jxr_linenumber" name="L20" href="#L20">20</a> +<a class="jxr_linenumber" name="L21" href="#L21">21</a> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L22" href="#L22">22</a> <em class="jxr_javadoccomment"> * Port from Marsaglia's <a href="<a href="http://www.cse.yorku.ca/~oz/marsaglia-rng.html" target="alexandria_uri">http://www.cse.yorku.ca/~oz/marsaglia-rng.html</a>"></em> +<a class="jxr_linenumber" name="L23" href="#L23">23</a> <em class="jxr_javadoccomment"> * "KISS" algorithm</a>.</em> +<a class="jxr_linenumber" name="L24" href="#L24">24</a> <em class="jxr_javadoccomment"> * This version contains the correction referred to</em> +<a class="jxr_linenumber" name="L25" href="#L25">25</a> <em class="jxr_javadoccomment"> * <a href="https://programmingpraxis.com/2010/10/05/george-marsaglias-random-number-generators/">here</a></em> +<a class="jxr_linenumber" name="L26" href="#L26">26</a> <em class="jxr_javadoccomment"> * in a reply to the original post.</em> +<a class="jxr_linenumber" name="L27" href="#L27">27</a> <em class="jxr_javadoccomment"> *</em> +<a class="jxr_linenumber" name="L28" href="#L28">28</a> <em class="jxr_javadoccomment"> * @since 1.0</em> +<a class="jxr_linenumber" name="L29" href="#L29">29</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L30" href="#L30">30</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">class</strong> <a href="../../../../../../org/apache/commons/rng/internal/source32/KISSRandom.html">KISSRandom</a> <strong class="jxr_keyword">extends</strong> <a href="../../../../../../org/apache/commons/rng/internal/source32/IntProvider.html">IntProvider</a> { +<a class="jxr_linenumber" name="L31" href="#L31">31</a> <em class="jxr_javadoccomment">/** Size of the seed. */</em> +<a class="jxr_linenumber" name="L32" href="#L32">32</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">static</strong> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> SEED_SIZE = 4; +<a class="jxr_linenumber" name="L33" href="#L33">33</a> <em class="jxr_javadoccomment">/** State variable. */</em> +<a class="jxr_linenumber" name="L34" href="#L34">34</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">int</strong> z; +<a class="jxr_linenumber" name="L35" href="#L35">35</a> <em class="jxr_javadoccomment">/** State variable. */</em> +<a class="jxr_linenumber" name="L36" href="#L36">36</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">int</strong> w; +<a class="jxr_linenumber" name="L37" href="#L37">37</a> <em class="jxr_javadoccomment">/** State variable. */</em> +<a class="jxr_linenumber" name="L38" href="#L38">38</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">int</strong> jsr; +<a class="jxr_linenumber" name="L39" href="#L39">39</a> <em class="jxr_javadoccomment">/** State variable. */</em> +<a class="jxr_linenumber" name="L40" href="#L40">40</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">int</strong> jcong; +<a class="jxr_linenumber" name="L41" href="#L41">41</a> +<a class="jxr_linenumber" name="L42" href="#L42">42</a> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L43" href="#L43">43</a> <em class="jxr_javadoccomment"> * Creates a new instance.</em> +<a class="jxr_linenumber" name="L44" href="#L44">44</a> <em class="jxr_javadoccomment"> *</em> +<a class="jxr_linenumber" name="L45" href="#L45">45</a> <em class="jxr_javadoccomment"> * @param seed Seed.</em> +<a class="jxr_linenumber" name="L46" href="#L46">46</a> <em class="jxr_javadoccomment"> * If the length is larger than 4, only the first 4 elements will</em> +<a class="jxr_linenumber" name="L47" href="#L47">47</a> <em class="jxr_javadoccomment"> * be used; if smaller, the remaining elements will be automatically</em> +<a class="jxr_linenumber" name="L48" href="#L48">48</a> <em class="jxr_javadoccomment"> * set.</em> +<a class="jxr_linenumber" name="L49" href="#L49">49</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L50" href="#L50">50</a> <strong class="jxr_keyword">public</strong> <a href="../../../../../../org/apache/commons/rng/internal/source32/KISSRandom.html">KISSRandom</a>(<strong class="jxr_keyword">int</strong>[] seed) { +<a class="jxr_linenumber" name="L51" href="#L51">51</a> setSeedInternal(seed); +<a class="jxr_linenumber" name="L52" href="#L52">52</a> } +<a class="jxr_linenumber" name="L53" href="#L53">53</a> +<a class="jxr_linenumber" name="L54" href="#L54">54</a> <em class="jxr_javadoccomment">/** {@inheritDoc} */</em> +<a class="jxr_linenumber" name="L55" href="#L55">55</a> @Override +<a class="jxr_linenumber" name="L56" href="#L56">56</a> <strong class="jxr_keyword">protected</strong> byte[] getStateInternal() { +<a class="jxr_linenumber" name="L57" href="#L57">57</a> <strong class="jxr_keyword">return</strong> NumberFactory.makeByteArray(<strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">int</strong>[] { z, w, jsr, jcong }); +<a class="jxr_linenumber" name="L58" href="#L58">58</a> } +<a class="jxr_linenumber" name="L59" href="#L59">59</a> +<a class="jxr_linenumber" name="L60" href="#L60">60</a> <em class="jxr_javadoccomment">/** {@inheritDoc} */</em> +<a class="jxr_linenumber" name="L61" href="#L61">61</a> @Override +<a class="jxr_linenumber" name="L62" href="#L62">62</a> <strong class="jxr_keyword">protected</strong> <strong class="jxr_keyword">void</strong> setStateInternal(byte[] s) { +<a class="jxr_linenumber" name="L63" href="#L63">63</a> checkStateSize(s, SEED_SIZE * 4); +<a class="jxr_linenumber" name="L64" href="#L64">64</a> +<a class="jxr_linenumber" name="L65" href="#L65">65</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong>[] tmp = NumberFactory.makeIntArray(s); +<a class="jxr_linenumber" name="L66" href="#L66">66</a> +<a class="jxr_linenumber" name="L67" href="#L67">67</a> z = tmp[0]; +<a class="jxr_linenumber" name="L68" href="#L68">68</a> w = tmp[1]; +<a class="jxr_linenumber" name="L69" href="#L69">69</a> jsr = tmp[2]; +<a class="jxr_linenumber" name="L70" href="#L70">70</a> jcong = tmp[3]; +<a class="jxr_linenumber" name="L71" href="#L71">71</a> } +<a class="jxr_linenumber" name="L72" href="#L72">72</a> +<a class="jxr_linenumber" name="L73" href="#L73">73</a> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L74" href="#L74">74</a> <em class="jxr_javadoccomment"> * Seeds the RNG.</em> +<a class="jxr_linenumber" name="L75" href="#L75">75</a> <em class="jxr_javadoccomment"> *</em> +<a class="jxr_linenumber" name="L76" href="#L76">76</a> <em class="jxr_javadoccomment"> * @param seed Seed.</em> +<a class="jxr_linenumber" name="L77" href="#L77">77</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L78" href="#L78">78</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">void</strong> setSeedInternal(<strong class="jxr_keyword">int</strong>[] seed) { +<a class="jxr_linenumber" name="L79" href="#L79">79</a> <em class="jxr_comment">// Reset the whole state of this RNG (i.e. the 4 state variables).</em> +<a class="jxr_linenumber" name="L80" href="#L80">80</a> <em class="jxr_comment">// Seeding procedure is not part of the reference code.</em> +<a class="jxr_linenumber" name="L81" href="#L81">81</a> +<a class="jxr_linenumber" name="L82" href="#L82">82</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong>[] tmp = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">int</strong>[SEED_SIZE]; +<a class="jxr_linenumber" name="L83" href="#L83">83</a> System.arraycopy(seed, 0, tmp, 0, Math.min(seed.length, tmp.length)); +<a class="jxr_linenumber" name="L84" href="#L84">84</a> +<a class="jxr_linenumber" name="L85" href="#L85">85</a> <strong class="jxr_keyword">if</strong> (seed.length < SEED_SIZE) { +<a class="jxr_linenumber" name="L86" href="#L86">86</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> i = seed.length; i < SEED_SIZE; i++) { +<a class="jxr_linenumber" name="L87" href="#L87">87</a> tmp[i] = 26021969 * i; +<a class="jxr_linenumber" name="L88" href="#L88">88</a> } +<a class="jxr_linenumber" name="L89" href="#L89">89</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> i = SEED_SIZE - 1; i > seed.length; i--) { +<a class="jxr_linenumber" name="L90" href="#L90">90</a> tmp[i] ^= tmp[SEED_SIZE - i - 1]; +<a class="jxr_linenumber" name="L91" href="#L91">91</a> } +<a class="jxr_linenumber" name="L92" href="#L92">92</a> +<a class="jxr_linenumber" name="L93" href="#L93">93</a> tmp[seed.length] = 0x80000000; <em class="jxr_comment">// Ensuring non-zero initial array.</em> +<a class="jxr_linenumber" name="L94" href="#L94">94</a> } +<a class="jxr_linenumber" name="L95" href="#L95">95</a> +<a class="jxr_linenumber" name="L96" href="#L96">96</a> z = tmp[0]; +<a class="jxr_linenumber" name="L97" href="#L97">97</a> w = tmp[1]; +<a class="jxr_linenumber" name="L98" href="#L98">98</a> jsr = tmp[2]; +<a class="jxr_linenumber" name="L99" href="#L99">99</a> jcong = tmp[3]; +<a class="jxr_linenumber" name="L100" href="#L100">100</a> } +<a class="jxr_linenumber" name="L101" href="#L101">101</a> +<a class="jxr_linenumber" name="L102" href="#L102">102</a> <em class="jxr_javadoccomment">/** {@inheritDoc} */</em> +<a class="jxr_linenumber" name="L103" href="#L103">103</a> @Override +<a class="jxr_linenumber" name="L104" href="#L104">104</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">int</strong> next() { +<a class="jxr_linenumber" name="L105" href="#L105">105</a> z = computeNew(36969, z); +<a class="jxr_linenumber" name="L106" href="#L106">106</a> w = computeNew(18000, w); +<a class="jxr_linenumber" name="L107" href="#L107">107</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">int</strong> mwc = (z << 16) + w; +<a class="jxr_linenumber" name="L108" href="#L108">108</a> +<a class="jxr_linenumber" name="L109" href="#L109">109</a> <em class="jxr_comment">// Cf. correction mentioned in the reply to the original post:</em> +<a class="jxr_linenumber" name="L110" href="#L110">110</a> <em class="jxr_comment">// https://programmingpraxis.com/2010/10/05/george-marsaglias-random-number-generators/</em> +<a class="jxr_linenumber" name="L111" href="#L111">111</a> jsr ^= jsr << 13; +<a class="jxr_linenumber" name="L112" href="#L112">112</a> jsr ^= jsr >>> 17; +<a class="jxr_linenumber" name="L113" href="#L113">113</a> jsr ^= jsr << 5; +<a class="jxr_linenumber" name="L114" href="#L114">114</a> +<a class="jxr_linenumber" name="L115" href="#L115">115</a> jcong = 69069 * jcong + 1234567; +<a class="jxr_linenumber" name="L116" href="#L116">116</a> +<a class="jxr_linenumber" name="L117" href="#L117">117</a> <strong class="jxr_keyword">return</strong> (mwc ^ jcong) + jsr; +<a class="jxr_linenumber" name="L118" href="#L118">118</a> } +<a class="jxr_linenumber" name="L119" href="#L119">119</a> +<a class="jxr_linenumber" name="L120" href="#L120">120</a> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L121" href="#L121">121</a> <em class="jxr_javadoccomment"> * Compute new value.</em> +<a class="jxr_linenumber" name="L122" href="#L122">122</a> <em class="jxr_javadoccomment"> *</em> +<a class="jxr_linenumber" name="L123" href="#L123">123</a> <em class="jxr_javadoccomment"> * @param mult Multiplier.</em> +<a class="jxr_linenumber" name="L124" href="#L124">124</a> <em class="jxr_javadoccomment"> * @param previous Previous value.</em> +<a class="jxr_linenumber" name="L125" href="#L125">125</a> <em class="jxr_javadoccomment"> * @return new value.</em> +<a class="jxr_linenumber" name="L126" href="#L126">126</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L127" href="#L127">127</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">int</strong> computeNew(<strong class="jxr_keyword">int</strong> mult, +<a class="jxr_linenumber" name="L128" href="#L128">128</a> <strong class="jxr_keyword">int</strong> previous) { +<a class="jxr_linenumber" name="L129" href="#L129">129</a> <strong class="jxr_keyword">return</strong> mult * (previous & 65535) + (previous >>> 16); +<a class="jxr_linenumber" name="L130" href="#L130">130</a> } +<a class="jxr_linenumber" name="L131" href="#L131">131</a> } +</pre> +<hr/> +<div id="footer">Copyright © 2016 <a href="https://www.apache.org/">The Apache Software Foundation</a>. All rights reserved.</div> +</body> +</html> \ No newline at end of file Propchange: websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/internal/source32/KISSRandom.html ------------------------------------------------------------------------------ svn:eol-style = native Modified: websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/internal/source32/MersenneTwister.html ============================================================================== --- websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/internal/source32/MersenneTwister.html (original) +++ websites/production/commons/content/proper/commons-rng/xref/org/apache/commons/rng/internal/source32/MersenneTwister.html Fri Sep 2 16:37:56 2016 @@ -143,95 +143,100 @@ <a class="jxr_linenumber" name="L135" href="#L135">135</a> } <a class="jxr_linenumber" name="L136" href="#L136">136</a> <a class="jxr_linenumber" name="L137" href="#L137">137</a> <em class="jxr_javadoccomment">/**</em> -<a class="jxr_linenumber" name="L138" href="#L138">138</a> <em class="jxr_javadoccomment"> * Reinitializes the generator as if just built with the given seed.</em> +<a class="jxr_linenumber" name="L138" href="#L138">138</a> <em class="jxr_javadoccomment"> * Initializes the generator with the given seed.</em> <a class="jxr_linenumber" name="L139" href="#L139">139</a> <em class="jxr_javadoccomment"> *</em> <a class="jxr_linenumber" name="L140" href="#L140">140</a> <em class="jxr_javadoccomment"> * @param seed Initial seed.</em> <a class="jxr_linenumber" name="L141" href="#L141">141</a> <em class="jxr_javadoccomment"> */</em> <a class="jxr_linenumber" name="L142" href="#L142">142</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">void</strong> setSeedInternal(<strong class="jxr_keyword">int</strong>[] seed) { -<a class="jxr_linenumber" name="L143" href="#L143">143</a> initState(19650218); -<a class="jxr_linenumber" name="L144" href="#L144">144</a> <strong class="jxr_keyword">int</strong> i = 1; -<a class="jxr_linenumber" name="L145" href="#L145">145</a> <strong class="jxr_keyword">int</strong> j = 0; -<a class="jxr_linenumber" name="L146" href="#L146">146</a> -<a class="jxr_linenumber" name="L147" href="#L147">147</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> k = Math.max(N, seed.length); k != 0; k--) { -<a class="jxr_linenumber" name="L148" href="#L148">148</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">long</strong> l0 = (mt[i] & LOWER_MASK_LONG) | ((mt[i] < 0) ? UPPER_MASK_LONG : 0); -<a class="jxr_linenumber" name="L149" href="#L149">149</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">long</strong> l1 = (mt[i - 1] & LOWER_MASK_LONG) | ((mt[i - 1] < 0) ? UPPER_MASK_LONG : 0); -<a class="jxr_linenumber" name="L150" href="#L150">150</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">long</strong> l = (l0 ^ ((l1 ^ (l1 >> 30)) * 1664525l)) + seed[j] + j; <em class="jxr_comment">// non linear</em> -<a class="jxr_linenumber" name="L151" href="#L151">151</a> mt[i] = (<strong class="jxr_keyword">int</strong>) (l & INT_MASK_LONG); -<a class="jxr_linenumber" name="L152" href="#L152">152</a> i++; j++; -<a class="jxr_linenumber" name="L153" href="#L153">153</a> <strong class="jxr_keyword">if</strong> (i >= N) { -<a class="jxr_linenumber" name="L154" href="#L154">154</a> mt[0] = mt[N - 1]; -<a class="jxr_linenumber" name="L155" href="#L155">155</a> i = 1; -<a class="jxr_linenumber" name="L156" href="#L156">156</a> } -<a class="jxr_linenumber" name="L157" href="#L157">157</a> <strong class="jxr_keyword">if</strong> (j >= seed.length) { -<a class="jxr_linenumber" name="L158" href="#L158">158</a> j = 0; -<a class="jxr_linenumber" name="L159" href="#L159">159</a> } -<a class="jxr_linenumber" name="L160" href="#L160">160</a> } -<a class="jxr_linenumber" name="L161" href="#L161">161</a> -<a class="jxr_linenumber" name="L162" href="#L162">162</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> k = N - 1; k != 0; k--) { -<a class="jxr_linenumber" name="L163" href="#L163">163</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">long</strong> l0 = (mt[i] & LOWER_MASK_LONG) | ((mt[i] < 0) ? UPPER_MASK_LONG : 0); -<a class="jxr_linenumber" name="L164" href="#L164">164</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">long</strong> l1 = (mt[i - 1] & LOWER_MASK_LONG) | ((mt[i - 1] < 0) ? UPPER_MASK_LONG : 0); -<a class="jxr_linenumber" name="L165" href="#L165">165</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">long</strong> l = (l0 ^ ((l1 ^ (l1 >> 30)) * 1566083941l)) - i; <em class="jxr_comment">// non linear</em> -<a class="jxr_linenumber" name="L166" href="#L166">166</a> mt[i] = (<strong class="jxr_keyword">int</strong>) (l & INT_MASK_LONG); -<a class="jxr_linenumber" name="L167" href="#L167">167</a> i++; -<a class="jxr_linenumber" name="L168" href="#L168">168</a> <strong class="jxr_keyword">if</strong> (i >= N) { -<a class="jxr_linenumber" name="L169" href="#L169">169</a> mt[0] = mt[N - 1]; -<a class="jxr_linenumber" name="L170" href="#L170">170</a> i = 1; -<a class="jxr_linenumber" name="L171" href="#L171">171</a> } -<a class="jxr_linenumber" name="L172" href="#L172">172</a> } -<a class="jxr_linenumber" name="L173" href="#L173">173</a> -<a class="jxr_linenumber" name="L174" href="#L174">174</a> mt[0] = UPPER_MASK; <em class="jxr_comment">// MSB is 1; assuring non-zero initial array</em> -<a class="jxr_linenumber" name="L175" href="#L175">175</a> } -<a class="jxr_linenumber" name="L176" href="#L176">176</a> -<a class="jxr_linenumber" name="L177" href="#L177">177</a> <em class="jxr_javadoccomment">/**</em> -<a class="jxr_linenumber" name="L178" href="#L178">178</a> <em class="jxr_javadoccomment"> * Initialize the internal state of this instance.</em> -<a class="jxr_linenumber" name="L179" href="#L179">179</a> <em class="jxr_javadoccomment"> *</em> -<a class="jxr_linenumber" name="L180" href="#L180">180</a> <em class="jxr_javadoccomment"> * @param seed Seed.</em> -<a class="jxr_linenumber" name="L181" href="#L181">181</a> <em class="jxr_javadoccomment"> */</em> -<a class="jxr_linenumber" name="L182" href="#L182">182</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">void</strong> initState(<strong class="jxr_keyword">int</strong> seed) { -<a class="jxr_linenumber" name="L183" href="#L183">183</a> <strong class="jxr_keyword">long</strong> longMT = seed & INT_MASK_LONG; -<a class="jxr_linenumber" name="L184" href="#L184">184</a> mt[0] = (<strong class="jxr_keyword">int</strong>) longMT; -<a class="jxr_linenumber" name="L185" href="#L185">185</a> <strong class="jxr_keyword">for</strong> (mti = 1; mti < N; ++mti) { -<a class="jxr_linenumber" name="L186" href="#L186">186</a> longMT = (1812433253L * (longMT ^ (longMT >> 30)) + mti) & INT_MASK_LONG; -<a class="jxr_linenumber" name="L187" href="#L187">187</a> mt[mti] = (<strong class="jxr_keyword">int</strong>) longMT; -<a class="jxr_linenumber" name="L188" href="#L188">188</a> } -<a class="jxr_linenumber" name="L189" href="#L189">189</a> } -<a class="jxr_linenumber" name="L190" href="#L190">190</a> -<a class="jxr_linenumber" name="L191" href="#L191">191</a> <em class="jxr_javadoccomment">/** {@inheritDoc} */</em> -<a class="jxr_linenumber" name="L192" href="#L192">192</a> @Override -<a class="jxr_linenumber" name="L193" href="#L193">193</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">int</strong> next() { -<a class="jxr_linenumber" name="L194" href="#L194">194</a> <strong class="jxr_keyword">int</strong> y; +<a class="jxr_linenumber" name="L143" href="#L143">143</a> <strong class="jxr_keyword">if</strong> (seed.length == 0) { +<a class="jxr_linenumber" name="L144" href="#L144">144</a> <em class="jxr_comment">// Accept empty seed.</em> +<a class="jxr_linenumber" name="L145" href="#L145">145</a> seed = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">int</strong>[1]; +<a class="jxr_linenumber" name="L146" href="#L146">146</a> } +<a class="jxr_linenumber" name="L147" href="#L147">147</a> +<a class="jxr_linenumber" name="L148" href="#L148">148</a> initState(19650218); +<a class="jxr_linenumber" name="L149" href="#L149">149</a> <strong class="jxr_keyword">int</strong> i = 1; +<a class="jxr_linenumber" name="L150" href="#L150">150</a> <strong class="jxr_keyword">int</strong> j = 0; +<a class="jxr_linenumber" name="L151" href="#L151">151</a> +<a class="jxr_linenumber" name="L152" href="#L152">152</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> k = Math.max(N, seed.length); k != 0; k--) { +<a class="jxr_linenumber" name="L153" href="#L153">153</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">long</strong> l0 = (mt[i] & LOWER_MASK_LONG) | ((mt[i] < 0) ? UPPER_MASK_LONG : 0); +<a class="jxr_linenumber" name="L154" href="#L154">154</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">long</strong> l1 = (mt[i - 1] & LOWER_MASK_LONG) | ((mt[i - 1] < 0) ? UPPER_MASK_LONG : 0); +<a class="jxr_linenumber" name="L155" href="#L155">155</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">long</strong> l = (l0 ^ ((l1 ^ (l1 >> 30)) * 1664525l)) + seed[j] + j; <em class="jxr_comment">// non linear</em> +<a class="jxr_linenumber" name="L156" href="#L156">156</a> mt[i] = (<strong class="jxr_keyword">int</strong>) (l & INT_MASK_LONG); +<a class="jxr_linenumber" name="L157" href="#L157">157</a> i++; j++; +<a class="jxr_linenumber" name="L158" href="#L158">158</a> <strong class="jxr_keyword">if</strong> (i >= N) { +<a class="jxr_linenumber" name="L159" href="#L159">159</a> mt[0] = mt[N - 1]; +<a class="jxr_linenumber" name="L160" href="#L160">160</a> i = 1; +<a class="jxr_linenumber" name="L161" href="#L161">161</a> } +<a class="jxr_linenumber" name="L162" href="#L162">162</a> <strong class="jxr_keyword">if</strong> (j >= seed.length) { +<a class="jxr_linenumber" name="L163" href="#L163">163</a> j = 0; +<a class="jxr_linenumber" name="L164" href="#L164">164</a> } +<a class="jxr_linenumber" name="L165" href="#L165">165</a> } +<a class="jxr_linenumber" name="L166" href="#L166">166</a> +<a class="jxr_linenumber" name="L167" href="#L167">167</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> k = N - 1; k != 0; k--) { +<a class="jxr_linenumber" name="L168" href="#L168">168</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">long</strong> l0 = (mt[i] & LOWER_MASK_LONG) | ((mt[i] < 0) ? UPPER_MASK_LONG : 0); +<a class="jxr_linenumber" name="L169" href="#L169">169</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">long</strong> l1 = (mt[i - 1] & LOWER_MASK_LONG) | ((mt[i - 1] < 0) ? UPPER_MASK_LONG : 0); +<a class="jxr_linenumber" name="L170" href="#L170">170</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">long</strong> l = (l0 ^ ((l1 ^ (l1 >> 30)) * 1566083941l)) - i; <em class="jxr_comment">// non linear</em> +<a class="jxr_linenumber" name="L171" href="#L171">171</a> mt[i] = (<strong class="jxr_keyword">int</strong>) (l & INT_MASK_LONG); +<a class="jxr_linenumber" name="L172" href="#L172">172</a> i++; +<a class="jxr_linenumber" name="L173" href="#L173">173</a> <strong class="jxr_keyword">if</strong> (i >= N) { +<a class="jxr_linenumber" name="L174" href="#L174">174</a> mt[0] = mt[N - 1]; +<a class="jxr_linenumber" name="L175" href="#L175">175</a> i = 1; +<a class="jxr_linenumber" name="L176" href="#L176">176</a> } +<a class="jxr_linenumber" name="L177" href="#L177">177</a> } +<a class="jxr_linenumber" name="L178" href="#L178">178</a> +<a class="jxr_linenumber" name="L179" href="#L179">179</a> mt[0] = UPPER_MASK; <em class="jxr_comment">// MSB is 1; assuring non-zero initial array</em> +<a class="jxr_linenumber" name="L180" href="#L180">180</a> } +<a class="jxr_linenumber" name="L181" href="#L181">181</a> +<a class="jxr_linenumber" name="L182" href="#L182">182</a> <em class="jxr_javadoccomment">/**</em> +<a class="jxr_linenumber" name="L183" href="#L183">183</a> <em class="jxr_javadoccomment"> * Initialize the internal state of this instance.</em> +<a class="jxr_linenumber" name="L184" href="#L184">184</a> <em class="jxr_javadoccomment"> *</em> +<a class="jxr_linenumber" name="L185" href="#L185">185</a> <em class="jxr_javadoccomment"> * @param seed Seed.</em> +<a class="jxr_linenumber" name="L186" href="#L186">186</a> <em class="jxr_javadoccomment"> */</em> +<a class="jxr_linenumber" name="L187" href="#L187">187</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">void</strong> initState(<strong class="jxr_keyword">int</strong> seed) { +<a class="jxr_linenumber" name="L188" href="#L188">188</a> <strong class="jxr_keyword">long</strong> longMT = seed & INT_MASK_LONG; +<a class="jxr_linenumber" name="L189" href="#L189">189</a> mt[0] = (<strong class="jxr_keyword">int</strong>) longMT; +<a class="jxr_linenumber" name="L190" href="#L190">190</a> <strong class="jxr_keyword">for</strong> (mti = 1; mti < N; ++mti) { +<a class="jxr_linenumber" name="L191" href="#L191">191</a> longMT = (1812433253L * (longMT ^ (longMT >> 30)) + mti) & INT_MASK_LONG; +<a class="jxr_linenumber" name="L192" href="#L192">192</a> mt[mti] = (<strong class="jxr_keyword">int</strong>) longMT; +<a class="jxr_linenumber" name="L193" href="#L193">193</a> } +<a class="jxr_linenumber" name="L194" href="#L194">194</a> } <a class="jxr_linenumber" name="L195" href="#L195">195</a> -<a class="jxr_linenumber" name="L196" href="#L196">196</a> <strong class="jxr_keyword">if</strong> (mti >= N) { <em class="jxr_comment">// Generate N words at one time.</em> -<a class="jxr_linenumber" name="L197" href="#L197">197</a> <strong class="jxr_keyword">int</strong> mtNext = mt[0]; -<a class="jxr_linenumber" name="L198" href="#L198">198</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> k = 0; k < N - M; ++k) { -<a class="jxr_linenumber" name="L199" href="#L199">199</a> <strong class="jxr_keyword">int</strong> mtCurr = mtNext; -<a class="jxr_linenumber" name="L200" href="#L200">200</a> mtNext = mt[k + 1]; -<a class="jxr_linenumber" name="L201" href="#L201">201</a> y = (mtCurr & UPPER_MASK) | (mtNext & LOWER_MASK); -<a class="jxr_linenumber" name="L202" href="#L202">202</a> mt[k] = mt[k + M] ^ (y >>> 1) ^ MAG01[y & 1]; -<a class="jxr_linenumber" name="L203" href="#L203">203</a> } -<a class="jxr_linenumber" name="L204" href="#L204">204</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> k = N - M; k < N - 1; ++k) { -<a class="jxr_linenumber" name="L205" href="#L205">205</a> <strong class="jxr_keyword">int</strong> mtCurr = mtNext; -<a class="jxr_linenumber" name="L206" href="#L206">206</a> mtNext = mt[k + 1]; -<a class="jxr_linenumber" name="L207" href="#L207">207</a> y = (mtCurr & UPPER_MASK) | (mtNext & LOWER_MASK); -<a class="jxr_linenumber" name="L208" href="#L208">208</a> mt[k] = mt[k + (M - N)] ^ (y >>> 1) ^ MAG01[y & 1]; -<a class="jxr_linenumber" name="L209" href="#L209">209</a> } -<a class="jxr_linenumber" name="L210" href="#L210">210</a> y = (mtNext & UPPER_MASK) | (mt[0] & LOWER_MASK); -<a class="jxr_linenumber" name="L211" href="#L211">211</a> mt[N - 1] = mt[M - 1] ^ (y >>> 1) ^ MAG01[y & 1]; -<a class="jxr_linenumber" name="L212" href="#L212">212</a> -<a class="jxr_linenumber" name="L213" href="#L213">213</a> mti = 0; -<a class="jxr_linenumber" name="L214" href="#L214">214</a> } -<a class="jxr_linenumber" name="L215" href="#L215">215</a> -<a class="jxr_linenumber" name="L216" href="#L216">216</a> y = mt[mti++]; +<a class="jxr_linenumber" name="L196" href="#L196">196</a> <em class="jxr_javadoccomment">/** {@inheritDoc} */</em> +<a class="jxr_linenumber" name="L197" href="#L197">197</a> @Override +<a class="jxr_linenumber" name="L198" href="#L198">198</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">int</strong> next() { +<a class="jxr_linenumber" name="L199" href="#L199">199</a> <strong class="jxr_keyword">int</strong> y; +<a class="jxr_linenumber" name="L200" href="#L200">200</a> +<a class="jxr_linenumber" name="L201" href="#L201">201</a> <strong class="jxr_keyword">if</strong> (mti >= N) { <em class="jxr_comment">// Generate N words at one time.</em> +<a class="jxr_linenumber" name="L202" href="#L202">202</a> <strong class="jxr_keyword">int</strong> mtNext = mt[0]; +<a class="jxr_linenumber" name="L203" href="#L203">203</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> k = 0; k < N - M; ++k) { +<a class="jxr_linenumber" name="L204" href="#L204">204</a> <strong class="jxr_keyword">int</strong> mtCurr = mtNext; +<a class="jxr_linenumber" name="L205" href="#L205">205</a> mtNext = mt[k + 1]; +<a class="jxr_linenumber" name="L206" href="#L206">206</a> y = (mtCurr & UPPER_MASK) | (mtNext & LOWER_MASK); +<a class="jxr_linenumber" name="L207" href="#L207">207</a> mt[k] = mt[k + M] ^ (y >>> 1) ^ MAG01[y & 1]; +<a class="jxr_linenumber" name="L208" href="#L208">208</a> } +<a class="jxr_linenumber" name="L209" href="#L209">209</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> k = N - M; k < N - 1; ++k) { +<a class="jxr_linenumber" name="L210" href="#L210">210</a> <strong class="jxr_keyword">int</strong> mtCurr = mtNext; +<a class="jxr_linenumber" name="L211" href="#L211">211</a> mtNext = mt[k + 1]; +<a class="jxr_linenumber" name="L212" href="#L212">212</a> y = (mtCurr & UPPER_MASK) | (mtNext & LOWER_MASK); +<a class="jxr_linenumber" name="L213" href="#L213">213</a> mt[k] = mt[k + (M - N)] ^ (y >>> 1) ^ MAG01[y & 1]; +<a class="jxr_linenumber" name="L214" href="#L214">214</a> } +<a class="jxr_linenumber" name="L215" href="#L215">215</a> y = (mtNext & UPPER_MASK) | (mt[0] & LOWER_MASK); +<a class="jxr_linenumber" name="L216" href="#L216">216</a> mt[N - 1] = mt[M - 1] ^ (y >>> 1) ^ MAG01[y & 1]; <a class="jxr_linenumber" name="L217" href="#L217">217</a> -<a class="jxr_linenumber" name="L218" href="#L218">218</a> <em class="jxr_comment">// Tempering.</em> -<a class="jxr_linenumber" name="L219" href="#L219">219</a> y ^= y >>> 11; -<a class="jxr_linenumber" name="L220" href="#L220">220</a> y ^= (y << 7) & 0x9d2c5680; -<a class="jxr_linenumber" name="L221" href="#L221">221</a> y ^= (y << 15) & 0xefc60000; -<a class="jxr_linenumber" name="L222" href="#L222">222</a> y ^= y >>> 18; -<a class="jxr_linenumber" name="L223" href="#L223">223</a> -<a class="jxr_linenumber" name="L224" href="#L224">224</a> <strong class="jxr_keyword">return</strong> y; -<a class="jxr_linenumber" name="L225" href="#L225">225</a> } -<a class="jxr_linenumber" name="L226" href="#L226">226</a> } +<a class="jxr_linenumber" name="L218" href="#L218">218</a> mti = 0; +<a class="jxr_linenumber" name="L219" href="#L219">219</a> } +<a class="jxr_linenumber" name="L220" href="#L220">220</a> +<a class="jxr_linenumber" name="L221" href="#L221">221</a> y = mt[mti++]; +<a class="jxr_linenumber" name="L222" href="#L222">222</a> +<a class="jxr_linenumber" name="L223" href="#L223">223</a> <em class="jxr_comment">// Tempering.</em> +<a class="jxr_linenumber" name="L224" href="#L224">224</a> y ^= y >>> 11; +<a class="jxr_linenumber" name="L225" href="#L225">225</a> y ^= (y << 7) & 0x9d2c5680; +<a class="jxr_linenumber" name="L226" href="#L226">226</a> y ^= (y << 15) & 0xefc60000; +<a class="jxr_linenumber" name="L227" href="#L227">227</a> y ^= y >>> 18; +<a class="jxr_linenumber" name="L228" href="#L228">228</a> +<a class="jxr_linenumber" name="L229" href="#L229">229</a> <strong class="jxr_keyword">return</strong> y; +<a class="jxr_linenumber" name="L230" href="#L230">230</a> } +<a class="jxr_linenumber" name="L231" href="#L231">231</a> } </pre> <hr/> <div id="footer">Copyright © 2016 <a href="https://www.apache.org/">The Apache Software Foundation</a>. All rights reserved.</div>