| <!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/rulinalg-309246e5a94bf5cf/1ed8b93/src/matrix/base/mod.rs`."><meta name="keywords" content="rust, rustlang, rust-lang"><title>mod.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="alternate icon" type="image/png" href="../../../../favicon-16x16.png"><link rel="alternate icon" type="image/png" href="../../../../favicon-32x32.png"><link rel="icon" type="image/svg+xml" href="../../../../favicon.svg"></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="../../../../rulinalg/index.html"><div class="logo-container"><img class="rust-logo" src="../../../../rust-logo.svg" alt="logo"></div></a></nav><main><div class="width-limiter"><nav class="sub"><a class="sub-logo-container" href="../../../../rulinalg/index.html"><img class="rust-logo" src="../../../../rust-logo.svg" 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> |
| <span id="1659">1659</span> |
| <span id="1660">1660</span> |
| <span id="1661">1661</span> |
| <span id="1662">1662</span> |
| <span id="1663">1663</span> |
| <span id="1664">1664</span> |
| <span id="1665">1665</span> |
| <span id="1666">1666</span> |
| <span id="1667">1667</span> |
| <span id="1668">1668</span> |
| <span id="1669">1669</span> |
| <span id="1670">1670</span> |
| <span id="1671">1671</span> |
| <span id="1672">1672</span> |
| <span id="1673">1673</span> |
| <span id="1674">1674</span> |
| <span id="1675">1675</span> |
| <span id="1676">1676</span> |
| <span id="1677">1677</span> |
| <span id="1678">1678</span> |
| <span id="1679">1679</span> |
| <span id="1680">1680</span> |
| <span id="1681">1681</span> |
| <span id="1682">1682</span> |
| <span id="1683">1683</span> |
| <span id="1684">1684</span> |
| <span id="1685">1685</span> |
| <span id="1686">1686</span> |
| <span id="1687">1687</span> |
| <span id="1688">1688</span> |
| <span id="1689">1689</span> |
| <span id="1690">1690</span> |
| <span id="1691">1691</span> |
| <span id="1692">1692</span> |
| <span id="1693">1693</span> |
| <span id="1694">1694</span> |
| <span id="1695">1695</span> |
| <span id="1696">1696</span> |
| <span id="1697">1697</span> |
| <span id="1698">1698</span> |
| <span id="1699">1699</span> |
| <span id="1700">1700</span> |
| <span id="1701">1701</span> |
| <span id="1702">1702</span> |
| <span id="1703">1703</span> |
| <span id="1704">1704</span> |
| <span id="1705">1705</span> |
| <span id="1706">1706</span> |
| <span id="1707">1707</span> |
| <span id="1708">1708</span> |
| <span id="1709">1709</span> |
| <span id="1710">1710</span> |
| <span id="1711">1711</span> |
| <span id="1712">1712</span> |
| <span id="1713">1713</span> |
| <span id="1714">1714</span> |
| <span id="1715">1715</span> |
| <span id="1716">1716</span> |
| <span id="1717">1717</span> |
| <span id="1718">1718</span> |
| <span id="1719">1719</span> |
| <span id="1720">1720</span> |
| <span id="1721">1721</span> |
| <span id="1722">1722</span> |
| <span id="1723">1723</span> |
| <span id="1724">1724</span> |
| <span id="1725">1725</span> |
| <span id="1726">1726</span> |
| <span id="1727">1727</span> |
| <span id="1728">1728</span> |
| <span id="1729">1729</span> |
| <span id="1730">1730</span> |
| <span id="1731">1731</span> |
| </pre><pre class="rust"><code><span class="doccomment">//! Traits for matrices operations. |
| //! |
| //! These traits defines operations for structs representing matrices arranged in row-major order. |
| //! |
| //! Implementations are provided for |
| //! - `Matrix`: an owned matrix |
| //! - `MatrixSlice`: a borrowed immutable block of `Matrix` |
| //! - `MatrixSliceMut`: a borrowed mutable block of `Matrix` |
| //! |
| //! ``` |
| //! use rulinalg::matrix::{Matrix, BaseMatrix}; |
| //! |
| //! let a = Matrix::new(3,3, (0..9).collect::<Vec<usize>>()); |
| //! |
| //! // Manually create our slice - [[4,5],[7,8]]. |
| //! let mat_slice = a.sub_slice([0,1], 3, 2); |
| //! |
| //! // We can perform arithmetic with mixing owned and borrowed versions |
| //! let _new_mat = &mat_slice.transpose() * &a; |
| //! ``` |
| |
| </span><span class="kw">use </span>matrix::{Matrix, MatrixSlice, MatrixSliceMut}; |
| <span class="kw">use </span>matrix::{Cols, ColsMut, Row, RowMut, Column, ColumnMut, Rows, RowsMut, Axes}; |
| <span class="kw">use </span>matrix::{DiagOffset, Diagonal, DiagonalMut}; |
| <span class="kw">use </span>matrix::{back_substitution, forward_substitution}; |
| <span class="kw">use </span>matrix::{SliceIter, SliceIterMut}; |
| <span class="kw">use </span>norm::{MatrixNorm, MatrixMetric}; |
| <span class="kw">use </span>vector::Vector; |
| <span class="kw">use </span>utils; |
| <span class="kw">use </span>libnum::{Zero, Float}; |
| <span class="kw">use </span>error::Error; |
| |
| <span class="kw">use </span>std::any::Any; |
| <span class="kw">use </span>std::cmp::min; |
| <span class="kw">use </span>std::marker::PhantomData; |
| <span class="kw">use </span>std::mem; |
| <span class="kw">use </span>std::ops::{Add, Mul, Div}; |
| <span class="kw">use </span>std::ptr; |
| <span class="kw">use </span>std::slice; |
| |
| <span class="kw">mod </span>impl_base; |
| |
| <span class="doccomment">/// Trait for immutable matrix structs. |
| </span><span class="kw">pub trait </span>BaseMatrix<T>: Sized { |
| <span class="doccomment">/// Rows in the matrix. |
| </span><span class="kw">fn </span>rows(<span class="kw-2">&</span><span class="self">self</span>) -> usize; |
| |
| <span class="doccomment">/// Columns in the matrix. |
| </span><span class="kw">fn </span>cols(<span class="kw-2">&</span><span class="self">self</span>) -> usize; |
| |
| <span class="doccomment">/// Row stride in the matrix. |
| </span><span class="kw">fn </span>row_stride(<span class="kw-2">&</span><span class="self">self</span>) -> usize; |
| |
| <span class="doccomment">/// Returns true if the matrix contais no elements |
| </span><span class="kw">fn </span>is_empty(<span class="kw-2">&</span><span class="self">self</span>) -> bool { |
| <span class="self">self</span>.rows() == <span class="number">0 </span>|| <span class="self">self</span>.cols() == <span class="number">0 |
| </span>} |
| |
| <span class="doccomment">/// Top left index of the matrix. |
| </span><span class="kw">fn </span>as_ptr(<span class="kw-2">&</span><span class="self">self</span>) -> <span class="kw-2">*const </span>T; |
| |
| <span class="doccomment">/// Returns a `MatrixSlice` over the whole matrix. |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// use rulinalg::matrix::{Matrix, BaseMatrix}; |
| /// |
| /// let a = Matrix::new(3, 3, vec![2.0; 9]); |
| /// let b = a.as_slice(); |
| /// ``` |
| </span><span class="kw">fn </span>as_slice(<span class="kw-2">&</span><span class="self">self</span>) -> MatrixSlice<T> { |
| <span class="kw">unsafe </span>{ |
| MatrixSlice::from_raw_parts(<span class="self">self</span>.as_ptr(), <span class="self">self</span>.rows(), <span class="self">self</span>.cols(), <span class="self">self</span>.row_stride()) |
| } |
| } |
| |
| <span class="doccomment">/// Get a reference to an element in the matrix without bounds checking. |
| </span><span class="kw">unsafe fn </span>get_unchecked(<span class="kw-2">&</span><span class="self">self</span>, index: [usize; <span class="number">2</span>]) -> <span class="kw-2">&</span>T { |
| <span class="kw-2">&*</span>(<span class="self">self</span>.as_ptr().offset((index[<span class="number">0</span>] * <span class="self">self</span>.row_stride() + index[<span class="number">1</span>]) <span class="kw">as </span>isize)) |
| } |
| |
| <span class="doccomment">/// Get a reference to an element in the matrix. |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// # #[macro_use] extern crate rulinalg; fn main() { |
| /// use rulinalg::matrix::{Matrix, BaseMatrix}; |
| /// |
| /// let mat = matrix![0, 1; |
| /// 3, 4; |
| /// 6, 7]; |
| /// |
| /// assert_eq!(mat.get([0, 2]), None); |
| /// assert_eq!(mat.get([3, 0]), None); |
| /// |
| /// assert_eq!( *mat.get([0, 0]).unwrap(), 0) |
| /// # } |
| /// ``` |
| </span><span class="kw">fn </span>get(<span class="kw-2">&</span><span class="self">self</span>, index: [usize; <span class="number">2</span>]) -> <span class="prelude-ty">Option</span><<span class="kw-2">&</span>T> { |
| <span class="kw">let </span>row_ind = index[<span class="number">0</span>]; |
| <span class="kw">let </span>col_ind = index[<span class="number">1</span>]; |
| |
| <span class="kw">if </span>row_ind >= <span class="self">self</span>.rows() || col_ind >= <span class="self">self</span>.cols() { |
| <span class="prelude-val">None |
| </span>} <span class="kw">else </span>{ |
| <span class="kw">unsafe </span>{ <span class="prelude-val">Some</span>(<span class="self">self</span>.get_unchecked(index)) } |
| } |
| } |
| |
| <span class="doccomment">/// Returns the column of a matrix at the given index. |
| /// `None` if the index is out of bounds. |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// # #[macro_use] extern crate rulinalg; fn main() { |
| /// use rulinalg::matrix::{Matrix, BaseMatrix}; |
| /// |
| /// let mat = matrix![0, 1, 2; |
| /// 3, 4, 5; |
| /// 6, 7, 8]; |
| /// let col = mat.col(1); |
| /// let expected = matrix![1usize; 4; 7]; |
| /// assert_matrix_eq!(*col, expected); |
| /// # } |
| /// ``` |
| /// |
| /// # Panics |
| /// |
| /// Will panic if the column index is out of bounds. |
| </span><span class="kw">fn </span>col(<span class="kw-2">&</span><span class="self">self</span>, index: usize) -> Column<T> { |
| <span class="kw">if </span>index < <span class="self">self</span>.cols() { |
| <span class="kw">unsafe </span>{ <span class="self">self</span>.col_unchecked(index) } |
| } <span class="kw">else </span>{ |
| <span class="macro">panic!</span>(<span class="string">"Column index out of bounds."</span>) |
| } |
| } |
| |
| <span class="doccomment">/// Returns the column of a matrix at the given |
| /// index without doing a bounds check. |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// # #[macro_use] extern crate rulinalg; fn main() { |
| /// use rulinalg::matrix::{Matrix, BaseMatrix}; |
| /// |
| /// let mat = matrix![0, 1, 2; |
| /// 3, 4, 5; |
| /// 6, 7, 8]; |
| /// let col = unsafe { mat.col_unchecked(2) }; |
| /// let expected = matrix![2usize; 5; 8]; |
| /// assert_matrix_eq!(*col, expected); |
| /// # } |
| /// ``` |
| </span><span class="kw">unsafe fn </span>col_unchecked(<span class="kw-2">&</span><span class="self">self</span>, index: usize) -> Column<T> { |
| <span class="kw">let </span>ptr = <span class="self">self</span>.as_ptr().offset(index <span class="kw">as </span>isize); |
| Column { col: MatrixSlice::from_raw_parts(ptr, <span class="self">self</span>.rows(), <span class="number">1</span>, <span class="self">self</span>.row_stride()) } |
| } |
| |
| <span class="doccomment">/// Returns the row of a matrix at the given index. |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// # #[macro_use] extern crate rulinalg; fn main() { |
| /// use rulinalg::matrix::{Matrix, BaseMatrix}; |
| /// |
| /// let mat = matrix![0, 1, 2; |
| /// 3, 4, 5; |
| /// 6, 7, 8]; |
| /// let row = mat.row(1); |
| /// let expected = matrix![3usize, 4, 5]; |
| /// assert_matrix_eq!(*row, expected); |
| /// # } |
| /// ``` |
| /// |
| /// # Panics |
| /// |
| /// Will panic if the row index is out of bounds. |
| </span><span class="kw">fn </span>row(<span class="kw-2">&</span><span class="self">self</span>, index: usize) -> Row<T> { |
| <span class="kw">if </span>index < <span class="self">self</span>.rows() { |
| <span class="kw">unsafe </span>{ <span class="self">self</span>.row_unchecked(index) } |
| } <span class="kw">else </span>{ |
| <span class="macro">panic!</span>(<span class="string">"Row index out of bounds."</span>) |
| } |
| } |
| |
| <span class="doccomment">/// Returns the row of a matrix at the given index without doing unbounds checking |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// # #[macro_use] extern crate rulinalg; fn main() { |
| /// use rulinalg::matrix::{Matrix, BaseMatrix}; |
| /// |
| /// let mat = matrix![0, 1, 2; |
| /// 3, 4, 5; |
| /// 6, 7, 8]; |
| /// let row = unsafe { mat.row_unchecked(2) }; |
| /// let expected = matrix![6usize, 7, 8]; |
| /// assert_matrix_eq!(*row, expected); |
| /// # } |
| /// ``` |
| </span><span class="kw">unsafe fn </span>row_unchecked(<span class="kw-2">&</span><span class="self">self</span>, index: usize) -> Row<T> { |
| <span class="kw">let </span>ptr = <span class="self">self</span>.as_ptr().offset((<span class="self">self</span>.row_stride() * index) <span class="kw">as </span>isize); |
| Row { row: MatrixSlice::from_raw_parts(ptr, <span class="number">1</span>, <span class="self">self</span>.cols(), <span class="self">self</span>.row_stride()) } |
| } |
| |
| <span class="doccomment">/// Returns an iterator over the matrix data. |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// # #[macro_use] extern crate rulinalg; fn main() { |
| /// use rulinalg::matrix::{Matrix, BaseMatrix}; |
| /// |
| /// let mat = matrix![0, 1, 2; |
| /// 3, 4, 5; |
| /// 6, 7, 8]; |
| /// let slice = mat.sub_slice([1, 1], 2, 2); |
| /// |
| /// let slice_data = slice.iter().map(|v| *v).collect::<Vec<usize>>(); |
| /// assert_eq!(slice_data, vec![4, 5, 7, 8]); |
| /// # } |
| /// ``` |
| </span><span class="kw">fn </span>iter<<span class="lifetime">'a</span>>(<span class="kw-2">&</span><span class="self">self</span>) -> SliceIter<<span class="lifetime">'a</span>, T> |
| <span class="kw">where </span>T: <span class="lifetime">'a |
| </span>{ |
| SliceIter { |
| slice_start: <span class="self">self</span>.as_ptr(), |
| row_pos: <span class="number">0</span>, |
| col_pos: <span class="number">0</span>, |
| slice_rows: <span class="self">self</span>.rows(), |
| slice_cols: <span class="self">self</span>.cols(), |
| row_stride: <span class="self">self</span>.row_stride(), |
| _marker: PhantomData::<<span class="kw-2">&</span>T>, |
| } |
| } |
| |
| <span class="doccomment">/// Iterate over the columns of the matrix. |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// # #[macro_use] extern crate rulinalg; fn main() { |
| /// use rulinalg::matrix::{Matrix, BaseMatrix}; |
| /// |
| /// let a = matrix![0, 1; |
| /// 2, 3; |
| /// 4, 5]; |
| /// |
| /// let mut iter = a.col_iter(); |
| /// |
| /// assert_matrix_eq!(*iter.next().unwrap(), matrix![ 0; 2; 4 ]); |
| /// assert_matrix_eq!(*iter.next().unwrap(), matrix![ 1; 3; 5 ]); |
| /// assert!(iter.next().is_none()); |
| /// # } |
| /// ``` |
| </span><span class="kw">fn </span>col_iter(<span class="kw-2">&</span><span class="self">self</span>) -> Cols<T> { |
| Cols { |
| _marker: PhantomData::<<span class="kw-2">&</span>T>, |
| col_pos: <span class="number">0</span>, |
| row_stride: <span class="self">self</span>.row_stride() <span class="kw">as </span>isize, |
| slice_cols: <span class="self">self</span>.cols(), |
| slice_rows: <span class="self">self</span>.rows(), |
| slice_start: <span class="self">self</span>.as_ptr(), |
| } |
| } |
| |
| <span class="doccomment">/// Iterate over the rows of the matrix. |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// # #[macro_use] extern crate rulinalg; fn main() { |
| /// use rulinalg::matrix::{Matrix, BaseMatrix}; |
| /// let a = matrix![0, 1, 2; |
| /// 3, 4, 5; |
| /// 6, 7, 8]; |
| /// |
| /// let mut iter = a.row_iter(); |
| /// |
| /// assert_matrix_eq!(*iter.next().unwrap(), matrix![ 0, 1, 2 ]); |
| /// assert_matrix_eq!(*iter.next().unwrap(), matrix![ 3, 4, 5 ]); |
| /// assert_matrix_eq!(*iter.next().unwrap(), matrix![ 6, 7, 8 ]); |
| /// assert!(iter.next().is_none()); |
| /// # } |
| /// ``` |
| </span><span class="kw">fn </span>row_iter(<span class="kw-2">&</span><span class="self">self</span>) -> Rows<T> { |
| Rows { |
| slice_start: <span class="self">self</span>.as_ptr(), |
| row_pos: <span class="number">0</span>, |
| slice_rows: <span class="self">self</span>.rows(), |
| slice_cols: <span class="self">self</span>.cols(), |
| row_stride: <span class="self">self</span>.row_stride() <span class="kw">as </span>isize, |
| _marker: PhantomData::<<span class="kw-2">&</span>T>, |
| } |
| } |
| |
| <span class="doccomment">/// Iterate over diagonal entries |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// # #[macro_use] extern crate rulinalg; |
| /// |
| /// # fn main() { |
| /// use rulinalg::matrix::{DiagOffset, Matrix, BaseMatrix}; |
| /// |
| /// let a = matrix![0, 1, 2; |
| /// 3, 4, 5; |
| /// 6, 7, 8]; |
| /// // Print super diag [1, 5] |
| /// for d in a.diag_iter(DiagOffset::Above(1)) { |
| /// println!("{}", d); |
| /// } |
| /// |
| /// // Print sub diag [3, 7] |
| /// // Equivalent to `diag_iter(DiagOffset::Below(1))` |
| /// for d in a.diag_iter(DiagOffset::from(-1)) { |
| /// println!("{}", d); |
| /// } |
| /// # } |
| /// ``` |
| /// |
| /// # Panics |
| /// |
| /// If using an `Above` or `Below` offset which is |
| /// out-of-bounds this function will panic. |
| /// |
| /// This function will never panic if the `Main` diagonal |
| /// offset is used. |
| </span><span class="kw">fn </span>diag_iter(<span class="kw-2">&</span><span class="self">self</span>, k: DiagOffset) -> Diagonal<T, <span class="self">Self</span>> { |
| <span class="kw">let </span>(diag_len, diag_start) = <span class="kw">match </span>k.into() { |
| DiagOffset::Main => (min(<span class="self">self</span>.rows(), <span class="self">self</span>.cols()), <span class="number">0</span>), |
| DiagOffset::Above(m) => { |
| <span class="macro">assert!</span>(m < <span class="self">self</span>.cols(), |
| <span class="string">"Offset diagonal is not within matrix dimensions."</span>); |
| (min(<span class="self">self</span>.rows(), <span class="self">self</span>.cols() - m), m) |
| } |
| DiagOffset::Below(m) => { |
| <span class="macro">assert!</span>(m < <span class="self">self</span>.rows(), |
| <span class="string">"Offset diagonal is not within matrix dimensions."</span>); |
| (min(<span class="self">self</span>.rows() - m, <span class="self">self</span>.cols()), m * <span class="self">self</span>.row_stride()) |
| } |
| }; |
| |
| Diagonal { |
| matrix: <span class="self">self</span>, |
| diag_pos: diag_start, |
| diag_end: diag_start + diag_len.saturating_sub(<span class="number">1</span>) * <span class="self">self</span>.row_stride() + diag_len, |
| _marker: PhantomData::<<span class="kw-2">&</span>T>, |
| } |
| } |
| |
| <span class="doccomment">/// The sum of the rows of the matrix. |
| /// |
| /// Returns a Vector equal to the sums of elements over the matrices rows. |
| /// |
| /// Note that the resulting vector is identical to the sums of |
| /// elements along each column of the matrix. |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// # #[macro_use] extern crate rulinalg; fn main() { |
| /// use rulinalg::matrix::{Matrix, BaseMatrix}; |
| /// |
| /// let a = matrix![1.0, 2.0; |
| /// 3.0, 4.0]; |
| /// |
| /// let c = a.sum_rows(); |
| /// assert_eq!(c, vector![4.0, 6.0]); |
| /// # } |
| /// ``` |
| </span><span class="kw">fn </span>sum_rows(<span class="kw-2">&</span><span class="self">self</span>) -> Vector<T> |
| <span class="kw">where </span>T: Copy + Zero + Add<T, Output = T> |
| { |
| <span class="kw">let </span><span class="kw-2">mut </span>sum_rows = <span class="macro">vec!</span>[T::zero(); <span class="self">self</span>.cols()]; |
| <span class="kw">for </span>row <span class="kw">in </span><span class="self">self</span>.row_iter() { |
| utils::in_place_vec_bin_op(<span class="kw-2">&mut </span>sum_rows, row.raw_slice(), |sum, <span class="kw-2">&</span>r| <span class="kw-2">*</span>sum = <span class="kw-2">*</span>sum + r); |
| } |
| Vector::new(sum_rows) |
| } |
| |
| <span class="doccomment">/// The sum of the columns of the matrix. |
| /// |
| /// Returns a Vector equal to the sums of elements over the matrices columns. |
| /// |
| /// Note that the resulting vector is identical to the sums of |
| /// elements along each row of the matrix. |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// # #[macro_use] extern crate rulinalg; fn main() { |
| /// use rulinalg::matrix::{Matrix, BaseMatrix}; |
| /// |
| /// let a = matrix![1.0, 2.0; |
| /// 3.0, 4.0]; |
| /// |
| /// let c = a.sum_cols(); |
| /// assert_eq!(c, vector![3.0, 7.0]); |
| /// # } |
| /// ``` |
| </span><span class="kw">fn </span>sum_cols(<span class="kw-2">&</span><span class="self">self</span>) -> Vector<T> |
| <span class="kw">where </span>T: Copy + Zero + Add<T, Output = T> |
| { |
| <span class="kw">let </span><span class="kw-2">mut </span>col_sum = Vec::with_capacity(<span class="self">self</span>.rows()); |
| col_sum.extend(<span class="self">self</span>.row_iter().map(|row| utils::unrolled_sum(row.raw_slice()))); |
| Vector::new(col_sum) |
| } |
| |
| <span class="doccomment">/// Compute given matrix norm for matrix. |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// # #[macro_use] extern crate rulinalg; fn main() { |
| /// use rulinalg::matrix::BaseMatrix; |
| /// use rulinalg::norm::Euclidean; |
| /// |
| /// let a = matrix![3.0, 4.0]; |
| /// let c = a.norm(Euclidean); |
| /// |
| /// assert_eq!(c, 5.0); |
| /// # } |
| /// ``` |
| </span><span class="kw">fn </span>norm<N: MatrixNorm<T, <span class="self">Self</span>>>(<span class="kw-2">&</span><span class="self">self</span>, norm: N) -> T |
| <span class="kw">where </span>T: Float |
| { |
| norm.norm(<span class="self">self</span>) |
| } |
| |
| <span class="doccomment">/// Compute the metric distance between two matrices. |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// # #[macro_use] extern crate rulinalg; fn main() { |
| /// use rulinalg::matrix::BaseMatrix; |
| /// use rulinalg::norm::Euclidean; |
| /// |
| /// let a = matrix![3.0, 4.0; |
| /// 1.0, 2.0]; |
| /// let b = matrix![2.0, 5.0; |
| /// 0.0, 3.0]; |
| /// |
| /// // Compute the square root of the sum of |
| /// // elementwise squared-differences |
| /// let c = a.metric(&b, Euclidean); |
| /// |
| /// assert_eq!(c, 2.0); |
| /// # } |
| /// ``` |
| </span><span class="kw">fn </span>metric<<span class="lifetime">'a</span>, <span class="lifetime">'b</span>, B, M>(<span class="kw-2">&</span><span class="lifetime">'a </span><span class="self">self</span>, mat: <span class="kw-2">&</span><span class="lifetime">'b </span>B, metric: M) -> T |
| <span class="kw">where </span>B: <span class="lifetime">'b </span>+ BaseMatrix<T>, |
| M: MatrixMetric<<span class="lifetime">'a</span>, <span class="lifetime">'b</span>, T, <span class="self">Self</span>, B> |
| { |
| metric.metric(<span class="self">self</span>, mat) |
| } |
| |
| <span class="doccomment">/// The sum of all elements in the matrix |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// # #[macro_use] extern crate rulinalg; fn main() { |
| /// use rulinalg::matrix::BaseMatrix; |
| /// |
| /// let a = matrix![1.0, 2.0; |
| /// 3.0, 4.0]; |
| /// |
| /// let c = a.sum(); |
| /// assert_eq!(c, 10.0); |
| /// # } |
| /// ``` |
| </span><span class="kw">fn </span>sum(<span class="kw-2">&</span><span class="self">self</span>) -> T |
| <span class="kw">where </span>T: Copy + Zero + Add<T, Output = T> |
| { |
| <span class="self">self</span>.row_iter() |
| .fold(T::zero(), |
| |sum, row| sum + utils::unrolled_sum(row.raw_slice())) |
| } |
| |
| <span class="doccomment">/// The min of the specified axis of the matrix. |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// # #[macro_use] extern crate rulinalg; fn main() { |
| /// use rulinalg::matrix::{Matrix, BaseMatrix, Axes}; |
| /// |
| /// let a = matrix![1.0, 2.0; |
| /// 3.0, 4.0]; |
| /// |
| /// let cmin = a.min(Axes::Col); |
| /// assert_eq!(cmin, vector![1.0, 3.0]); |
| /// |
| /// let rmin = a.min(Axes::Row); |
| /// assert_eq!(rmin, vector![1.0, 2.0]); |
| /// # } |
| /// ``` |
| </span><span class="kw">fn </span>min(<span class="kw-2">&</span><span class="self">self</span>, axis: Axes) -> Vector<T> |
| <span class="kw">where </span>T: Copy + PartialOrd |
| { |
| <span class="kw">match </span>axis { |
| Axes::Col => { |
| <span class="kw">let </span><span class="kw-2">mut </span>mins: Vec<T> = Vec::with_capacity(<span class="self">self</span>.rows()); |
| <span class="kw">for </span>row <span class="kw">in </span><span class="self">self</span>.row_iter() { |
| <span class="kw">let </span>min = row.iter() |
| .skip(<span class="number">1</span>) |
| .fold(row[<span class="number">0</span>], |m, <span class="kw-2">&</span>v| <span class="kw">if </span>v < m { v } <span class="kw">else </span>{ m }); |
| mins.push(min); |
| } |
| Vector::new(mins) |
| } |
| Axes::Row => { |
| <span class="kw">let </span><span class="kw-2">mut </span>mins: Vec<T> = <span class="self">self</span>.row(<span class="number">0</span>).raw_slice().into(); |
| <span class="kw">for </span>row <span class="kw">in </span><span class="self">self</span>.row_iter().skip(<span class="number">1</span>) { |
| utils::in_place_vec_bin_op(<span class="kw-2">&mut </span>mins, row.raw_slice(), |min, <span class="kw-2">&</span>r| <span class="kw">if </span>r < <span class="kw-2">*</span>min { |
| <span class="kw-2">*</span>min = r; |
| }); |
| } |
| Vector::new(mins) |
| } |
| } |
| } |
| |
| <span class="doccomment">/// The max of the specified axis of the matrix. |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// # #[macro_use] extern crate rulinalg; fn main() { |
| /// use rulinalg::matrix::{BaseMatrix, Axes}; |
| /// |
| /// let a = matrix![1.0, 2.0; |
| /// 3.0, 4.0]; |
| /// |
| /// let cmax = a.max(Axes::Col); |
| /// assert_eq!(cmax, vector![2.0, 4.0]); |
| /// |
| /// let rmax = a.max(Axes::Row); |
| /// assert_eq!(rmax, vector![3.0, 4.0]); |
| /// # } |
| /// ``` |
| </span><span class="kw">fn </span>max(<span class="kw-2">&</span><span class="self">self</span>, axis: Axes) -> Vector<T> |
| <span class="kw">where </span>T: Copy + PartialOrd |
| { |
| <span class="kw">match </span>axis { |
| Axes::Col => { |
| <span class="kw">let </span><span class="kw-2">mut </span>maxs: Vec<T> = Vec::with_capacity(<span class="self">self</span>.rows()); |
| <span class="kw">for </span>row <span class="kw">in </span><span class="self">self</span>.row_iter() { |
| <span class="kw">let </span>max = row.iter() |
| .skip(<span class="number">1</span>) |
| .fold(row[<span class="number">0</span>], |m, <span class="kw-2">&</span>v| <span class="kw">if </span>v > m { v } <span class="kw">else </span>{ m }); |
| maxs.push(max); |
| } |
| Vector::new(maxs) |
| } |
| Axes::Row => { |
| <span class="kw">let </span><span class="kw-2">mut </span>maxs: Vec<T> = <span class="self">self</span>.row(<span class="number">0</span>).raw_slice().into(); |
| <span class="kw">for </span>row <span class="kw">in </span><span class="self">self</span>.row_iter().skip(<span class="number">1</span>) { |
| utils::in_place_vec_bin_op(<span class="kw-2">&mut </span>maxs, row.raw_slice(), |max, <span class="kw-2">&</span>r| <span class="kw">if </span>r > <span class="kw-2">*</span>max { |
| <span class="kw-2">*</span>max = r; |
| }); |
| } |
| Vector::new(maxs) |
| } |
| } |
| } |
| |
| <span class="doccomment">/// Convert the matrix struct into a owned Matrix. |
| </span><span class="kw">fn </span>into_matrix(<span class="self">self</span>) -> Matrix<T> |
| <span class="kw">where </span>T: Copy |
| { |
| <span class="self">self</span>.row_iter().collect() |
| } |
| |
| <span class="doccomment">/// Select rows from matrix |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// use rulinalg::matrix::{Matrix, BaseMatrix}; |
| /// |
| /// let a = Matrix::<f64>::ones(3,3); |
| /// |
| /// let b = &a.select_rows(&[2]); |
| /// assert_eq!(b.rows(), 1); |
| /// assert_eq!(b.cols(), 3); |
| /// |
| /// let c = &a.select_rows(&[1,2]); |
| /// assert_eq!(c.rows(), 2); |
| /// assert_eq!(c.cols(), 3); |
| /// ``` |
| /// |
| /// # Panics |
| /// |
| /// - Panics if row indices exceed the matrix dimensions. |
| </span><span class="kw">fn </span>select_rows<<span class="lifetime">'a</span>, I>(<span class="kw-2">&</span><span class="self">self</span>, rows: I) -> Matrix<T> |
| <span class="kw">where </span>T: Copy, |
| I: IntoIterator<Item = <span class="kw-2">&</span><span class="lifetime">'a </span>usize>, |
| I::IntoIter: ExactSizeIterator + Clone |
| { |
| <span class="kw">let </span>row_iter = rows.into_iter(); |
| <span class="kw">let </span><span class="kw-2">mut </span>mat_vec = Vec::with_capacity(row_iter.len() * <span class="self">self</span>.cols()); |
| |
| <span class="kw">for </span>row <span class="kw">in </span>row_iter.clone() { |
| <span class="macro">assert!</span>(<span class="kw-2">*</span>row < <span class="self">self</span>.rows(), |
| <span class="string">"Row index is greater than number of rows."</span>); |
| } |
| |
| <span class="kw">for </span>row_idx <span class="kw">in </span>row_iter.clone() { |
| <span class="kw">unsafe </span>{ |
| <span class="kw">let </span>row = <span class="self">self</span>.row_unchecked(<span class="kw-2">*</span>row_idx); |
| mat_vec.extend_from_slice(row.raw_slice()); |
| } |
| } |
| |
| Matrix { |
| cols: <span class="self">self</span>.cols(), |
| rows: row_iter.len(), |
| data: mat_vec, |
| } |
| } |
| |
| <span class="doccomment">/// Select columns from matrix |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// use rulinalg::matrix::{Matrix, BaseMatrix}; |
| /// |
| /// let a = Matrix::<f64>::ones(3,3); |
| /// let b = &a.select_cols(&[2]); |
| /// assert_eq!(b.rows(), 3); |
| /// assert_eq!(b.cols(), 1); |
| /// |
| /// let c = &a.select_cols(&[1,2]); |
| /// assert_eq!(c.rows(), 3); |
| /// assert_eq!(c.cols(), 2); |
| /// ``` |
| /// |
| /// # Panics |
| /// |
| /// - Panics if column indices exceed the matrix dimensions. |
| </span><span class="kw">fn </span>select_cols<<span class="lifetime">'a</span>, I>(<span class="kw-2">&</span><span class="self">self</span>, cols: I) -> Matrix<T> |
| <span class="kw">where </span>T: Copy, |
| I: IntoIterator<Item = <span class="kw-2">&</span><span class="lifetime">'a </span>usize>, |
| I::IntoIter: ExactSizeIterator + Clone |
| { |
| <span class="kw">let </span>col_iter = cols.into_iter(); |
| <span class="kw">let </span><span class="kw-2">mut </span>mat_vec = Vec::with_capacity(col_iter.len() * <span class="self">self</span>.rows()); |
| |
| <span class="kw">for </span>col <span class="kw">in </span>col_iter.clone() { |
| <span class="macro">assert!</span>(<span class="kw-2">*</span>col < <span class="self">self</span>.cols(), |
| <span class="string">"Column index is greater than number of columns."</span>); |
| } |
| |
| <span class="kw">unsafe </span>{ |
| <span class="kw">for </span>i <span class="kw">in </span><span class="number">0</span>..<span class="self">self</span>.rows() { |
| <span class="kw">for </span>col <span class="kw">in </span>col_iter.clone() { |
| mat_vec.push(<span class="kw-2">*</span><span class="self">self</span>.get_unchecked([i, <span class="kw-2">*</span>col])); |
| } |
| } |
| } |
| |
| Matrix { |
| cols: col_iter.len(), |
| rows: <span class="self">self</span>.rows(), |
| data: mat_vec, |
| } |
| } |
| |
| <span class="doccomment">/// The elementwise product of two matrices. |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// # #[macro_use] extern crate rulinalg; fn main() { |
| /// use rulinalg::matrix::{Matrix, BaseMatrix}; |
| /// |
| /// let a = matrix![1.0, 2.0; |
| /// 3.0, 4.0]; |
| /// let b = matrix![1.0, 2.0; |
| /// 3.0, 4.0]; |
| /// |
| /// let c = &a.elemul(&b); |
| /// assert_matrix_eq!(c, &matrix![1.0, 4.0; 9.0, 16.0]); |
| /// } |
| /// ``` |
| /// |
| /// # Panics |
| /// |
| /// - The matrices have different row counts. |
| /// - The matrices have different column counts. |
| </span><span class="kw">fn </span>elemul(<span class="kw-2">&</span><span class="self">self</span>, m: <span class="kw-2">&</span><span class="self">Self</span>) -> Matrix<T> |
| <span class="kw">where </span>T: Copy + Mul<T, Output = T> |
| { |
| <span class="macro">assert!</span>(<span class="self">self</span>.rows() == m.rows(), <span class="string">"Matrix row counts not equal."</span>); |
| <span class="macro">assert!</span>(<span class="self">self</span>.cols() == m.cols(), <span class="string">"Matrix column counts not equal."</span>); |
| |
| <span class="kw">let </span><span class="kw-2">mut </span>data = Vec::with_capacity(<span class="self">self</span>.rows() * <span class="self">self</span>.cols()); |
| <span class="kw">for </span>(self_r, m_r) <span class="kw">in </span><span class="self">self</span>.row_iter().zip(m.row_iter()) { |
| data.extend_from_slice(<span class="kw-2">&</span>utils::vec_bin_op(self_r.raw_slice(), m_r.raw_slice(), T::mul)); |
| } |
| Matrix::new(<span class="self">self</span>.rows(), <span class="self">self</span>.cols(), data) |
| } |
| |
| <span class="doccomment">/// The elementwise division of two matrices. |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// # #[macro_use] extern crate rulinalg; fn main() { |
| /// use rulinalg::matrix::{Matrix, BaseMatrix}; |
| /// |
| /// let a = matrix![1.0, 2.0; |
| /// 3.0, 4.0]; |
| /// let b = matrix![1.0, 2.0; |
| /// 3.0, 4.0]; |
| /// |
| /// let c = &a.elediv(&b); |
| /// assert_matrix_eq!(c, &matrix![1.0, 1.0; 1.0, 1.0]); |
| /// # } |
| /// ``` |
| /// |
| /// # Panics |
| /// |
| /// - The matrices have different row counts. |
| /// - The matrices have different column counts. |
| </span><span class="kw">fn </span>elediv(<span class="kw-2">&</span><span class="self">self</span>, m: <span class="kw-2">&</span><span class="self">Self</span>) -> Matrix<T> |
| <span class="kw">where </span>T: Copy + Div<T, Output = T> |
| { |
| <span class="macro">assert!</span>(<span class="self">self</span>.rows() == m.rows(), <span class="string">"Matrix row counts not equal."</span>); |
| <span class="macro">assert!</span>(<span class="self">self</span>.cols() == m.cols(), <span class="string">"Matrix column counts not equal."</span>); |
| |
| <span class="kw">let </span><span class="kw-2">mut </span>data = Vec::with_capacity(<span class="self">self</span>.rows() * <span class="self">self</span>.cols()); |
| <span class="kw">for </span>(self_r, m_r) <span class="kw">in </span><span class="self">self</span>.row_iter().zip(m.row_iter()) { |
| data.extend_from_slice(<span class="kw-2">&</span>utils::vec_bin_op(self_r.raw_slice(), m_r.raw_slice(), T::div)); |
| } |
| Matrix::new(<span class="self">self</span>.rows(), <span class="self">self</span>.cols(), data) |
| } |
| |
| <span class="doccomment">/// Select block matrix from matrix |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// use rulinalg::matrix::{Matrix, BaseMatrix}; |
| /// |
| /// let a = Matrix::<f64>::identity(3); |
| /// let b = &a.select(&[0,1], &[1,2]); |
| /// |
| /// // We get the 2x2 block matrix in the upper right corner. |
| /// assert_eq!(b.rows(), 2); |
| /// assert_eq!(b.cols(), 2); |
| /// |
| /// // Prints [0,0, 1,0] |
| /// println!("{:?}", b.data()); |
| /// ``` |
| /// |
| /// # Panics |
| /// |
| /// - Panics if row or column indices exceed the matrix dimensions. |
| </span><span class="kw">fn </span>select(<span class="kw-2">&</span><span class="self">self</span>, rows: <span class="kw-2">&</span>[usize], cols: <span class="kw-2">&</span>[usize]) -> Matrix<T> |
| <span class="kw">where </span>T: Copy |
| { |
| |
| <span class="kw">let </span><span class="kw-2">mut </span>mat_vec = Vec::with_capacity(cols.len() * rows.len()); |
| |
| <span class="kw">for </span>col <span class="kw">in </span>cols { |
| <span class="macro">assert!</span>(<span class="kw-2">*</span>col < <span class="self">self</span>.cols(), |
| <span class="string">"Column index is greater than number of columns."</span>); |
| } |
| |
| <span class="kw">for </span>row <span class="kw">in </span>rows { |
| <span class="macro">assert!</span>(<span class="kw-2">*</span>row < <span class="self">self</span>.rows(), |
| <span class="string">"Row index is greater than number of columns."</span>); |
| } |
| |
| <span class="kw">unsafe </span>{ |
| <span class="kw">for </span>row <span class="kw">in </span>rows { |
| <span class="kw">for </span>col <span class="kw">in </span>cols { |
| mat_vec.push(<span class="kw-2">*</span><span class="self">self</span>.get_unchecked([<span class="kw-2">*</span>row, <span class="kw-2">*</span>col])); |
| } |
| } |
| } |
| |
| Matrix { |
| cols: cols.len(), |
| rows: rows.len(), |
| data: mat_vec, |
| } |
| } |
| |
| <span class="doccomment">/// Horizontally concatenates two matrices. With self on the left. |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// # #[macro_use] extern crate rulinalg; fn main() { |
| /// use rulinalg::matrix::{Matrix, BaseMatrix}; |
| /// |
| /// let a = matrix![1.0, 2.0; |
| /// 3.0, 4.0; |
| /// 5.0, 6.0]; |
| /// let b = matrix![4.0; |
| /// 5.0; |
| /// 6.0]; |
| /// |
| /// let c = &a.hcat(&b); |
| /// assert_eq!(c.cols(), a.cols() + b.cols()); |
| /// assert_eq!(c[[1, 2]], 5.0); |
| /// # } |
| /// ``` |
| /// |
| /// # Panics |
| /// |
| /// - Self and m have different row counts. |
| </span><span class="kw">fn </span>hcat<S>(<span class="kw-2">&</span><span class="self">self</span>, m: <span class="kw-2">&</span>S) -> Matrix<T> |
| <span class="kw">where </span>T: Copy, |
| S: BaseMatrix<T> |
| { |
| <span class="macro">assert!</span>(<span class="self">self</span>.rows() == m.rows(), <span class="string">"Matrix row counts are not equal."</span>); |
| |
| <span class="kw">let </span><span class="kw-2">mut </span>new_data = Vec::with_capacity((<span class="self">self</span>.cols() + m.cols()) * <span class="self">self</span>.rows()); |
| |
| <span class="kw">for </span>(self_row, m_row) <span class="kw">in </span><span class="self">self</span>.row_iter().zip(m.row_iter()) { |
| new_data.extend_from_slice(self_row.raw_slice()); |
| new_data.extend_from_slice(m_row.raw_slice()); |
| } |
| |
| Matrix { |
| cols: (<span class="self">self</span>.cols() + m.cols()), |
| rows: <span class="self">self</span>.rows(), |
| data: new_data, |
| } |
| } |
| |
| <span class="doccomment">/// Vertically concatenates two matrices. With self on top. |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// # #[macro_use] extern crate rulinalg; fn main() { |
| /// use rulinalg::matrix::{Matrix, BaseMatrix}; |
| /// |
| /// let a = matrix![1.0, 2.0, 3.0; |
| /// 4.0, 5.0, 6.0]; |
| /// let b = matrix![4.0, 5.0, 6.0];; |
| /// |
| /// let c = &a.vcat(&b); |
| /// assert_eq!(c.rows(), a.rows() + b.rows()); |
| /// assert_eq!(c[[2, 2]], 6.0); |
| /// # } |
| /// ``` |
| /// |
| /// # Panics |
| /// |
| /// - Self and m have different column counts. |
| </span><span class="kw">fn </span>vcat<S>(<span class="kw-2">&</span><span class="self">self</span>, m: <span class="kw-2">&</span>S) -> Matrix<T> |
| <span class="kw">where </span>T: Copy, |
| S: BaseMatrix<T> |
| { |
| <span class="macro">assert!</span>(<span class="self">self</span>.cols() == m.cols(), |
| <span class="string">"Matrix column counts are not equal."</span>); |
| |
| <span class="kw">let </span><span class="kw-2">mut </span>new_data = Vec::with_capacity((<span class="self">self</span>.rows() + m.rows()) * <span class="self">self</span>.cols()); |
| |
| <span class="kw">for </span>row <span class="kw">in </span><span class="self">self</span>.row_iter().chain(m.row_iter()) { |
| new_data.extend_from_slice(row.raw_slice()); |
| } |
| |
| Matrix { |
| cols: <span class="self">self</span>.cols(), |
| rows: (<span class="self">self</span>.rows() + m.rows()), |
| data: new_data, |
| } |
| } |
| |
| <span class="doccomment">/// Extract the diagonal of the matrix |
| /// |
| /// Examples |
| /// |
| /// ``` |
| /// # #[macro_use] |
| /// # extern crate rulinalg; |
| /// |
| /// use rulinalg::matrix::BaseMatrix; |
| /// |
| /// # fn main() { |
| /// let a = matrix![1, 2, 3; |
| /// 4, 5, 6; |
| /// 7, 8, 9].diag().cloned().collect::<Vec<_>>(); |
| /// let b = matrix![1, 2; |
| /// 3, 4; |
| /// 5, 6].diag().cloned().collect::<Vec<_>>(); |
| /// |
| /// assert_eq!(a, vec![1, 5, 9]); |
| /// assert_eq!(b, vec![1, 4]); |
| /// # } |
| /// ``` |
| </span><span class="kw">fn </span>diag(<span class="kw-2">&</span><span class="self">self</span>) -> Diagonal<T, <span class="self">Self</span>> { |
| <span class="self">self</span>.diag_iter(DiagOffset::Main) |
| } |
| |
| <span class="doccomment">/// Tranposes the given matrix |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// # #[macro_use] extern crate rulinalg; fn main() { |
| /// use rulinalg::matrix::{Matrix, BaseMatrix}; |
| /// |
| /// let mat = matrix![1.0, 2.0, 3.0; |
| /// 4.0, 5.0, 6.0]; |
| /// |
| /// let expected = matrix![1.0, 4.0; |
| /// 2.0, 5.0; |
| /// 3.0, 6.0]; |
| /// assert_matrix_eq!(mat.transpose(), expected); |
| /// # } |
| /// ``` |
| </span><span class="kw">fn </span>transpose(<span class="kw-2">&</span><span class="self">self</span>) -> Matrix<T> |
| <span class="kw">where </span>T: Copy |
| { |
| <span class="kw">let </span><span class="kw-2">mut </span>new_data = Vec::with_capacity(<span class="self">self</span>.rows() * <span class="self">self</span>.cols()); |
| |
| <span class="kw">unsafe </span>{ |
| new_data.set_len(<span class="self">self</span>.rows() * <span class="self">self</span>.cols()); |
| <span class="kw">for </span>i <span class="kw">in </span><span class="number">0</span>..<span class="self">self</span>.cols() { |
| <span class="kw">for </span>j <span class="kw">in </span><span class="number">0</span>..<span class="self">self</span>.rows() { |
| <span class="kw-2">*</span>new_data.get_unchecked_mut(i * <span class="self">self</span>.rows() + j) = <span class="kw-2">*</span><span class="self">self</span>.get_unchecked([j, i]); |
| } |
| } |
| } |
| |
| Matrix { |
| cols: <span class="self">self</span>.rows(), |
| rows: <span class="self">self</span>.cols(), |
| data: new_data, |
| } |
| } |
| |
| <span class="doccomment">/// Checks if matrix is diagonal. |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// # #[macro_use] extern crate rulinalg; fn main() { |
| /// use rulinalg::matrix::{Matrix, BaseMatrix}; |
| /// |
| /// let a = matrix![1.0, 0.0; |
| /// 0.0, 1.0]; |
| /// let a_diag = a.is_diag(); |
| /// |
| /// assert_eq!(a_diag, true); |
| /// |
| /// let b = matrix![1.0, 0.0; |
| /// 1.0, 0.0]; |
| /// let b_diag = b.is_diag(); |
| /// |
| /// assert_eq!(b_diag, false); |
| /// # } |
| /// ``` |
| </span><span class="kw">fn </span>is_diag(<span class="kw-2">&</span><span class="self">self</span>) -> bool |
| <span class="kw">where </span>T: Zero + PartialEq |
| { |
| <span class="kw">let </span><span class="kw-2">mut </span>next_diag = <span class="number">0usize</span>; |
| <span class="self">self</span>.iter().enumerate().all(|(i, data)| <span class="kw">if </span>i == next_diag { |
| next_diag += <span class="self">self</span>.cols() + <span class="number">1</span>; |
| <span class="bool-val">true |
| </span>} <span class="kw">else </span>{ |
| data == <span class="kw-2">&</span>T::zero() |
| }) |
| } |
| |
| <span class="doccomment">/// Solves an upper triangular linear system. |
| /// |
| /// Given a matrix `A` and a vector `b`, this function returns the |
| /// solution of the upper triangular system `Ux = b`, where `U` is |
| /// the upper triangular part of `A`. |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// # #[macro_use] extern crate rulinalg; fn main() { |
| /// use rulinalg::matrix::BaseMatrix; |
| /// use std::f32; |
| /// |
| /// let u = matrix![1.0, 2.0; |
| /// 0.0, 1.0]; |
| /// let y = vector![3.0, 1.0]; |
| /// |
| /// let x = u.solve_u_triangular(y).expect("A solution should exist!"); |
| /// assert!((x[0] - 1.0) < f32::EPSILON); |
| /// assert!((x[1] - 1.0) < f32::EPSILON); |
| /// # } |
| /// ``` |
| /// |
| /// # Panics |
| /// |
| /// - Vector size and matrix column count are not equal. |
| /// |
| /// # Failures |
| /// |
| /// - There is no valid solution to the system (matrix is singular). |
| /// - The matrix is empty. |
| </span><span class="kw">fn </span>solve_u_triangular(<span class="kw-2">&</span><span class="self">self</span>, y: Vector<T>) -> <span class="prelude-ty">Result</span><Vector<T>, Error> |
| <span class="kw">where </span>T: Any + Float |
| { |
| <span class="macro">assert!</span>(<span class="self">self</span>.cols() == y.size(), |
| <span class="macro">format!</span>(<span class="string">"Vector size {0} != {1} Matrix column count."</span>, |
| y.size(), |
| <span class="self">self</span>.cols())); |
| |
| back_substitution(<span class="self">self</span>, y) |
| } |
| |
| <span class="doccomment">/// Solves a lower triangular linear system. |
| /// |
| /// Given a matrix `A` and a vector `b`, this function returns the |
| /// solution of the lower triangular system `Lx = b`, where `L` is |
| /// the lower triangular part of `A`. |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// # #[macro_use] extern crate rulinalg; fn main() { |
| /// use rulinalg::matrix::BaseMatrix; |
| /// use std::f32; |
| /// |
| /// let l = matrix![1.0, 0.0; |
| /// 2.0, 1.0]; |
| /// let y = vector![1.0, 3.0]; |
| /// |
| /// let x = l.solve_l_triangular(y).expect("A solution should exist!"); |
| /// println!("{:?}", x); |
| /// assert!((x[0] - 1.0) < f32::EPSILON); |
| /// assert!((x[1] - 1.0) < f32::EPSILON); |
| /// # } |
| /// ``` |
| /// |
| /// # Panics |
| /// |
| /// - Vector size and matrix column count are not equal. |
| /// |
| /// # Failures |
| /// |
| /// - There is no valid solution to the system (matrix is singular). |
| /// - The matrix is empty. |
| </span><span class="kw">fn </span>solve_l_triangular(<span class="kw-2">&</span><span class="self">self</span>, y: Vector<T>) -> <span class="prelude-ty">Result</span><Vector<T>, Error> |
| <span class="kw">where </span>T: Any + Float |
| { |
| <span class="macro">assert!</span>(<span class="self">self</span>.cols() == y.size(), |
| <span class="macro">format!</span>(<span class="string">"Vector size {0} != {1} Matrix column count."</span>, |
| y.size(), |
| <span class="self">self</span>.cols())); |
| |
| forward_substitution(<span class="self">self</span>, y) |
| } |
| |
| <span class="doccomment">/// Split the matrix at the specified axis returning two `MatrixSlice`s. |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// use rulinalg::matrix::{Axes, Matrix, BaseMatrix}; |
| /// |
| /// let a = Matrix::new(3,3, vec![2.0; 9]); |
| /// let (b,c) = a.split_at(1, Axes::Row); |
| /// ``` |
| </span><span class="kw">fn </span>split_at(<span class="kw-2">&</span><span class="self">self</span>, mid: usize, axis: Axes) -> (MatrixSlice<T>, MatrixSlice<T>) { |
| <span class="kw">let </span>slice_1: MatrixSlice<T>; |
| <span class="kw">let </span>slice_2: MatrixSlice<T>; |
| |
| <span class="kw">match </span>axis { |
| Axes::Row => { |
| <span class="macro">assert!</span>(mid < <span class="self">self</span>.rows()); |
| <span class="kw">unsafe </span>{ |
| slice_1 = MatrixSlice::from_raw_parts(<span class="self">self</span>.as_ptr(), |
| mid, |
| <span class="self">self</span>.cols(), |
| <span class="self">self</span>.row_stride()); |
| slice_2 = MatrixSlice::from_raw_parts(<span class="self">self</span>.as_ptr() |
| .offset((mid * <span class="self">self</span>.row_stride()) <span class="kw">as |
| </span>isize), |
| <span class="self">self</span>.rows() - mid, |
| <span class="self">self</span>.cols(), |
| <span class="self">self</span>.row_stride()); |
| } |
| } |
| Axes::Col => { |
| <span class="macro">assert!</span>(mid < <span class="self">self</span>.cols()); |
| <span class="kw">unsafe </span>{ |
| slice_1 = MatrixSlice::from_raw_parts(<span class="self">self</span>.as_ptr(), |
| <span class="self">self</span>.rows(), |
| mid, |
| <span class="self">self</span>.row_stride()); |
| slice_2 = MatrixSlice::from_raw_parts(<span class="self">self</span>.as_ptr().offset(mid <span class="kw">as </span>isize), |
| <span class="self">self</span>.rows(), |
| <span class="self">self</span>.cols() - mid, |
| <span class="self">self</span>.row_stride()); |
| } |
| } |
| } |
| |
| (slice_1, slice_2) |
| } |
| |
| <span class="doccomment">/// Produce a `MatrixSlice` from an existing matrix. |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// use rulinalg::matrix::{Matrix, BaseMatrix, MatrixSlice}; |
| /// |
| /// let a = Matrix::new(3,3, (0..9).collect::<Vec<usize>>()); |
| /// let slice = MatrixSlice::from_matrix(&a, [1,1], 2, 2); |
| /// let new_slice = slice.sub_slice([0,0], 1, 1); |
| /// ``` |
| </span><span class="kw">fn </span>sub_slice<<span class="lifetime">'a</span>>(<span class="kw-2">&</span><span class="self">self</span>, start: [usize; <span class="number">2</span>], rows: usize, cols: usize) -> MatrixSlice<<span class="lifetime">'a</span>, T> |
| <span class="kw">where </span>T: <span class="lifetime">'a |
| </span>{ |
| <span class="macro">assert!</span>(start[<span class="number">0</span>] + rows <= <span class="self">self</span>.rows(), |
| <span class="string">"View dimensions exceed matrix dimensions."</span>); |
| <span class="macro">assert!</span>(start[<span class="number">1</span>] + cols <= <span class="self">self</span>.cols(), |
| <span class="string">"View dimensions exceed matrix dimensions."</span>); |
| |
| <span class="kw">unsafe </span>{ |
| MatrixSlice::from_raw_parts(<span class="self">self</span>.as_ptr() |
| .offset((start[<span class="number">0</span>] * <span class="self">self</span>.row_stride() + start[<span class="number">1</span>]) <span class="kw">as |
| </span>isize), |
| rows, |
| cols, |
| <span class="self">self</span>.row_stride()) |
| } |
| } |
| } |
| |
| <span class="doccomment">/// Trait for mutable matrices. |
| </span><span class="kw">pub trait </span>BaseMatrixMut<T>: BaseMatrix<T> { |
| <span class="doccomment">/// Top left index of the slice. |
| </span><span class="kw">fn </span>as_mut_ptr(<span class="kw-2">&mut </span><span class="self">self</span>) -> <span class="kw-2">*mut </span>T; |
| |
| <span class="doccomment">/// Returns a `MatrixSliceMut` over the whole matrix. |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// use rulinalg::matrix::{Matrix, BaseMatrixMut}; |
| /// |
| /// let mut a = Matrix::new(3, 3, vec![2.0; 9]); |
| /// let b = a.as_mut_slice(); |
| /// ``` |
| </span><span class="kw">fn </span>as_mut_slice(<span class="kw-2">&mut </span><span class="self">self</span>) -> MatrixSliceMut<T> { |
| <span class="kw">unsafe </span>{ |
| MatrixSliceMut::from_raw_parts(<span class="self">self</span>.as_mut_ptr(), |
| <span class="self">self</span>.rows(), |
| <span class="self">self</span>.cols(), |
| <span class="self">self</span>.row_stride()) |
| } |
| } |
| |
| <span class="doccomment">/// Get a mutable reference to an element in the matrix without bounds checks. |
| </span><span class="kw">unsafe fn </span>get_unchecked_mut(<span class="kw-2">&mut </span><span class="self">self</span>, index: [usize; <span class="number">2</span>]) -> <span class="kw-2">&mut </span>T { |
| <span class="kw-2">&mut *</span>(<span class="self">self</span>.as_mut_ptr().offset((index[<span class="number">0</span>] * <span class="self">self</span>.row_stride() + index[<span class="number">1</span>]) <span class="kw">as </span>isize)) |
| } |
| |
| <span class="doccomment">/// Get a mutable reference to an element in the matrix. |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// # #[macro_use] extern crate rulinalg; fn main() { |
| /// use rulinalg::matrix::{Matrix, BaseMatrix, BaseMatrixMut}; |
| /// |
| /// let mut mat = matrix![0, 1; |
| /// 3, 4; |
| /// 6, 7]; |
| /// |
| /// assert_eq!(mat.get_mut([0, 2]), None); |
| /// assert_eq!(mat.get_mut([3, 0]), None); |
| /// |
| /// assert_eq!(*mat.get_mut([0, 0]).unwrap(), 0); |
| /// *mat.get_mut([0,0]).unwrap() = 2; |
| /// assert_eq!(*mat.get_mut([0, 0]).unwrap(), 2); |
| /// # } |
| /// ``` |
| </span><span class="kw">fn </span>get_mut(<span class="kw-2">&mut </span><span class="self">self</span>, index: [usize; <span class="number">2</span>]) -> <span class="prelude-ty">Option</span><<span class="kw-2">&mut </span>T> { |
| <span class="kw">let </span>row_ind = index[<span class="number">0</span>]; |
| <span class="kw">let </span>col_ind = index[<span class="number">1</span>]; |
| |
| <span class="kw">if </span>row_ind >= <span class="self">self</span>.rows() || col_ind >= <span class="self">self</span>.cols() { |
| <span class="prelude-val">None |
| </span>} <span class="kw">else </span>{ |
| <span class="kw">unsafe </span>{ <span class="prelude-val">Some</span>(<span class="self">self</span>.get_unchecked_mut(index)) } |
| } |
| } |
| |
| <span class="doccomment">/// Returns a mutable iterator over the matrix. |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// # #[macro_use] extern crate rulinalg; fn main() { |
| /// use rulinalg::matrix::{Matrix, BaseMatrixMut}; |
| /// |
| /// let mut a = Matrix::new(3,3, (0..9).collect::<Vec<usize>>()); |
| /// |
| /// { |
| /// let mut slice = a.sub_slice_mut([1,1], 2, 2); |
| /// |
| /// for d in slice.iter_mut() { |
| /// *d = *d + 2; |
| /// } |
| /// } |
| /// |
| /// // Only the matrix slice is updated. |
| /// assert_matrix_eq!(a, matrix![0, 1, 2; 3, 6, 7; 6, 9, 10]); |
| /// # } |
| /// ``` |
| </span><span class="kw">fn </span>iter_mut<<span class="lifetime">'a</span>>(<span class="kw-2">&mut </span><span class="self">self</span>) -> SliceIterMut<<span class="lifetime">'a</span>, T> |
| <span class="kw">where </span>T: <span class="lifetime">'a |
| </span>{ |
| SliceIterMut { |
| slice_start: <span class="self">self</span>.as_mut_ptr(), |
| row_pos: <span class="number">0</span>, |
| col_pos: <span class="number">0</span>, |
| slice_rows: <span class="self">self</span>.rows(), |
| slice_cols: <span class="self">self</span>.cols(), |
| row_stride: <span class="self">self</span>.row_stride(), |
| _marker: PhantomData::<<span class="kw-2">&mut </span>T>, |
| } |
| } |
| |
| <span class="doccomment">/// Returns a mutable reference to the column of a matrix at the given index. |
| /// `None` if the index is out of bounds. |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// # #[macro_use] |
| /// # extern crate rulinalg; |
| /// |
| /// # fn main() { |
| /// use rulinalg::matrix::{Matrix, BaseMatrixMut}; |
| /// |
| /// let mut mat = matrix![0, 1, 2; |
| /// 3, 4, 5; |
| /// 6, 7, 8]; |
| /// let mut slice = mat.sub_slice_mut([1,1], 2, 2); |
| /// { |
| /// let col = slice.col_mut(1); |
| /// let mut expected = matrix![5usize; 8]; |
| /// assert_matrix_eq!(*col, expected); |
| /// } |
| /// # } |
| /// ``` |
| /// |
| /// # Panics |
| /// |
| /// Will panic if the column index is out of bounds. |
| </span><span class="kw">fn </span>col_mut(<span class="kw-2">&mut </span><span class="self">self</span>, index: usize) -> ColumnMut<T> { |
| <span class="kw">if </span>index < <span class="self">self</span>.cols() { |
| <span class="kw">unsafe </span>{ <span class="self">self</span>.col_unchecked_mut(index) } |
| } <span class="kw">else </span>{ |
| <span class="macro">panic!</span>(<span class="string">"Column index out of bounds."</span>) |
| } |
| } |
| |
| <span class="doccomment">/// Returns a mutable reference to the column of a matrix at the given index |
| /// without doing a bounds check. |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// # #[macro_use] |
| /// # extern crate rulinalg; |
| /// |
| /// # fn main() { |
| /// use rulinalg::matrix::{Matrix, BaseMatrixMut}; |
| /// |
| /// let mut mat = matrix![0, 1, 2; |
| /// 3, 4, 5; |
| /// 6, 7, 8]; |
| /// let mut slice = mat.sub_slice_mut([1,1], 2, 2); |
| /// let col = unsafe { slice.col_unchecked_mut(1) }; |
| /// let mut expected = matrix![5usize; 8]; |
| /// assert_matrix_eq!(*col, expected); |
| /// # } |
| /// ``` |
| </span><span class="kw">unsafe fn </span>col_unchecked_mut(<span class="kw-2">&mut </span><span class="self">self</span>, index: usize) -> ColumnMut<T> { |
| <span class="kw">let </span>ptr = <span class="self">self</span>.as_mut_ptr().offset(index <span class="kw">as </span>isize); |
| ColumnMut { col: MatrixSliceMut::from_raw_parts(ptr, <span class="self">self</span>.rows(), <span class="number">1</span>, <span class="self">self</span>.row_stride()) } |
| } |
| |
| <span class="doccomment">/// Returns a mutable reference to the row of a matrix at the given index. |
| /// `None` if the index is out of bounds. |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// # #[macro_use] |
| /// # extern crate rulinalg; |
| /// |
| /// # fn main() { |
| /// use rulinalg::matrix::{Matrix, BaseMatrixMut}; |
| /// |
| /// let mut mat = matrix![0, 1, 2; |
| /// 3, 4, 5; |
| /// 6, 7, 8]; |
| /// let mut slice = mat.sub_slice_mut([1,1], 2, 2); |
| /// { |
| /// let row = slice.row_mut(1); |
| /// let mut expected = matrix![7usize, 8]; |
| /// assert_matrix_eq!(*row, expected); |
| /// } |
| /// # } |
| /// ``` |
| /// |
| /// # Panics |
| /// |
| /// Will panic if the row index is out of bounds. |
| </span><span class="kw">fn </span>row_mut(<span class="kw-2">&mut </span><span class="self">self</span>, index: usize) -> RowMut<T> { |
| <span class="kw">if </span>index < <span class="self">self</span>.rows() { |
| <span class="kw">unsafe </span>{ <span class="self">self</span>.row_unchecked_mut(index) } |
| } <span class="kw">else </span>{ |
| <span class="macro">panic!</span>(<span class="string">"Row index out of bounds."</span>) |
| } |
| } |
| |
| <span class="doccomment">/// Returns a mutable reference to the row of a matrix at the given index |
| /// without doing a bounds check. |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// # #[macro_use] |
| /// # extern crate rulinalg; |
| /// |
| /// # fn main() { |
| /// use rulinalg::matrix::{Matrix, BaseMatrixMut}; |
| /// |
| /// let mut mat = matrix![0, 1, 2; |
| /// 3, 4, 5; |
| /// 6, 7, 8]; |
| /// let mut slice = mat.sub_slice_mut([1,1], 2, 2); |
| /// let row = unsafe { slice.row_unchecked_mut(1) }; |
| /// let mut expected = matrix![7usize, 8]; |
| /// assert_matrix_eq!(*row, expected); |
| /// # } |
| /// ``` |
| </span><span class="kw">unsafe fn </span>row_unchecked_mut(<span class="kw-2">&mut </span><span class="self">self</span>, index: usize) -> RowMut<T> { |
| <span class="kw">let </span>ptr = <span class="self">self</span>.as_mut_ptr().offset((<span class="self">self</span>.row_stride() * index) <span class="kw">as </span>isize); |
| RowMut { row: MatrixSliceMut::from_raw_parts(ptr, <span class="number">1</span>, <span class="self">self</span>.cols(), <span class="self">self</span>.row_stride()) } |
| } |
| |
| <span class="doccomment">/// Swaps two rows in a matrix. |
| /// |
| /// If `a == b`, this method does nothing. |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// # #[macro_use] extern crate rulinalg; |
| /// # fn main() { |
| /// use rulinalg::matrix::{Matrix, BaseMatrixMut}; |
| /// |
| /// let mut x = matrix![0, 1; |
| /// 2, 3; |
| /// 4, 5; |
| /// 6, 7]; |
| /// |
| /// x.swap_rows(1, 3); |
| /// let expected = matrix![0, 1; |
| /// 6, 7; |
| /// 4, 5; |
| /// 2, 3]; |
| /// |
| /// assert_matrix_eq!(x, expected); |
| /// # } |
| /// ``` |
| /// |
| /// # Panics |
| /// |
| /// Panics if `a` or `b` are out of bounds. |
| </span><span class="kw">fn </span>swap_rows(<span class="kw-2">&mut </span><span class="self">self</span>, a: usize, b: usize) { |
| <span class="macro">assert!</span>(a < <span class="self">self</span>.rows(), |
| <span class="macro">format!</span>(<span class="string">"Row index {0} larger than row count {1}"</span>, a, <span class="self">self</span>.rows())); |
| <span class="macro">assert!</span>(b < <span class="self">self</span>.rows(), |
| <span class="macro">format!</span>(<span class="string">"Row index {0} larger than row count {1}"</span>, b, <span class="self">self</span>.rows())); |
| |
| <span class="kw">if </span>a != b { |
| <span class="kw">unsafe </span>{ |
| <span class="kw">let </span>row_a = slice::from_raw_parts_mut(<span class="self">self</span>.as_mut_ptr() |
| .offset((<span class="self">self</span>.row_stride() * a) <span class="kw">as |
| </span>isize), |
| <span class="self">self</span>.cols()); |
| <span class="kw">let </span>row_b = slice::from_raw_parts_mut(<span class="self">self</span>.as_mut_ptr() |
| .offset((<span class="self">self</span>.row_stride() * b) <span class="kw">as |
| </span>isize), |
| <span class="self">self</span>.cols()); |
| |
| <span class="kw">for </span>(x, y) <span class="kw">in </span>row_a.into_iter().zip(row_b.into_iter()) { |
| mem::swap(x, y); |
| } |
| } |
| } |
| |
| } |
| |
| <span class="doccomment">/// Swaps two columns in a matrix. |
| /// |
| /// If `a == b`, this method does nothing. |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// # #[macro_use] extern crate rulinalg; |
| /// # fn main() { |
| /// use rulinalg::matrix::{Matrix, BaseMatrixMut}; |
| /// |
| /// let mut x = matrix![0, 1; |
| /// 2, 3; |
| /// 4, 5]; |
| /// |
| /// x.swap_cols(0, 1); |
| /// let expected = matrix![1, 0; |
| /// 3, 2; |
| /// 5, 4]; |
| /// |
| /// assert_matrix_eq!(x, expected); |
| /// # } |
| /// ``` |
| /// |
| /// # Panics |
| /// |
| /// Panics if `a` or `b` are out of bounds. |
| </span><span class="kw">fn </span>swap_cols(<span class="kw-2">&mut </span><span class="self">self</span>, a: usize, b: usize) { |
| <span class="macro">assert!</span>(a < <span class="self">self</span>.cols(), |
| <span class="macro">format!</span>(<span class="string">"Row index {0} larger than row count {1}"</span>, a, <span class="self">self</span>.rows())); |
| <span class="macro">assert!</span>(b < <span class="self">self</span>.cols(), |
| <span class="macro">format!</span>(<span class="string">"Row index {0} larger than row count {1}"</span>, b, <span class="self">self</span>.rows())); |
| |
| <span class="kw">if </span>a != b { |
| <span class="kw">unsafe </span>{ |
| <span class="kw">for </span>i <span class="kw">in </span><span class="number">0</span>..<span class="self">self</span>.rows() { |
| <span class="kw">let </span>a_ptr: <span class="kw-2">*mut </span>T = <span class="self">self</span>.get_unchecked_mut([i, a]); |
| <span class="kw">let </span>b_ptr: <span class="kw-2">*mut </span>T = <span class="self">self</span>.get_unchecked_mut([i, b]); |
| ptr::swap(a_ptr, b_ptr); |
| } |
| } |
| } |
| |
| } |
| |
| <span class="doccomment">/// Iterate over the mutable columns of the matrix. |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// # #[macro_use] extern crate rulinalg; fn main() { |
| /// use rulinalg::matrix::{Matrix, BaseMatrixMut}; |
| /// |
| /// let mut a = matrix![0, 1; |
| /// 2, 3; |
| /// 4, 5]; |
| /// |
| /// for mut col in a.col_iter_mut() { |
| /// *col += 1; |
| /// } |
| /// |
| /// // Now contains the range 1..7 |
| /// println!("{}", a); |
| /// # } |
| /// ``` |
| </span><span class="kw">fn </span>col_iter_mut(<span class="kw-2">&mut </span><span class="self">self</span>) -> ColsMut<T> { |
| ColsMut { |
| _marker: PhantomData::<<span class="kw-2">&mut </span>T>, |
| col_pos: <span class="number">0</span>, |
| row_stride: <span class="self">self</span>.row_stride() <span class="kw">as </span>isize, |
| slice_cols: <span class="self">self</span>.cols(), |
| slice_rows: <span class="self">self</span>.rows(), |
| slice_start: <span class="self">self</span>.as_mut_ptr(), |
| } |
| } |
| |
| <span class="doccomment">/// Iterate over the mutable rows of the matrix. |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// # #[macro_use] extern crate rulinalg; fn main() { |
| /// use rulinalg::matrix::{Matrix, BaseMatrixMut}; |
| /// |
| /// let mut a = matrix![0, 1; |
| /// 2, 3; |
| /// 4, 5]; |
| /// |
| /// for mut row in a.row_iter_mut() { |
| /// *row += 1; |
| /// } |
| /// |
| /// // Now contains the range 1..7 |
| /// println!("{}", a); |
| /// # } |
| /// ``` |
| </span><span class="kw">fn </span>row_iter_mut(<span class="kw-2">&mut </span><span class="self">self</span>) -> RowsMut<T> { |
| RowsMut { |
| slice_start: <span class="self">self</span>.as_mut_ptr(), |
| row_pos: <span class="number">0</span>, |
| slice_rows: <span class="self">self</span>.rows(), |
| slice_cols: <span class="self">self</span>.cols(), |
| row_stride: <span class="self">self</span>.row_stride() <span class="kw">as </span>isize, |
| _marker: PhantomData::<<span class="kw-2">&mut </span>T>, |
| } |
| } |
| |
| <span class="doccomment">/// Iterate over diagonal entries mutably |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// # #[macro_use] extern crate rulinalg; |
| /// |
| /// # fn main() { |
| /// use rulinalg::matrix::{Matrix, BaseMatrixMut, DiagOffset}; |
| /// |
| /// let mut a = matrix![0, 1, 2; |
| /// 3, 4, 5; |
| /// 6, 7, 8]; |
| /// |
| /// // Increment super diag |
| /// for d in a.diag_iter_mut(DiagOffset::Above(1)) { |
| /// *d = *d + 1; |
| /// } |
| /// |
| /// // Zero the sub-diagonal (sets 3 and 7 to 0) |
| /// // Equivalent to `diag_iter(DiagOffset::Below(1))` |
| /// for sub_d in a.diag_iter_mut(DiagOffset::from(-1)) { |
| /// *sub_d = 0; |
| /// } |
| /// |
| /// println!("{}", a); |
| /// # } |
| /// ``` |
| /// |
| /// # Panics |
| /// |
| /// If using an `Above` or `Below` offset which is |
| /// out-of-bounds this function will panic. |
| /// |
| /// This function will never panic if the `Main` diagonal |
| /// offset is used. |
| </span><span class="kw">fn </span>diag_iter_mut(<span class="kw-2">&mut </span><span class="self">self</span>, k: DiagOffset) -> DiagonalMut<T, <span class="self">Self</span>> { |
| <span class="kw">let </span>(diag_len, diag_start) = <span class="kw">match </span>k.into() { |
| DiagOffset::Main => (min(<span class="self">self</span>.rows(), <span class="self">self</span>.cols()), <span class="number">0</span>), |
| DiagOffset::Above(m) => { |
| <span class="macro">assert!</span>(m < <span class="self">self</span>.cols(), |
| <span class="string">"Offset diagonal is not within matrix dimensions."</span>); |
| (min(<span class="self">self</span>.rows(), <span class="self">self</span>.cols() - m), m) |
| } |
| DiagOffset::Below(m) => { |
| <span class="macro">assert!</span>(m < <span class="self">self</span>.rows(), |
| <span class="string">"Offset diagonal is not within matrix dimensions."</span>); |
| (min(<span class="self">self</span>.rows() - m, <span class="self">self</span>.cols()), m * <span class="self">self</span>.row_stride()) |
| } |
| }; |
| |
| |
| <span class="kw">let </span>diag_end = diag_start + (diag_len - <span class="number">1</span>) * <span class="self">self</span>.row_stride() + diag_len; |
| DiagonalMut { |
| matrix: <span class="self">self</span>, |
| diag_pos: diag_start, |
| diag_end: diag_end, |
| _marker: PhantomData::<<span class="kw-2">&mut </span>T>, |
| } |
| } |
| |
| <span class="doccomment">/// Sets the underlying matrix data to the target data. |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// use rulinalg::matrix::{Matrix, BaseMatrixMut}; |
| /// |
| /// let mut mat = Matrix::<f32>::zeros(4,4); |
| /// let one_block = Matrix::<f32>::ones(2,2); |
| /// |
| /// // Get a mutable slice of the upper left 2x2 block. |
| /// let mat_block = mat.sub_slice_mut([0,0], 2, 2); |
| /// |
| /// // Set the upper left 2x2 block to be ones. |
| /// mat_block.set_to(one_block); |
| /// ``` |
| /// |
| /// # Panics |
| /// |
| /// Panics if the dimensions of `self` and `target` are not the same. |
| </span><span class="kw">fn </span>set_to<M: BaseMatrix<T>>(<span class="kw-2">mut </span><span class="self">self</span>, target: M) |
| <span class="kw">where </span>T: Copy |
| { |
| <span class="macro">assert!</span>(<span class="self">self</span>.rows() == target.rows(), |
| <span class="string">"Target has different row count to self."</span>); |
| <span class="macro">assert!</span>(<span class="self">self</span>.cols() == target.cols(), |
| <span class="string">"Target has different column count to self."</span>); |
| <span class="kw">for </span>(<span class="kw-2">mut </span>s, t) <span class="kw">in </span><span class="self">self</span>.row_iter_mut().zip(target.row_iter()) { |
| <span class="comment">// Vectorized assignment per row. |
| </span>utils::in_place_vec_bin_op(s.raw_slice_mut(), t.raw_slice(), |x, <span class="kw-2">&</span>y| <span class="kw-2">*</span>x = y); |
| } |
| } |
| |
| <span class="doccomment">/// Applies a function to each element in the matrix. |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// # #[macro_use] extern crate rulinalg; fn main() { |
| /// use rulinalg::matrix::{Matrix, BaseMatrixMut}; |
| /// |
| /// fn add_two(a: f64) -> f64 { |
| /// a + 2f64 |
| /// } |
| /// |
| /// let a = Matrix::new(2, 2, vec![0.;4]); |
| /// |
| /// let b = a.apply(&add_two); |
| /// |
| /// assert_eq!(b, matrix![2.0, 2.0; 2.0, 2.0]); |
| /// # } |
| /// ``` |
| </span><span class="kw">fn </span>apply(<span class="kw-2">mut </span><span class="self">self</span>, f: <span class="kw-2">&</span>Fn(T) -> T) -> <span class="self">Self |
| </span><span class="kw">where </span>T: Copy |
| { |
| <span class="kw">for </span>val <span class="kw">in </span><span class="self">self</span>.iter_mut() { |
| <span class="kw-2">*</span>val = f(<span class="kw-2">*</span>val); |
| } |
| <span class="self">self |
| </span>} |
| |
| <span class="doccomment">/// Split the matrix at the specified axis returning two `MatrixSliceMut`s. |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// use rulinalg::matrix::{Axes, Matrix, BaseMatrixMut}; |
| /// |
| /// let mut a = Matrix::new(3,3, vec![2.0; 9]); |
| /// let (b, c) = a.split_at_mut(1, Axes::Col); |
| /// ``` |
| </span><span class="kw">fn </span>split_at_mut(<span class="kw-2">&mut </span><span class="self">self</span>, mid: usize, axis: Axes) -> (MatrixSliceMut<T>, MatrixSliceMut<T>) { |
| |
| <span class="kw">let </span>slice_1: MatrixSliceMut<T>; |
| <span class="kw">let </span>slice_2: MatrixSliceMut<T>; |
| |
| <span class="kw">match </span>axis { |
| Axes::Row => { |
| <span class="macro">assert!</span>(mid < <span class="self">self</span>.rows()); |
| <span class="kw">unsafe </span>{ |
| slice_1 = MatrixSliceMut::from_raw_parts(<span class="self">self</span>.as_mut_ptr(), |
| mid, |
| <span class="self">self</span>.cols(), |
| <span class="self">self</span>.row_stride()); |
| slice_2 = MatrixSliceMut::from_raw_parts(<span class="self">self</span>.as_mut_ptr() |
| .offset((mid * |
| <span class="self">self</span>.row_stride()) <span class="kw">as |
| </span>isize), |
| <span class="self">self</span>.rows() - mid, |
| <span class="self">self</span>.cols(), |
| <span class="self">self</span>.row_stride()); |
| } |
| } |
| Axes::Col => { |
| <span class="macro">assert!</span>(mid < <span class="self">self</span>.cols()); |
| <span class="kw">unsafe </span>{ |
| slice_1 = MatrixSliceMut::from_raw_parts(<span class="self">self</span>.as_mut_ptr(), |
| <span class="self">self</span>.rows(), |
| mid, |
| <span class="self">self</span>.row_stride()); |
| slice_2 = MatrixSliceMut::from_raw_parts(<span class="self">self</span>.as_mut_ptr() |
| .offset(mid <span class="kw">as </span>isize), |
| <span class="self">self</span>.rows(), |
| <span class="self">self</span>.cols() - mid, |
| <span class="self">self</span>.row_stride()); |
| } |
| } |
| } |
| |
| (slice_1, slice_2) |
| } |
| |
| <span class="doccomment">/// Produce a `MatrixSliceMut` from an existing matrix. |
| /// |
| /// # Examples |
| /// |
| /// ``` |
| /// use rulinalg::matrix::{Matrix, MatrixSliceMut, BaseMatrixMut}; |
| /// |
| /// let mut a = Matrix::new(3,3, (0..9).collect::<Vec<usize>>()); |
| /// let mut slice = MatrixSliceMut::from_matrix(&mut a, [1,1], 2, 2); |
| /// let new_slice = slice.sub_slice_mut([0,0], 1, 1); |
| /// ``` |
| </span><span class="kw">fn </span>sub_slice_mut<<span class="lifetime">'a</span>>(<span class="kw-2">&mut </span><span class="self">self</span>, |
| start: [usize; <span class="number">2</span>], |
| rows: usize, |
| cols: usize) |
| -> MatrixSliceMut<<span class="lifetime">'a</span>, T> |
| <span class="kw">where </span>T: <span class="lifetime">'a |
| </span>{ |
| <span class="macro">assert!</span>(start[<span class="number">0</span>] + rows <= <span class="self">self</span>.rows(), |
| <span class="string">"View dimensions exceed matrix dimensions."</span>); |
| <span class="macro">assert!</span>(start[<span class="number">1</span>] + cols <= <span class="self">self</span>.cols(), |
| <span class="string">"View dimensions exceed matrix dimensions."</span>); |
| |
| <span class="kw">unsafe </span>{ |
| MatrixSliceMut::from_raw_parts(<span class="self">self</span>.as_mut_ptr() |
| .offset((start[<span class="number">0</span>] * <span class="self">self</span>.row_stride() + start[<span class="number">1</span>]) <span class="kw">as |
| </span>isize), |
| rows, |
| cols, |
| <span class="self">self</span>.row_stride()) |
| } |
| } |
| } |
| </code></pre></div> |
| </section></div></main><div id="rustdoc-vars" data-root-path="../../../../" data-current-crate="rulinalg" data-themes="ayu,dark,light" data-resource-suffix="" data-rustdoc-version="1.66.0-nightly (5c8bff74b 2022-10-21)" ></div></body></html> |