作者j0958322080 (Tidus)
看板MATLAB
标题[问题] 回圈数减少
时间Mon Dec 14 22:54:26 2015
大家好，最近遇到一个二维的问题想用MATLAB跑，
虽然C++已经OK，但是转到MATLAB好像会因为矩阵过大而运算过久，
目前已经把一些主要回圈用colon表示，但速度还是过慢，
https://drive.google.com/file/d/0B-jMD726BbgualpvS3ZMazc3Wjg/view
C++(最上面少一个 J 的回圈)
https://drive.google.com/file/d/0B-jMD726BbguMXFtOWNnb19XUTA/view
matlab，
执行中断好像主要是在计算H场那边，希望大家可以给点意见，感谢。
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 140.115.223.6
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/MATLAB/M.1450104871.A.561.html
1F：→ celestialgod: 回圈停止条件？ 12/14 23:11
3F：→ celestialgod: ez(1:IE,1:JE) = ga(1:IE,1:JE).*dz(1:IE,1:JE); 12/14 23:18
4F：→ celestialgod: 可以直接改成 ez = ga .* dz; 12/14 23:19
5F：→ celestialgod: 剩下看ic, jc是什麽... 12/14 23:19
6F：→ celestialgod: 没看到全貌也不知道可不可以直接vectorize计算 12/14 23:21
7F：→ j0958322080: 主要回圈就是上面那个C++档案的 12/14 23:56
※ 编辑: j0958322080 (140.115.223.6), 12/15/2015 00:05:51
8F：→ j0958322080: 另外nsteps是要算到50，只是我跑两个就无法了QQ 12/15 00:07
9F：→ j0958322080: 我觉得我应该是有那里陷入无穷回圈....... 12/15 00:33
10F：→ celestialgod: while nsteps <= 50 ? 12/15 01:34
11F：→ celestialgod: ic jc还是看不懂是啥 12/15 01:37
12F：→ celestialgod: ic, jc看到了 你的while停止条件没有设定阿= = 12/15 01:41
13F：→ celestialgod: nsteps = nsteps + 1; ?? 12/15 01:41
14F：→ j0958322080: 对应该就是跑到2就结束，只是这样好像结束不了 12/15 08:50
15F：→ j0958322080: 嗯我应该知道哪里错了 12/15 08:53
16F：→ j0958322080: 已经解决了，感谢各位 12/15 10:53