u2 = exp(i.*zeros(8192,8192));
for ite1 = 1:800
   for ite2 = 1:900
      u = exp(i.*sqrt(xy_2d- z_2d(ite1,ite2)));
      u2 = u2.*u;
   end
end

这个程序越来越慢，请问是什么原因？请问有什么办法加速或者让他匀速运行吗？
※ 修改:·novicer 于 Sep  9 10:01:02 2021 修改本文·[FROM: 143.198.64.*]
※ 来源:·水木社区 http://www.mysmth.net·[FROM: 49.66.217.*]

修改:novicer FROM 143.198.64.*
FROM 49.66.217.*

能稍微详细说下怎么该吗？多谢！

u2 = exp(i.*zeros(8192,8192));
for ite1 = 1:800
  for ite2 = 1:900
     u = exp(i.*sqrt(xy_2d- z_2d(ite1,ite2)));
     u2 = u2.*u;
  end
end


【 在 zszqzzzf 的大作中提到: 】
: prod
--
修改:novicer FROM 143.198.64.*
FROM 143.198.64.*

整个u2提到循环外，怎么保存u呢？u本身是一个800*900的矩阵，如果要保存，循环之后好像是8192*8192*800*900个元素


【 在 zszqzzzf 的大作中提到: 】
: u2整个提取到循环外，可以用prod，对矩阵求累积。
: 双重循环计算也是慢的。如果你能做到消灭两个for，改成矩阵运算，速度可以提高百倍。
: cellfun翻一翻看看能不能就取代ite1和ite2，……
: ...................
--
FROM 143.198.231.*

temp = zeros(8192,8192);
for ite = 1:720000
  temp = temp + i*sqrt(xy_2d - z_2d(ite));
end
u2 = exp(temp);

【 在 novicer 的大作中提到: 】
: u2 = exp(i.*zeros(8192,8192));
: for ite1 = 1:800
:    for ite2 = 1:900
: ...................
--
修改:shaolimin FROM 218.104.71.*
FROM 218.104.71.*

【 在 shaolimin (魏晋风度) 的大作中提到: 】
: 标  题: Re: 请问这个程序为什么越来越慢?
: 发信站: 水木社区 (Thu Sep  9 13:06:52 2021), 站内
:
: temp = zeros(8192,8192);
: for ite = 1:720000
:   temp = temp + i*sqrt(xy_2d - z_2d(ite));
: end
temp =arrayfun(@(ite) sum(i*sqrt(xy_2d-z2d(ite)),[1:720000]);

: u2 = exp(temp);
:
: 【 在 novicer 的大作中提到: 】
: : u2 = exp(i.*zeros(8192,8192));
: : for ite1 = 1:800
: :    for ite2 = 1:900
: : ...................
:
: --
: ※ 修改:·shaolimin 于 Sep  9 13:07:45 2021 修改本文·[FROM: 218.104.71.*]
: ※ 来源:·水木社区 http://www.mysmth.net·[FROM: 218.104.71.*]
--
修改:shaolimin FROM 218.104.71.*
FROM 112.47.163.*