EDA365电子论坛网

标题: 相同计算公式，为什么某一句耗时特别长 [打印本页]
作者: Lionoon    时间: 2020-2-26 13:26
标题: 相同计算公式，为什么某一句耗时特别长
是一个用来计算颗粒集体内应力分布的程序，其中出问题的是这个循环，用来计算内应力矢量在XY方向上的分布
: v0 p) o* ^% r! f* ^9 P. X
; H, e7 _; b* O
  • for i=1:sandnum
  •     for j=1:sandnum
  •         h=((spp(i,1)-spp(j,1))^2+(spp(i,2)-spp(j,2))^2)^0.5;
  •         r1=max(spp(i,3),spp(j,3));
  •         r2=min(spp(i,3),spp(j,3));
  •         a=asin(r2*sin(angel/2)/r1)*2;
  •         h=h-r2*cos(angel/2)-r1*cos(a/2);
  •         h=(h*2*r2*sin(angel/2)-(0.5*angel*r2^2-0.5*sin(angel)*r2^2)-(0.5*a*r1^2-0.5*sin(a)*r1^2))/(2*r2*sin(angel/2));
  •         rr=r2*sin(angel/2)*kcem;
  •         %A=2*rr/10^3;
  •         A=pi*rr^2/10^6;
  •         a(i,j)=e1(i,j,1)*E*A/h*10^6;%[N]x
  •         b(i,j)=e1(i,j,2)*E*A/h*10^6;%[N]y
  •         c(i,j)=e2(i,j,1)*G*A/h*10^6;%[T]x
  •         d(i,j)=e2(i,j,2)*G*A/h*10^6;%[T]y
  •     end
  • end
  • : J  l  ^: ]7 g: b7 t1 c, o
! Z* Q7 v5 L8 Y0 \6 L9 K
复制代码
- v/ M# ?. ]9 e' U% g4 b
' K0 \3 {- j1 T0 Y" `$ w# i6 m" I
用探查器优化的时候发现一个奇怪的现象：下图中的37行，也就是上面代码13行，a(i,j)那行耗时远大于其他。但是37-40行的计算公式明明是一样的，循环次数也一样。如果增加计算颗粒数量这个现象更明显，37行耗时能超过总运行时间的95%。各参数我都点开看了，不是因为有虚数。
# @6 B+ D5 V5 X# y3 N3 Z到现在也没搞清楚怎么回事 各位看看有没有知道啥情况的。我用的matlab2016a，电脑系统是64位win10
: `* S" m8 }4 E0 ~0 d
+ O6 C. ]) E0 E+ A* q) z5 y! d0 c2 O7 m) e

4 b/ L3 U' ^5 T* K7 V9 b. Y7 q4 f# c& O

# j1 c& y3 b( g5 C5 g2 h5 ~: K
/ I4 ]0 m8 s- B3 @5 f. S7 Y' F9 T
作者: ExxNEN    时间: 2020-2-27 16:03
后面是啥啊？



欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/) Powered by Discuz! X3.2