看板 [ Fortran ]
討論串[問題] 迴圈拆開計算後與原本值有微小誤差
共 2 篇文章
內容預覽： 開啟 | 關閉 | 只限未讀
首頁
上一頁
1
下一頁
尾頁
#1
[問題] 迴圈拆開計算後與原本值有微小誤差
推噓2(2推 0噓 6→)留言8則，0人參與, 最新作者boa85391 ( )時間9年前 (2015/08/12 21:27), 編輯資訊
1
0
0
內容預覽:
do i=1, 3999, 1. do j=i+1, 4000, 1. rx=bead(j)%x_p - bead(i)%x_p. ry=bead(j)%y_p - bead(i)%y_p. rz=bead(j)%z_p - bead(i)%z_p. rpq= dsqrt(rx**2+ry**2+r
(還有3088個字)
#2
Re: [問題] 迴圈拆開計算後與原本值有微小誤差
推噓2(2推 0噓 1→)留言3則，0人參與, 最新作者latinboy (阿賢)時間9年前 (2015/08/13 10:51), 編輯資訊
0
0
0
內容預覽:
上述浮點數的運算為了避免誤差累積以及增加效率. 可以改成. r2 = 1.d0 / (rx*rx + ry*ry + rz*rz). r6 = r2 * r2 * r2. a2 = sigma * sigma. a6 = a2 * a2 * a2. ar6 = a6 * r6. rep_term =
(還有74個字)
首頁
上一頁
1
下一頁
尾頁