Octave的奇妙功能之一是它如何處理稀疏矩陣。我已經(jīng)完成了將一些Octave代碼復(fù)制到fortran中的任務(wù)（速度至關(guān)重要），具體來(lái)說(shuō)，鑒于稀疏矩陣Sparse（節(jié)點(diǎn)，節(jié)點(diǎn)），我需要提出子矩陣SubSparse（afew，afew）：r=4;c=6;k=5;nz=r*k;Sparse=spalloc(r,c,nz);Sparse(1,1)=1;Sparse(1,2)=2;Sparse(3,5)=4;Sparse(4,6)=6;SubSparse=Sparse([1,4],[1,6])完整的矩陣為：Sparse1.000 2.000 0.000 0.000 0.000 0.0000.000 0.000 0.000 0.000 0.000 0.0000.000 0.000 0.000 0.000 4.000 0.0000.000 0.000 0.000 0.000 0.000 6.000SubSparse1.000 0.0000.000 6.000我已經(jīng)讀過(guò)了，這個(gè)我不能做的一件事是重新創(chuàng)建完整的矩陣，因?yàn)樗鼈兲蠖鵁o(wú)法放入機(jī)器中。在稀疏矩陣的fortran實(shí)現(xiàn)中，我使用Hartwell-Boeing存儲(chǔ)，這是由我對(duì)SuperLU求解線性系統(tǒng)的要求所決定的。我相信八度的壓縮列格式非常相似。我的問(wèn)題是，在哪里可以找到（a）執(zhí)行MWE的第4行時(shí)發(fā)生的邏輯，以及（b）代碼