api-docs/crates-enclave/src/rulinalg/matrix/base/mod.rs.html - incubator-teaclave-website - Git at Google

 <!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::&lt;Vec&lt;usize&gt;&gt;());
 //!
 //! // 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 = &amp;mat_slice.transpose() * &amp;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&lt;T&gt;: Sized {
     <span class="doccomment">/// Rows in the matrix.
     </span><span class="kw">fn </span>rows(<span class="kw-2">&amp;</span><span class="self">self</span>) -&gt; usize;

     <span class="doccomment">/// Columns in the matrix.
     </span><span class="kw">fn </span>cols(<span class="kw-2">&amp;</span><span class="self">self</span>) -&gt; usize;

     <span class="doccomment">/// Row stride in the matrix.
     </span><span class="kw">fn </span>row_stride(<span class="kw-2">&amp;</span><span class="self">self</span>) -&gt; usize;

     <span class="doccomment">/// Returns true if the matrix contais no elements
     </span><span class="kw">fn </span>is_empty(<span class="kw-2">&amp;</span><span class="self">self</span>) -&gt; 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">&amp;</span><span class="self">self</span>) -&gt; <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">&amp;</span><span class="self">self</span>) -&gt; MatrixSlice&lt;T&gt; {
         <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">&amp;</span><span class="self">self</span>, index: [usize; <span class="number">2</span>]) -&gt; <span class="kw-2">&amp;</span>T {
         <span class="kw-2">&amp;*</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">&amp;</span><span class="self">self</span>, index: [usize; <span class="number">2</span>]) -&gt; <span class="prelude-ty">Option</span>&lt;<span class="kw-2">&amp;</span>T&gt; {
         <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 &gt;= <span class="self">self</span>.rows() || col_ind &gt;= <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">&amp;</span><span class="self">self</span>, index: usize) -&gt; Column&lt;T&gt; {
         <span class="kw">if </span>index &lt; <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">&quot;Column index out of bounds.&quot;</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">&amp;</span><span class="self">self</span>, index: usize) -&gt; Column&lt;T&gt; {
         <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">&amp;</span><span class="self">self</span>, index: usize) -&gt; Row&lt;T&gt; {
         <span class="kw">if </span>index &lt; <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">&quot;Row index out of bounds.&quot;</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">&amp;</span><span class="self">self</span>, index: usize) -&gt; Row&lt;T&gt; {
         <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::&lt;Vec&lt;usize&gt;&gt;();
     /// assert_eq!(slice_data, vec![4, 5, 7, 8]);
     /// # }
     /// ```
     </span><span class="kw">fn </span>iter&lt;<span class="lifetime">&#39;a</span>&gt;(<span class="kw-2">&amp;</span><span class="self">self</span>) -&gt; SliceIter&lt;<span class="lifetime">&#39;a</span>, T&gt;
         <span class="kw">where </span>T: <span class="lifetime">&#39;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::&lt;<span class="kw-2">&amp;</span>T&gt;,
         }
     }

     <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">&amp;</span><span class="self">self</span>) -&gt; Cols&lt;T&gt; {
         Cols {
             _marker: PhantomData::&lt;<span class="kw-2">&amp;</span>T&gt;,
             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">&amp;</span><span class="self">self</span>) -&gt; Rows&lt;T&gt; {
         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::&lt;<span class="kw-2">&amp;</span>T&gt;,
         }
     }

     <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!(&quot;{}&quot;, d);
     /// }
     ///
     /// // Print sub diag [3, 7]
     /// // Equivalent to `diag_iter(DiagOffset::Below(1))`
     /// for d in a.diag_iter(DiagOffset::from(-1)) {
     ///     println!(&quot;{}&quot;, 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">&amp;</span><span class="self">self</span>, k: DiagOffset) -&gt; Diagonal&lt;T, <span class="self">Self</span>&gt; {
         <span class="kw">let </span>(diag_len, diag_start) = <span class="kw">match </span>k.into() {
             DiagOffset::Main =&gt; (min(<span class="self">self</span>.rows(), <span class="self">self</span>.cols()), <span class="number">0</span>),
             DiagOffset::Above(m) =&gt; {
                 <span class="macro">assert!</span>(m &lt; <span class="self">self</span>.cols(),
                         <span class="string">&quot;Offset diagonal is not within matrix dimensions.&quot;</span>);
                 (min(<span class="self">self</span>.rows(), <span class="self">self</span>.cols() - m), m)
             }
             DiagOffset::Below(m) =&gt; {
                 <span class="macro">assert!</span>(m &lt; <span class="self">self</span>.rows(),
                         <span class="string">&quot;Offset diagonal is not within matrix dimensions.&quot;</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::&lt;<span class="kw-2">&amp;</span>T&gt;,
         }
     }

     <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">&amp;</span><span class="self">self</span>) -&gt; Vector&lt;T&gt;
         <span class="kw">where </span>T: Copy + Zero + Add&lt;T, Output = T&gt;
     {
         <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">&amp;mut </span>sum_rows, row.raw_slice(), |sum, <span class="kw-2">&amp;</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">&amp;</span><span class="self">self</span>) -&gt; Vector&lt;T&gt;
         <span class="kw">where </span>T: Copy + Zero + Add&lt;T, Output = T&gt;
     {
         <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&lt;N: MatrixNorm&lt;T, <span class="self">Self</span>&gt;&gt;(<span class="kw-2">&amp;</span><span class="self">self</span>, norm: N) -&gt; 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(&amp;b, Euclidean);
     ///
     /// assert_eq!(c, 2.0);
     /// # }
     /// ```
     </span><span class="kw">fn </span>metric&lt;<span class="lifetime">&#39;a</span>, <span class="lifetime">&#39;b</span>, B, M&gt;(<span class="kw-2">&amp;</span><span class="lifetime">&#39;a </span><span class="self">self</span>, mat: <span class="kw-2">&amp;</span><span class="lifetime">&#39;b </span>B, metric: M) -&gt; T
         <span class="kw">where </span>B: <span class="lifetime">&#39;b </span>+ BaseMatrix&lt;T&gt;,
               M: MatrixMetric&lt;<span class="lifetime">&#39;a</span>, <span class="lifetime">&#39;b</span>, T, <span class="self">Self</span>, B&gt;
     {
         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">&amp;</span><span class="self">self</span>) -&gt; T
         <span class="kw">where </span>T: Copy + Zero + Add&lt;T, Output = T&gt;
     {
         <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">&amp;</span><span class="self">self</span>, axis: Axes) -&gt; Vector&lt;T&gt;
         <span class="kw">where </span>T: Copy + PartialOrd
     {
         <span class="kw">match </span>axis {
             Axes::Col =&gt; {
                 <span class="kw">let </span><span class="kw-2">mut </span>mins: Vec&lt;T&gt; = 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">&amp;</span>v| <span class="kw">if </span>v &lt; m { v } <span class="kw">else </span>{ m });
                     mins.push(min);
                 }
                 Vector::new(mins)
             }
             Axes::Row =&gt; {
                 <span class="kw">let </span><span class="kw-2">mut </span>mins: Vec&lt;T&gt; = <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">&amp;mut </span>mins, row.raw_slice(), |min, <span class="kw-2">&amp;</span>r| <span class="kw">if </span>r &lt; <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">&amp;</span><span class="self">self</span>, axis: Axes) -&gt; Vector&lt;T&gt;
         <span class="kw">where </span>T: Copy + PartialOrd
     {
         <span class="kw">match </span>axis {
             Axes::Col =&gt; {
                 <span class="kw">let </span><span class="kw-2">mut </span>maxs: Vec&lt;T&gt; = 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">&amp;</span>v| <span class="kw">if </span>v &gt; m { v } <span class="kw">else </span>{ m });
                     maxs.push(max);
                 }
                 Vector::new(maxs)
             }
             Axes::Row =&gt; {
                 <span class="kw">let </span><span class="kw-2">mut </span>maxs: Vec&lt;T&gt; = <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">&amp;mut </span>maxs, row.raw_slice(), |max, <span class="kw-2">&amp;</span>r| <span class="kw">if </span>r &gt; <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>) -&gt; Matrix&lt;T&gt;
         <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::&lt;f64&gt;::ones(3,3);
     ///
     /// let b = &amp;a.select_rows(&amp;[2]);
     /// assert_eq!(b.rows(), 1);
     /// assert_eq!(b.cols(), 3);
     ///
     /// let c = &amp;a.select_rows(&amp;[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&lt;<span class="lifetime">&#39;a</span>, I&gt;(<span class="kw-2">&amp;</span><span class="self">self</span>, rows: I) -&gt; Matrix&lt;T&gt;
         <span class="kw">where </span>T: Copy,
               I: IntoIterator&lt;Item = <span class="kw-2">&amp;</span><span class="lifetime">&#39;a </span>usize&gt;,
               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 &lt; <span class="self">self</span>.rows(),
                     <span class="string">&quot;Row index is greater than number of rows.&quot;</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::&lt;f64&gt;::ones(3,3);
     /// let b = &amp;a.select_cols(&amp;[2]);
     /// assert_eq!(b.rows(), 3);
     /// assert_eq!(b.cols(), 1);
     ///
     /// let c = &amp;a.select_cols(&amp;[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&lt;<span class="lifetime">&#39;a</span>, I&gt;(<span class="kw-2">&amp;</span><span class="self">self</span>, cols: I) -&gt; Matrix&lt;T&gt;
         <span class="kw">where </span>T: Copy,
               I: IntoIterator&lt;Item = <span class="kw-2">&amp;</span><span class="lifetime">&#39;a </span>usize&gt;,
               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 &lt; <span class="self">self</span>.cols(),
                     <span class="string">&quot;Column index is greater than number of columns.&quot;</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 = &amp;a.elemul(&amp;b);
     /// assert_matrix_eq!(c, &amp;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">&amp;</span><span class="self">self</span>, m: <span class="kw-2">&amp;</span><span class="self">Self</span>) -&gt; Matrix&lt;T&gt;
         <span class="kw">where </span>T: Copy + Mul&lt;T, Output = T&gt;
     {
         <span class="macro">assert!</span>(<span class="self">self</span>.rows() == m.rows(), <span class="string">&quot;Matrix row counts not equal.&quot;</span>);
         <span class="macro">assert!</span>(<span class="self">self</span>.cols() == m.cols(), <span class="string">&quot;Matrix column counts not equal.&quot;</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">&amp;</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 = &amp;a.elediv(&amp;b);
     /// assert_matrix_eq!(c, &amp;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">&amp;</span><span class="self">self</span>, m: <span class="kw-2">&amp;</span><span class="self">Self</span>) -&gt; Matrix&lt;T&gt;
         <span class="kw">where </span>T: Copy + Div&lt;T, Output = T&gt;
     {
         <span class="macro">assert!</span>(<span class="self">self</span>.rows() == m.rows(), <span class="string">&quot;Matrix row counts not equal.&quot;</span>);
         <span class="macro">assert!</span>(<span class="self">self</span>.cols() == m.cols(), <span class="string">&quot;Matrix column counts not equal.&quot;</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">&amp;</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::&lt;f64&gt;::identity(3);
     /// let b = &amp;a.select(&amp;[0,1], &amp;[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!(&quot;{:?}&quot;, b.data());
     /// ```
     ///
     /// # Panics
     ///
     /// - Panics if row or column indices exceed the matrix dimensions.
     </span><span class="kw">fn </span>select(<span class="kw-2">&amp;</span><span class="self">self</span>, rows: <span class="kw-2">&amp;</span>[usize], cols: <span class="kw-2">&amp;</span>[usize]) -&gt; Matrix&lt;T&gt;
         <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 &lt; <span class="self">self</span>.cols(),
                     <span class="string">&quot;Column index is greater than number of columns.&quot;</span>);
         }

         <span class="kw">for </span>row <span class="kw">in </span>rows {
             <span class="macro">assert!</span>(<span class="kw-2">*</span>row &lt; <span class="self">self</span>.rows(),
                     <span class="string">&quot;Row index is greater than number of columns.&quot;</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 = &amp;a.hcat(&amp;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&lt;S&gt;(<span class="kw-2">&amp;</span><span class="self">self</span>, m: <span class="kw-2">&amp;</span>S) -&gt; Matrix&lt;T&gt;
         <span class="kw">where </span>T: Copy,
               S: BaseMatrix&lt;T&gt;
     {
         <span class="macro">assert!</span>(<span class="self">self</span>.rows() == m.rows(), <span class="string">&quot;Matrix row counts are not equal.&quot;</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 = &amp;a.vcat(&amp;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&lt;S&gt;(<span class="kw-2">&amp;</span><span class="self">self</span>, m: <span class="kw-2">&amp;</span>S) -&gt; Matrix&lt;T&gt;
         <span class="kw">where </span>T: Copy,
               S: BaseMatrix&lt;T&gt;
     {
         <span class="macro">assert!</span>(<span class="self">self</span>.cols() == m.cols(),
                 <span class="string">&quot;Matrix column counts are not equal.&quot;</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::&lt;Vec&lt;_&gt;&gt;();
     /// let b = matrix![1, 2;
     ///                 3, 4;
     ///                 5, 6].diag().cloned().collect::&lt;Vec&lt;_&gt;&gt;();
     ///
     /// assert_eq!(a, vec![1, 5, 9]);
     /// assert_eq!(b, vec![1, 4]);
     /// # }
     /// ```
     </span><span class="kw">fn </span>diag(<span class="kw-2">&amp;</span><span class="self">self</span>) -&gt; Diagonal&lt;T, <span class="self">Self</span>&gt; {
         <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">&amp;</span><span class="self">self</span>) -&gt; Matrix&lt;T&gt;
         <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">&amp;</span><span class="self">self</span>) -&gt; 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">&amp;</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(&quot;A solution should exist!&quot;);
     /// assert!((x[0] - 1.0) &lt; f32::EPSILON);
     /// assert!((x[1] - 1.0) &lt; 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">&amp;</span><span class="self">self</span>, y: Vector&lt;T&gt;) -&gt; <span class="prelude-ty">Result</span>&lt;Vector&lt;T&gt;, Error&gt;
         <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">&quot;Vector size {0} != {1} Matrix column count.&quot;</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(&quot;A solution should exist!&quot;);
     /// println!(&quot;{:?}&quot;, x);
     /// assert!((x[0] - 1.0) &lt; f32::EPSILON);
     /// assert!((x[1] - 1.0) &lt; 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">&amp;</span><span class="self">self</span>, y: Vector&lt;T&gt;) -&gt; <span class="prelude-ty">Result</span>&lt;Vector&lt;T&gt;, Error&gt;
         <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">&quot;Vector size {0} != {1} Matrix column count.&quot;</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">&amp;</span><span class="self">self</span>, mid: usize, axis: Axes) -&gt; (MatrixSlice&lt;T&gt;, MatrixSlice&lt;T&gt;) {
         <span class="kw">let </span>slice_1: MatrixSlice&lt;T&gt;;
         <span class="kw">let </span>slice_2: MatrixSlice&lt;T&gt;;

         <span class="kw">match </span>axis {
             Axes::Row =&gt; {
                 <span class="macro">assert!</span>(mid &lt; <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 =&gt; {
                 <span class="macro">assert!</span>(mid &lt; <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::&lt;Vec&lt;usize&gt;&gt;());
     /// let slice = MatrixSlice::from_matrix(&amp;a, [1,1], 2, 2);
     /// let new_slice = slice.sub_slice([0,0], 1, 1);
     /// ```
     </span><span class="kw">fn </span>sub_slice&lt;<span class="lifetime">&#39;a</span>&gt;(<span class="kw-2">&amp;</span><span class="self">self</span>, start: [usize; <span class="number">2</span>], rows: usize, cols: usize) -&gt; MatrixSlice&lt;<span class="lifetime">&#39;a</span>, T&gt;
         <span class="kw">where </span>T: <span class="lifetime">&#39;a
     </span>{
         <span class="macro">assert!</span>(start[<span class="number">0</span>] + rows &lt;= <span class="self">self</span>.rows(),
                 <span class="string">&quot;View dimensions exceed matrix dimensions.&quot;</span>);
         <span class="macro">assert!</span>(start[<span class="number">1</span>] + cols &lt;= <span class="self">self</span>.cols(),
                 <span class="string">&quot;View dimensions exceed matrix dimensions.&quot;</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&lt;T&gt;: BaseMatrix&lt;T&gt; {
     <span class="doccomment">/// Top left index of the slice.
     </span><span class="kw">fn </span>as_mut_ptr(<span class="kw-2">&amp;mut </span><span class="self">self</span>) -&gt; <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">&amp;mut </span><span class="self">self</span>) -&gt; MatrixSliceMut&lt;T&gt; {
         <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">&amp;mut </span><span class="self">self</span>, index: [usize; <span class="number">2</span>]) -&gt; <span class="kw-2">&amp;mut </span>T {
         <span class="kw-2">&amp;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">&amp;mut </span><span class="self">self</span>, index: [usize; <span class="number">2</span>]) -&gt; <span class="prelude-ty">Option</span>&lt;<span class="kw-2">&amp;mut </span>T&gt; {
         <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 &gt;= <span class="self">self</span>.rows() || col_ind &gt;= <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::&lt;Vec&lt;usize&gt;&gt;());
     ///
     /// {
     ///     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&lt;<span class="lifetime">&#39;a</span>&gt;(<span class="kw-2">&amp;mut </span><span class="self">self</span>) -&gt; SliceIterMut&lt;<span class="lifetime">&#39;a</span>, T&gt;
         <span class="kw">where </span>T: <span class="lifetime">&#39;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::&lt;<span class="kw-2">&amp;mut </span>T&gt;,
         }
     }

     <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">&amp;mut </span><span class="self">self</span>, index: usize) -&gt; ColumnMut&lt;T&gt; {
         <span class="kw">if </span>index &lt; <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">&quot;Column index out of bounds.&quot;</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">&amp;mut </span><span class="self">self</span>, index: usize) -&gt; ColumnMut&lt;T&gt; {
         <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">&amp;mut </span><span class="self">self</span>, index: usize) -&gt; RowMut&lt;T&gt; {
         <span class="kw">if </span>index &lt; <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">&quot;Row index out of bounds.&quot;</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">&amp;mut </span><span class="self">self</span>, index: usize) -&gt; RowMut&lt;T&gt; {
         <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">&amp;mut </span><span class="self">self</span>, a: usize, b: usize) {
         <span class="macro">assert!</span>(a &lt; <span class="self">self</span>.rows(),
                 <span class="macro">format!</span>(<span class="string">&quot;Row index {0} larger than row count {1}&quot;</span>, a, <span class="self">self</span>.rows()));
         <span class="macro">assert!</span>(b &lt; <span class="self">self</span>.rows(),
                 <span class="macro">format!</span>(<span class="string">&quot;Row index {0} larger than row count {1}&quot;</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">&amp;mut </span><span class="self">self</span>, a: usize, b: usize) {
         <span class="macro">assert!</span>(a &lt; <span class="self">self</span>.cols(),
                 <span class="macro">format!</span>(<span class="string">&quot;Row index {0} larger than row count {1}&quot;</span>, a, <span class="self">self</span>.rows()));
         <span class="macro">assert!</span>(b &lt; <span class="self">self</span>.cols(),
                 <span class="macro">format!</span>(<span class="string">&quot;Row index {0} larger than row count {1}&quot;</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!(&quot;{}&quot;, a);
     /// # }
     /// ```
     </span><span class="kw">fn </span>col_iter_mut(<span class="kw-2">&amp;mut </span><span class="self">self</span>) -&gt; ColsMut&lt;T&gt; {
         ColsMut {
             _marker: PhantomData::&lt;<span class="kw-2">&amp;mut </span>T&gt;,
             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!(&quot;{}&quot;, a);
     /// # }
     /// ```
     </span><span class="kw">fn </span>row_iter_mut(<span class="kw-2">&amp;mut </span><span class="self">self</span>) -&gt; RowsMut&lt;T&gt; {
         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::&lt;<span class="kw-2">&amp;mut </span>T&gt;,
         }
     }

     <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!(&quot;{}&quot;, 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">&amp;mut </span><span class="self">self</span>, k: DiagOffset) -&gt; DiagonalMut&lt;T, <span class="self">Self</span>&gt; {
         <span class="kw">let </span>(diag_len, diag_start) = <span class="kw">match </span>k.into() {
             DiagOffset::Main =&gt; (min(<span class="self">self</span>.rows(), <span class="self">self</span>.cols()), <span class="number">0</span>),
             DiagOffset::Above(m) =&gt; {
                 <span class="macro">assert!</span>(m &lt; <span class="self">self</span>.cols(),
                         <span class="string">&quot;Offset diagonal is not within matrix dimensions.&quot;</span>);
                 (min(<span class="self">self</span>.rows(), <span class="self">self</span>.cols() - m), m)
             }
             DiagOffset::Below(m) =&gt; {
                 <span class="macro">assert!</span>(m &lt; <span class="self">self</span>.rows(),
                         <span class="string">&quot;Offset diagonal is not within matrix dimensions.&quot;</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::&lt;<span class="kw-2">&amp;mut </span>T&gt;,
         }
     }

     <span class="doccomment">/// Sets the underlying matrix data to the target data.
     ///
     /// # Examples
     ///
     /// ```
     /// use rulinalg::matrix::{Matrix, BaseMatrixMut};
     ///
     /// let mut mat = Matrix::&lt;f32&gt;::zeros(4,4);
     /// let one_block = Matrix::&lt;f32&gt;::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&lt;M: BaseMatrix&lt;T&gt;&gt;(<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">&quot;Target has different row count to self.&quot;</span>);
         <span class="macro">assert!</span>(<span class="self">self</span>.cols() == target.cols(),
                 <span class="string">&quot;Target has different column count to self.&quot;</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">&amp;</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) -&gt; f64 {
     ///     a + 2f64
     /// }
     ///
     /// let a = Matrix::new(2, 2, vec![0.;4]);
     ///
     /// let b = a.apply(&amp;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">&amp;</span>Fn(T) -&gt; T) -&gt; <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">&amp;mut </span><span class="self">self</span>, mid: usize, axis: Axes) -&gt; (MatrixSliceMut&lt;T&gt;, MatrixSliceMut&lt;T&gt;) {

         <span class="kw">let </span>slice_1: MatrixSliceMut&lt;T&gt;;
         <span class="kw">let </span>slice_2: MatrixSliceMut&lt;T&gt;;

         <span class="kw">match </span>axis {
             Axes::Row =&gt; {
                 <span class="macro">assert!</span>(mid &lt; <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 =&gt; {
                 <span class="macro">assert!</span>(mid &lt; <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::&lt;Vec&lt;usize&gt;&gt;());
     /// let mut slice = MatrixSliceMut::from_matrix(&amp;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&lt;<span class="lifetime">&#39;a</span>&gt;(<span class="kw-2">&amp;mut </span><span class="self">self</span>,
                          start: [usize; <span class="number">2</span>],
                          rows: usize,
                          cols: usize)
                          -&gt; MatrixSliceMut&lt;<span class="lifetime">&#39;a</span>, T&gt;
         <span class="kw">where </span>T: <span class="lifetime">&#39;a
     </span>{
         <span class="macro">assert!</span>(start[<span class="number">0</span>] + rows &lt;= <span class="self">self</span>.rows(),
                 <span class="string">&quot;View dimensions exceed matrix dimensions.&quot;</span>);
         <span class="macro">assert!</span>(start[<span class="number">1</span>] + cols &lt;= <span class="self">self</span>.cols(),
                 <span class="string">&quot;View dimensions exceed matrix dimensions.&quot;</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>