- 1 名前：デフォルトの名無しさん [2006/11/22(水) 00:00:36 ]
- このスレッドは、他のスレッドでは書き込めない超低レベル、
もしくは質問者自身何が何だか分からない質問を勇気を持って書き込むスレッドです。
FORTRAN使いが優しくコメントを返しますが、
お礼はFORTRANの布教と初心者の救済と次期Fortran2008規格でのCOMEFROM文採用をお願いします。
●注意事項
・質問する前にGoogle等の検索サイトで検索しましょう。
・回答者にわかりやすい様に、質問内容はできる限り詳しく書きましょう。
・エラーの場合は起きた状況、環境（OS・コンパイラ）、エラーメッセージも詳しく書きましょう。
●前スレ
くだすれFORTRAN（超初心者用）
pc8.2ch.net/test/read.cgi/tech/1138063703/
●関連スレ
FORTRAN Ⅳ
pc8.2ch.net/test/read.cgi/tech/1163319215/
- 293 名前：デフォルトの名無しさん [2007/02/25(日) 14:25:10 ]
- s = a**b
の計算が必要なときは
s = exp(b*log(a))
と書いた方が計算が速いと聞いたのですが，本当でしょうか？
- 294 名前：デフォルトの名無しさん mailto:sage [2007/02/25(日) 17:37:12 ]
- >>293
自分でベンチマーク作って確かめろｗ
どっちにしろ一般のべき乗はexp log で計算しているはずなので、
８０８７のようにハードウェア的に実装されたｅｘｐを利用したかしないかの違いだろう。
- 295 名前：293 mailto:sage [2007/02/27(火) 10:33:43 ]
- >>294
うちのPCではどうやら後者の方が2割程度早いようです．
ｺﾝﾊﾟｲﾗはIntel Visual Fortran 9.1です．
それでこれからは全てexp logで計算しようと思ったのですが，
問題が発生しました．
次のような，モジュールで変数を共有する仕様の
プログラムを書いたのですが，コメントしているように，
共有用モジュールの変数を介して値を渡した場合と
直接値を代入した場合で結果が異なってしまいます．
ちなみにb**bで計算した場合はどちらも同じ結果を出力します．
b=0の場合を除いて演算してもいいのですが，
NaNになる原因が分からないので何とも怖くて使えません．
何かアドバイスあったらお願いします．．
- 296 名前：293:プログラムリスト mailto:sage [2007/02/27(火) 10:34:24 ]
- module commons
implicit none
real(8) a(1:10)
end module commons
program main
use commons
implicit none
a(1)=0.d0
call sub
end program main
subroutine sub
use commons
implicit none
real(8) b
!こちらを使うと出力は NaN
b=a(1)
!こちらを使うと出力は 1
b=0.d0
write(6,*) exp(b*log(b))
return
end subroutine sub
- 297 名前：デフォルトの名無しさん mailto:sage [2007/02/27(火) 21:19:52 ]
- >>295
組み込み関数 LOG(X) はXが実数の場合、Xは正(X>0)でなければならない。
IVFでは -∞ を返しているが、b * log(b) = 0 * -∞ = NaN となるので、
exp(b * log(b)) = exp(NaN) = NaN となる。
b=0.d0 の時に 1 と出力されるのは、コンパイル時に bの値が既知の定数(=0)なので
最適化によりプログラム中で b*log(b) が計算されず、直接 exp0.0D0) を求めてるため。
最適化なし(-Od)でコンパイルすると両方とも NaN になる。
- 298 名前：293 mailto:sage [2007/02/27(火) 21:35:35 ]
- >>297
おおー！確かにそうなります！
実にすっきりしました．
どうもありがとうございます．
- 299 名前：293 mailto:sage [2007/02/27(火) 21:39:28 ]
- 途中で書き込んでしまいました．
あまりに訳が分からないので
さっきVinu Linuxをわざわざインストールしてgfortranでも
同じことをしたんですが，こちらは b=0.d0 と書いても b=a(1) で渡しても
NaNとなりました．
今まで最適化で問題が発生したことが無かったこともあって，
原因がコンパイラの最適化にあったとは全然気がつきませんでした．
- 300 名前：デフォルトの名無しさん mailto:sage [2007/02/27(火) 23:39:40 ]
- Xが0以下の実数の場合、LOG(X) の結果は処理系依存になるので
実行時エラーを出してプログラムを止める処理系もある。
- 301 名前：デフォルトの名無しさん mailto:sage [2007/02/28(水) 00:44:46 ]
- 最近の処理系はIEEE準拠になって途中で実行時エラーで止まることが少なくなってありがたい。
せっかくIF文で例外を排除しているのに、コンパイラが最適化して勝手に投機的実行をして、
挙句に０割とかで実行時エラーを起こしてとまるとか、ウンコのようなことが昔はよくあったｗｗｗｗ
- 302 名前：デフォルトの名無しさん [2007/02/28(水) 04:01:39 ]
- 学校で今週までの課題で
『ｎ元連立一次方程式をガウス・ザイデル法を利用して解くプログラムを作成し、テストを行え』
という課題がだされたのですが、全くわかりません(;_;)
頑張って本を読んでみたのですが理解できなくて(;_;)
最終手段として、ここに来ました(>_<)
どなたか、ご教授下されば幸いです(;O;)
初めての書き込みなので粗相があったらすみませんm(__)m
ちなみに学校で使っているのはフォートラン77です。
- 303 名前：デフォルトの名無しさん [2007/02/28(水) 12:05:08 ]
- またガウス・ザイデル法かｗ
こんなの教科書にも載りまくりだろう。
誰かFAQにでも入れておけ。
|
 |
