关于稀疏矩阵第二次赋值当机的问题

swinder

矩阵赋值
要得到一个大矩阵tem，他的行按某种条件可以分为两种epl与opl，列做同样的事有epr与opr，这个大矩阵是分块的。只有行属于epl同时列属于epr位置的元素，与行opl同时列opr 的元素非零，
计算得到对应两个小矩阵etem与otem,这时要把值赋给大矩阵
code如下：

tem=sparse(q^3,q^5*d^2);
tic
tem(opl,opr)=otem;clear otem opl opr;
toc
tic
tem(epl,epr)=etem;clear etem epl epr;
toc
测试时,第一个小矩阵赋值可以顺利完成,第二个就会死在那里,我等了半个小时也赋不完,但是如果把第二个小矩阵与第一个换位,它本身也可以顺利赋完.所以情况是：第一次赋值顺利，第二次赋值就会卡在那里。请问这是为什么。我应该怎样改。谢谢

lengyunfeng

1# swinder

问题可能出在你的otem 和etem上，一般来说它应该是和opl、opr以及epl、epr同维的向量，而在这里你的数据却是两个矩阵。