| 查看: 1105 | 回复: 5 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
破碎东瓜金虫 (初入文坛)
|
[求助]
谁能帮我把一段matlab程序转为fortran啊！！！（很短的一段）
|
|
|
while err>errlim s=find(P>0); sn=find(P<=0); dd=real(ifft2(fA.*fft2(P,n*2,n*2))); u=dd(1:n,1:n)*2.0*dx/pi^2; rk=h+u; do=mean(rk(s)); rk=rk-do; G=sum(sum(rk(s).*rk(s))); pk(s)=rk(s)+G/Gold*pk(s); pk(sn)=0; Gold=G; dd=real(ifft2(fA.*fft2(pk,2*n,2*n))); qk=dd(1:n,1:n)*2.0*dx/pi^2; rb=mean(qk(s)); qk=qk-rb; dp=sum(sum(rk(s).*pk(s)))/sum(sum(qk(s).*pk(s))); P(s)=P(s)-dp*pk(s); s=find(P<=0);P(s)=0; sol=find((P==0)&(rk<0)); P(sol)=P(sol)-dp*rk(sol); W=sum(sum(P))*dx*dx; P=charge/W*P; err=sqrt(Gold*dx*dx); end |
回复此楼» 猜你喜欢
调剂
已经有6人回复
-
280求调剂
已经有6人回复
-
302求调剂一志愿华中师范大学
已经有5人回复
-
346求调剂
已经有4人回复
-
材料专硕322分
已经有9人回复
-
372分材料与化工（085600）一志愿湖南大学求调剂
已经有5人回复
-
318求调剂，计算材料方向
已经有10人回复
-
310求调剂
已经有14人回复
-
315求调剂
已经有3人回复
-
265求调剂
已经有13人回复
破碎东瓜
金虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 1520.2
- 散金: 53
- 帖子: 46
- 在线: 32.1小时
- 虫号: 2231948
- 注册: 2013-01-09
- 性别: GG
- 专业: 机械摩擦学与表面技术
6楼2015-12-13 14:49:55
wxws.2008
木虫 (正式写手)
- 应助: 5 (幼儿园)
- 金币: 5015.5
- 散金: 256
- 红花: 9
- 帖子: 354
- 在线: 147.9小时
- 虫号: 639701
- 注册: 2008-10-29
- 性别: GG
- 专业: 荒漠化与水土保持
2楼2015-12-12 09:57:31
破碎东瓜
金虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 1520.2
- 散金: 53
- 帖子: 46
- 在线: 32.1小时
- 虫号: 2231948
- 注册: 2013-01-09
- 性别: GG
- 专业: 机械摩擦学与表面技术
3楼2015-12-12 15:36:58
wxws.2008
木虫 (正式写手)
- 应助: 5 (幼儿园)
- 金币: 5015.5
- 散金: 256
- 红花: 9
- 帖子: 354
- 在线: 147.9小时
- 虫号: 639701
- 注册: 2008-10-29
- 性别: GG
- 专业: 荒漠化与水土保持
4楼2015-12-12 18:07:54
