只是preallocating Rup应该加快:
第一行:
Rup = 0(大小(A));
我如何删除for循环?
1视图(30天)
显示旧的评论
我想删除的循环代码我写了很久以前的事了。
函数(Rup) =更新器(A、E S)
[imax, jmax zmax] =大小(一个);
为z = 1: zmax
Stmp = S (:,:, setdiff (1: zmax, z));
Atmp = (:,:, setdiff (1: zmax, z));
Etmp = E (::,:,:, setdiff (1: zmax, z));
为我= 1:imax
为j = 1: jmax
sumtmp =排列(和(和(Etmp (i, j,::,:)、3)、4), (1 2 3 4 5));
β= (Stmp (i, j) + Atmp (i, j:) + sumtmp);
n = - max(β);
EEtmp = E (i, j,:,:, z);
Rup (i, j, z) = S (i, j, z) + n +总和(EEtmp (:));
结束
结束
结束
R, S是[N M x P]矩阵。
我自己试过很多,但是我没有自信与技巧有助于消除循环。谁能帮我吗?
马特奥
