| <!DOCTYPE html><html lang="en"><head><meta charset="utf-8"><meta name="viewport" content="width=device-width, initial-scale=1.0"><meta name="generator" content="rustdoc"><meta name="description" content="Source of the Rust file `/root/.cargo/git/checkouts/incubator-teaclave-crates-c8106113f74feefc/ede1f68/rand/src/distributions/uniform.rs`."><meta name="keywords" content="rust, rustlang, rust-lang"><title>uniform.rs - source</title><link rel="preload" as="font" type="font/woff2" crossorigin href="../../../SourceSerif4-Regular.ttf.woff2"><link rel="preload" as="font" type="font/woff2" crossorigin href="../../../FiraSans-Regular.woff2"><link rel="preload" as="font" type="font/woff2" crossorigin href="../../../FiraSans-Medium.woff2"><link rel="preload" as="font" type="font/woff2" crossorigin href="../../../SourceCodePro-Regular.ttf.woff2"><link rel="preload" as="font" type="font/woff2" crossorigin href="../../../SourceSerif4-Bold.ttf.woff2"><link rel="preload" as="font" type="font/woff2" crossorigin href="../../../SourceCodePro-Semibold.ttf.woff2"><link rel="stylesheet" href="../../../normalize.css"><link rel="stylesheet" href="../../../rustdoc.css" id="mainThemeStyle"><link rel="stylesheet" href="../../../ayu.css" disabled><link rel="stylesheet" href="../../../dark.css" disabled><link rel="stylesheet" href="../../../light.css" id="themeStyle"><script id="default-settings" ></script><script src="../../../storage.js"></script><script defer src="../../../source-script.js"></script><script defer src="../../../source-files.js"></script><script defer src="../../../main.js"></script><noscript><link rel="stylesheet" href="../../../noscript.css"></noscript><link rel="icon" href="https://www.rust-lang.org/favicon.ico"></head><body class="rustdoc source"><!--[if lte IE 11]><div class="warning">This old browser is unsupported and will most likely display funky things.</div><![endif]--><nav class="sidebar"><a class="sidebar-logo" href="../../../rand/index.html"><div class="logo-container"> |
| <img src="https://www.rust-lang.org/logos/rust-logo-128x128-blk.png" alt="logo"></div></a></nav><main><div class="width-limiter"><nav class="sub"><a class="sub-logo-container" href="../../../rand/index.html"> |
| <img src="https://www.rust-lang.org/logos/rust-logo-128x128-blk.png" alt="logo"></a><form class="search-form"><div class="search-container"><span></span><input class="search-input" name="search" autocomplete="off" spellcheck="false" placeholder="Click or press ‘S’ to search, ‘?’ for more options…" type="search"><div id="help-button" title="help" tabindex="-1"><a href="../../../help.html">?</a></div><div id="settings-menu" tabindex="-1"><a href="../../../settings.html" title="settings"><img width="22" height="22" alt="Change settings" src="../../../wheel.svg"></a></div></div></form></nav><section id="main-content" class="content"><div class="example-wrap"><pre class="src-line-numbers"><span id="1">1</span> |
| <span id="2">2</span> |
| <span id="3">3</span> |
| <span id="4">4</span> |
| <span id="5">5</span> |
| <span id="6">6</span> |
| <span id="7">7</span> |
| <span id="8">8</span> |
| <span id="9">9</span> |
| <span id="10">10</span> |
| <span id="11">11</span> |
| <span id="12">12</span> |
| <span id="13">13</span> |
| <span id="14">14</span> |
| <span id="15">15</span> |
| <span id="16">16</span> |
| <span id="17">17</span> |
| <span id="18">18</span> |
| <span id="19">19</span> |
| <span id="20">20</span> |
| <span id="21">21</span> |
| <span id="22">22</span> |
| <span id="23">23</span> |
| <span id="24">24</span> |
| <span id="25">25</span> |
| <span id="26">26</span> |
| <span id="27">27</span> |
| <span id="28">28</span> |
| <span id="29">29</span> |
| <span id="30">30</span> |
| <span id="31">31</span> |
| <span id="32">32</span> |
| <span id="33">33</span> |
| <span id="34">34</span> |
| <span id="35">35</span> |
| <span id="36">36</span> |
| <span id="37">37</span> |
| <span id="38">38</span> |
| <span id="39">39</span> |
| <span id="40">40</span> |
| <span id="41">41</span> |
| <span id="42">42</span> |
| <span id="43">43</span> |
| <span id="44">44</span> |
| <span id="45">45</span> |
| <span id="46">46</span> |
| <span id="47">47</span> |
| <span id="48">48</span> |
| <span id="49">49</span> |
| <span id="50">50</span> |
| <span id="51">51</span> |
| <span id="52">52</span> |
| <span id="53">53</span> |
| <span id="54">54</span> |
| <span id="55">55</span> |
| <span id="56">56</span> |
| <span id="57">57</span> |
| <span id="58">58</span> |
| <span id="59">59</span> |
| <span id="60">60</span> |
| <span id="61">61</span> |
| <span id="62">62</span> |
| <span id="63">63</span> |
| <span id="64">64</span> |
| <span id="65">65</span> |
| <span id="66">66</span> |
| <span id="67">67</span> |
| <span id="68">68</span> |
| <span id="69">69</span> |
| <span id="70">70</span> |
| <span id="71">71</span> |
| <span id="72">72</span> |
| <span id="73">73</span> |
| <span id="74">74</span> |
| <span id="75">75</span> |
| <span id="76">76</span> |
| <span id="77">77</span> |
| <span id="78">78</span> |
| <span id="79">79</span> |
| <span id="80">80</span> |
| <span id="81">81</span> |
| <span id="82">82</span> |
| <span id="83">83</span> |
| <span id="84">84</span> |
| <span id="85">85</span> |
| <span id="86">86</span> |
| <span id="87">87</span> |
| <span id="88">88</span> |
| <span id="89">89</span> |
| <span id="90">90</span> |
| <span id="91">91</span> |
| <span id="92">92</span> |
| <span id="93">93</span> |
| <span id="94">94</span> |
| <span id="95">95</span> |
| <span id="96">96</span> |
| <span id="97">97</span> |
| <span id="98">98</span> |
| <span id="99">99</span> |
| <span id="100">100</span> |
| <span id="101">101</span> |
| <span id="102">102</span> |
| <span id="103">103</span> |
| <span id="104">104</span> |
| <span id="105">105</span> |
| <span id="106">106</span> |
| <span id="107">107</span> |
| <span id="108">108</span> |
| <span id="109">109</span> |
| <span id="110">110</span> |
| <span id="111">111</span> |
| <span id="112">112</span> |
| <span id="113">113</span> |
| <span id="114">114</span> |
| <span id="115">115</span> |
| <span id="116">116</span> |
| <span id="117">117</span> |
| <span id="118">118</span> |
| <span id="119">119</span> |
| <span id="120">120</span> |
| <span id="121">121</span> |
| <span id="122">122</span> |
| <span id="123">123</span> |
| <span id="124">124</span> |
| <span id="125">125</span> |
| <span id="126">126</span> |
| <span id="127">127</span> |
| <span id="128">128</span> |
| <span id="129">129</span> |
| <span id="130">130</span> |
| <span id="131">131</span> |
| <span id="132">132</span> |
| <span id="133">133</span> |
| <span id="134">134</span> |
| <span id="135">135</span> |
| <span id="136">136</span> |
| <span id="137">137</span> |
| <span id="138">138</span> |
| <span id="139">139</span> |
| <span id="140">140</span> |
| <span id="141">141</span> |
| <span id="142">142</span> |
| <span id="143">143</span> |
| <span id="144">144</span> |
| <span id="145">145</span> |
| <span id="146">146</span> |
| <span id="147">147</span> |
| <span id="148">148</span> |
| <span id="149">149</span> |
| <span id="150">150</span> |
| <span id="151">151</span> |
| <span id="152">152</span> |
| <span id="153">153</span> |
| <span id="154">154</span> |
| <span id="155">155</span> |
| <span id="156">156</span> |
| <span id="157">157</span> |
| <span id="158">158</span> |
| <span id="159">159</span> |
| <span id="160">160</span> |
| <span id="161">161</span> |
| <span id="162">162</span> |
| <span id="163">163</span> |
| <span id="164">164</span> |
| <span id="165">165</span> |
| <span id="166">166</span> |
| <span id="167">167</span> |
| <span id="168">168</span> |
| <span id="169">169</span> |
| <span id="170">170</span> |
| <span id="171">171</span> |
| <span id="172">172</span> |
| <span id="173">173</span> |
| <span id="174">174</span> |
| <span id="175">175</span> |
| <span id="176">176</span> |
| <span id="177">177</span> |
| <span id="178">178</span> |
| <span id="179">179</span> |
| <span id="180">180</span> |
| <span id="181">181</span> |
| <span id="182">182</span> |
| <span id="183">183</span> |
| <span id="184">184</span> |
| <span id="185">185</span> |
| <span id="186">186</span> |
| <span id="187">187</span> |
| <span id="188">188</span> |
| <span id="189">189</span> |
| <span id="190">190</span> |
| <span id="191">191</span> |
| <span id="192">192</span> |
| <span id="193">193</span> |
| <span id="194">194</span> |
| <span id="195">195</span> |
| <span id="196">196</span> |
| <span id="197">197</span> |
| <span id="198">198</span> |
| <span id="199">199</span> |
| <span id="200">200</span> |
| <span id="201">201</span> |
| <span id="202">202</span> |
| <span id="203">203</span> |
| <span id="204">204</span> |
| <span id="205">205</span> |
| <span id="206">206</span> |
| <span id="207">207</span> |
| <span id="208">208</span> |
| <span id="209">209</span> |
| <span id="210">210</span> |
| <span id="211">211</span> |
| <span id="212">212</span> |
| <span id="213">213</span> |
| <span id="214">214</span> |
| <span id="215">215</span> |
| <span id="216">216</span> |
| <span id="217">217</span> |
| <span id="218">218</span> |
| <span id="219">219</span> |
| <span id="220">220</span> |
| <span id="221">221</span> |
| <span id="222">222</span> |
| <span id="223">223</span> |
| <span id="224">224</span> |
| <span id="225">225</span> |
| <span id="226">226</span> |
| <span id="227">227</span> |
| <span id="228">228</span> |
| <span id="229">229</span> |
| <span id="230">230</span> |
| <span id="231">231</span> |
| <span id="232">232</span> |
| <span id="233">233</span> |
| <span id="234">234</span> |
| <span id="235">235</span> |
| <span id="236">236</span> |
| <span id="237">237</span> |
| <span id="238">238</span> |
| <span id="239">239</span> |
| <span id="240">240</span> |
| <span id="241">241</span> |
| <span id="242">242</span> |
| <span id="243">243</span> |
| <span id="244">244</span> |
| <span id="245">245</span> |
| <span id="246">246</span> |
| <span id="247">247</span> |
| <span id="248">248</span> |
| <span id="249">249</span> |
| <span id="250">250</span> |
| <span id="251">251</span> |
| <span id="252">252</span> |
| <span id="253">253</span> |
| <span id="254">254</span> |
| <span id="255">255</span> |
| <span id="256">256</span> |
| <span id="257">257</span> |
| <span id="258">258</span> |
| <span id="259">259</span> |
| <span id="260">260</span> |
| <span id="261">261</span> |
| <span id="262">262</span> |
| <span id="263">263</span> |
| <span id="264">264</span> |
| <span id="265">265</span> |
| <span id="266">266</span> |
| <span id="267">267</span> |
| <span id="268">268</span> |
| <span id="269">269</span> |
| <span id="270">270</span> |
| <span id="271">271</span> |
| <span id="272">272</span> |
| <span id="273">273</span> |
| <span id="274">274</span> |
| <span id="275">275</span> |
| <span id="276">276</span> |
| <span id="277">277</span> |
| <span id="278">278</span> |
| <span id="279">279</span> |
| <span id="280">280</span> |
| <span id="281">281</span> |
| <span id="282">282</span> |
| <span id="283">283</span> |
| <span id="284">284</span> |
| <span id="285">285</span> |
| <span id="286">286</span> |
| <span id="287">287</span> |
| <span id="288">288</span> |
| <span id="289">289</span> |
| <span id="290">290</span> |
| <span id="291">291</span> |
| <span id="292">292</span> |
| <span id="293">293</span> |
| <span id="294">294</span> |
| <span id="295">295</span> |
| <span id="296">296</span> |
| <span id="297">297</span> |
| <span id="298">298</span> |
| <span id="299">299</span> |
| <span id="300">300</span> |
| <span id="301">301</span> |
| <span id="302">302</span> |
| <span id="303">303</span> |
| <span id="304">304</span> |
| <span id="305">305</span> |
| <span id="306">306</span> |
| <span id="307">307</span> |
| <span id="308">308</span> |
| <span id="309">309</span> |
| <span id="310">310</span> |
| <span id="311">311</span> |
| <span id="312">312</span> |
| <span id="313">313</span> |
| <span id="314">314</span> |
| <span id="315">315</span> |
| <span id="316">316</span> |
| <span id="317">317</span> |
| <span id="318">318</span> |
| <span id="319">319</span> |
| <span id="320">320</span> |
| <span id="321">321</span> |
| <span id="322">322</span> |
| <span id="323">323</span> |
| <span id="324">324</span> |
| <span id="325">325</span> |
| <span id="326">326</span> |
| <span id="327">327</span> |
| <span id="328">328</span> |
| <span id="329">329</span> |
| <span id="330">330</span> |
| <span id="331">331</span> |
| <span id="332">332</span> |
| <span id="333">333</span> |
| <span id="334">334</span> |
| <span id="335">335</span> |
| <span id="336">336</span> |
| <span id="337">337</span> |
| <span id="338">338</span> |
| <span id="339">339</span> |
| <span id="340">340</span> |
| <span id="341">341</span> |
| <span id="342">342</span> |
| <span id="343">343</span> |
| <span id="344">344</span> |
| <span id="345">345</span> |
| <span id="346">346</span> |
| <span id="347">347</span> |
| <span id="348">348</span> |
| <span id="349">349</span> |
| <span id="350">350</span> |
| <span id="351">351</span> |
| <span id="352">352</span> |
| <span id="353">353</span> |
| <span id="354">354</span> |
| <span id="355">355</span> |
| <span id="356">356</span> |
| <span id="357">357</span> |
| <span id="358">358</span> |
| <span id="359">359</span> |
| <span id="360">360</span> |
| <span id="361">361</span> |
| <span id="362">362</span> |
| <span id="363">363</span> |
| <span id="364">364</span> |
| <span id="365">365</span> |
| <span id="366">366</span> |
| <span id="367">367</span> |
| <span id="368">368</span> |
| <span id="369">369</span> |
| <span id="370">370</span> |
| <span id="371">371</span> |
| <span id="372">372</span> |
| <span id="373">373</span> |
| <span id="374">374</span> |
| <span id="375">375</span> |
| <span id="376">376</span> |
| <span id="377">377</span> |
| <span id="378">378</span> |
| <span id="379">379</span> |
| <span id="380">380</span> |
| <span id="381">381</span> |
| <span id="382">382</span> |
| <span id="383">383</span> |
| <span id="384">384</span> |
| <span id="385">385</span> |
| <span id="386">386</span> |
| <span id="387">387</span> |
| <span id="388">388</span> |
| <span id="389">389</span> |
| <span id="390">390</span> |
| <span id="391">391</span> |
| <span id="392">392</span> |
| <span id="393">393</span> |
| <span id="394">394</span> |
| <span id="395">395</span> |
| <span id="396">396</span> |
| <span id="397">397</span> |
| <span id="398">398</span> |
| <span id="399">399</span> |
| <span id="400">400</span> |
| <span id="401">401</span> |
| <span id="402">402</span> |
| <span id="403">403</span> |
| <span id="404">404</span> |
| <span id="405">405</span> |
| <span id="406">406</span> |
| <span id="407">407</span> |
| <span id="408">408</span> |
| <span id="409">409</span> |
| <span id="410">410</span> |
| <span id="411">411</span> |
| <span id="412">412</span> |
| <span id="413">413</span> |
| <span id="414">414</span> |
| <span id="415">415</span> |
| <span id="416">416</span> |
| <span id="417">417</span> |
| <span id="418">418</span> |
| <span id="419">419</span> |
| <span id="420">420</span> |
| <span id="421">421</span> |
| <span id="422">422</span> |
| <span id="423">423</span> |
| <span id="424">424</span> |
| <span id="425">425</span> |
| <span id="426">426</span> |
| <span id="427">427</span> |
| <span id="428">428</span> |
| <span id="429">429</span> |
| <span id="430">430</span> |
| <span id="431">431</span> |
| <span id="432">432</span> |
| <span id="433">433</span> |
| <span id="434">434</span> |
| <span id="435">435</span> |
| <span id="436">436</span> |
| <span id="437">437</span> |
| <span id="438">438</span> |
| <span id="439">439</span> |
| <span id="440">440</span> |
| <span id="441">441</span> |
| <span id="442">442</span> |
| <span id="443">443</span> |
| <span id="444">444</span> |
| <span id="445">445</span> |
| <span id="446">446</span> |
| <span id="447">447</span> |
| <span id="448">448</span> |
| <span id="449">449</span> |
| <span id="450">450</span> |
| <span id="451">451</span> |
| <span id="452">452</span> |
| <span id="453">453</span> |
| <span id="454">454</span> |
| <span id="455">455</span> |
| <span id="456">456</span> |
| <span id="457">457</span> |
| <span id="458">458</span> |
| <span id="459">459</span> |
| <span id="460">460</span> |
| <span id="461">461</span> |
| <span id="462">462</span> |
| <span id="463">463</span> |
| <span id="464">464</span> |
| <span id="465">465</span> |
| <span id="466">466</span> |
| <span id="467">467</span> |
| <span id="468">468</span> |
| <span id="469">469</span> |
| <span id="470">470</span> |
| <span id="471">471</span> |
| <span id="472">472</span> |
| <span id="473">473</span> |
| <span id="474">474</span> |
| <span id="475">475</span> |
| <span id="476">476</span> |
| <span id="477">477</span> |
| <span id="478">478</span> |
| <span id="479">479</span> |
| <span id="480">480</span> |
| <span id="481">481</span> |
| <span id="482">482</span> |
| <span id="483">483</span> |
| <span id="484">484</span> |
| <span id="485">485</span> |
| <span id="486">486</span> |
| <span id="487">487</span> |
| <span id="488">488</span> |
| <span id="489">489</span> |
| <span id="490">490</span> |
| <span id="491">491</span> |
| <span id="492">492</span> |
| <span id="493">493</span> |
| <span id="494">494</span> |
| <span id="495">495</span> |
| <span id="496">496</span> |
| <span id="497">497</span> |
| <span id="498">498</span> |
| <span id="499">499</span> |
| <span id="500">500</span> |
| <span id="501">501</span> |
| <span id="502">502</span> |
| <span id="503">503</span> |
| <span id="504">504</span> |
| <span id="505">505</span> |
| <span id="506">506</span> |
| <span id="507">507</span> |
| <span id="508">508</span> |
| <span id="509">509</span> |
| <span id="510">510</span> |
| <span id="511">511</span> |
| <span id="512">512</span> |
| <span id="513">513</span> |
| <span id="514">514</span> |
| <span id="515">515</span> |
| <span id="516">516</span> |
| <span id="517">517</span> |
| <span id="518">518</span> |
| <span id="519">519</span> |
| <span id="520">520</span> |
| <span id="521">521</span> |
| <span id="522">522</span> |
| <span id="523">523</span> |
| <span id="524">524</span> |
| <span id="525">525</span> |
| <span id="526">526</span> |
| <span id="527">527</span> |
| <span id="528">528</span> |
| <span id="529">529</span> |
| <span id="530">530</span> |
| <span id="531">531</span> |
| <span id="532">532</span> |
| <span id="533">533</span> |
| <span id="534">534</span> |
| <span id="535">535</span> |
| <span id="536">536</span> |
| <span id="537">537</span> |
| <span id="538">538</span> |
| <span id="539">539</span> |
| <span id="540">540</span> |
| <span id="541">541</span> |
| <span id="542">542</span> |
| <span id="543">543</span> |
| <span id="544">544</span> |
| <span id="545">545</span> |
| <span id="546">546</span> |
| <span id="547">547</span> |
| <span id="548">548</span> |
| <span id="549">549</span> |
| <span id="550">550</span> |
| <span id="551">551</span> |
| <span id="552">552</span> |
| <span id="553">553</span> |
| <span id="554">554</span> |
| <span id="555">555</span> |
| <span id="556">556</span> |
| <span id="557">557</span> |
| <span id="558">558</span> |
| <span id="559">559</span> |
| <span id="560">560</span> |
| <span id="561">561</span> |
| <span id="562">562</span> |
| <span id="563">563</span> |
| <span id="564">564</span> |
| <span id="565">565</span> |
| <span id="566">566</span> |
| <span id="567">567</span> |
| <span id="568">568</span> |
| <span id="569">569</span> |
| <span id="570">570</span> |
| <span id="571">571</span> |
| <span id="572">572</span> |
| <span id="573">573</span> |
| <span id="574">574</span> |
| <span id="575">575</span> |
| <span id="576">576</span> |
| <span id="577">577</span> |
| <span id="578">578</span> |
| <span id="579">579</span> |
| <span id="580">580</span> |
| <span id="581">581</span> |
| <span id="582">582</span> |
| <span id="583">583</span> |
| <span id="584">584</span> |
| <span id="585">585</span> |
| <span id="586">586</span> |
| <span id="587">587</span> |
| <span id="588">588</span> |
| <span id="589">589</span> |
| <span id="590">590</span> |
| <span id="591">591</span> |
| <span id="592">592</span> |
| <span id="593">593</span> |
| <span id="594">594</span> |
| <span id="595">595</span> |
| <span id="596">596</span> |
| <span id="597">597</span> |
| <span id="598">598</span> |
| <span id="599">599</span> |
| <span id="600">600</span> |
| <span id="601">601</span> |
| <span id="602">602</span> |
| <span id="603">603</span> |
| <span id="604">604</span> |
| <span id="605">605</span> |
| <span id="606">606</span> |
| <span id="607">607</span> |
| <span id="608">608</span> |
| <span id="609">609</span> |
| <span id="610">610</span> |
| <span id="611">611</span> |
| <span id="612">612</span> |
| <span id="613">613</span> |
| <span id="614">614</span> |
| <span id="615">615</span> |
| <span id="616">616</span> |
| <span id="617">617</span> |
| <span id="618">618</span> |
| <span id="619">619</span> |
| <span id="620">620</span> |
| <span id="621">621</span> |
| <span id="622">622</span> |
| <span id="623">623</span> |
| <span id="624">624</span> |
| <span id="625">625</span> |
| <span id="626">626</span> |
| <span id="627">627</span> |
| <span id="628">628</span> |
| <span id="629">629</span> |
| <span id="630">630</span> |
| <span id="631">631</span> |
| <span id="632">632</span> |
| <span id="633">633</span> |
| <span id="634">634</span> |
| <span id="635">635</span> |
| <span id="636">636</span> |
| <span id="637">637</span> |
| <span id="638">638</span> |
| <span id="639">639</span> |
| <span id="640">640</span> |
| <span id="641">641</span> |
| <span id="642">642</span> |
| <span id="643">643</span> |
| <span id="644">644</span> |
| <span id="645">645</span> |
| <span id="646">646</span> |
| <span id="647">647</span> |
| <span id="648">648</span> |
| <span id="649">649</span> |
| <span id="650">650</span> |
| <span id="651">651</span> |
| <span id="652">652</span> |
| <span id="653">653</span> |
| <span id="654">654</span> |
| <span id="655">655</span> |
| <span id="656">656</span> |
| <span id="657">657</span> |
| <span id="658">658</span> |
| <span id="659">659</span> |
| <span id="660">660</span> |
| <span id="661">661</span> |
| <span id="662">662</span> |
| <span id="663">663</span> |
| <span id="664">664</span> |
| <span id="665">665</span> |
| <span id="666">666</span> |
| <span id="667">667</span> |
| <span id="668">668</span> |
| <span id="669">669</span> |
| <span id="670">670</span> |
| <span id="671">671</span> |
| <span id="672">672</span> |
| <span id="673">673</span> |
| <span id="674">674</span> |
| <span id="675">675</span> |
| <span id="676">676</span> |
| <span id="677">677</span> |
| <span id="678">678</span> |
| <span id="679">679</span> |
| <span id="680">680</span> |
| <span id="681">681</span> |
| <span id="682">682</span> |
| <span id="683">683</span> |
| <span id="684">684</span> |
| <span id="685">685</span> |
| <span id="686">686</span> |
| <span id="687">687</span> |
| <span id="688">688</span> |
| <span id="689">689</span> |
| <span id="690">690</span> |
| <span id="691">691</span> |
| <span id="692">692</span> |
| <span id="693">693</span> |
| <span id="694">694</span> |
| <span id="695">695</span> |
| <span id="696">696</span> |
| <span id="697">697</span> |
| <span id="698">698</span> |
| <span id="699">699</span> |
| <span id="700">700</span> |
| <span id="701">701</span> |
| <span id="702">702</span> |
| <span id="703">703</span> |
| <span id="704">704</span> |
| <span id="705">705</span> |
| <span id="706">706</span> |
| <span id="707">707</span> |
| <span id="708">708</span> |
| <span id="709">709</span> |
| <span id="710">710</span> |
| <span id="711">711</span> |
| <span id="712">712</span> |
| <span id="713">713</span> |
| <span id="714">714</span> |
| <span id="715">715</span> |
| <span id="716">716</span> |
| <span id="717">717</span> |
| <span id="718">718</span> |
| <span id="719">719</span> |
| <span id="720">720</span> |
| <span id="721">721</span> |
| <span id="722">722</span> |
| <span id="723">723</span> |
| <span id="724">724</span> |
| <span id="725">725</span> |
| <span id="726">726</span> |
| <span id="727">727</span> |
| <span id="728">728</span> |
| <span id="729">729</span> |
| <span id="730">730</span> |
| <span id="731">731</span> |
| <span id="732">732</span> |
| <span id="733">733</span> |
| <span id="734">734</span> |
| <span id="735">735</span> |
| <span id="736">736</span> |
| <span id="737">737</span> |
| <span id="738">738</span> |
| <span id="739">739</span> |
| <span id="740">740</span> |
| <span id="741">741</span> |
| <span id="742">742</span> |
| <span id="743">743</span> |
| <span id="744">744</span> |
| <span id="745">745</span> |
| <span id="746">746</span> |
| <span id="747">747</span> |
| <span id="748">748</span> |
| <span id="749">749</span> |
| <span id="750">750</span> |
| <span id="751">751</span> |
| <span id="752">752</span> |
| <span id="753">753</span> |
| <span id="754">754</span> |
| <span id="755">755</span> |
| <span id="756">756</span> |
| <span id="757">757</span> |
| <span id="758">758</span> |
| <span id="759">759</span> |
| <span id="760">760</span> |
| <span id="761">761</span> |
| <span id="762">762</span> |
| <span id="763">763</span> |
| <span id="764">764</span> |
| <span id="765">765</span> |
| <span id="766">766</span> |
| <span id="767">767</span> |
| <span id="768">768</span> |
| <span id="769">769</span> |
| <span id="770">770</span> |
| <span id="771">771</span> |
| <span id="772">772</span> |
| <span id="773">773</span> |
| <span id="774">774</span> |
| <span id="775">775</span> |
| <span id="776">776</span> |
| <span id="777">777</span> |
| <span id="778">778</span> |
| <span id="779">779</span> |
| <span id="780">780</span> |
| <span id="781">781</span> |
| <span id="782">782</span> |
| <span id="783">783</span> |
| <span id="784">784</span> |
| <span id="785">785</span> |
| <span id="786">786</span> |
| <span id="787">787</span> |
| <span id="788">788</span> |
| <span id="789">789</span> |
| <span id="790">790</span> |
| <span id="791">791</span> |
| <span id="792">792</span> |
| <span id="793">793</span> |
| <span id="794">794</span> |
| <span id="795">795</span> |
| <span id="796">796</span> |
| <span id="797">797</span> |
| <span id="798">798</span> |
| <span id="799">799</span> |
| <span id="800">800</span> |
| <span id="801">801</span> |
| <span id="802">802</span> |
| <span id="803">803</span> |
| <span id="804">804</span> |
| <span id="805">805</span> |
| <span id="806">806</span> |
| <span id="807">807</span> |
| <span id="808">808</span> |
| <span id="809">809</span> |
| <span id="810">810</span> |
| <span id="811">811</span> |
| <span id="812">812</span> |
| <span id="813">813</span> |
| <span id="814">814</span> |
| <span id="815">815</span> |
| <span id="816">816</span> |
| <span id="817">817</span> |
| <span id="818">818</span> |
| <span id="819">819</span> |
| <span id="820">820</span> |
| <span id="821">821</span> |
| <span id="822">822</span> |
| <span id="823">823</span> |
| <span id="824">824</span> |
| <span id="825">825</span> |
| <span id="826">826</span> |
| <span id="827">827</span> |
| <span id="828">828</span> |
| <span id="829">829</span> |
| <span id="830">830</span> |
| <span id="831">831</span> |
| <span id="832">832</span> |
| <span id="833">833</span> |
| <span id="834">834</span> |
| <span id="835">835</span> |
| <span id="836">836</span> |
| <span id="837">837</span> |
| <span id="838">838</span> |
| <span id="839">839</span> |
| <span id="840">840</span> |
| <span id="841">841</span> |
| <span id="842">842</span> |
| <span id="843">843</span> |
| <span id="844">844</span> |
| <span id="845">845</span> |
| <span id="846">846</span> |
| <span id="847">847</span> |
| <span id="848">848</span> |
| <span id="849">849</span> |
| <span id="850">850</span> |
| <span id="851">851</span> |
| <span id="852">852</span> |
| <span id="853">853</span> |
| <span id="854">854</span> |
| <span id="855">855</span> |
| <span id="856">856</span> |
| <span id="857">857</span> |
| <span id="858">858</span> |
| <span id="859">859</span> |
| <span id="860">860</span> |
| <span id="861">861</span> |
| <span id="862">862</span> |
| <span id="863">863</span> |
| <span id="864">864</span> |
| <span id="865">865</span> |
| <span id="866">866</span> |
| <span id="867">867</span> |
| <span id="868">868</span> |
| <span id="869">869</span> |
| <span id="870">870</span> |
| <span id="871">871</span> |
| <span id="872">872</span> |
| <span id="873">873</span> |
| <span id="874">874</span> |
| <span id="875">875</span> |
| <span id="876">876</span> |
| <span id="877">877</span> |
| <span id="878">878</span> |
| <span id="879">879</span> |
| <span id="880">880</span> |
| <span id="881">881</span> |
| <span id="882">882</span> |
| <span id="883">883</span> |
| <span id="884">884</span> |
| <span id="885">885</span> |
| <span id="886">886</span> |
| <span id="887">887</span> |
| <span id="888">888</span> |
| <span id="889">889</span> |
| <span id="890">890</span> |
| <span id="891">891</span> |
| <span id="892">892</span> |
| <span id="893">893</span> |
| <span id="894">894</span> |
| <span id="895">895</span> |
| <span id="896">896</span> |
| <span id="897">897</span> |
| <span id="898">898</span> |
| <span id="899">899</span> |
| <span id="900">900</span> |
| <span id="901">901</span> |
| <span id="902">902</span> |
| <span id="903">903</span> |
| <span id="904">904</span> |
| <span id="905">905</span> |
| <span id="906">906</span> |
| <span id="907">907</span> |
| <span id="908">908</span> |
| <span id="909">909</span> |
| <span id="910">910</span> |
| <span id="911">911</span> |
| <span id="912">912</span> |
| <span id="913">913</span> |
| <span id="914">914</span> |
| <span id="915">915</span> |
| <span id="916">916</span> |
| <span id="917">917</span> |
| <span id="918">918</span> |
| <span id="919">919</span> |
| <span id="920">920</span> |
| <span id="921">921</span> |
| <span id="922">922</span> |
| <span id="923">923</span> |
| <span id="924">924</span> |
| <span id="925">925</span> |
| <span id="926">926</span> |
| <span id="927">927</span> |
| <span id="928">928</span> |
| <span id="929">929</span> |
| <span id="930">930</span> |
| <span id="931">931</span> |
| <span id="932">932</span> |
| <span id="933">933</span> |
| <span id="934">934</span> |
| <span id="935">935</span> |
| <span id="936">936</span> |
| <span id="937">937</span> |
| <span id="938">938</span> |
| <span id="939">939</span> |
| <span id="940">940</span> |
| <span id="941">941</span> |
| <span id="942">942</span> |
| <span id="943">943</span> |
| <span id="944">944</span> |
| <span id="945">945</span> |
| <span id="946">946</span> |
| <span id="947">947</span> |
| <span id="948">948</span> |
| <span id="949">949</span> |
| <span id="950">950</span> |
| <span id="951">951</span> |
| <span id="952">952</span> |
| <span id="953">953</span> |
| <span id="954">954</span> |
| <span id="955">955</span> |
| <span id="956">956</span> |
| <span id="957">957</span> |
| <span id="958">958</span> |
| <span id="959">959</span> |
| <span id="960">960</span> |
| <span id="961">961</span> |
| <span id="962">962</span> |
| <span id="963">963</span> |
| <span id="964">964</span> |
| <span id="965">965</span> |
| <span id="966">966</span> |
| <span id="967">967</span> |
| <span id="968">968</span> |
| <span id="969">969</span> |
| <span id="970">970</span> |
| <span id="971">971</span> |
| <span id="972">972</span> |
| <span id="973">973</span> |
| <span id="974">974</span> |
| <span id="975">975</span> |
| <span id="976">976</span> |
| <span id="977">977</span> |
| <span id="978">978</span> |
| <span id="979">979</span> |
| <span id="980">980</span> |
| <span id="981">981</span> |
| <span id="982">982</span> |
| <span id="983">983</span> |
| <span id="984">984</span> |
| <span id="985">985</span> |
| <span id="986">986</span> |
| <span id="987">987</span> |
| <span id="988">988</span> |
| <span id="989">989</span> |
| <span id="990">990</span> |
| <span id="991">991</span> |
| <span id="992">992</span> |
| <span id="993">993</span> |
| <span id="994">994</span> |
| <span id="995">995</span> |
| <span id="996">996</span> |
| <span id="997">997</span> |
| <span id="998">998</span> |
| <span id="999">999</span> |
| <span id="1000">1000</span> |
| <span id="1001">1001</span> |
| <span id="1002">1002</span> |
| <span id="1003">1003</span> |
| <span id="1004">1004</span> |
| <span id="1005">1005</span> |
| <span id="1006">1006</span> |
| <span id="1007">1007</span> |
| <span id="1008">1008</span> |
| <span id="1009">1009</span> |
| <span id="1010">1010</span> |
| <span id="1011">1011</span> |
| <span id="1012">1012</span> |
| <span id="1013">1013</span> |
| <span id="1014">1014</span> |
| <span id="1015">1015</span> |
| <span id="1016">1016</span> |
| <span id="1017">1017</span> |
| <span id="1018">1018</span> |
| <span id="1019">1019</span> |
| <span id="1020">1020</span> |
| <span id="1021">1021</span> |
| <span id="1022">1022</span> |
| <span id="1023">1023</span> |
| <span id="1024">1024</span> |
| <span id="1025">1025</span> |
| <span id="1026">1026</span> |
| <span id="1027">1027</span> |
| <span id="1028">1028</span> |
| <span id="1029">1029</span> |
| <span id="1030">1030</span> |
| <span id="1031">1031</span> |
| <span id="1032">1032</span> |
| <span id="1033">1033</span> |
| <span id="1034">1034</span> |
| <span id="1035">1035</span> |
| <span id="1036">1036</span> |
| <span id="1037">1037</span> |
| <span id="1038">1038</span> |
| <span id="1039">1039</span> |
| <span id="1040">1040</span> |
| <span id="1041">1041</span> |
| <span id="1042">1042</span> |
| <span id="1043">1043</span> |
| <span id="1044">1044</span> |
| <span id="1045">1045</span> |
| <span id="1046">1046</span> |
| <span id="1047">1047</span> |
| <span id="1048">1048</span> |
| <span id="1049">1049</span> |
| <span id="1050">1050</span> |
| <span id="1051">1051</span> |
| <span id="1052">1052</span> |
| <span id="1053">1053</span> |
| <span id="1054">1054</span> |
| <span id="1055">1055</span> |
| <span id="1056">1056</span> |
| <span id="1057">1057</span> |
| <span id="1058">1058</span> |
| <span id="1059">1059</span> |
| <span id="1060">1060</span> |
| <span id="1061">1061</span> |
| <span id="1062">1062</span> |
| <span id="1063">1063</span> |
| <span id="1064">1064</span> |
| <span id="1065">1065</span> |
| <span id="1066">1066</span> |
| <span id="1067">1067</span> |
| <span id="1068">1068</span> |
| <span id="1069">1069</span> |
| <span id="1070">1070</span> |
| <span id="1071">1071</span> |
| <span id="1072">1072</span> |
| <span id="1073">1073</span> |
| <span id="1074">1074</span> |
| <span id="1075">1075</span> |
| <span id="1076">1076</span> |
| <span id="1077">1077</span> |
| <span id="1078">1078</span> |
| <span id="1079">1079</span> |
| <span id="1080">1080</span> |
| <span id="1081">1081</span> |
| <span id="1082">1082</span> |
| <span id="1083">1083</span> |
| <span id="1084">1084</span> |
| <span id="1085">1085</span> |
| <span id="1086">1086</span> |
| <span id="1087">1087</span> |
| <span id="1088">1088</span> |
| <span id="1089">1089</span> |
| <span id="1090">1090</span> |
| <span id="1091">1091</span> |
| <span id="1092">1092</span> |
| <span id="1093">1093</span> |
| <span id="1094">1094</span> |
| <span id="1095">1095</span> |
| <span id="1096">1096</span> |
| <span id="1097">1097</span> |
| <span id="1098">1098</span> |
| <span id="1099">1099</span> |
| <span id="1100">1100</span> |
| <span id="1101">1101</span> |
| <span id="1102">1102</span> |
| <span id="1103">1103</span> |
| <span id="1104">1104</span> |
| <span id="1105">1105</span> |
| <span id="1106">1106</span> |
| <span id="1107">1107</span> |
| <span id="1108">1108</span> |
| <span id="1109">1109</span> |
| <span id="1110">1110</span> |
| <span id="1111">1111</span> |
| <span id="1112">1112</span> |
| <span id="1113">1113</span> |
| <span id="1114">1114</span> |
| <span id="1115">1115</span> |
| <span id="1116">1116</span> |
| <span id="1117">1117</span> |
| <span id="1118">1118</span> |
| <span id="1119">1119</span> |
| <span id="1120">1120</span> |
| <span id="1121">1121</span> |
| <span id="1122">1122</span> |
| <span id="1123">1123</span> |
| <span id="1124">1124</span> |
| <span id="1125">1125</span> |
| <span id="1126">1126</span> |
| <span id="1127">1127</span> |
| <span id="1128">1128</span> |
| <span id="1129">1129</span> |
| <span id="1130">1130</span> |
| <span id="1131">1131</span> |
| <span id="1132">1132</span> |
| <span id="1133">1133</span> |
| <span id="1134">1134</span> |
| <span id="1135">1135</span> |
| <span id="1136">1136</span> |
| <span id="1137">1137</span> |
| <span id="1138">1138</span> |
| <span id="1139">1139</span> |
| <span id="1140">1140</span> |
| <span id="1141">1141</span> |
| <span id="1142">1142</span> |
| <span id="1143">1143</span> |
| <span id="1144">1144</span> |
| <span id="1145">1145</span> |
| <span id="1146">1146</span> |
| <span id="1147">1147</span> |
| <span id="1148">1148</span> |
| <span id="1149">1149</span> |
| <span id="1150">1150</span> |
| <span id="1151">1151</span> |
| <span id="1152">1152</span> |
| <span id="1153">1153</span> |
| <span id="1154">1154</span> |
| <span id="1155">1155</span> |
| <span id="1156">1156</span> |
| <span id="1157">1157</span> |
| <span id="1158">1158</span> |
| <span id="1159">1159</span> |
| <span id="1160">1160</span> |
| <span id="1161">1161</span> |
| <span id="1162">1162</span> |
| <span id="1163">1163</span> |
| <span id="1164">1164</span> |
| <span id="1165">1165</span> |
| <span id="1166">1166</span> |
| <span id="1167">1167</span> |
| <span id="1168">1168</span> |
| <span id="1169">1169</span> |
| <span id="1170">1170</span> |
| <span id="1171">1171</span> |
| <span id="1172">1172</span> |
| <span id="1173">1173</span> |
| <span id="1174">1174</span> |
| <span id="1175">1175</span> |
| <span id="1176">1176</span> |
| <span id="1177">1177</span> |
| <span id="1178">1178</span> |
| <span id="1179">1179</span> |
| <span id="1180">1180</span> |
| <span id="1181">1181</span> |
| <span id="1182">1182</span> |
| <span id="1183">1183</span> |
| <span id="1184">1184</span> |
| <span id="1185">1185</span> |
| <span id="1186">1186</span> |
| <span id="1187">1187</span> |
| <span id="1188">1188</span> |
| <span id="1189">1189</span> |
| <span id="1190">1190</span> |
| <span id="1191">1191</span> |
| <span id="1192">1192</span> |
| <span id="1193">1193</span> |
| <span id="1194">1194</span> |
| <span id="1195">1195</span> |
| <span id="1196">1196</span> |
| <span id="1197">1197</span> |
| <span id="1198">1198</span> |
| <span id="1199">1199</span> |
| <span id="1200">1200</span> |
| <span id="1201">1201</span> |
| <span id="1202">1202</span> |
| <span id="1203">1203</span> |
| <span id="1204">1204</span> |
| <span id="1205">1205</span> |
| <span id="1206">1206</span> |
| <span id="1207">1207</span> |
| <span id="1208">1208</span> |
| <span id="1209">1209</span> |
| <span id="1210">1210</span> |
| <span id="1211">1211</span> |
| <span id="1212">1212</span> |
| <span id="1213">1213</span> |
| <span id="1214">1214</span> |
| <span id="1215">1215</span> |
| <span id="1216">1216</span> |
| <span id="1217">1217</span> |
| <span id="1218">1218</span> |
| <span id="1219">1219</span> |
| <span id="1220">1220</span> |
| <span id="1221">1221</span> |
| <span id="1222">1222</span> |
| <span id="1223">1223</span> |
| <span id="1224">1224</span> |
| <span id="1225">1225</span> |
| <span id="1226">1226</span> |
| <span id="1227">1227</span> |
| <span id="1228">1228</span> |
| <span id="1229">1229</span> |
| <span id="1230">1230</span> |
| <span id="1231">1231</span> |
| <span id="1232">1232</span> |
| <span id="1233">1233</span> |
| <span id="1234">1234</span> |
| <span id="1235">1235</span> |
| <span id="1236">1236</span> |
| <span id="1237">1237</span> |
| <span id="1238">1238</span> |
| <span id="1239">1239</span> |
| <span id="1240">1240</span> |
| <span id="1241">1241</span> |
| <span id="1242">1242</span> |
| <span id="1243">1243</span> |
| <span id="1244">1244</span> |
| <span id="1245">1245</span> |
| <span id="1246">1246</span> |
| <span id="1247">1247</span> |
| <span id="1248">1248</span> |
| <span id="1249">1249</span> |
| <span id="1250">1250</span> |
| <span id="1251">1251</span> |
| <span id="1252">1252</span> |
| <span id="1253">1253</span> |
| <span id="1254">1254</span> |
| <span id="1255">1255</span> |
| <span id="1256">1256</span> |
| <span id="1257">1257</span> |
| <span id="1258">1258</span> |
| <span id="1259">1259</span> |
| <span id="1260">1260</span> |
| <span id="1261">1261</span> |
| <span id="1262">1262</span> |
| <span id="1263">1263</span> |
| <span id="1264">1264</span> |
| <span id="1265">1265</span> |
| <span id="1266">1266</span> |
| <span id="1267">1267</span> |
| <span id="1268">1268</span> |
| <span id="1269">1269</span> |
| <span id="1270">1270</span> |
| <span id="1271">1271</span> |
| <span id="1272">1272</span> |
| <span id="1273">1273</span> |
| <span id="1274">1274</span> |
| <span id="1275">1275</span> |
| <span id="1276">1276</span> |
| <span id="1277">1277</span> |
| <span id="1278">1278</span> |
| <span id="1279">1279</span> |
| <span id="1280">1280</span> |
| <span id="1281">1281</span> |
| <span id="1282">1282</span> |
| <span id="1283">1283</span> |
| <span id="1284">1284</span> |
| <span id="1285">1285</span> |
| <span id="1286">1286</span> |
| <span id="1287">1287</span> |
| <span id="1288">1288</span> |
| <span id="1289">1289</span> |
| <span id="1290">1290</span> |
| <span id="1291">1291</span> |
| <span id="1292">1292</span> |
| <span id="1293">1293</span> |
| <span id="1294">1294</span> |
| <span id="1295">1295</span> |
| <span id="1296">1296</span> |
| <span id="1297">1297</span> |
| <span id="1298">1298</span> |
| <span id="1299">1299</span> |
| <span id="1300">1300</span> |
| <span id="1301">1301</span> |
| <span id="1302">1302</span> |
| <span id="1303">1303</span> |
| <span id="1304">1304</span> |
| <span id="1305">1305</span> |
| <span id="1306">1306</span> |
| <span id="1307">1307</span> |
| <span id="1308">1308</span> |
| <span id="1309">1309</span> |
| <span id="1310">1310</span> |
| <span id="1311">1311</span> |
| <span id="1312">1312</span> |
| <span id="1313">1313</span> |
| <span id="1314">1314</span> |
| <span id="1315">1315</span> |
| <span id="1316">1316</span> |
| <span id="1317">1317</span> |
| <span id="1318">1318</span> |
| <span id="1319">1319</span> |
| <span id="1320">1320</span> |
| <span id="1321">1321</span> |
| <span id="1322">1322</span> |
| <span id="1323">1323</span> |
| <span id="1324">1324</span> |
| <span id="1325">1325</span> |
| <span id="1326">1326</span> |
| <span id="1327">1327</span> |
| <span id="1328">1328</span> |
| <span id="1329">1329</span> |
| <span id="1330">1330</span> |
| <span id="1331">1331</span> |
| <span id="1332">1332</span> |
| <span id="1333">1333</span> |
| <span id="1334">1334</span> |
| <span id="1335">1335</span> |
| <span id="1336">1336</span> |
| <span id="1337">1337</span> |
| <span id="1338">1338</span> |
| <span id="1339">1339</span> |
| <span id="1340">1340</span> |
| <span id="1341">1341</span> |
| <span id="1342">1342</span> |
| <span id="1343">1343</span> |
| <span id="1344">1344</span> |
| <span id="1345">1345</span> |
| <span id="1346">1346</span> |
| <span id="1347">1347</span> |
| <span id="1348">1348</span> |
| <span id="1349">1349</span> |
| <span id="1350">1350</span> |
| <span id="1351">1351</span> |
| <span id="1352">1352</span> |
| <span id="1353">1353</span> |
| <span id="1354">1354</span> |
| <span id="1355">1355</span> |
| <span id="1356">1356</span> |
| <span id="1357">1357</span> |
| <span id="1358">1358</span> |
| <span id="1359">1359</span> |
| <span id="1360">1360</span> |
| <span id="1361">1361</span> |
| <span id="1362">1362</span> |
| <span id="1363">1363</span> |
| <span id="1364">1364</span> |
| <span id="1365">1365</span> |
| <span id="1366">1366</span> |
| <span id="1367">1367</span> |
| <span id="1368">1368</span> |
| <span id="1369">1369</span> |
| <span id="1370">1370</span> |
| <span id="1371">1371</span> |
| <span id="1372">1372</span> |
| <span id="1373">1373</span> |
| <span id="1374">1374</span> |
| <span id="1375">1375</span> |
| <span id="1376">1376</span> |
| <span id="1377">1377</span> |
| <span id="1378">1378</span> |
| <span id="1379">1379</span> |
| <span id="1380">1380</span> |
| <span id="1381">1381</span> |
| <span id="1382">1382</span> |
| <span id="1383">1383</span> |
| <span id="1384">1384</span> |
| <span id="1385">1385</span> |
| <span id="1386">1386</span> |
| <span id="1387">1387</span> |
| <span id="1388">1388</span> |
| <span id="1389">1389</span> |
| <span id="1390">1390</span> |
| <span id="1391">1391</span> |
| <span id="1392">1392</span> |
| <span id="1393">1393</span> |
| <span id="1394">1394</span> |
| <span id="1395">1395</span> |
| <span id="1396">1396</span> |
| <span id="1397">1397</span> |
| <span id="1398">1398</span> |
| <span id="1399">1399</span> |
| <span id="1400">1400</span> |
| <span id="1401">1401</span> |
| <span id="1402">1402</span> |
| <span id="1403">1403</span> |
| <span id="1404">1404</span> |
| <span id="1405">1405</span> |
| <span id="1406">1406</span> |
| <span id="1407">1407</span> |
| <span id="1408">1408</span> |
| <span id="1409">1409</span> |
| <span id="1410">1410</span> |
| <span id="1411">1411</span> |
| <span id="1412">1412</span> |
| <span id="1413">1413</span> |
| <span id="1414">1414</span> |
| <span id="1415">1415</span> |
| <span id="1416">1416</span> |
| <span id="1417">1417</span> |
| <span id="1418">1418</span> |
| <span id="1419">1419</span> |
| <span id="1420">1420</span> |
| <span id="1421">1421</span> |
| <span id="1422">1422</span> |
| <span id="1423">1423</span> |
| <span id="1424">1424</span> |
| <span id="1425">1425</span> |
| <span id="1426">1426</span> |
| <span id="1427">1427</span> |
| <span id="1428">1428</span> |
| <span id="1429">1429</span> |
| <span id="1430">1430</span> |
| <span id="1431">1431</span> |
| <span id="1432">1432</span> |
| <span id="1433">1433</span> |
| <span id="1434">1434</span> |
| <span id="1435">1435</span> |
| <span id="1436">1436</span> |
| <span id="1437">1437</span> |
| <span id="1438">1438</span> |
| <span id="1439">1439</span> |
| <span id="1440">1440</span> |
| <span id="1441">1441</span> |
| <span id="1442">1442</span> |
| <span id="1443">1443</span> |
| <span id="1444">1444</span> |
| <span id="1445">1445</span> |
| <span id="1446">1446</span> |
| <span id="1447">1447</span> |
| <span id="1448">1448</span> |
| <span id="1449">1449</span> |
| <span id="1450">1450</span> |
| <span id="1451">1451</span> |
| <span id="1452">1452</span> |
| <span id="1453">1453</span> |
| <span id="1454">1454</span> |
| <span id="1455">1455</span> |
| <span id="1456">1456</span> |
| <span id="1457">1457</span> |
| <span id="1458">1458</span> |
| <span id="1459">1459</span> |
| <span id="1460">1460</span> |
| <span id="1461">1461</span> |
| <span id="1462">1462</span> |
| <span id="1463">1463</span> |
| <span id="1464">1464</span> |
| <span id="1465">1465</span> |
| <span id="1466">1466</span> |
| <span id="1467">1467</span> |
| <span id="1468">1468</span> |
| <span id="1469">1469</span> |
| <span id="1470">1470</span> |
| <span id="1471">1471</span> |
| <span id="1472">1472</span> |
| <span id="1473">1473</span> |
| <span id="1474">1474</span> |
| <span id="1475">1475</span> |
| <span id="1476">1476</span> |
| <span id="1477">1477</span> |
| <span id="1478">1478</span> |
| <span id="1479">1479</span> |
| <span id="1480">1480</span> |
| <span id="1481">1481</span> |
| <span id="1482">1482</span> |
| <span id="1483">1483</span> |
| <span id="1484">1484</span> |
| <span id="1485">1485</span> |
| <span id="1486">1486</span> |
| <span id="1487">1487</span> |
| <span id="1488">1488</span> |
| <span id="1489">1489</span> |
| <span id="1490">1490</span> |
| <span id="1491">1491</span> |
| <span id="1492">1492</span> |
| <span id="1493">1493</span> |
| <span id="1494">1494</span> |
| <span id="1495">1495</span> |
| <span id="1496">1496</span> |
| <span id="1497">1497</span> |
| <span id="1498">1498</span> |
| <span id="1499">1499</span> |
| <span id="1500">1500</span> |
| <span id="1501">1501</span> |
| <span id="1502">1502</span> |
| <span id="1503">1503</span> |
| <span id="1504">1504</span> |
| <span id="1505">1505</span> |
| <span id="1506">1506</span> |
| <span id="1507">1507</span> |
| <span id="1508">1508</span> |
| <span id="1509">1509</span> |
| <span id="1510">1510</span> |
| <span id="1511">1511</span> |
| <span id="1512">1512</span> |
| <span id="1513">1513</span> |
| <span id="1514">1514</span> |
| <span id="1515">1515</span> |
| <span id="1516">1516</span> |
| <span id="1517">1517</span> |
| <span id="1518">1518</span> |
| <span id="1519">1519</span> |
| <span id="1520">1520</span> |
| <span id="1521">1521</span> |
| <span id="1522">1522</span> |
| <span id="1523">1523</span> |
| <span id="1524">1524</span> |
| <span id="1525">1525</span> |
| <span id="1526">1526</span> |
| <span id="1527">1527</span> |
| <span id="1528">1528</span> |
| <span id="1529">1529</span> |
| <span id="1530">1530</span> |
| <span id="1531">1531</span> |
| <span id="1532">1532</span> |
| <span id="1533">1533</span> |
| <span id="1534">1534</span> |
| <span id="1535">1535</span> |
| <span id="1536">1536</span> |
| <span id="1537">1537</span> |
| <span id="1538">1538</span> |
| <span id="1539">1539</span> |
| <span id="1540">1540</span> |
| <span id="1541">1541</span> |
| <span id="1542">1542</span> |
| <span id="1543">1543</span> |
| <span id="1544">1544</span> |
| <span id="1545">1545</span> |
| <span id="1546">1546</span> |
| <span id="1547">1547</span> |
| <span id="1548">1548</span> |
| <span id="1549">1549</span> |
| <span id="1550">1550</span> |
| <span id="1551">1551</span> |
| <span id="1552">1552</span> |
| <span id="1553">1553</span> |
| <span id="1554">1554</span> |
| <span id="1555">1555</span> |
| <span id="1556">1556</span> |
| <span id="1557">1557</span> |
| <span id="1558">1558</span> |
| <span id="1559">1559</span> |
| <span id="1560">1560</span> |
| <span id="1561">1561</span> |
| <span id="1562">1562</span> |
| <span id="1563">1563</span> |
| <span id="1564">1564</span> |
| <span id="1565">1565</span> |
| <span id="1566">1566</span> |
| <span id="1567">1567</span> |
| <span id="1568">1568</span> |
| <span id="1569">1569</span> |
| <span id="1570">1570</span> |
| <span id="1571">1571</span> |
| <span id="1572">1572</span> |
| <span id="1573">1573</span> |
| <span id="1574">1574</span> |
| <span id="1575">1575</span> |
| <span id="1576">1576</span> |
| <span id="1577">1577</span> |
| <span id="1578">1578</span> |
| <span id="1579">1579</span> |
| <span id="1580">1580</span> |
| <span id="1581">1581</span> |
| <span id="1582">1582</span> |
| <span id="1583">1583</span> |
| <span id="1584">1584</span> |
| <span id="1585">1585</span> |
| <span id="1586">1586</span> |
| <span id="1587">1587</span> |
| <span id="1588">1588</span> |
| <span id="1589">1589</span> |
| <span id="1590">1590</span> |
| <span id="1591">1591</span> |
| <span id="1592">1592</span> |
| <span id="1593">1593</span> |
| <span id="1594">1594</span> |
| <span id="1595">1595</span> |
| <span id="1596">1596</span> |
| <span id="1597">1597</span> |
| <span id="1598">1598</span> |
| <span id="1599">1599</span> |
| <span id="1600">1600</span> |
| <span id="1601">1601</span> |
| <span id="1602">1602</span> |
| <span id="1603">1603</span> |
| <span id="1604">1604</span> |
| <span id="1605">1605</span> |
| <span id="1606">1606</span> |
| <span id="1607">1607</span> |
| <span id="1608">1608</span> |
| <span id="1609">1609</span> |
| <span id="1610">1610</span> |
| <span id="1611">1611</span> |
| <span id="1612">1612</span> |
| <span id="1613">1613</span> |
| <span id="1614">1614</span> |
| <span id="1615">1615</span> |
| <span id="1616">1616</span> |
| <span id="1617">1617</span> |
| <span id="1618">1618</span> |
| <span id="1619">1619</span> |
| <span id="1620">1620</span> |
| <span id="1621">1621</span> |
| <span id="1622">1622</span> |
| <span id="1623">1623</span> |
| <span id="1624">1624</span> |
| <span id="1625">1625</span> |
| <span id="1626">1626</span> |
| <span id="1627">1627</span> |
| <span id="1628">1628</span> |
| <span id="1629">1629</span> |
| <span id="1630">1630</span> |
| <span id="1631">1631</span> |
| <span id="1632">1632</span> |
| <span id="1633">1633</span> |
| <span id="1634">1634</span> |
| <span id="1635">1635</span> |
| <span id="1636">1636</span> |
| <span id="1637">1637</span> |
| <span id="1638">1638</span> |
| <span id="1639">1639</span> |
| <span id="1640">1640</span> |
| <span id="1641">1641</span> |
| <span id="1642">1642</span> |
| <span id="1643">1643</span> |
| <span id="1644">1644</span> |
| <span id="1645">1645</span> |
| <span id="1646">1646</span> |
| <span id="1647">1647</span> |
| <span id="1648">1648</span> |
| <span id="1649">1649</span> |
| <span id="1650">1650</span> |
| <span id="1651">1651</span> |
| <span id="1652">1652</span> |
| <span id="1653">1653</span> |
| <span id="1654">1654</span> |
| <span id="1655">1655</span> |
| <span id="1656">1656</span> |
| <span id="1657">1657</span> |
| <span id="1658">1658</span> |
| </pre><pre class="rust"><code><span class="comment">// Copyright 2018-2020 Developers of the Rand project. |
| // Copyright 2017 The Rust Project Developers. |
| // |
| // Licensed under the Apache License, Version 2.0 <LICENSE-APACHE or |
| // https://www.apache.org/licenses/LICENSE-2.0> or the MIT license |
| // <LICENSE-MIT or https://opensource.org/licenses/MIT>, at your |
| // option. This file may not be copied, modified, or distributed |
| // except according to those terms. |
| |
| </span><span class="doccomment">//! A distribution uniformly sampling numbers within a given range. |
| //! |
| //! [`Uniform`] is the standard distribution to sample uniformly from a range; |
| //! e.g. `Uniform::new_inclusive(1, 6)` can sample integers from 1 to 6, like a |
| //! standard die. [`Rng::gen_range`] supports any type supported by |
| //! [`Uniform`]. |
| //! |
| //! This distribution is provided with support for several primitive types |
| //! (all integer and floating-point types) as well as [`std::time::Duration`], |
| //! and supports extension to user-defined types via a type-specific *back-end* |
| //! implementation. |
| //! |
| //! The types [`UniformInt`], [`UniformFloat`] and [`UniformDuration`] are the |
| //! back-ends supporting sampling from primitive integer and floating-point |
| //! ranges as well as from [`std::time::Duration`]; these types do not normally |
| //! need to be used directly (unless implementing a derived back-end). |
| //! |
| //! # Example usage |
| //! |
| //! ``` |
| //! use rand::{Rng, thread_rng}; |
| //! use rand::distributions::Uniform; |
| //! |
| //! let mut rng = thread_rng(); |
| //! let side = Uniform::new(-10.0, 10.0); |
| //! |
| //! // sample between 1 and 10 points |
| //! for _ in 0..rng.gen_range(1..=10) { |
| //! // sample a point from the square with sides -10 - 10 in two dimensions |
| //! let (x, y) = (rng.sample(side), rng.sample(side)); |
| //! println!("Point: {}, {}", x, y); |
| //! } |
| //! ``` |
| //! |
| //! # Extending `Uniform` to support a custom type |
| //! |
| //! To extend [`Uniform`] to support your own types, write a back-end which |
| //! implements the [`UniformSampler`] trait, then implement the [`SampleUniform`] |
| //! helper trait to "register" your back-end. See the `MyF32` example below. |
| //! |
| //! At a minimum, the back-end needs to store any parameters needed for sampling |
| //! (e.g. the target range) and implement `new`, `new_inclusive` and `sample`. |
| //! Those methods should include an assert to check the range is valid (i.e. |
| //! `low < high`). The example below merely wraps another back-end. |
| //! |
| //! The `new`, `new_inclusive` and `sample_single` functions use arguments of |
| //! type SampleBorrow<X> in order to support passing in values by reference or |
| //! by value. In the implementation of these functions, you can choose to |
| //! simply use the reference returned by [`SampleBorrow::borrow`], or you can choose |
| //! to copy or clone the value, whatever is appropriate for your type. |
| //! |
| //! ``` |
| //! use rand::prelude::*; |
| //! use rand::distributions::uniform::{Uniform, SampleUniform, |
| //! UniformSampler, UniformFloat, SampleBorrow}; |
| //! |
| //! struct MyF32(f32); |
| //! |
| //! #[derive(Clone, Copy, Debug)] |
| //! struct UniformMyF32(UniformFloat<f32>); |
| //! |
| //! impl UniformSampler for UniformMyF32 { |
| //! type X = MyF32; |
| //! fn new<B1, B2>(low: B1, high: B2) -> Self |
| //! where B1: SampleBorrow<Self::X> + Sized, |
| //! B2: SampleBorrow<Self::X> + Sized |
| //! { |
| //! UniformMyF32(UniformFloat::<f32>::new(low.borrow().0, high.borrow().0)) |
| //! } |
| //! fn new_inclusive<B1, B2>(low: B1, high: B2) -> Self |
| //! where B1: SampleBorrow<Self::X> + Sized, |
| //! B2: SampleBorrow<Self::X> + Sized |
| //! { |
| //! UniformMyF32(UniformFloat::<f32>::new_inclusive( |
| //! low.borrow().0, |
| //! high.borrow().0, |
| //! )) |
| //! } |
| //! fn sample<R: Rng + ?Sized>(&self, rng: &mut R) -> Self::X { |
| //! MyF32(self.0.sample(rng)) |
| //! } |
| //! } |
| //! |
| //! impl SampleUniform for MyF32 { |
| //! type Sampler = UniformMyF32; |
| //! } |
| //! |
| //! let (low, high) = (MyF32(17.0f32), MyF32(22.0f32)); |
| //! let uniform = Uniform::new(low, high); |
| //! let x = uniform.sample(&mut thread_rng()); |
| //! ``` |
| //! |
| //! [`SampleUniform`]: crate::distributions::uniform::SampleUniform |
| //! [`UniformSampler`]: crate::distributions::uniform::UniformSampler |
| //! [`UniformInt`]: crate::distributions::uniform::UniformInt |
| //! [`UniformFloat`]: crate::distributions::uniform::UniformFloat |
| //! [`UniformDuration`]: crate::distributions::uniform::UniformDuration |
| //! [`SampleBorrow::borrow`]: crate::distributions::uniform::SampleBorrow::borrow |
| |
| </span><span class="kw">use </span>core::time::Duration; |
| <span class="kw">use </span>core::ops::{Range, RangeInclusive}; |
| |
| <span class="kw">use </span><span class="kw">crate</span>::distributions::float::IntoFloat; |
| <span class="kw">use </span><span class="kw">crate</span>::distributions::utils::{BoolAsSIMD, FloatAsSIMD, FloatSIMDUtils, WideningMultiply}; |
| <span class="kw">use </span><span class="kw">crate</span>::distributions::Distribution; |
| <span class="kw">use crate</span>::{Rng, RngCore}; |
| |
| <span class="attribute">#[cfg(not(feature = <span class="string">"std"</span>))] |
| #[allow(unused_imports)] </span><span class="comment">// rustc doesn't detect that this is actually used |
| </span><span class="kw">use </span><span class="kw">crate</span>::distributions::utils::Float; |
| |
| <span class="attribute">#[cfg(feature = <span class="string">"simd_support"</span>)] </span><span class="kw">use </span>packed_simd::<span class="kw-2">*</span>; |
| |
| <span class="attribute">#[cfg(feature = <span class="string">"serde1"</span>)] |
| </span><span class="kw">use </span>serde::{Serialize, Deserialize}; |
| |
| <span class="doccomment">/// Sample values uniformly between two bounds. |
| /// |
| /// [`Uniform::new`] and [`Uniform::new_inclusive`] construct a uniform |
| /// distribution sampling from the given range; these functions may do extra |
| /// work up front to make sampling of multiple values faster. If only one sample |
| /// from the range is required, [`Rng::gen_range`] can be more efficient. |
| /// |
| /// When sampling from a constant range, many calculations can happen at |
| /// compile-time and all methods should be fast; for floating-point ranges and |
| /// the full range of integer types this should have comparable performance to |
| /// the `Standard` distribution. |
| /// |
| /// Steps are taken to avoid bias which might be present in naive |
| /// implementations; for example `rng.gen::<u8>() % 170` samples from the range |
| /// `[0, 169]` but is twice as likely to select numbers less than 85 than other |
| /// values. Further, the implementations here give more weight to the high-bits |
| /// generated by the RNG than the low bits, since with some RNGs the low-bits |
| /// are of lower quality than the high bits. |
| /// |
| /// Implementations must sample in `[low, high)` range for |
| /// `Uniform::new(low, high)`, i.e., excluding `high`. In particular, care must |
| /// be taken to ensure that rounding never results values `< low` or `>= high`. |
| /// |
| /// # Example |
| /// |
| /// ``` |
| /// use rand::distributions::{Distribution, Uniform}; |
| /// |
| /// let between = Uniform::from(10..10000); |
| /// let mut rng = rand::thread_rng(); |
| /// let mut sum = 0; |
| /// for _ in 0..1000 { |
| /// sum += between.sample(&mut rng); |
| /// } |
| /// println!("{}", sum); |
| /// ``` |
| /// |
| /// For a single sample, [`Rng::gen_range`] may be preferred: |
| /// |
| /// ``` |
| /// use rand::Rng; |
| /// |
| /// let mut rng = rand::thread_rng(); |
| /// println!("{}", rng.gen_range(0..10)); |
| /// ``` |
| /// |
| /// [`new`]: Uniform::new |
| /// [`new_inclusive`]: Uniform::new_inclusive |
| /// [`Rng::gen_range`]: Rng::gen_range |
| </span><span class="attribute">#[derive(Clone, Copy, Debug, PartialEq)] |
| #[cfg_attr(feature = <span class="string">"serde1"</span>, derive(Serialize, Deserialize))] |
| #[cfg_attr(feature = <span class="string">"serde1"</span>, serde(bound(serialize = <span class="string">"X::Sampler: Serialize"</span>)))] |
| #[cfg_attr(feature = <span class="string">"serde1"</span>, serde(bound(deserialize = <span class="string">"X::Sampler: Deserialize<'de>"</span>)))] |
| </span><span class="kw">pub struct </span>Uniform<X: SampleUniform>(X::Sampler); |
| |
| <span class="kw">impl</span><X: SampleUniform> Uniform<X> { |
| <span class="doccomment">/// Create a new `Uniform` instance which samples uniformly from the half |
| /// open range `[low, high)` (excluding `high`). Panics if `low >= high`. |
| </span><span class="kw">pub fn </span>new<B1, B2>(low: B1, high: B2) -> Uniform<X> |
| <span class="kw">where |
| </span>B1: SampleBorrow<X> + Sized, |
| B2: SampleBorrow<X> + Sized, |
| { |
| Uniform(X::Sampler::new(low, high)) |
| } |
| |
| <span class="doccomment">/// Create a new `Uniform` instance which samples uniformly from the closed |
| /// range `[low, high]` (inclusive). Panics if `low > high`. |
| </span><span class="kw">pub fn </span>new_inclusive<B1, B2>(low: B1, high: B2) -> Uniform<X> |
| <span class="kw">where |
| </span>B1: SampleBorrow<X> + Sized, |
| B2: SampleBorrow<X> + Sized, |
| { |
| Uniform(X::Sampler::new_inclusive(low, high)) |
| } |
| } |
| |
| <span class="kw">impl</span><X: SampleUniform> Distribution<X> <span class="kw">for </span>Uniform<X> { |
| <span class="kw">fn </span>sample<R: Rng + <span class="question-mark">?</span>Sized>(<span class="kw-2">&</span><span class="self">self</span>, rng: <span class="kw-2">&mut </span>R) -> X { |
| <span class="self">self</span>.<span class="number">0</span>.sample(rng) |
| } |
| } |
| |
| <span class="doccomment">/// Helper trait for creating objects using the correct implementation of |
| /// [`UniformSampler`] for the sampling type. |
| /// |
| /// See the [module documentation] on how to implement [`Uniform`] range |
| /// sampling for a custom type. |
| /// |
| /// [module documentation]: crate::distributions::uniform |
| </span><span class="kw">pub trait </span>SampleUniform: Sized { |
| <span class="doccomment">/// The `UniformSampler` implementation supporting type `X`. |
| </span><span class="kw">type </span>Sampler: UniformSampler<X = <span class="self">Self</span>>; |
| } |
| |
| <span class="doccomment">/// Helper trait handling actual uniform sampling. |
| /// |
| /// See the [module documentation] on how to implement [`Uniform`] range |
| /// sampling for a custom type. |
| /// |
| /// Implementation of [`sample_single`] is optional, and is only useful when |
| /// the implementation can be faster than `Self::new(low, high).sample(rng)`. |
| /// |
| /// [module documentation]: crate::distributions::uniform |
| /// [`sample_single`]: UniformSampler::sample_single |
| </span><span class="kw">pub trait </span>UniformSampler: Sized { |
| <span class="doccomment">/// The type sampled by this implementation. |
| </span><span class="kw">type </span>X; |
| |
| <span class="doccomment">/// Construct self, with inclusive lower bound and exclusive upper bound |
| /// `[low, high)`. |
| /// |
| /// Usually users should not call this directly but instead use |
| /// `Uniform::new`, which asserts that `low < high` before calling this. |
| </span><span class="kw">fn </span>new<B1, B2>(low: B1, high: B2) -> <span class="self">Self |
| </span><span class="kw">where |
| </span>B1: SampleBorrow<<span class="self">Self</span>::X> + Sized, |
| B2: SampleBorrow<<span class="self">Self</span>::X> + Sized; |
| |
| <span class="doccomment">/// Construct self, with inclusive bounds `[low, high]`. |
| /// |
| /// Usually users should not call this directly but instead use |
| /// `Uniform::new_inclusive`, which asserts that `low <= high` before |
| /// calling this. |
| </span><span class="kw">fn </span>new_inclusive<B1, B2>(low: B1, high: B2) -> <span class="self">Self |
| </span><span class="kw">where |
| </span>B1: SampleBorrow<<span class="self">Self</span>::X> + Sized, |
| B2: SampleBorrow<<span class="self">Self</span>::X> + Sized; |
| |
| <span class="doccomment">/// Sample a value. |
| </span><span class="kw">fn </span>sample<R: Rng + <span class="question-mark">?</span>Sized>(<span class="kw-2">&</span><span class="self">self</span>, rng: <span class="kw-2">&mut </span>R) -> <span class="self">Self</span>::X; |
| |
| <span class="doccomment">/// Sample a single value uniformly from a range with inclusive lower bound |
| /// and exclusive upper bound `[low, high)`. |
| /// |
| /// By default this is implemented using |
| /// `UniformSampler::new(low, high).sample(rng)`. However, for some types |
| /// more optimal implementations for single usage may be provided via this |
| /// method (which is the case for integers and floats). |
| /// Results may not be identical. |
| /// |
| /// Note that to use this method in a generic context, the type needs to be |
| /// retrieved via `SampleUniform::Sampler` as follows: |
| /// ``` |
| /// use rand::{thread_rng, distributions::uniform::{SampleUniform, UniformSampler}}; |
| /// # #[allow(unused)] |
| /// fn sample_from_range<T: SampleUniform>(lb: T, ub: T) -> T { |
| /// let mut rng = thread_rng(); |
| /// <T as SampleUniform>::Sampler::sample_single(lb, ub, &mut rng) |
| /// } |
| /// ``` |
| </span><span class="kw">fn </span>sample_single<R: Rng + <span class="question-mark">?</span>Sized, B1, B2>(low: B1, high: B2, rng: <span class="kw-2">&mut </span>R) -> <span class="self">Self</span>::X |
| <span class="kw">where |
| </span>B1: SampleBorrow<<span class="self">Self</span>::X> + Sized, |
| B2: SampleBorrow<<span class="self">Self</span>::X> + Sized, |
| { |
| <span class="kw">let </span>uniform: <span class="self">Self </span>= UniformSampler::new(low, high); |
| uniform.sample(rng) |
| } |
| |
| <span class="doccomment">/// Sample a single value uniformly from a range with inclusive lower bound |
| /// and inclusive upper bound `[low, high]`. |
| /// |
| /// By default this is implemented using |
| /// `UniformSampler::new_inclusive(low, high).sample(rng)`. However, for |
| /// some types more optimal implementations for single usage may be provided |
| /// via this method. |
| /// Results may not be identical. |
| </span><span class="kw">fn </span>sample_single_inclusive<R: Rng + <span class="question-mark">?</span>Sized, B1, B2>(low: B1, high: B2, rng: <span class="kw-2">&mut </span>R) |
| -> <span class="self">Self</span>::X |
| <span class="kw">where </span>B1: SampleBorrow<<span class="self">Self</span>::X> + Sized, |
| B2: SampleBorrow<<span class="self">Self</span>::X> + Sized |
| { |
| <span class="kw">let </span>uniform: <span class="self">Self </span>= UniformSampler::new_inclusive(low, high); |
| uniform.sample(rng) |
| } |
| } |
| |
| <span class="kw">impl</span><X: SampleUniform> From<Range<X>> <span class="kw">for </span>Uniform<X> { |
| <span class="kw">fn </span>from(r: ::core::ops::Range<X>) -> Uniform<X> { |
| Uniform::new(r.start, r.end) |
| } |
| } |
| |
| <span class="kw">impl</span><X: SampleUniform> From<RangeInclusive<X>> <span class="kw">for </span>Uniform<X> { |
| <span class="kw">fn </span>from(r: ::core::ops::RangeInclusive<X>) -> Uniform<X> { |
| Uniform::new_inclusive(r.start(), r.end()) |
| } |
| } |
| |
| |
| <span class="doccomment">/// Helper trait similar to [`Borrow`] but implemented |
| /// only for SampleUniform and references to SampleUniform in |
| /// order to resolve ambiguity issues. |
| /// |
| /// [`Borrow`]: std::borrow::Borrow |
| </span><span class="kw">pub trait </span>SampleBorrow<Borrowed> { |
| <span class="doccomment">/// Immutably borrows from an owned value. See [`Borrow::borrow`] |
| /// |
| /// [`Borrow::borrow`]: std::borrow::Borrow::borrow |
| </span><span class="kw">fn </span>borrow(<span class="kw-2">&</span><span class="self">self</span>) -> <span class="kw-2">&</span>Borrowed; |
| } |
| <span class="kw">impl</span><Borrowed> SampleBorrow<Borrowed> <span class="kw">for </span>Borrowed |
| <span class="kw">where </span>Borrowed: SampleUniform |
| { |
| <span class="attribute">#[inline(always)] |
| </span><span class="kw">fn </span>borrow(<span class="kw-2">&</span><span class="self">self</span>) -> <span class="kw-2">&</span>Borrowed { |
| <span class="self">self |
| </span>} |
| } |
| <span class="kw">impl</span><<span class="lifetime">'a</span>, Borrowed> SampleBorrow<Borrowed> <span class="kw">for </span><span class="kw-2">&</span><span class="lifetime">'a </span>Borrowed |
| <span class="kw">where </span>Borrowed: SampleUniform |
| { |
| <span class="attribute">#[inline(always)] |
| </span><span class="kw">fn </span>borrow(<span class="kw-2">&</span><span class="self">self</span>) -> <span class="kw-2">&</span>Borrowed { |
| <span class="kw-2">*</span><span class="self">self |
| </span>} |
| } |
| |
| <span class="doccomment">/// Range that supports generating a single sample efficiently. |
| /// |
| /// Any type implementing this trait can be used to specify the sampled range |
| /// for `Rng::gen_range`. |
| </span><span class="kw">pub trait </span>SampleRange<T> { |
| <span class="doccomment">/// Generate a sample from the given range. |
| </span><span class="kw">fn </span>sample_single<R: RngCore + <span class="question-mark">?</span>Sized>(<span class="self">self</span>, rng: <span class="kw-2">&mut </span>R) -> T; |
| |
| <span class="doccomment">/// Check whether the range is empty. |
| </span><span class="kw">fn </span>is_empty(<span class="kw-2">&</span><span class="self">self</span>) -> bool; |
| } |
| |
| <span class="kw">impl</span><T: SampleUniform + PartialOrd> SampleRange<T> <span class="kw">for </span>Range<T> { |
| <span class="attribute">#[inline] |
| </span><span class="kw">fn </span>sample_single<R: RngCore + <span class="question-mark">?</span>Sized>(<span class="self">self</span>, rng: <span class="kw-2">&mut </span>R) -> T { |
| T::Sampler::sample_single(<span class="self">self</span>.start, <span class="self">self</span>.end, rng) |
| } |
| |
| <span class="attribute">#[inline] |
| </span><span class="kw">fn </span>is_empty(<span class="kw-2">&</span><span class="self">self</span>) -> bool { |
| !(<span class="self">self</span>.start < <span class="self">self</span>.end) |
| } |
| } |
| |
| <span class="kw">impl</span><T: SampleUniform + PartialOrd> SampleRange<T> <span class="kw">for </span>RangeInclusive<T> { |
| <span class="attribute">#[inline] |
| </span><span class="kw">fn </span>sample_single<R: RngCore + <span class="question-mark">?</span>Sized>(<span class="self">self</span>, rng: <span class="kw-2">&mut </span>R) -> T { |
| T::Sampler::sample_single_inclusive(<span class="self">self</span>.start(), <span class="self">self</span>.end(), rng) |
| } |
| |
| <span class="attribute">#[inline] |
| </span><span class="kw">fn </span>is_empty(<span class="kw-2">&</span><span class="self">self</span>) -> bool { |
| !(<span class="self">self</span>.start() <= <span class="self">self</span>.end()) |
| } |
| } |
| |
| |
| <span class="comment">//////////////////////////////////////////////////////////////////////////////// |
| |
| // What follows are all back-ends. |
| |
| |
| </span><span class="doccomment">/// The back-end implementing [`UniformSampler`] for integer types. |
| /// |
| /// Unless you are implementing [`UniformSampler`] for your own type, this type |
| /// should not be used directly, use [`Uniform`] instead. |
| /// |
| /// # Implementation notes |
| /// |
| /// For simplicity, we use the same generic struct `UniformInt<X>` for all |
| /// integer types `X`. This gives us only one field type, `X`; to store unsigned |
| /// values of this size, we take use the fact that these conversions are no-ops. |
| /// |
| /// For a closed range, the number of possible numbers we should generate is |
| /// `range = (high - low + 1)`. To avoid bias, we must ensure that the size of |
| /// our sample space, `zone`, is a multiple of `range`; other values must be |
| /// rejected (by replacing with a new random sample). |
| /// |
| /// As a special case, we use `range = 0` to represent the full range of the |
| /// result type (i.e. for `new_inclusive($ty::MIN, $ty::MAX)`). |
| /// |
| /// The optimum `zone` is the largest product of `range` which fits in our |
| /// (unsigned) target type. We calculate this by calculating how many numbers we |
| /// must reject: `reject = (MAX + 1) % range = (MAX - range + 1) % range`. Any (large) |
| /// product of `range` will suffice, thus in `sample_single` we multiply by a |
| /// power of 2 via bit-shifting (faster but may cause more rejections). |
| /// |
| /// The smallest integer PRNGs generate is `u32`. For 8- and 16-bit outputs we |
| /// use `u32` for our `zone` and samples (because it's not slower and because |
| /// it reduces the chance of having to reject a sample). In this case we cannot |
| /// store `zone` in the target type since it is too large, however we know |
| /// `ints_to_reject < range <= $unsigned::MAX`. |
| /// |
| /// An alternative to using a modulus is widening multiply: After a widening |
| /// multiply by `range`, the result is in the high word. Then comparing the low |
| /// word against `zone` makes sure our distribution is uniform. |
| </span><span class="attribute">#[derive(Clone, Copy, Debug, PartialEq)] |
| #[cfg_attr(feature = <span class="string">"serde1"</span>, derive(Serialize, Deserialize))] |
| </span><span class="kw">pub struct </span>UniformInt<X> { |
| low: X, |
| range: X, |
| z: X, <span class="comment">// either ints_to_reject or zone depending on implementation |
| </span>} |
| |
| <span class="macro">macro_rules! </span>uniform_int_impl { |
| (<span class="macro-nonterminal">$ty</span>:ty, <span class="macro-nonterminal">$unsigned</span>:ident, <span class="macro-nonterminal">$u_large</span>:ident) => { |
| <span class="kw">impl </span>SampleUniform <span class="kw">for </span><span class="macro-nonterminal">$ty </span>{ |
| <span class="kw">type </span>Sampler = UniformInt<<span class="macro-nonterminal">$ty</span>>; |
| } |
| |
| <span class="kw">impl </span>UniformSampler <span class="kw">for </span>UniformInt<<span class="macro-nonterminal">$ty</span>> { |
| <span class="comment">// We play free and fast with unsigned vs signed here |
| // (when $ty is signed), but that's fine, since the |
| // contract of this macro is for $ty and $unsigned to be |
| // "bit-equal", so casting between them is a no-op. |
| |
| </span><span class="kw">type </span>X = <span class="macro-nonterminal">$ty</span>; |
| |
| <span class="attribute">#[inline] </span><span class="comment">// if the range is constant, this helps LLVM to do the |
| // calculations at compile-time. |
| </span><span class="kw">fn </span>new<B1, B2>(low_b: B1, high_b: B2) -> <span class="self">Self |
| </span><span class="kw">where |
| </span>B1: SampleBorrow<<span class="self">Self</span>::X> + Sized, |
| B2: SampleBorrow<<span class="self">Self</span>::X> + Sized, |
| { |
| <span class="kw">let </span>low = <span class="kw-2">*</span>low_b.borrow(); |
| <span class="kw">let </span>high = <span class="kw-2">*</span>high_b.borrow(); |
| <span class="macro">assert!</span>(low < high, <span class="string">"Uniform::new called with `low >= high`"</span>); |
| UniformSampler::new_inclusive(low, high - <span class="number">1</span>) |
| } |
| |
| <span class="attribute">#[inline] </span><span class="comment">// if the range is constant, this helps LLVM to do the |
| // calculations at compile-time. |
| </span><span class="kw">fn </span>new_inclusive<B1, B2>(low_b: B1, high_b: B2) -> <span class="self">Self |
| </span><span class="kw">where |
| </span>B1: SampleBorrow<<span class="self">Self</span>::X> + Sized, |
| B2: SampleBorrow<<span class="self">Self</span>::X> + Sized, |
| { |
| <span class="kw">let </span>low = <span class="kw-2">*</span>low_b.borrow(); |
| <span class="kw">let </span>high = <span class="kw-2">*</span>high_b.borrow(); |
| <span class="macro">assert!</span>( |
| low <= high, |
| <span class="string">"Uniform::new_inclusive called with `low > high`" |
| </span>); |
| <span class="kw">let </span>unsigned_max = ::core::<span class="macro-nonterminal">$u_large::MAX</span>; |
| |
| <span class="kw">let </span>range = high.wrapping_sub(low).wrapping_add(<span class="number">1</span>) <span class="kw">as </span><span class="macro-nonterminal">$unsigned</span>; |
| <span class="kw">let </span>ints_to_reject = <span class="kw">if </span>range > <span class="number">0 </span>{ |
| <span class="kw">let </span>range = <span class="macro-nonterminal">$u_large::from</span>(range); |
| (unsigned_max - range + <span class="number">1</span>) % range |
| } <span class="kw">else </span>{ |
| <span class="number">0 |
| </span>}; |
| |
| UniformInt { |
| low, |
| <span class="comment">// These are really $unsigned values, but store as $ty: |
| </span>range: range <span class="kw">as </span><span class="macro-nonterminal">$ty</span>, |
| z: ints_to_reject <span class="kw">as </span><span class="macro-nonterminal">$unsigned </span><span class="kw">as </span><span class="macro-nonterminal">$ty</span>, |
| } |
| } |
| |
| <span class="attribute">#[inline] |
| </span><span class="kw">fn </span>sample<R: Rng + <span class="question-mark">?</span>Sized>(<span class="kw-2">&</span><span class="self">self</span>, rng: <span class="kw-2">&mut </span>R) -> <span class="self">Self</span>::X { |
| <span class="kw">let </span>range = <span class="self">self</span>.range <span class="kw">as </span><span class="macro-nonterminal">$unsigned </span><span class="kw">as </span><span class="macro-nonterminal">$u_large</span>; |
| <span class="kw">if </span>range > <span class="number">0 </span>{ |
| <span class="kw">let </span>unsigned_max = ::core::<span class="macro-nonterminal">$u_large::MAX</span>; |
| <span class="kw">let </span>zone = unsigned_max - (<span class="self">self</span>.z <span class="kw">as </span><span class="macro-nonterminal">$unsigned </span><span class="kw">as </span><span class="macro-nonterminal">$u_large</span>); |
| <span class="kw">loop </span>{ |
| <span class="kw">let </span>v: <span class="macro-nonterminal">$u_large </span>= rng.gen(); |
| <span class="kw">let </span>(hi, lo) = v.wmul(range); |
| <span class="kw">if </span>lo <= zone { |
| <span class="kw">return </span><span class="self">self</span>.low.wrapping_add(hi <span class="kw">as </span><span class="macro-nonterminal">$ty</span>); |
| } |
| } |
| } <span class="kw">else </span>{ |
| <span class="comment">// Sample from the entire integer range. |
| </span>rng.gen() |
| } |
| } |
| |
| <span class="attribute">#[inline] |
| </span><span class="kw">fn </span>sample_single<R: Rng + <span class="question-mark">?</span>Sized, B1, B2>(low_b: B1, high_b: B2, rng: <span class="kw-2">&mut </span>R) -> <span class="self">Self</span>::X |
| <span class="kw">where |
| </span>B1: SampleBorrow<<span class="self">Self</span>::X> + Sized, |
| B2: SampleBorrow<<span class="self">Self</span>::X> + Sized, |
| { |
| <span class="kw">let </span>low = <span class="kw-2">*</span>low_b.borrow(); |
| <span class="kw">let </span>high = <span class="kw-2">*</span>high_b.borrow(); |
| <span class="macro">assert!</span>(low < high, <span class="string">"UniformSampler::sample_single: low >= high"</span>); |
| <span class="self">Self</span>::sample_single_inclusive(low, high - <span class="number">1</span>, rng) |
| } |
| |
| <span class="attribute">#[inline] |
| </span><span class="kw">fn </span>sample_single_inclusive<R: Rng + <span class="question-mark">?</span>Sized, B1, B2>(low_b: B1, high_b: B2, rng: <span class="kw-2">&mut </span>R) -> <span class="self">Self</span>::X |
| <span class="kw">where |
| </span>B1: SampleBorrow<<span class="self">Self</span>::X> + Sized, |
| B2: SampleBorrow<<span class="self">Self</span>::X> + Sized, |
| { |
| <span class="kw">let </span>low = <span class="kw-2">*</span>low_b.borrow(); |
| <span class="kw">let </span>high = <span class="kw-2">*</span>high_b.borrow(); |
| <span class="macro">assert!</span>(low <= high, <span class="string">"UniformSampler::sample_single_inclusive: low > high"</span>); |
| <span class="kw">let </span>range = high.wrapping_sub(low).wrapping_add(<span class="number">1</span>) <span class="kw">as </span><span class="macro-nonterminal">$unsigned </span><span class="kw">as </span><span class="macro-nonterminal">$u_large</span>; |
| <span class="comment">// If the above resulted in wrap-around to 0, the range is $ty::MIN..=$ty::MAX, |
| // and any integer will do. |
| </span><span class="kw">if </span>range == <span class="number">0 </span>{ |
| <span class="kw">return </span>rng.gen(); |
| } |
| |
| <span class="kw">let </span>zone = <span class="kw">if </span>::core::<span class="macro-nonterminal">$unsigned::MAX </span><= ::core::u16::MAX <span class="kw">as </span><span class="macro-nonterminal">$unsigned </span>{ |
| <span class="comment">// Using a modulus is faster than the approximation for |
| // i8 and i16. I suppose we trade the cost of one |
| // modulus for near-perfect branch prediction. |
| </span><span class="kw">let </span>unsigned_max: <span class="macro-nonterminal">$u_large </span>= ::core::<span class="macro-nonterminal">$u_large::MAX</span>; |
| <span class="kw">let </span>ints_to_reject = (unsigned_max - range + <span class="number">1</span>) % range; |
| unsigned_max - ints_to_reject |
| } <span class="kw">else </span>{ |
| <span class="comment">// conservative but fast approximation. `- 1` is necessary to allow the |
| // same comparison without bias. |
| </span>(range << range.leading_zeros()).wrapping_sub(<span class="number">1</span>) |
| }; |
| |
| <span class="kw">loop </span>{ |
| <span class="kw">let </span>v: <span class="macro-nonterminal">$u_large </span>= rng.gen(); |
| <span class="kw">let </span>(hi, lo) = v.wmul(range); |
| <span class="kw">if </span>lo <= zone { |
| <span class="kw">return </span>low.wrapping_add(hi <span class="kw">as </span><span class="macro-nonterminal">$ty</span>); |
| } |
| } |
| } |
| } |
| }; |
| } |
| |
| <span class="macro">uniform_int_impl! </span>{ i8, u8, u32 } |
| <span class="macro">uniform_int_impl! </span>{ i16, u16, u32 } |
| <span class="macro">uniform_int_impl! </span>{ i32, u32, u32 } |
| <span class="macro">uniform_int_impl! </span>{ i64, u64, u64 } |
| <span class="macro">uniform_int_impl! </span>{ i128, u128, u128 } |
| <span class="macro">uniform_int_impl! </span>{ isize, usize, usize } |
| <span class="macro">uniform_int_impl! </span>{ u8, u8, u32 } |
| <span class="macro">uniform_int_impl! </span>{ u16, u16, u32 } |
| <span class="macro">uniform_int_impl! </span>{ u32, u32, u32 } |
| <span class="macro">uniform_int_impl! </span>{ u64, u64, u64 } |
| <span class="macro">uniform_int_impl! </span>{ usize, usize, usize } |
| <span class="macro">uniform_int_impl! </span>{ u128, u128, u128 } |
| |
| <span class="attribute">#[cfg(feature = <span class="string">"simd_support"</span>)] |
| </span><span class="macro">macro_rules! </span>uniform_simd_int_impl { |
| (<span class="macro-nonterminal">$ty</span>:ident, <span class="macro-nonterminal">$unsigned</span>:ident, <span class="macro-nonterminal">$u_scalar</span>:ident) => { |
| <span class="comment">// The "pick the largest zone that can fit in an `u32`" optimization |
| // is less useful here. Multiple lanes complicate things, we don't |
| // know the PRNG's minimal output size, and casting to a larger vector |
| // is generally a bad idea for SIMD performance. The user can still |
| // implement it manually. |
| |
| // TODO: look into `Uniform::<u32x4>::new(0u32, 100)` functionality |
| // perhaps `impl SampleUniform for $u_scalar`? |
| </span><span class="kw">impl </span>SampleUniform <span class="kw">for </span><span class="macro-nonterminal">$ty </span>{ |
| <span class="kw">type </span>Sampler = UniformInt<<span class="macro-nonterminal">$ty</span>>; |
| } |
| |
| <span class="kw">impl </span>UniformSampler <span class="kw">for </span>UniformInt<<span class="macro-nonterminal">$ty</span>> { |
| <span class="kw">type </span>X = <span class="macro-nonterminal">$ty</span>; |
| |
| <span class="attribute">#[inline] </span><span class="comment">// if the range is constant, this helps LLVM to do the |
| // calculations at compile-time. |
| </span><span class="kw">fn </span>new<B1, B2>(low_b: B1, high_b: B2) -> <span class="self">Self |
| </span><span class="kw">where </span>B1: SampleBorrow<<span class="self">Self</span>::X> + Sized, |
| B2: SampleBorrow<<span class="self">Self</span>::X> + Sized |
| { |
| <span class="kw">let </span>low = <span class="kw-2">*</span>low_b.borrow(); |
| <span class="kw">let </span>high = <span class="kw-2">*</span>high_b.borrow(); |
| <span class="macro">assert!</span>(low.lt(high).all(), <span class="string">"Uniform::new called with `low >= high`"</span>); |
| UniformSampler::new_inclusive(low, high - <span class="number">1</span>) |
| } |
| |
| <span class="attribute">#[inline] </span><span class="comment">// if the range is constant, this helps LLVM to do the |
| // calculations at compile-time. |
| </span><span class="kw">fn </span>new_inclusive<B1, B2>(low_b: B1, high_b: B2) -> <span class="self">Self |
| </span><span class="kw">where </span>B1: SampleBorrow<<span class="self">Self</span>::X> + Sized, |
| B2: SampleBorrow<<span class="self">Self</span>::X> + Sized |
| { |
| <span class="kw">let </span>low = <span class="kw-2">*</span>low_b.borrow(); |
| <span class="kw">let </span>high = <span class="kw-2">*</span>high_b.borrow(); |
| <span class="macro">assert!</span>(low.le(high).all(), |
| <span class="string">"Uniform::new_inclusive called with `low > high`"</span>); |
| <span class="kw">let </span>unsigned_max = ::core::<span class="macro-nonterminal">$u_scalar::MAX</span>; |
| |
| <span class="comment">// NOTE: these may need to be replaced with explicitly |
| // wrapping operations if `packed_simd` changes |
| </span><span class="kw">let </span>range: <span class="macro-nonterminal">$unsigned </span>= ((high - low) + <span class="number">1</span>).cast(); |
| <span class="comment">// `% 0` will panic at runtime. |
| </span><span class="kw">let </span>not_full_range = range.gt(<span class="macro-nonterminal">$unsigned::splat</span>(<span class="number">0</span>)); |
| <span class="comment">// replacing 0 with `unsigned_max` allows a faster `select` |
| // with bitwise OR |
| </span><span class="kw">let </span>modulo = not_full_range.select(range, <span class="macro-nonterminal">$unsigned::splat</span>(unsigned_max)); |
| <span class="comment">// wrapping addition |
| </span><span class="kw">let </span>ints_to_reject = (unsigned_max - range + <span class="number">1</span>) % modulo; |
| <span class="comment">// When `range` is 0, `lo` of `v.wmul(range)` will always be |
| // zero which means only one sample is needed. |
| </span><span class="kw">let </span>zone = unsigned_max - ints_to_reject; |
| |
| UniformInt { |
| low, |
| <span class="comment">// These are really $unsigned values, but store as $ty: |
| </span>range: range.cast(), |
| z: zone.cast(), |
| } |
| } |
| |
| <span class="kw">fn </span>sample<R: Rng + <span class="question-mark">?</span>Sized>(<span class="kw-2">&</span><span class="self">self</span>, rng: <span class="kw-2">&mut </span>R) -> <span class="self">Self</span>::X { |
| <span class="kw">let </span>range: <span class="macro-nonterminal">$unsigned </span>= <span class="self">self</span>.range.cast(); |
| <span class="kw">let </span>zone: <span class="macro-nonterminal">$unsigned </span>= <span class="self">self</span>.z.cast(); |
| |
| <span class="comment">// This might seem very slow, generating a whole new |
| // SIMD vector for every sample rejection. For most uses |
| // though, the chance of rejection is small and provides good |
| // general performance. With multiple lanes, that chance is |
| // multiplied. To mitigate this, we replace only the lanes of |
| // the vector which fail, iteratively reducing the chance of |
| // rejection. The replacement method does however add a little |
| // overhead. Benchmarking or calculating probabilities might |
| // reveal contexts where this replacement method is slower. |
| </span><span class="kw">let </span><span class="kw-2">mut </span>v: <span class="macro-nonterminal">$unsigned </span>= rng.gen(); |
| <span class="kw">loop </span>{ |
| <span class="kw">let </span>(hi, lo) = v.wmul(range); |
| <span class="kw">let </span>mask = lo.le(zone); |
| <span class="kw">if </span>mask.all() { |
| <span class="kw">let </span>hi: <span class="macro-nonterminal">$ty </span>= hi.cast(); |
| <span class="comment">// wrapping addition |
| </span><span class="kw">let </span>result = <span class="self">self</span>.low + hi; |
| <span class="comment">// `select` here compiles to a blend operation |
| // When `range.eq(0).none()` the compare and blend |
| // operations are avoided. |
| </span><span class="kw">let </span>v: <span class="macro-nonterminal">$ty </span>= v.cast(); |
| <span class="kw">return </span>range.gt(<span class="macro-nonterminal">$unsigned::splat</span>(<span class="number">0</span>)).select(result, v); |
| } |
| <span class="comment">// Replace only the failing lanes |
| </span>v = mask.select(v, rng.gen()); |
| } |
| } |
| } |
| }; |
| |
| <span class="comment">// bulk implementation |
| </span>($((<span class="macro-nonterminal">$unsigned</span>:ident, <span class="macro-nonterminal">$signed</span>:ident),)+ <span class="macro-nonterminal">$u_scalar</span>:ident) => { |
| $( |
| <span class="macro">uniform_simd_int_impl!</span>(<span class="macro-nonterminal">$unsigned</span>, <span class="macro-nonterminal">$unsigned</span>, <span class="macro-nonterminal">$u_scalar</span>); |
| <span class="macro">uniform_simd_int_impl!</span>(<span class="macro-nonterminal">$signed</span>, <span class="macro-nonterminal">$unsigned</span>, <span class="macro-nonterminal">$u_scalar</span>); |
| )+ |
| }; |
| } |
| |
| <span class="attribute">#[cfg(feature = <span class="string">"simd_support"</span>)] |
| </span><span class="macro">uniform_simd_int_impl! </span>{ |
| (u64x2, i64x2), |
| (u64x4, i64x4), |
| (u64x8, i64x8), |
| u64 |
| } |
| |
| <span class="attribute">#[cfg(feature = <span class="string">"simd_support"</span>)] |
| </span><span class="macro">uniform_simd_int_impl! </span>{ |
| (u32x2, i32x2), |
| (u32x4, i32x4), |
| (u32x8, i32x8), |
| (u32x16, i32x16), |
| u32 |
| } |
| |
| <span class="attribute">#[cfg(feature = <span class="string">"simd_support"</span>)] |
| </span><span class="macro">uniform_simd_int_impl! </span>{ |
| (u16x2, i16x2), |
| (u16x4, i16x4), |
| (u16x8, i16x8), |
| (u16x16, i16x16), |
| (u16x32, i16x32), |
| u16 |
| } |
| |
| <span class="attribute">#[cfg(feature = <span class="string">"simd_support"</span>)] |
| </span><span class="macro">uniform_simd_int_impl! </span>{ |
| (u8x2, i8x2), |
| (u8x4, i8x4), |
| (u8x8, i8x8), |
| (u8x16, i8x16), |
| (u8x32, i8x32), |
| (u8x64, i8x64), |
| u8 |
| } |
| |
| <span class="kw">impl </span>SampleUniform <span class="kw">for </span>char { |
| <span class="kw">type </span>Sampler = UniformChar; |
| } |
| |
| <span class="doccomment">/// The back-end implementing [`UniformSampler`] for `char`. |
| /// |
| /// Unless you are implementing [`UniformSampler`] for your own type, this type |
| /// should not be used directly, use [`Uniform`] instead. |
| /// |
| /// This differs from integer range sampling since the range `0xD800..=0xDFFF` |
| /// are used for surrogate pairs in UCS and UTF-16, and consequently are not |
| /// valid Unicode code points. We must therefore avoid sampling values in this |
| /// range. |
| </span><span class="attribute">#[derive(Clone, Copy, Debug)] |
| #[cfg_attr(feature = <span class="string">"serde1"</span>, derive(Serialize, Deserialize))] |
| </span><span class="kw">pub struct </span>UniformChar { |
| sampler: UniformInt<u32>, |
| } |
| |
| <span class="doccomment">/// UTF-16 surrogate range start |
| </span><span class="kw">const </span>CHAR_SURROGATE_START: u32 = <span class="number">0xD800</span>; |
| <span class="doccomment">/// UTF-16 surrogate range size |
| </span><span class="kw">const </span>CHAR_SURROGATE_LEN: u32 = <span class="number">0xE000 </span>- CHAR_SURROGATE_START; |
| |
| <span class="doccomment">/// Convert `char` to compressed `u32` |
| </span><span class="kw">fn </span>char_to_comp_u32(c: char) -> u32 { |
| <span class="kw">match </span>c <span class="kw">as </span>u32 { |
| c <span class="kw">if </span>c >= CHAR_SURROGATE_START => c - CHAR_SURROGATE_LEN, |
| c => c, |
| } |
| } |
| |
| <span class="kw">impl </span>UniformSampler <span class="kw">for </span>UniformChar { |
| <span class="kw">type </span>X = char; |
| |
| <span class="attribute">#[inline] </span><span class="comment">// if the range is constant, this helps LLVM to do the |
| // calculations at compile-time. |
| </span><span class="kw">fn </span>new<B1, B2>(low_b: B1, high_b: B2) -> <span class="self">Self |
| </span><span class="kw">where |
| </span>B1: SampleBorrow<<span class="self">Self</span>::X> + Sized, |
| B2: SampleBorrow<<span class="self">Self</span>::X> + Sized, |
| { |
| <span class="kw">let </span>low = char_to_comp_u32(<span class="kw-2">*</span>low_b.borrow()); |
| <span class="kw">let </span>high = char_to_comp_u32(<span class="kw-2">*</span>high_b.borrow()); |
| <span class="kw">let </span>sampler = UniformInt::<u32>::new(low, high); |
| UniformChar { sampler } |
| } |
| |
| <span class="attribute">#[inline] </span><span class="comment">// if the range is constant, this helps LLVM to do the |
| // calculations at compile-time. |
| </span><span class="kw">fn </span>new_inclusive<B1, B2>(low_b: B1, high_b: B2) -> <span class="self">Self |
| </span><span class="kw">where |
| </span>B1: SampleBorrow<<span class="self">Self</span>::X> + Sized, |
| B2: SampleBorrow<<span class="self">Self</span>::X> + Sized, |
| { |
| <span class="kw">let </span>low = char_to_comp_u32(<span class="kw-2">*</span>low_b.borrow()); |
| <span class="kw">let </span>high = char_to_comp_u32(<span class="kw-2">*</span>high_b.borrow()); |
| <span class="kw">let </span>sampler = UniformInt::<u32>::new_inclusive(low, high); |
| UniformChar { sampler } |
| } |
| |
| <span class="kw">fn </span>sample<R: Rng + <span class="question-mark">?</span>Sized>(<span class="kw-2">&</span><span class="self">self</span>, rng: <span class="kw-2">&mut </span>R) -> <span class="self">Self</span>::X { |
| <span class="kw">let </span><span class="kw-2">mut </span>x = <span class="self">self</span>.sampler.sample(rng); |
| <span class="kw">if </span>x >= CHAR_SURROGATE_START { |
| x += CHAR_SURROGATE_LEN; |
| } |
| <span class="comment">// SAFETY: x must not be in surrogate range or greater than char::MAX. |
| // This relies on range constructors which accept char arguments. |
| // Validity of input char values is assumed. |
| </span><span class="kw">unsafe </span>{ core::char::from_u32_unchecked(x) } |
| } |
| } |
| |
| <span class="doccomment">/// The back-end implementing [`UniformSampler`] for floating-point types. |
| /// |
| /// Unless you are implementing [`UniformSampler`] for your own type, this type |
| /// should not be used directly, use [`Uniform`] instead. |
| /// |
| /// # Implementation notes |
| /// |
| /// Instead of generating a float in the `[0, 1)` range using [`Standard`], the |
| /// `UniformFloat` implementation converts the output of an PRNG itself. This |
| /// way one or two steps can be optimized out. |
| /// |
| /// The floats are first converted to a value in the `[1, 2)` interval using a |
| /// transmute-based method, and then mapped to the expected range with a |
| /// multiply and addition. Values produced this way have what equals 23 bits of |
| /// random digits for an `f32`, and 52 for an `f64`. |
| /// |
| /// [`new`]: UniformSampler::new |
| /// [`new_inclusive`]: UniformSampler::new_inclusive |
| /// [`Standard`]: crate::distributions::Standard |
| </span><span class="attribute">#[derive(Clone, Copy, Debug, PartialEq)] |
| #[cfg_attr(feature = <span class="string">"serde1"</span>, derive(Serialize, Deserialize))] |
| </span><span class="kw">pub struct </span>UniformFloat<X> { |
| low: X, |
| scale: X, |
| } |
| |
| <span class="macro">macro_rules! </span>uniform_float_impl { |
| (<span class="macro-nonterminal">$ty</span>:ty, <span class="macro-nonterminal">$uty</span>:ident, <span class="macro-nonterminal">$f_scalar</span>:ident, <span class="macro-nonterminal">$u_scalar</span>:ident, <span class="macro-nonterminal">$bits_to_discard</span>:expr) => { |
| <span class="kw">impl </span>SampleUniform <span class="kw">for </span><span class="macro-nonterminal">$ty </span>{ |
| <span class="kw">type </span>Sampler = UniformFloat<<span class="macro-nonterminal">$ty</span>>; |
| } |
| |
| <span class="kw">impl </span>UniformSampler <span class="kw">for </span>UniformFloat<<span class="macro-nonterminal">$ty</span>> { |
| <span class="kw">type </span>X = <span class="macro-nonterminal">$ty</span>; |
| |
| <span class="kw">fn </span>new<B1, B2>(low_b: B1, high_b: B2) -> <span class="self">Self |
| </span><span class="kw">where |
| </span>B1: SampleBorrow<<span class="self">Self</span>::X> + Sized, |
| B2: SampleBorrow<<span class="self">Self</span>::X> + Sized, |
| { |
| <span class="kw">let </span>low = <span class="kw-2">*</span>low_b.borrow(); |
| <span class="kw">let </span>high = <span class="kw-2">*</span>high_b.borrow(); |
| <span class="macro">debug_assert!</span>( |
| low.all_finite(), |
| <span class="string">"Uniform::new called with `low` non-finite." |
| </span>); |
| <span class="macro">debug_assert!</span>( |
| high.all_finite(), |
| <span class="string">"Uniform::new called with `high` non-finite." |
| </span>); |
| <span class="macro">assert!</span>(low.all_lt(high), <span class="string">"Uniform::new called with `low >= high`"</span>); |
| <span class="kw">let </span>max_rand = <<span class="macro-nonterminal">$ty</span>>::splat( |
| (::core::<span class="macro-nonterminal">$u_scalar::MAX </span>>> <span class="macro-nonterminal">$bits_to_discard</span>).into_float_with_exponent(<span class="number">0</span>) - <span class="number">1.0</span>, |
| ); |
| |
| <span class="kw">let </span><span class="kw-2">mut </span>scale = high - low; |
| <span class="macro">assert!</span>(scale.all_finite(), <span class="string">"Uniform::new: range overflow"</span>); |
| |
| <span class="kw">loop </span>{ |
| <span class="kw">let </span>mask = (scale * max_rand + low).ge_mask(high); |
| <span class="kw">if </span>mask.none() { |
| <span class="kw">break</span>; |
| } |
| scale = scale.decrease_masked(mask); |
| } |
| |
| <span class="macro">debug_assert!</span>(<<span class="macro-nonterminal">$ty</span>>::splat(<span class="number">0.0</span>).all_le(scale)); |
| |
| UniformFloat { low, scale } |
| } |
| |
| <span class="kw">fn </span>new_inclusive<B1, B2>(low_b: B1, high_b: B2) -> <span class="self">Self |
| </span><span class="kw">where |
| </span>B1: SampleBorrow<<span class="self">Self</span>::X> + Sized, |
| B2: SampleBorrow<<span class="self">Self</span>::X> + Sized, |
| { |
| <span class="kw">let </span>low = <span class="kw-2">*</span>low_b.borrow(); |
| <span class="kw">let </span>high = <span class="kw-2">*</span>high_b.borrow(); |
| <span class="macro">debug_assert!</span>( |
| low.all_finite(), |
| <span class="string">"Uniform::new_inclusive called with `low` non-finite." |
| </span>); |
| <span class="macro">debug_assert!</span>( |
| high.all_finite(), |
| <span class="string">"Uniform::new_inclusive called with `high` non-finite." |
| </span>); |
| <span class="macro">assert!</span>( |
| low.all_le(high), |
| <span class="string">"Uniform::new_inclusive called with `low > high`" |
| </span>); |
| <span class="kw">let </span>max_rand = <<span class="macro-nonterminal">$ty</span>>::splat( |
| (::core::<span class="macro-nonterminal">$u_scalar::MAX </span>>> <span class="macro-nonterminal">$bits_to_discard</span>).into_float_with_exponent(<span class="number">0</span>) - <span class="number">1.0</span>, |
| ); |
| |
| <span class="kw">let </span><span class="kw-2">mut </span>scale = (high - low) / max_rand; |
| <span class="macro">assert!</span>(scale.all_finite(), <span class="string">"Uniform::new_inclusive: range overflow"</span>); |
| |
| <span class="kw">loop </span>{ |
| <span class="kw">let </span>mask = (scale * max_rand + low).gt_mask(high); |
| <span class="kw">if </span>mask.none() { |
| <span class="kw">break</span>; |
| } |
| scale = scale.decrease_masked(mask); |
| } |
| |
| <span class="macro">debug_assert!</span>(<<span class="macro-nonterminal">$ty</span>>::splat(<span class="number">0.0</span>).all_le(scale)); |
| |
| UniformFloat { low, scale } |
| } |
| |
| <span class="kw">fn </span>sample<R: Rng + <span class="question-mark">?</span>Sized>(<span class="kw-2">&</span><span class="self">self</span>, rng: <span class="kw-2">&mut </span>R) -> <span class="self">Self</span>::X { |
| <span class="comment">// Generate a value in the range [1, 2) |
| </span><span class="kw">let </span>value1_2 = (rng.gen::<<span class="macro-nonterminal">$uty</span>>() >> <span class="macro-nonterminal">$bits_to_discard</span>).into_float_with_exponent(<span class="number">0</span>); |
| |
| <span class="comment">// Get a value in the range [0, 1) in order to avoid |
| // overflowing into infinity when multiplying with scale |
| </span><span class="kw">let </span>value0_1 = value1_2 - <span class="number">1.0</span>; |
| |
| <span class="comment">// We don't use `f64::mul_add`, because it is not available with |
| // `no_std`. Furthermore, it is slower for some targets (but |
| // faster for others). However, the order of multiplication and |
| // addition is important, because on some platforms (e.g. ARM) |
| // it will be optimized to a single (non-FMA) instruction. |
| </span>value0_1 * <span class="self">self</span>.scale + <span class="self">self</span>.low |
| } |
| |
| <span class="attribute">#[inline] |
| </span><span class="kw">fn </span>sample_single<R: Rng + <span class="question-mark">?</span>Sized, B1, B2>(low_b: B1, high_b: B2, rng: <span class="kw-2">&mut </span>R) -> <span class="self">Self</span>::X |
| <span class="kw">where |
| </span>B1: SampleBorrow<<span class="self">Self</span>::X> + Sized, |
| B2: SampleBorrow<<span class="self">Self</span>::X> + Sized, |
| { |
| <span class="kw">let </span>low = <span class="kw-2">*</span>low_b.borrow(); |
| <span class="kw">let </span>high = <span class="kw-2">*</span>high_b.borrow(); |
| <span class="macro">debug_assert!</span>( |
| low.all_finite(), |
| <span class="string">"UniformSampler::sample_single called with `low` non-finite." |
| </span>); |
| <span class="macro">debug_assert!</span>( |
| high.all_finite(), |
| <span class="string">"UniformSampler::sample_single called with `high` non-finite." |
| </span>); |
| <span class="macro">assert!</span>( |
| low.all_lt(high), |
| <span class="string">"UniformSampler::sample_single: low >= high" |
| </span>); |
| <span class="kw">let </span><span class="kw-2">mut </span>scale = high - low; |
| <span class="macro">assert!</span>(scale.all_finite(), <span class="string">"UniformSampler::sample_single: range overflow"</span>); |
| |
| <span class="kw">loop </span>{ |
| <span class="comment">// Generate a value in the range [1, 2) |
| </span><span class="kw">let </span>value1_2 = |
| (rng.gen::<<span class="macro-nonterminal">$uty</span>>() >> <span class="macro-nonterminal">$bits_to_discard</span>).into_float_with_exponent(<span class="number">0</span>); |
| |
| <span class="comment">// Get a value in the range [0, 1) in order to avoid |
| // overflowing into infinity when multiplying with scale |
| </span><span class="kw">let </span>value0_1 = value1_2 - <span class="number">1.0</span>; |
| |
| <span class="comment">// Doing multiply before addition allows some architectures |
| // to use a single instruction. |
| </span><span class="kw">let </span>res = value0_1 * scale + low; |
| |
| <span class="macro">debug_assert!</span>(low.all_le(res) || !scale.all_finite()); |
| <span class="kw">if </span>res.all_lt(high) { |
| <span class="kw">return </span>res; |
| } |
| |
| <span class="comment">// This handles a number of edge cases. |
| // * `low` or `high` is NaN. In this case `scale` and |
| // `res` are going to end up as NaN. |
| // * `low` is negative infinity and `high` is finite. |
| // `scale` is going to be infinite and `res` will be |
| // NaN. |
| // * `high` is positive infinity and `low` is finite. |
| // `scale` is going to be infinite and `res` will |
| // be infinite or NaN (if value0_1 is 0). |
| // * `low` is negative infinity and `high` is positive |
| // infinity. `scale` will be infinite and `res` will |
| // be NaN. |
| // * `low` and `high` are finite, but `high - low` |
| // overflows to infinite. `scale` will be infinite |
| // and `res` will be infinite or NaN (if value0_1 is 0). |
| // So if `high` or `low` are non-finite, we are guaranteed |
| // to fail the `res < high` check above and end up here. |
| // |
| // While we technically should check for non-finite `low` |
| // and `high` before entering the loop, by doing the checks |
| // here instead, we allow the common case to avoid these |
| // checks. But we are still guaranteed that if `low` or |
| // `high` are non-finite we'll end up here and can do the |
| // appropriate checks. |
| // |
| // Likewise `high - low` overflowing to infinity is also |
| // rare, so handle it here after the common case. |
| </span><span class="kw">let </span>mask = !scale.finite_mask(); |
| <span class="kw">if </span>mask.any() { |
| <span class="macro">assert!</span>( |
| low.all_finite() && high.all_finite(), |
| <span class="string">"Uniform::sample_single: low and high must be finite" |
| </span>); |
| scale = scale.decrease_masked(mask); |
| } |
| } |
| } |
| } |
| }; |
| } |
| |
| <span class="macro">uniform_float_impl! </span>{ f32, u32, f32, u32, <span class="number">32 </span>- <span class="number">23 </span>} |
| <span class="macro">uniform_float_impl! </span>{ f64, u64, f64, u64, <span class="number">64 </span>- <span class="number">52 </span>} |
| |
| <span class="attribute">#[cfg(feature = <span class="string">"simd_support"</span>)] |
| </span><span class="macro">uniform_float_impl! </span>{ f32x2, u32x2, f32, u32, <span class="number">32 </span>- <span class="number">23 </span>} |
| <span class="attribute">#[cfg(feature = <span class="string">"simd_support"</span>)] |
| </span><span class="macro">uniform_float_impl! </span>{ f32x4, u32x4, f32, u32, <span class="number">32 </span>- <span class="number">23 </span>} |
| <span class="attribute">#[cfg(feature = <span class="string">"simd_support"</span>)] |
| </span><span class="macro">uniform_float_impl! </span>{ f32x8, u32x8, f32, u32, <span class="number">32 </span>- <span class="number">23 </span>} |
| <span class="attribute">#[cfg(feature = <span class="string">"simd_support"</span>)] |
| </span><span class="macro">uniform_float_impl! </span>{ f32x16, u32x16, f32, u32, <span class="number">32 </span>- <span class="number">23 </span>} |
| |
| <span class="attribute">#[cfg(feature = <span class="string">"simd_support"</span>)] |
| </span><span class="macro">uniform_float_impl! </span>{ f64x2, u64x2, f64, u64, <span class="number">64 </span>- <span class="number">52 </span>} |
| <span class="attribute">#[cfg(feature = <span class="string">"simd_support"</span>)] |
| </span><span class="macro">uniform_float_impl! </span>{ f64x4, u64x4, f64, u64, <span class="number">64 </span>- <span class="number">52 </span>} |
| <span class="attribute">#[cfg(feature = <span class="string">"simd_support"</span>)] |
| </span><span class="macro">uniform_float_impl! </span>{ f64x8, u64x8, f64, u64, <span class="number">64 </span>- <span class="number">52 </span>} |
| |
| |
| <span class="doccomment">/// The back-end implementing [`UniformSampler`] for `Duration`. |
| /// |
| /// Unless you are implementing [`UniformSampler`] for your own types, this type |
| /// should not be used directly, use [`Uniform`] instead. |
| </span><span class="attribute">#[derive(Clone, Copy, Debug)] |
| #[cfg_attr(feature = <span class="string">"serde1"</span>, derive(Serialize, Deserialize))] |
| </span><span class="kw">pub struct </span>UniformDuration { |
| mode: UniformDurationMode, |
| offset: u32, |
| } |
| |
| <span class="attribute">#[derive(Debug, Copy, Clone)] |
| #[cfg_attr(feature = <span class="string">"serde1"</span>, derive(Serialize, Deserialize))] |
| </span><span class="kw">enum </span>UniformDurationMode { |
| Small { |
| secs: u64, |
| nanos: Uniform<u32>, |
| }, |
| Medium { |
| nanos: Uniform<u64>, |
| }, |
| Large { |
| max_secs: u64, |
| max_nanos: u32, |
| secs: Uniform<u64>, |
| }, |
| } |
| |
| <span class="kw">impl </span>SampleUniform <span class="kw">for </span>Duration { |
| <span class="kw">type </span>Sampler = UniformDuration; |
| } |
| |
| <span class="kw">impl </span>UniformSampler <span class="kw">for </span>UniformDuration { |
| <span class="kw">type </span>X = Duration; |
| |
| <span class="attribute">#[inline] |
| </span><span class="kw">fn </span>new<B1, B2>(low_b: B1, high_b: B2) -> <span class="self">Self |
| </span><span class="kw">where |
| </span>B1: SampleBorrow<<span class="self">Self</span>::X> + Sized, |
| B2: SampleBorrow<<span class="self">Self</span>::X> + Sized, |
| { |
| <span class="kw">let </span>low = <span class="kw-2">*</span>low_b.borrow(); |
| <span class="kw">let </span>high = <span class="kw-2">*</span>high_b.borrow(); |
| <span class="macro">assert!</span>(low < high, <span class="string">"Uniform::new called with `low >= high`"</span>); |
| UniformDuration::new_inclusive(low, high - Duration::new(<span class="number">0</span>, <span class="number">1</span>)) |
| } |
| |
| <span class="attribute">#[inline] |
| </span><span class="kw">fn </span>new_inclusive<B1, B2>(low_b: B1, high_b: B2) -> <span class="self">Self |
| </span><span class="kw">where |
| </span>B1: SampleBorrow<<span class="self">Self</span>::X> + Sized, |
| B2: SampleBorrow<<span class="self">Self</span>::X> + Sized, |
| { |
| <span class="kw">let </span>low = <span class="kw-2">*</span>low_b.borrow(); |
| <span class="kw">let </span>high = <span class="kw-2">*</span>high_b.borrow(); |
| <span class="macro">assert!</span>( |
| low <= high, |
| <span class="string">"Uniform::new_inclusive called with `low > high`" |
| </span>); |
| |
| <span class="kw">let </span>low_s = low.as_secs(); |
| <span class="kw">let </span>low_n = low.subsec_nanos(); |
| <span class="kw">let </span><span class="kw-2">mut </span>high_s = high.as_secs(); |
| <span class="kw">let </span><span class="kw-2">mut </span>high_n = high.subsec_nanos(); |
| |
| <span class="kw">if </span>high_n < low_n { |
| high_s -= <span class="number">1</span>; |
| high_n += <span class="number">1_000_000_000</span>; |
| } |
| |
| <span class="kw">let </span>mode = <span class="kw">if </span>low_s == high_s { |
| UniformDurationMode::Small { |
| secs: low_s, |
| nanos: Uniform::new_inclusive(low_n, high_n), |
| } |
| } <span class="kw">else </span>{ |
| <span class="kw">let </span>max = high_s |
| .checked_mul(<span class="number">1_000_000_000</span>) |
| .and_then(|n| n.checked_add(u64::from(high_n))); |
| |
| <span class="kw">if let </span><span class="prelude-val">Some</span>(higher_bound) = max { |
| <span class="kw">let </span>lower_bound = low_s * <span class="number">1_000_000_000 </span>+ u64::from(low_n); |
| UniformDurationMode::Medium { |
| nanos: Uniform::new_inclusive(lower_bound, higher_bound), |
| } |
| } <span class="kw">else </span>{ |
| <span class="comment">// An offset is applied to simplify generation of nanoseconds |
| </span><span class="kw">let </span>max_nanos = high_n - low_n; |
| UniformDurationMode::Large { |
| max_secs: high_s, |
| max_nanos, |
| secs: Uniform::new_inclusive(low_s, high_s), |
| } |
| } |
| }; |
| UniformDuration { |
| mode, |
| offset: low_n, |
| } |
| } |
| |
| <span class="attribute">#[inline] |
| </span><span class="kw">fn </span>sample<R: Rng + <span class="question-mark">?</span>Sized>(<span class="kw-2">&</span><span class="self">self</span>, rng: <span class="kw-2">&mut </span>R) -> Duration { |
| <span class="kw">match </span><span class="self">self</span>.mode { |
| UniformDurationMode::Small { secs, nanos } => { |
| <span class="kw">let </span>n = nanos.sample(rng); |
| Duration::new(secs, n) |
| } |
| UniformDurationMode::Medium { nanos } => { |
| <span class="kw">let </span>nanos = nanos.sample(rng); |
| Duration::new(nanos / <span class="number">1_000_000_000</span>, (nanos % <span class="number">1_000_000_000</span>) <span class="kw">as </span>u32) |
| } |
| UniformDurationMode::Large { |
| max_secs, |
| max_nanos, |
| secs, |
| } => { |
| <span class="comment">// constant folding means this is at least as fast as `Rng::sample(Range)` |
| </span><span class="kw">let </span>nano_range = Uniform::new(<span class="number">0</span>, <span class="number">1_000_000_000</span>); |
| <span class="kw">loop </span>{ |
| <span class="kw">let </span>s = secs.sample(rng); |
| <span class="kw">let </span>n = nano_range.sample(rng); |
| <span class="kw">if </span>!(s == max_secs && n > max_nanos) { |
| <span class="kw">let </span>sum = n + <span class="self">self</span>.offset; |
| <span class="kw">break </span>Duration::new(s, sum); |
| } |
| } |
| } |
| } |
| } |
| } |
| |
| <span class="attribute">#[cfg(test)] |
| </span><span class="kw">mod </span>tests { |
| <span class="kw">use super</span>::<span class="kw-2">*</span>; |
| <span class="kw">use </span><span class="kw">crate</span>::rngs::mock::StepRng; |
| |
| <span class="attribute">#[test] |
| #[cfg(feature = <span class="string">"serde1"</span>)] |
| </span><span class="kw">fn </span>test_serialization_uniform_duration() { |
| <span class="kw">let </span>distr = UniformDuration::new(Duration::from_secs(<span class="number">10</span>), Duration::from_secs(<span class="number">60</span>)); |
| <span class="kw">let </span>de_distr: UniformDuration = bincode::deserialize(<span class="kw-2">&</span>bincode::serialize(<span class="kw-2">&</span>distr).unwrap()).unwrap(); |
| <span class="macro">assert_eq!</span>( |
| distr.offset, de_distr.offset |
| ); |
| <span class="kw">match </span>(distr.mode, de_distr.mode) { |
| (UniformDurationMode::Small {secs: a_secs, nanos: a_nanos}, UniformDurationMode::Small {secs, nanos}) => { |
| <span class="macro">assert_eq!</span>(a_secs, secs); |
| |
| <span class="macro">assert_eq!</span>(a_nanos.<span class="number">0</span>.low, nanos.<span class="number">0</span>.low); |
| <span class="macro">assert_eq!</span>(a_nanos.<span class="number">0</span>.range, nanos.<span class="number">0</span>.range); |
| <span class="macro">assert_eq!</span>(a_nanos.<span class="number">0</span>.z, nanos.<span class="number">0</span>.z); |
| } |
| (UniformDurationMode::Medium {nanos: a_nanos} , UniformDurationMode::Medium {nanos}) => { |
| <span class="macro">assert_eq!</span>(a_nanos.<span class="number">0</span>.low, nanos.<span class="number">0</span>.low); |
| <span class="macro">assert_eq!</span>(a_nanos.<span class="number">0</span>.range, nanos.<span class="number">0</span>.range); |
| <span class="macro">assert_eq!</span>(a_nanos.<span class="number">0</span>.z, nanos.<span class="number">0</span>.z); |
| } |
| (UniformDurationMode::Large {max_secs:a_max_secs, max_nanos:a_max_nanos, secs:a_secs}, UniformDurationMode::Large {max_secs, max_nanos, secs} ) => { |
| <span class="macro">assert_eq!</span>(a_max_secs, max_secs); |
| <span class="macro">assert_eq!</span>(a_max_nanos, max_nanos); |
| |
| <span class="macro">assert_eq!</span>(a_secs.<span class="number">0</span>.low, secs.<span class="number">0</span>.low); |
| <span class="macro">assert_eq!</span>(a_secs.<span class="number">0</span>.range, secs.<span class="number">0</span>.range); |
| <span class="macro">assert_eq!</span>(a_secs.<span class="number">0</span>.z, secs.<span class="number">0</span>.z); |
| } |
| <span class="kw">_ </span>=> <span class="macro">panic!</span>(<span class="string">"`UniformDurationMode` was not serialized/deserialized correctly"</span>) |
| } |
| } |
| |
| <span class="attribute">#[test] |
| #[cfg(feature = <span class="string">"serde1"</span>)] |
| </span><span class="kw">fn </span>test_uniform_serialization() { |
| <span class="kw">let </span>unit_box: Uniform<i32> = Uniform::new(-<span class="number">1</span>, <span class="number">1</span>); |
| <span class="kw">let </span>de_unit_box: Uniform<i32> = bincode::deserialize(<span class="kw-2">&</span>bincode::serialize(<span class="kw-2">&</span>unit_box).unwrap()).unwrap(); |
| |
| <span class="macro">assert_eq!</span>(unit_box.<span class="number">0</span>.low, de_unit_box.<span class="number">0</span>.low); |
| <span class="macro">assert_eq!</span>(unit_box.<span class="number">0</span>.range, de_unit_box.<span class="number">0</span>.range); |
| <span class="macro">assert_eq!</span>(unit_box.<span class="number">0</span>.z, de_unit_box.<span class="number">0</span>.z); |
| |
| <span class="kw">let </span>unit_box: Uniform<f32> = Uniform::new(-<span class="number">1.</span>, <span class="number">1.</span>); |
| <span class="kw">let </span>de_unit_box: Uniform<f32> = bincode::deserialize(<span class="kw-2">&</span>bincode::serialize(<span class="kw-2">&</span>unit_box).unwrap()).unwrap(); |
| |
| <span class="macro">assert_eq!</span>(unit_box.<span class="number">0</span>.low, de_unit_box.<span class="number">0</span>.low); |
| <span class="macro">assert_eq!</span>(unit_box.<span class="number">0</span>.scale, de_unit_box.<span class="number">0</span>.scale); |
| } |
| |
| <span class="attribute">#[should_panic] |
| #[test] |
| </span><span class="kw">fn </span>test_uniform_bad_limits_equal_int() { |
| Uniform::new(<span class="number">10</span>, <span class="number">10</span>); |
| } |
| |
| <span class="attribute">#[test] |
| </span><span class="kw">fn </span>test_uniform_good_limits_equal_int() { |
| <span class="kw">let </span><span class="kw-2">mut </span>rng = <span class="kw">crate</span>::test::rng(<span class="number">804</span>); |
| <span class="kw">let </span>dist = Uniform::new_inclusive(<span class="number">10</span>, <span class="number">10</span>); |
| <span class="kw">for _ in </span><span class="number">0</span>..<span class="number">20 </span>{ |
| <span class="macro">assert_eq!</span>(rng.sample(dist), <span class="number">10</span>); |
| } |
| } |
| |
| <span class="attribute">#[should_panic] |
| #[test] |
| </span><span class="kw">fn </span>test_uniform_bad_limits_flipped_int() { |
| Uniform::new(<span class="number">10</span>, <span class="number">5</span>); |
| } |
| |
| <span class="attribute">#[test] |
| #[cfg_attr(miri, ignore)] </span><span class="comment">// Miri is too slow |
| </span><span class="kw">fn </span>test_integers() { |
| <span class="kw">use </span>core::{i128, u128}; |
| <span class="kw">use </span>core::{i16, i32, i64, i8, isize}; |
| <span class="kw">use </span>core::{u16, u32, u64, u8, usize}; |
| |
| <span class="kw">let </span><span class="kw-2">mut </span>rng = <span class="kw">crate</span>::test::rng(<span class="number">251</span>); |
| <span class="macro">macro_rules! </span>t { |
| (<span class="macro-nonterminal">$ty</span>:ident, <span class="macro-nonterminal">$v</span>:expr, <span class="macro-nonterminal">$le</span>:expr, <span class="macro-nonterminal">$lt</span>:expr) => {{ |
| <span class="kw">for </span><span class="kw-2">&</span>(low, high) <span class="kw">in </span><span class="macro-nonterminal">$v</span>.iter() { |
| <span class="kw">let </span>my_uniform = Uniform::new(low, high); |
| <span class="kw">for _ in </span><span class="number">0</span>..<span class="number">1000 </span>{ |
| <span class="kw">let </span>v: <span class="macro-nonterminal">$ty </span>= rng.sample(my_uniform); |
| <span class="macro">assert!</span>(<span class="macro-nonterminal">$le</span>(low, v) && <span class="macro-nonterminal">$lt</span>(v, high)); |
| } |
| |
| <span class="kw">let </span>my_uniform = Uniform::new_inclusive(low, high); |
| <span class="kw">for _ in </span><span class="number">0</span>..<span class="number">1000 </span>{ |
| <span class="kw">let </span>v: <span class="macro-nonterminal">$ty </span>= rng.sample(my_uniform); |
| <span class="macro">assert!</span>(<span class="macro-nonterminal">$le</span>(low, v) && <span class="macro-nonterminal">$le</span>(v, high)); |
| } |
| |
| <span class="kw">let </span>my_uniform = Uniform::new(<span class="kw-2">&</span>low, high); |
| <span class="kw">for _ in </span><span class="number">0</span>..<span class="number">1000 </span>{ |
| <span class="kw">let </span>v: <span class="macro-nonterminal">$ty </span>= rng.sample(my_uniform); |
| <span class="macro">assert!</span>(<span class="macro-nonterminal">$le</span>(low, v) && <span class="macro-nonterminal">$lt</span>(v, high)); |
| } |
| |
| <span class="kw">let </span>my_uniform = Uniform::new_inclusive(<span class="kw-2">&</span>low, <span class="kw-2">&</span>high); |
| <span class="kw">for _ in </span><span class="number">0</span>..<span class="number">1000 </span>{ |
| <span class="kw">let </span>v: <span class="macro-nonterminal">$ty </span>= rng.sample(my_uniform); |
| <span class="macro">assert!</span>(<span class="macro-nonterminal">$le</span>(low, v) && <span class="macro-nonterminal">$le</span>(v, high)); |
| } |
| |
| <span class="kw">for _ in </span><span class="number">0</span>..<span class="number">1000 </span>{ |
| <span class="kw">let </span>v = <<span class="macro-nonterminal">$ty </span><span class="kw">as </span>SampleUniform>::Sampler::sample_single(low, high, <span class="kw-2">&mut </span>rng); |
| <span class="macro">assert!</span>(<span class="macro-nonterminal">$le</span>(low, v) && <span class="macro-nonterminal">$lt</span>(v, high)); |
| } |
| |
| <span class="kw">for _ in </span><span class="number">0</span>..<span class="number">1000 </span>{ |
| <span class="kw">let </span>v = <<span class="macro-nonterminal">$ty </span><span class="kw">as </span>SampleUniform>::Sampler::sample_single_inclusive(low, high, <span class="kw-2">&mut </span>rng); |
| <span class="macro">assert!</span>(<span class="macro-nonterminal">$le</span>(low, v) && <span class="macro-nonterminal">$le</span>(v, high)); |
| } |
| } |
| }}; |
| |
| <span class="comment">// scalar bulk |
| </span>($(<span class="macro-nonterminal">$ty</span>:ident),<span class="kw-2">*</span>) => {{ |
| $(<span class="macro">t!</span>( |
| <span class="macro-nonterminal">$ty</span>, |
| [(<span class="number">0</span>, <span class="number">10</span>), (<span class="number">10</span>, <span class="number">127</span>), (<span class="macro-nonterminal">$ty::MIN</span>, <span class="macro-nonterminal">$ty::MAX</span>)], |
| |x, y| x <= y, |
| |x, y| x < y |
| );)* |
| }}; |
| |
| <span class="comment">// simd bulk |
| </span>($(<span class="macro-nonterminal">$ty</span>:ident),* => <span class="macro-nonterminal">$scalar</span>:ident) => {{ |
| $(<span class="macro">t!</span>( |
| <span class="macro-nonterminal">$ty</span>, |
| [ |
| (<span class="macro-nonterminal">$ty::splat</span>(<span class="number">0</span>), <span class="macro-nonterminal">$ty::splat</span>(<span class="number">10</span>)), |
| (<span class="macro-nonterminal">$ty::splat</span>(<span class="number">10</span>), <span class="macro-nonterminal">$ty::splat</span>(<span class="number">127</span>)), |
| (<span class="macro-nonterminal">$ty::splat</span>(<span class="macro-nonterminal">$scalar::MIN</span>), <span class="macro-nonterminal">$ty::splat</span>(<span class="macro-nonterminal">$scalar::MAX</span>)), |
| ], |
| |x: <span class="macro-nonterminal">$ty</span>, y| x.le(y).all(), |
| |x: <span class="macro-nonterminal">$ty</span>, y| x.lt(y).all() |
| );)* |
| }}; |
| } |
| <span class="macro">t!</span>(i8, i16, i32, i64, isize, u8, u16, u32, u64, usize, i128, u128); |
| |
| <span class="attribute">#[cfg(feature = <span class="string">"simd_support"</span>)] |
| </span>{ |
| <span class="macro">t!</span>(u8x2, u8x4, u8x8, u8x16, u8x32, u8x64 => u8); |
| <span class="macro">t!</span>(i8x2, i8x4, i8x8, i8x16, i8x32, i8x64 => i8); |
| <span class="macro">t!</span>(u16x2, u16x4, u16x8, u16x16, u16x32 => u16); |
| <span class="macro">t!</span>(i16x2, i16x4, i16x8, i16x16, i16x32 => i16); |
| <span class="macro">t!</span>(u32x2, u32x4, u32x8, u32x16 => u32); |
| <span class="macro">t!</span>(i32x2, i32x4, i32x8, i32x16 => i32); |
| <span class="macro">t!</span>(u64x2, u64x4, u64x8 => u64); |
| <span class="macro">t!</span>(i64x2, i64x4, i64x8 => i64); |
| } |
| } |
| |
| <span class="attribute">#[test] |
| #[cfg_attr(miri, ignore)] </span><span class="comment">// Miri is too slow |
| </span><span class="kw">fn </span>test_char() { |
| <span class="kw">let </span><span class="kw-2">mut </span>rng = <span class="kw">crate</span>::test::rng(<span class="number">891</span>); |
| <span class="kw">let </span><span class="kw-2">mut </span>max = core::char::from_u32(<span class="number">0</span>).unwrap(); |
| <span class="kw">for _ in </span><span class="number">0</span>..<span class="number">100 </span>{ |
| <span class="kw">let </span>c = rng.gen_range(<span class="string">'A'</span>..=<span class="string">'Z'</span>); |
| <span class="macro">assert!</span>((<span class="string">'A'</span>..=<span class="string">'Z'</span>).contains(<span class="kw-2">&</span>c)); |
| max = max.max(c); |
| } |
| <span class="macro">assert_eq!</span>(max, <span class="string">'Z'</span>); |
| <span class="kw">let </span>d = Uniform::new( |
| core::char::from_u32(<span class="number">0xD7F0</span>).unwrap(), |
| core::char::from_u32(<span class="number">0xE010</span>).unwrap(), |
| ); |
| <span class="kw">for _ in </span><span class="number">0</span>..<span class="number">100 </span>{ |
| <span class="kw">let </span>c = d.sample(<span class="kw-2">&mut </span>rng); |
| <span class="macro">assert!</span>((c <span class="kw">as </span>u32) < <span class="number">0xD800 </span>|| (c <span class="kw">as </span>u32) > <span class="number">0xDFFF</span>); |
| } |
| } |
| |
| <span class="attribute">#[test] |
| #[cfg_attr(miri, ignore)] </span><span class="comment">// Miri is too slow |
| </span><span class="kw">fn </span>test_floats() { |
| <span class="kw">let </span><span class="kw-2">mut </span>rng = <span class="kw">crate</span>::test::rng(<span class="number">252</span>); |
| <span class="kw">let </span><span class="kw-2">mut </span>zero_rng = StepRng::new(<span class="number">0</span>, <span class="number">0</span>); |
| <span class="kw">let </span><span class="kw-2">mut </span>max_rng = StepRng::new(<span class="number">0xffff_ffff_ffff_ffff</span>, <span class="number">0</span>); |
| <span class="macro">macro_rules! </span>t { |
| (<span class="macro-nonterminal">$ty</span>:ty, <span class="macro-nonterminal">$f_scalar</span>:ident, <span class="macro-nonterminal">$bits_shifted</span>:expr) => {{ |
| <span class="kw">let </span>v: <span class="kw-2">&</span>[(<span class="macro-nonterminal">$f_scalar</span>, <span class="macro-nonterminal">$f_scalar</span>)] = <span class="kw-2">&</span>[ |
| (<span class="number">0.0</span>, <span class="number">100.0</span>), |
| (-<span class="number">1e35</span>, -<span class="number">1e25</span>), |
| (<span class="number">1e-35</span>, <span class="number">1e-25</span>), |
| (-<span class="number">1e35</span>, <span class="number">1e35</span>), |
| (<<span class="macro-nonterminal">$f_scalar</span>>::from_bits(<span class="number">0</span>), <<span class="macro-nonterminal">$f_scalar</span>>::from_bits(<span class="number">3</span>)), |
| (-<<span class="macro-nonterminal">$f_scalar</span>>::from_bits(<span class="number">10</span>), -<<span class="macro-nonterminal">$f_scalar</span>>::from_bits(<span class="number">1</span>)), |
| (-<<span class="macro-nonterminal">$f_scalar</span>>::from_bits(<span class="number">5</span>), <span class="number">0.0</span>), |
| (-<<span class="macro-nonterminal">$f_scalar</span>>::from_bits(<span class="number">7</span>), -<span class="number">0.0</span>), |
| (<span class="number">0.1 </span>* ::core::<span class="macro-nonterminal">$f_scalar::MAX</span>, ::core::<span class="macro-nonterminal">$f_scalar::MAX</span>), |
| (-::core::<span class="macro-nonterminal">$f_scalar::MAX </span>* <span class="number">0.2</span>, ::core::<span class="macro-nonterminal">$f_scalar::MAX </span>* <span class="number">0.7</span>), |
| ]; |
| <span class="kw">for </span><span class="kw-2">&</span>(low_scalar, high_scalar) <span class="kw">in </span>v.iter() { |
| <span class="kw">for </span>lane <span class="kw">in </span><span class="number">0</span>..<<span class="macro-nonterminal">$ty</span>>::lanes() { |
| <span class="kw">let </span>low = <<span class="macro-nonterminal">$ty</span>>::splat(<span class="number">0.0 </span><span class="kw">as </span><span class="macro-nonterminal">$f_scalar</span>).replace(lane, low_scalar); |
| <span class="kw">let </span>high = <<span class="macro-nonterminal">$ty</span>>::splat(<span class="number">1.0 </span><span class="kw">as </span><span class="macro-nonterminal">$f_scalar</span>).replace(lane, high_scalar); |
| <span class="kw">let </span>my_uniform = Uniform::new(low, high); |
| <span class="kw">let </span>my_incl_uniform = Uniform::new_inclusive(low, high); |
| <span class="kw">for _ in </span><span class="number">0</span>..<span class="number">100 </span>{ |
| <span class="kw">let </span>v = rng.sample(my_uniform).extract(lane); |
| <span class="macro">assert!</span>(low_scalar <= v && v < high_scalar); |
| <span class="kw">let </span>v = rng.sample(my_incl_uniform).extract(lane); |
| <span class="macro">assert!</span>(low_scalar <= v && v <= high_scalar); |
| <span class="kw">let </span>v = <<span class="macro-nonterminal">$ty </span><span class="kw">as </span>SampleUniform>::Sampler |
| ::sample_single(low, high, <span class="kw-2">&mut </span>rng).extract(lane); |
| <span class="macro">assert!</span>(low_scalar <= v && v < high_scalar); |
| } |
| |
| <span class="macro">assert_eq!</span>( |
| rng.sample(Uniform::new_inclusive(low, low)).extract(lane), |
| low_scalar |
| ); |
| |
| <span class="macro">assert_eq!</span>(zero_rng.sample(my_uniform).extract(lane), low_scalar); |
| <span class="macro">assert_eq!</span>(zero_rng.sample(my_incl_uniform).extract(lane), low_scalar); |
| <span class="macro">assert_eq!</span>(<<span class="macro-nonterminal">$ty </span><span class="kw">as </span>SampleUniform>::Sampler |
| ::sample_single(low, high, <span class="kw-2">&mut </span>zero_rng) |
| .extract(lane), low_scalar); |
| <span class="macro">assert!</span>(max_rng.sample(my_uniform).extract(lane) < high_scalar); |
| <span class="macro">assert!</span>(max_rng.sample(my_incl_uniform).extract(lane) <= high_scalar); |
| |
| <span class="comment">// Don't run this test for really tiny differences between high and low |
| // since for those rounding might result in selecting high for a very |
| // long time. |
| </span><span class="kw">if </span>(high_scalar - low_scalar) > <span class="number">0.0001 </span>{ |
| <span class="kw">let </span><span class="kw-2">mut </span>lowering_max_rng = StepRng::new( |
| <span class="number">0xffff_ffff_ffff_ffff</span>, |
| (-<span class="number">1i64 </span><< <span class="macro-nonterminal">$bits_shifted</span>) <span class="kw">as </span>u64, |
| ); |
| <span class="macro">assert!</span>( |
| <<span class="macro-nonterminal">$ty </span><span class="kw">as </span>SampleUniform>::Sampler |
| ::sample_single(low, high, <span class="kw-2">&mut </span>lowering_max_rng) |
| .extract(lane) < high_scalar |
| ); |
| } |
| } |
| } |
| |
| <span class="macro">assert_eq!</span>( |
| rng.sample(Uniform::new_inclusive( |
| ::core::<span class="macro-nonterminal">$f_scalar::MAX</span>, |
| ::core::<span class="macro-nonterminal">$f_scalar::MAX |
| </span>)), |
| ::core::<span class="macro-nonterminal">$f_scalar::MAX |
| </span>); |
| <span class="macro">assert_eq!</span>( |
| rng.sample(Uniform::new_inclusive( |
| -::core::<span class="macro-nonterminal">$f_scalar::MAX</span>, |
| -::core::<span class="macro-nonterminal">$f_scalar::MAX |
| </span>)), |
| -::core::<span class="macro-nonterminal">$f_scalar::MAX |
| </span>); |
| }}; |
| } |
| |
| <span class="macro">t!</span>(f32, f32, <span class="number">32 </span>- <span class="number">23</span>); |
| <span class="macro">t!</span>(f64, f64, <span class="number">64 </span>- <span class="number">52</span>); |
| <span class="attribute">#[cfg(feature = <span class="string">"simd_support"</span>)] |
| </span>{ |
| <span class="macro">t!</span>(f32x2, f32, <span class="number">32 </span>- <span class="number">23</span>); |
| <span class="macro">t!</span>(f32x4, f32, <span class="number">32 </span>- <span class="number">23</span>); |
| <span class="macro">t!</span>(f32x8, f32, <span class="number">32 </span>- <span class="number">23</span>); |
| <span class="macro">t!</span>(f32x16, f32, <span class="number">32 </span>- <span class="number">23</span>); |
| <span class="macro">t!</span>(f64x2, f64, <span class="number">64 </span>- <span class="number">52</span>); |
| <span class="macro">t!</span>(f64x4, f64, <span class="number">64 </span>- <span class="number">52</span>); |
| <span class="macro">t!</span>(f64x8, f64, <span class="number">64 </span>- <span class="number">52</span>); |
| } |
| } |
| |
| <span class="attribute">#[test] |
| #[should_panic] |
| </span><span class="kw">fn </span>test_float_overflow() { |
| <span class="kw">let _ </span>= Uniform::from(::core::f64::MIN..::core::f64::MAX); |
| } |
| |
| <span class="attribute">#[test] |
| #[should_panic] |
| </span><span class="kw">fn </span>test_float_overflow_single() { |
| <span class="kw">let </span><span class="kw-2">mut </span>rng = <span class="kw">crate</span>::test::rng(<span class="number">252</span>); |
| rng.gen_range(::core::f64::MIN..::core::f64::MAX); |
| } |
| |
| <span class="attribute">#[test] |
| #[cfg(all( |
| feature = <span class="string">"std"</span>, |
| not(target_arch = <span class="string">"wasm32"</span>), |
| not(target_arch = <span class="string">"asmjs"</span>) |
| ))] |
| </span><span class="kw">fn </span>test_float_assertions() { |
| <span class="kw">use </span><span class="kw">super</span>::SampleUniform; |
| <span class="kw">use </span>std::panic::catch_unwind; |
| <span class="kw">fn </span>range<T: SampleUniform>(low: T, high: T) { |
| <span class="kw">let </span><span class="kw-2">mut </span>rng = <span class="kw">crate</span>::test::rng(<span class="number">253</span>); |
| T::Sampler::sample_single(low, high, <span class="kw-2">&mut </span>rng); |
| } |
| |
| <span class="macro">macro_rules! </span>t { |
| (<span class="macro-nonterminal">$ty</span>:ident, <span class="macro-nonterminal">$f_scalar</span>:ident) => {{ |
| <span class="kw">let </span>v: <span class="kw-2">&</span>[(<span class="macro-nonterminal">$f_scalar</span>, <span class="macro-nonterminal">$f_scalar</span>)] = <span class="kw-2">&</span>[ |
| (::std::<span class="macro-nonterminal">$f_scalar::NAN</span>, <span class="number">0.0</span>), |
| (<span class="number">1.0</span>, ::std::<span class="macro-nonterminal">$f_scalar::NAN</span>), |
| (::std::<span class="macro-nonterminal">$f_scalar::NAN</span>, ::std::<span class="macro-nonterminal">$f_scalar::NAN</span>), |
| (<span class="number">1.0</span>, <span class="number">0.5</span>), |
| (::std::<span class="macro-nonterminal">$f_scalar::MAX</span>, -::std::<span class="macro-nonterminal">$f_scalar::MAX</span>), |
| (::std::<span class="macro-nonterminal">$f_scalar::INFINITY</span>, ::std::<span class="macro-nonterminal">$f_scalar::INFINITY</span>), |
| ( |
| ::std::<span class="macro-nonterminal">$f_scalar::NEG_INFINITY</span>, |
| ::std::<span class="macro-nonterminal">$f_scalar::NEG_INFINITY</span>, |
| ), |
| (::std::<span class="macro-nonterminal">$f_scalar::NEG_INFINITY</span>, <span class="number">5.0</span>), |
| (<span class="number">5.0</span>, ::std::<span class="macro-nonterminal">$f_scalar::INFINITY</span>), |
| (::std::<span class="macro-nonterminal">$f_scalar::NAN</span>, ::std::<span class="macro-nonterminal">$f_scalar::INFINITY</span>), |
| (::std::<span class="macro-nonterminal">$f_scalar::NEG_INFINITY</span>, ::std::<span class="macro-nonterminal">$f_scalar::NAN</span>), |
| (::std::<span class="macro-nonterminal">$f_scalar::NEG_INFINITY</span>, ::std::<span class="macro-nonterminal">$f_scalar::INFINITY</span>), |
| ]; |
| <span class="kw">for </span><span class="kw-2">&</span>(low_scalar, high_scalar) <span class="kw">in </span>v.iter() { |
| <span class="kw">for </span>lane <span class="kw">in </span><span class="number">0</span>..<<span class="macro-nonterminal">$ty</span>>::lanes() { |
| <span class="kw">let </span>low = <<span class="macro-nonterminal">$ty</span>>::splat(<span class="number">0.0 </span><span class="kw">as </span><span class="macro-nonterminal">$f_scalar</span>).replace(lane, low_scalar); |
| <span class="kw">let </span>high = <<span class="macro-nonterminal">$ty</span>>::splat(<span class="number">1.0 </span><span class="kw">as </span><span class="macro-nonterminal">$f_scalar</span>).replace(lane, high_scalar); |
| <span class="macro">assert!</span>(catch_unwind(|| range(low, high)).is_err()); |
| <span class="macro">assert!</span>(catch_unwind(|| Uniform::new(low, high)).is_err()); |
| <span class="macro">assert!</span>(catch_unwind(|| Uniform::new_inclusive(low, high)).is_err()); |
| <span class="macro">assert!</span>(catch_unwind(|| range(low, low)).is_err()); |
| <span class="macro">assert!</span>(catch_unwind(|| Uniform::new(low, low)).is_err()); |
| } |
| } |
| }}; |
| } |
| |
| <span class="macro">t!</span>(f32, f32); |
| <span class="macro">t!</span>(f64, f64); |
| <span class="attribute">#[cfg(feature = <span class="string">"simd_support"</span>)] |
| </span>{ |
| <span class="macro">t!</span>(f32x2, f32); |
| <span class="macro">t!</span>(f32x4, f32); |
| <span class="macro">t!</span>(f32x8, f32); |
| <span class="macro">t!</span>(f32x16, f32); |
| <span class="macro">t!</span>(f64x2, f64); |
| <span class="macro">t!</span>(f64x4, f64); |
| <span class="macro">t!</span>(f64x8, f64); |
| } |
| } |
| |
| |
| <span class="attribute">#[test] |
| #[cfg_attr(miri, ignore)] </span><span class="comment">// Miri is too slow |
| </span><span class="kw">fn </span>test_durations() { |
| <span class="kw">let </span><span class="kw-2">mut </span>rng = <span class="kw">crate</span>::test::rng(<span class="number">253</span>); |
| |
| <span class="kw">let </span>v = <span class="kw-2">&</span>[ |
| (Duration::new(<span class="number">10</span>, <span class="number">50000</span>), Duration::new(<span class="number">100</span>, <span class="number">1234</span>)), |
| (Duration::new(<span class="number">0</span>, <span class="number">100</span>), Duration::new(<span class="number">1</span>, <span class="number">50</span>)), |
| ( |
| Duration::new(<span class="number">0</span>, <span class="number">0</span>), |
| Duration::new(u64::max_value(), <span class="number">999_999_999</span>), |
| ), |
| ]; |
| <span class="kw">for </span><span class="kw-2">&</span>(low, high) <span class="kw">in </span>v.iter() { |
| <span class="kw">let </span>my_uniform = Uniform::new(low, high); |
| <span class="kw">for _ in </span><span class="number">0</span>..<span class="number">1000 </span>{ |
| <span class="kw">let </span>v = rng.sample(my_uniform); |
| <span class="macro">assert!</span>(low <= v && v < high); |
| } |
| } |
| } |
| |
| <span class="attribute">#[test] |
| </span><span class="kw">fn </span>test_custom_uniform() { |
| <span class="kw">use </span><span class="kw">crate</span>::distributions::uniform::{ |
| SampleBorrow, SampleUniform, UniformFloat, UniformSampler, |
| }; |
| <span class="attribute">#[derive(Clone, Copy, PartialEq, PartialOrd)] |
| </span><span class="kw">struct </span>MyF32 { |
| x: f32, |
| } |
| <span class="attribute">#[derive(Clone, Copy, Debug)] |
| </span><span class="kw">struct </span>UniformMyF32(UniformFloat<f32>); |
| <span class="kw">impl </span>UniformSampler <span class="kw">for </span>UniformMyF32 { |
| <span class="kw">type </span>X = MyF32; |
| |
| <span class="kw">fn </span>new<B1, B2>(low: B1, high: B2) -> <span class="self">Self |
| </span><span class="kw">where |
| </span>B1: SampleBorrow<<span class="self">Self</span>::X> + Sized, |
| B2: SampleBorrow<<span class="self">Self</span>::X> + Sized, |
| { |
| UniformMyF32(UniformFloat::<f32>::new(low.borrow().x, high.borrow().x)) |
| } |
| |
| <span class="kw">fn </span>new_inclusive<B1, B2>(low: B1, high: B2) -> <span class="self">Self |
| </span><span class="kw">where |
| </span>B1: SampleBorrow<<span class="self">Self</span>::X> + Sized, |
| B2: SampleBorrow<<span class="self">Self</span>::X> + Sized, |
| { |
| UniformSampler::new(low, high) |
| } |
| |
| <span class="kw">fn </span>sample<R: Rng + <span class="question-mark">?</span>Sized>(<span class="kw-2">&</span><span class="self">self</span>, rng: <span class="kw-2">&mut </span>R) -> <span class="self">Self</span>::X { |
| MyF32 { |
| x: <span class="self">self</span>.<span class="number">0</span>.sample(rng), |
| } |
| } |
| } |
| <span class="kw">impl </span>SampleUniform <span class="kw">for </span>MyF32 { |
| <span class="kw">type </span>Sampler = UniformMyF32; |
| } |
| |
| <span class="kw">let </span>(low, high) = (MyF32 { x: <span class="number">17.0f32 </span>}, MyF32 { x: <span class="number">22.0f32 </span>}); |
| <span class="kw">let </span>uniform = Uniform::new(low, high); |
| <span class="kw">let </span><span class="kw-2">mut </span>rng = <span class="kw">crate</span>::test::rng(<span class="number">804</span>); |
| <span class="kw">for _ in </span><span class="number">0</span>..<span class="number">100 </span>{ |
| <span class="kw">let </span>x: MyF32 = rng.sample(uniform); |
| <span class="macro">assert!</span>(low <= x && x < high); |
| } |
| } |
| |
| <span class="attribute">#[test] |
| </span><span class="kw">fn </span>test_uniform_from_std_range() { |
| <span class="kw">let </span>r = Uniform::from(<span class="number">2u32</span>..<span class="number">7</span>); |
| <span class="macro">assert_eq!</span>(r.<span class="number">0</span>.low, <span class="number">2</span>); |
| <span class="macro">assert_eq!</span>(r.<span class="number">0</span>.range, <span class="number">5</span>); |
| <span class="kw">let </span>r = Uniform::from(<span class="number">2.0f64</span>..<span class="number">7.0</span>); |
| <span class="macro">assert_eq!</span>(r.<span class="number">0</span>.low, <span class="number">2.0</span>); |
| <span class="macro">assert_eq!</span>(r.<span class="number">0</span>.scale, <span class="number">5.0</span>); |
| } |
| |
| <span class="attribute">#[test] |
| </span><span class="kw">fn </span>test_uniform_from_std_range_inclusive() { |
| <span class="kw">let </span>r = Uniform::from(<span class="number">2u32</span>..=<span class="number">6</span>); |
| <span class="macro">assert_eq!</span>(r.<span class="number">0</span>.low, <span class="number">2</span>); |
| <span class="macro">assert_eq!</span>(r.<span class="number">0</span>.range, <span class="number">5</span>); |
| <span class="kw">let </span>r = Uniform::from(<span class="number">2.0f64</span>..=<span class="number">7.0</span>); |
| <span class="macro">assert_eq!</span>(r.<span class="number">0</span>.low, <span class="number">2.0</span>); |
| <span class="macro">assert!</span>(r.<span class="number">0</span>.scale > <span class="number">5.0</span>); |
| <span class="macro">assert!</span>(r.<span class="number">0</span>.scale < <span class="number">5.0 </span>+ <span class="number">1e-14</span>); |
| } |
| |
| <span class="attribute">#[test] |
| </span><span class="kw">fn </span>value_stability() { |
| <span class="kw">fn </span>test_samples<T: SampleUniform + Copy + core::fmt::Debug + PartialEq>( |
| lb: T, ub: T, expected_single: <span class="kw-2">&</span>[T], expected_multiple: <span class="kw-2">&</span>[T], |
| ) <span class="kw">where </span>Uniform<T>: Distribution<T> { |
| <span class="kw">let </span><span class="kw-2">mut </span>rng = <span class="kw">crate</span>::test::rng(<span class="number">897</span>); |
| <span class="kw">let </span><span class="kw-2">mut </span>buf = [lb; <span class="number">3</span>]; |
| |
| <span class="kw">for </span>x <span class="kw">in </span><span class="kw-2">&mut </span>buf { |
| <span class="kw-2">*</span>x = T::Sampler::sample_single(lb, ub, <span class="kw-2">&mut </span>rng); |
| } |
| <span class="macro">assert_eq!</span>(<span class="kw-2">&</span>buf, expected_single); |
| |
| <span class="kw">let </span>distr = Uniform::new(lb, ub); |
| <span class="kw">for </span>x <span class="kw">in </span><span class="kw-2">&mut </span>buf { |
| <span class="kw-2">*</span>x = rng.sample(<span class="kw-2">&</span>distr); |
| } |
| <span class="macro">assert_eq!</span>(<span class="kw-2">&</span>buf, expected_multiple); |
| } |
| |
| <span class="comment">// We test on a sub-set of types; possibly we should do more. |
| // TODO: SIMD types |
| |
| </span>test_samples(<span class="number">11u8</span>, <span class="number">219</span>, <span class="kw-2">&</span>[<span class="number">17</span>, <span class="number">66</span>, <span class="number">214</span>], <span class="kw-2">&</span>[<span class="number">181</span>, <span class="number">93</span>, <span class="number">165</span>]); |
| test_samples(<span class="number">11u32</span>, <span class="number">219</span>, <span class="kw-2">&</span>[<span class="number">17</span>, <span class="number">66</span>, <span class="number">214</span>], <span class="kw-2">&</span>[<span class="number">181</span>, <span class="number">93</span>, <span class="number">165</span>]); |
| |
| test_samples(<span class="number">0f32</span>, <span class="number">1e-2f32</span>, <span class="kw-2">&</span>[<span class="number">0.0003070104</span>, <span class="number">0.0026630748</span>, <span class="number">0.00979833</span>], <span class="kw-2">&</span>[ |
| <span class="number">0.008194133</span>, |
| <span class="number">0.00398172</span>, |
| <span class="number">0.007428536</span>, |
| ]); |
| test_samples( |
| -<span class="number">1e10f64</span>, |
| <span class="number">1e10f64</span>, |
| <span class="kw-2">&</span>[-<span class="number">4673848682.871551</span>, <span class="number">6388267422.932352</span>, <span class="number">4857075081.198343</span>], |
| <span class="kw-2">&</span>[<span class="number">1173375212.1808167</span>, <span class="number">1917642852.109581</span>, <span class="number">2365076174.3153973</span>], |
| ); |
| |
| test_samples( |
| Duration::new(<span class="number">2</span>, <span class="number">0</span>), |
| Duration::new(<span class="number">4</span>, <span class="number">0</span>), |
| <span class="kw-2">&</span>[ |
| Duration::new(<span class="number">2</span>, <span class="number">532615131</span>), |
| Duration::new(<span class="number">3</span>, <span class="number">638826742</span>), |
| Duration::new(<span class="number">3</span>, <span class="number">485707508</span>), |
| ], |
| <span class="kw-2">&</span>[ |
| Duration::new(<span class="number">3</span>, <span class="number">117337521</span>), |
| Duration::new(<span class="number">3</span>, <span class="number">191764285</span>), |
| Duration::new(<span class="number">3</span>, <span class="number">236507617</span>), |
| ], |
| ); |
| } |
| |
| <span class="attribute">#[test] |
| </span><span class="kw">fn </span>uniform_distributions_can_be_compared() { |
| <span class="macro">assert_eq!</span>(Uniform::new(<span class="number">1.0</span>, <span class="number">2.0</span>), Uniform::new(<span class="number">1.0</span>, <span class="number">2.0</span>)); |
| |
| <span class="comment">// To cover UniformInt |
| </span><span class="macro">assert_eq!</span>(Uniform::new(<span class="number">1 </span><span class="kw">as </span>u32, <span class="number">2 </span><span class="kw">as </span>u32), Uniform::new(<span class="number">1 </span><span class="kw">as </span>u32, <span class="number">2 </span><span class="kw">as </span>u32)); |
| } |
| } |
| </code></pre></div> |
| </section></div></main><div id="rustdoc-vars" data-root-path="../../../" data-current-crate="rand" data-themes="ayu,dark,light" data-resource-suffix="" data-rustdoc-version="1.66.0-nightly (5c8bff74b 2022-10-21)" ></div></body></html> |