| 查看: 1554 | 回复: 10 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
duwei789551新虫 (小有名气)
|
[求助]
求助，两个uds变量相除就显示发散是什么问题？ 已有1人参与
|
|
|
如题，需要计算保存为两个标量方程源项的比值。当只用这一段代码时 r=C_UDSI(c,c_thread,1)/C_UDSI(c,c_thread,0); 计算一开始就显示uds-1发散； 为避免发散出现，对分母做出限制 if (C_UDSI(c,c_thread,0)<=1.0) { r=0.00000001; } else { r=C_UDSI(c,c_thread,1)/C_UDSI(c,c_thread,0); } 这样的话算到途中，uds迭代曲线将要出现的时候显示uds-1发散； 不加r=C_UDSI(c,c_thread,1)/C_UDSI(c,c_thread,0);这句话就没有问题，不知是什么原因 谢谢解答！ |
回复此楼» 猜你喜欢
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有194人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
qc12345
至尊木虫 (职业作家)
- 应助: 146 (高中生)
- 金币: 13299.8
- 红花: 26
- 帖子: 3461
- 在线: 717.9小时
- 虫号: 3071487
- 注册: 2014-03-19
- 专业: 计算机科学的基础理论
|
real rr,ss,tt； ... rr=C_UDSI(c,c_thread,0)； ss=C_UDSI(c,c_thread,1)； if (rr<=1.0) { tt=0.00000001; } else { tt=ss/rr; } 這樣試試 |
7楼2016-10-13 15:11:43
qc12345
至尊木虫 (职业作家)
- 应助: 146 (高中生)
- 金币: 13299.8
- 红花: 26
- 帖子: 3461
- 在线: 717.9小时
- 虫号: 3071487
- 注册: 2014-03-19
- 专业: 计算机科学的基础理论
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
duwei789551: 金币+100, ★★★★★最佳答案, 解决问题 2016-10-13 20:05:11
感谢参与，应助指数 +1
duwei789551: 金币+100, ★★★★★最佳答案, 解决问题 2016-10-13 20:05:11
| 設s=C_UDSI(c,c_thread,0)再相除 |
» 本帖已获得的红花（最新10朵）
2楼2016-10-12 17:59:34
duwei789551
新虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 709.4
- 帖子: 65
- 在线: 39.5小时
- 虫号: 3854710
- 注册: 2015-05-07
- 性别: GG
- 专业: 多相流热物理学
送红花一朵 |
我的程序是这样的 标量方程1源相： real s； ... s=C_UDSI(c,c_thread,0； if (C_UDSI(c,c_thread,0)<=1.0) { r=0.00000001; } else { r=C_UDSI(c,c_thread,1)/s; } J=pow(rou,2.0)/roul； C_UDSI(c,c_thread,0)=J; C_UDMI(c,c_thread,0)=J; dS[eqn]=0.0; return J; 源相方程2 real s； ... s=C_UDSI(c,c_thread,0； if (C_UDSI(c,c_thread,0)<=1.0) { r=0.00000001; } else { r=C_UDSI(c,c_thread,1)/s; } ... yihao=C_UDMI(c,c_thread,0)*r+C_UDSI(c,c_thread,0)*rou； C_UDSI(c,c_thread,1)=yihao; C_UDMI(c,c_thread,5)=yihao; dS[eqn]=0.0; return yihao; 改写成这样还是算到迭代曲线出现时发散，不知道有什么问题 |
3楼2016-10-13 11:28:41
qc12345
至尊木虫 (职业作家)
- 应助: 146 (高中生)
- 金币: 13299.8
- 红花: 26
- 帖子: 3461
- 在线: 717.9小时
- 虫号: 3071487
- 注册: 2014-03-19
- 专业: 计算机科学的基础理论
4楼2016-10-13 14:00:53
投票:
10