くだすれFORTRAN（超初心者用）その4
■ このスレッドは過去ログ倉庫に格納されています
0001デフォルトの名無しさん
2009/01/24(土) 18:32:01もしくは質問者自身何が何だか分からない質問を勇気を持って書き込むスレッドです。
FORTRAN使いが優しくコメントを返しますが、
お礼はFORTRANの布教と初心者の救済と次期Fortran2008規格でのCOMEFROM文採用をお願いします。
●注意事項
・質問する前にGoogle等の検索サイトで検索しましょう。
・回答者にわかりやすい様に、質問内容はできる限り詳しく書きましょう。
・エラーの場合は起きた状況、環境（OS・コンパイラ・バージョン）、エラーメッセージも詳しく書きましょう。
●前スレ
くだすれFORTRAN（超初心者用）その3
http://pc11.2ch.net/test/read.cgi/tech/1196384126/
くだすれFORTRAN（超初心者用）その2
http://pc11.2ch.net/test/read.cgi/tech/1164121236/
くだすれFORTRAN（超初心者用）
http://pc8.2ch.net/test/read.cgi/tech/1138063703/
●関連スレ
FORTRAN Ⅳ
http://pc11.2ch.net/test/read.cgi/tech/1163319215/
0002デフォルトの名無しさん
2009/01/24(土) 21:22:49スレ立て乙です
0003デフォルトの名無しさん
2009/01/25(日) 21:38:41まとめがあった。
ttp://www.ifremer.fr/ditigo/molagnon/fortran90/f90col.html
たとえば、以下のような場合
CALL( A ）
CALL( A(:) ) －＞ ALLOCATE( TMP(SIZE(A)) ) ; TMP = A ; CALL (TMP) ; A = TMP ;DEALLOCATE(TMP)
ソースの意味的には同じだが、下の書き方だと、コンパイラが一度テンポラリに
Aの中身をコピーしてテンポラリを読んだりする実装がされたりする。
最近は少し賢くなったようだｗ
0004デフォルトの名無しさん
2009/01/26(月) 01:32:560005デフォルトの名無しさん
2009/01/26(月) 06:57:57allocateして動的にメモリの空きを探して確保してコピーして、って
素で考えると手間がかかるだけな気がするけど、わざわざそうするということは
何か別な利得が発生するのかな？
00065
2009/01/26(月) 06:59:41a(:)が添字範囲全部の場合も、一部だけの場合でも両方対応するようにしたけど、
範囲全部の場合だけ別にする事を忘れただけなのかな？
0007デフォルトの名無しさん
2009/01/26(月) 12:03:14出ているコンパイラやOSを見る限り。
0008デフォルトの名無しさん
2009/01/26(月) 16:23:180009デフォルトの名無しさん
2009/01/27(火) 21:09:52if(条件式) then
endif
done
と
do I = 1 , N
if( ! 条件式) cycle
done
は、やってることは等価だけど、速さも同じ？
0010デフォルトの名無しさん
2009/01/28(水) 00:00:36まあ、コンパイラにもよるから自分で試してみてくれ。
0011デフォルトの名無しさん
2009/01/28(水) 01:53:01「ifortl」というコマンドを使うように指示されたのですがうまくいきません
どなたかわかる方いらっしゃいませんでしょうか？
環境は
OS Linux バージョン不明
コンパイラ Intel Fortaran 10.0
MKL Intel MKL9.1
コンパイラの購入時に「ifortl」コマンドを使うと自動でMKLのリンクが張られる～
のような説明を受けたらしいのですが、使用すると引用したモジュールが見つかりません
という旨のエラーメッセージがでます。
リンクオプションをきちんと記述するとコンパイルできます
MKLへのリンクの環境変数が通っていないんだと思うのですが、
設定方法がわからない状態です
0013デフォルトの名無しさん
2009/01/28(水) 02:49:33ウチのIntel Fortran 11.0 x86_64版には"ifortl"なるコマンドはないし、
ちょっとググってもそれらしい物はヒットしない。
とりあえず
file `which ifortl`
で ifortl の位置と種類を確認して、シェルスクリプトか何かなら
cat `which ifortl`
で中身を確認、書き換えられそうなら自力で修正
という感じで。
0014デフォルトの名無しさん
2009/01/28(水) 03:00:18「ifort」の間違いってオチじゃないよね…
0015デフォルトの名無しさん
2009/01/28(水) 09:59:39１３の通りwhichしてその内容をさらしてみるべし。
0016デフォルトの名無しさん
2009/01/28(水) 10:00:09001711
2009/01/28(水) 12:33:23ありがとうございます、無事解決できました。
>>15
まったくその通りでした
うっかり「ifortl」という隠しコマンドでもあるのかと思ってました。
リンクを張ったつもり～でしたが、実際は一部しか張ってなかったです
それにしてもなんというゆとりコマンド・・・
0018デフォルトの名無しさん
2009/01/28(水) 12:45:500019デフォルトの名無しさん
2009/01/28(水) 17:55:31FW=CB(2)-CB(1) //計算領域の幅を算出
CW=FW/NCX //幅をセル数で割ってセルの幅を算出
となる部分で、FWが1.33、NCXが200で、CWが6.65E-03となるはずが、
なぜか6.6500003E-3となっています。ちなみにFWはここで突如出現、
それ以外は特に型の宣言無くCOMMON文の中で定義されています。
何となく型の宣言が問題のような気がするんですが、原因と解決法
（ぴったり6.65E-3にしたい）を教えていただけませんか？
002019
2009/01/28(水) 18:05:010021デフォルトの名無しさん
2009/01/29(木) 01:01:07単精度なら有効桁からいって７～８桁目にゴミが出てもおかしくない。
ゆえに気にするな。
フォーマットを与えずに＊で書かせれば、自動的に精度の範囲内で出力するので
ゴミは見えなくなるだろう。フォーマットは浮動小数点での表現を単に何桁でも
十進法に直して書くだけだから、単精度でもF30.20とか書けば小数点以下２０桁とか
書いてくれるが、有効数字的には６～７個しか意味ない。
002219
2009/01/29(木) 09:49:23回答ありがとうございます。なるほど、７～８桁目は有効数字外なんですね。
実害があってですね、粒子の位置の判定ををするとき、粒子の位置Xを
セルの幅CWで割って、
LCX=(X/CW)+0.999999 //Xが3.5、CWが1だとすると4.499999、整数に代入して4で4番目のセルに入る
みたいな処理があるんです。このとき、CWがずれてるとセルナンバーが大きいところだと、
想定外の変なセルに入っちゃう粒子が増えてきます。でも有効数字の問題だとすると、倍精度にしたとしても、
頻度は減るにせよ、根本的解決ではないってことになりますね？やっぱりこの処理自体に問題があるんでしょうか？
問題があるとしたらどういう方法にするべきでしょうか？
0023デフォルトの名無しさん
2009/01/29(木) 10:17:31結果が変わるようなモデルが悪い。そんなの窓から投げ捨てちゃえw
002423
2009/01/29(木) 10:25:56変わらないだろうな。同じ問題は倍精度にしてもたぶん起こるよ。いまのCPUは単精度と倍精度で
ほとんど計算コストが変わらないので気が向いたら倍精度でやってみよう。
0025デフォルトの名無しさん
2009/01/29(木) 11:05:29ほんとに？！
普通に考えれば倍精度になればデータが倍になるから、スピードは半分になるような
その辺考えてあえて倍精度じゃなく単精度を使ってるコードもあるじゃん
0026デフォルトの名無しさん
2009/01/29(木) 11:33:48いまの計算機ってメモリ・CPU間のデータ転送をクレバーに
さばくからあまり目立っては違わない。
ほとんど変わらないとは思わないけどコードにもよるけど
単精度でしても２割速くなれば御の字じゃないかな。
メモリ節約して粒子数なり格子数を増やす方にメリットがあると
当然単精度でするだろうけど。仮に２倍速になるにしても、
高々２倍のスピードを求めてあえて単精度で計算する、
というのは何か間違っていると思うよ。
0027デフォルトの名無しさん
2009/01/29(木) 12:04:26単精度ですることは良くある
たくさんの測定データを読み込んで処理するときとかは特に
測定値だと元から倍精度の精度はないし
0028デフォルトの名無しさん
2009/01/29(木) 13:25:25はじめに一言。
四捨五入をしたいだけなら、NINTという関数があるのでそれを使うとよい。
ただNINTを使っても、十進表記にこだわると、今回の問題は解決しないかもしれない。
その理由を簡単に言うと以下のようになる。
Fortranでいうところの実数は、数学的な意味での実数ではない。
いわゆる浮動小数であって、有限個の有理数点で実数を代表させている。
ある種の区間をあらわす代表の数である。
その区間は、二進数表記で区切りのいい点で切っているので、十進数表記では
切りのいい点になっていない。
十進数で四捨五入しようとしても、そもそも閾値に幅があるのだから、その幅分の誤差が
出てもしょうがない。単精度ではそれが、有効数字７～８桁目にくる。
数直線をかみそりで切り分けるイメージではなく、木材をのこぎりで切るイメージで
考えればおｋ。切るとき、おがくずが出て、切った材木にはのこぎりの歯の幅だけの誤差がでる。
それは避けられない。
0029デフォルトの名無しさん
2009/01/29(木) 15:08:31クレバーってどういう意味？データ転送のレイテンシ隠蔽をがんがってくれるって意味？
普通のアプリでは役に立つけど、数値計算のようなデータを単純にガバガバ読むプログラムじゃ
そういうテクの恩恵はあんまないでしょ。
メモリから演算器へ転送するデータ量が単純に倍になるんだから・・・
0030デフォルトの名無しさん
2009/01/29(木) 15:30:11ベクトル化をすごく気にしたコードとか
一つのループを短くしたコードだとデータ転送の頻度がたぶん
すごく高いから２９の通り単精度倍精度の違いは大きいだろう。
ループを細切れにせずに一つのループ内でながめな処理をするタイプだと
計算中に次に使うデータを取ってきたり更新の済んだ分をメモリに戻す暇が
出来そうだ。
0031デフォルトの名無しさん
2009/01/29(木) 15:56:06簡単に見えて意外に厄介
0032デフォルトの名無しさん
2009/01/29(木) 16:36:30んだ。超厄介
>>30
キャッシュ内ですっぽり収まる程度の長さのループなら
ループの頭でlocalな変数に大配列の一部を複写してからの方が
編集やバグ探しが楽になるからまあ良いんだけど。
いまはだいぶ楽ちんだけど、ちょっと前はある長さを境にして
ガクンと遅くなることがあったな。他人のコードってあんまり見た
こと無いから判らないけど,ループは長いor短い方が良い、
とかtemporalな変数にうつすor大きい配列をガンガン右辺に持ってこよう、
とか、今の流行というか推奨はどっちなんだろう。
003319
2009/01/29(木) 17:45:46見てみることにします。情報を小出しにしちゃって申し訳ないんですが、
一番問題なのは、上の処理をすると、領域ぎりぎりに粒子があったとき、
５０までしかセルがないのに５１セル目を見ようとしちゃうことですね。
このエラー自体が不可避ならとりあえずそこの例外処理だけは入れといたほうがよさそうですね。
いろいろありがとうございます。上記件について何かいいアイデアあればぜひ。
0034デフォルトの名無しさん
2009/01/29(木) 23:00:06５０セルがMaxなら、いったん１００００倍なり拡大して整数に入れて、その整数の世界で
切り分ければいい。
整数がオーバーフローすると問題だが、５０セルがMaxなら大丈夫だろう。
0035デフォルトの名無しさん
2009/01/30(金) 00:07:25拡大係数行列Pに対して、部分ぴぼっと選択と行基本変形を繰り返して、
ガウスの消去法をするプログラムを書きました。
ここでPのある縦列が全て0のとき、たとえぴぼっと選択していても、
(＊)のところで0の割算が出てしまいます。
これを回避するうまい方法or常套手段はありますか？
003635
2009/01/30(金) 00:07:53do j = 1,sizeN
ValMax = 0.0d0
ImaxHeld = 0
do k = j, sizeN
if( ValMax < abs( P(j,k) ) )then
ImaxHeld = k
ValMax = abs( P(j,k) )
end if
end do
write(*,*) 'ImaxHeld=', ImaxHeld, 'ValMax=', ValMax
! Two rows are exchanged so that the maximum factor is placed at a diagonal
tmp(:) = P(:,j)
P(:,j) = P(:,ImaxHeld)
P(:,ImaxHeld) = tmp(:)
! Now the diagonal factor P(j,j) is a pibot
P(:,j) = P(:,j) / P(j,j) ! <=(＊)
do i = 1, sizeN
if (i /= j) then
P(:,i) = P(:,i) - P(j,i)*P(:,j)
end if
end do
end do
0037デフォルトの名無しさん
2009/01/30(金) 01:55:43一列全部０なら、シンギュラーな行列なんだから、そもそも解が無いだろ！
線形代数の教科書を１００回読み直せ！ｗ
0038デフォルトの名無しさん
2009/01/30(金) 07:23:57003935
2009/01/30(金) 11:25:16『行列式も0になって、クラメールも使えなくてどうしよう』などと考えてしまいました。
やりたいことは最小二乗法の正規方程式を解くことなのですが、問題は
データの与えかたによって、解が有ったり無かったりする
(列が全部0になったり、ならなかったりする)事です。
こんなとき、システマティックにやる方法はありますか？
004035
2009/01/30(金) 11:31:51間違ってsin(0)を与えていたために係数行列の列が全て0になっていました。
>>37さんありがとうございました。
0041デフォルトの名無しさん
2009/02/07(土) 00:28:19ついでにFortran2008にも部分対応らしいな。
といっても、CoArrayはもともとCRAYがやっていたので当然なんだが（笑
0042デフォルトの名無しさん
2009/02/07(土) 02:48:210043デフォルトの名無しさん
2009/02/07(土) 11:42:24Crayはここ１０年くらいのアメリカの護主義のおかげで復活してる。
一時期の低迷は脱している。
そういえば、東工大のつばめに刺さっていたクリアスピード社が、イギリスの経済壊滅の
余波で株価急落、資金繰りにもつまってきたというニュースが先週流れてきてた。
ETA１０の夢ふたたびｗ
0044デフォルトの名無しさん
2009/02/07(土) 16:20:54つーかCrayってもう石は作ってないんだよな？
その他大勢の単なるシステム屋さんに成り下がっちまった
0045デフォルトの名無しさん
2009/02/07(土) 22:16:59アセンブルコードを吐かせてサイズを比較しろ
0046デフォルトの名無しさん
2009/02/08(日) 07:53:42１CPUなら後者が速いことが多い・・・とおもう。
並列だとたぶん前者・・・とおもう。並列処理すると代償でマスキング処理が抑制されることがあるから。
個人的にはif-endifで処理しない・する範囲を明示する方が好き。
0047デフォルトの名無しさん
2009/02/08(日) 16:14:46違ったらごめん
0048デフォルトの名無しさん
2009/02/08(日) 18:01:42先輩からいただいたプログラムのファイルは2つあり、ひとつが本プログラム、もうひとつがサブルーチンプログラムです。
この状況で、もうひとつのサブルーチンプログラムをworkspaceのところで作成し、compile、excecutionをしたのですが、
どうもerrorなどはないのですが、write文などの結果が既存のサブルーチンプログラムからの結果が出力され、
新しいプログラムでの結果が出力できません。どうやったら、新しいプログラムの結果が
出力できるのか、教えていただけますでしょうか。
大変、初心者の質問で申し訳ございません。
0049デフォルトの名無しさん
2009/02/08(日) 18:08:07無視してください。また、お世話になるかと思いますが、どうぞよろしくお願いいたします。
0050デフォルトの名無しさん
2009/02/08(日) 19:23:15あってるよ。
でも並列処理するときに誤作動防止で
ベクトル処理を抑制することがあるというのはありそうではある
0051デフォルトの名無しさん
2009/02/10(火) 15:42:39Matlabから移行してきたんですが、公式マニュアル見てもプログラミング関係の単語がわからずで。ラッパーとか。
数値計算に関しては、噂に聞いてたよりも使いやすいなあ、と思ってるんですが。
よろしくお願いします
0052デフォルトの名無しさん
2009/02/10(火) 17:32:15gfortran -i8 '-mcmodel=medium' -fdefault-integer-8 -O2 -mtune=opteron -m64 -c hoge.f
In file hoge.f
C(1) = LShift(I1,32).or.(I2.and.NM32)
1
Error: Operands of logical operator '.and.' at (1) are INTEGER(8)/INTEGER(8)
似たような問題を見つけましたが、対処方法がよく分かりません。
ttp://gcc.gnu.org/ml/fortran/2005-07/msg00416.html
どなたか教えて頂けないでしょうか。
ENV: RedHat 4, GCC 4.1.2
0053デフォルトの名無しさん
2009/02/11(水) 00:36:13何をしたいのか、それだけでは意味不明だが、『hoge』などというUNIX系の方言を
しゃべるところを見ると、C言語とかと混同して論理演算とビット操作を勘違いしていると邪推？
＞ C(1) = LShift(I1,32).or.(I2.and.NM32)
（LShift という関数は、Fortranの標準関数には存在しない。が、ベンダー拡張の
ビット操作演算子としてよく存在しているので、ビット操作演算子だと仮定する。）
これは、論理演算とビット操作と論理演算を混用しているので、コンパイラが怒っていると
思われる。
正しくは、
INTEGER(8) :: i1, i2, nm32
INTEGER(8) :: ic(整数)
ic(1) = IOR( ISHFT(i1, 32), IAND(i2, nm32) )
Fortranには、論理型と整数型の区別がある。 .and. や .or. は論理演算子であって
真偽の二値しかもたない論理型の変数（あるいは定数）に作用する。
一方、整数型の変数（や定数）に対するビット演算子も存在している。
それは二引数をとる関数で IAND や IOR、あるいはビットシフト演算子 ISHFT である。
ｃを整数として使うのは、暗黙の型に反するので、どうしてもと言う理由がなければ
やめたほうがいいと思う。
なお、Fortranの論理変数の内部表現は規格では不定であって、処理系によってことなる。
一般的にLSBが１か０かで真偽を区別していることが多い。
C言語のように０が偽で、その他が真というわけでない。
よって、論理演算を行ってもLSBしか変化していなかったりする。
しりたければEQUIVALENCEをかけるか、TRANSFER関数で内部表現を見れば分かる。
0054デフォルトの名無しさん
2009/02/11(水) 10:15:01的確かつ詳細な回答有り難う御座います。大変参考になりました。
お察しの通りC言語の経験者ですが、Fortranは皆無です。
ご指摘を参考に勉強致します。有り難う御座いました:-)
0055デフォルトの名無しさん
2009/02/15(日) 15:04:380056デフォルトの名無しさん
2009/02/16(月) 22:01:19Vistaに対応した無料FORTRANが出来るようになるサイトを教えてください。
申し訳ありませんが、よろしくお願いします。
0057デフォルトの名無しさん
2009/02/16(月) 23:42:15http://pc11.2ch.net/test/read.cgi/tech/1163319215/2
0058デフォルトの名無しさん
2009/02/17(火) 11:11:13俺的にこれがお勧め
ttp://www.download.com/Silverfrost-FTN95/3000-2069_4-10491439.html
Plato3と言うソフトが一緒にインストールされるから、それを起動して
File->Newもしくは左上の白いアイコンをクリックするとNew Fileと言うウィンドウが出る
FORTRAN77のコードを書きたいならFixed Format Fortran fileを
Fortran90/95のコードを書くならFree Fomat Forrtran fileを選択しOpenボタンを押す。
あとはコードを書いて、左上に青い再生ボタンがあるのでクリック、変なダイアログが出ても
OKおせばおｋ。あと、黒いウインドウを閉じる時は、右上の×じゃなくて、キーボードの
Enterを押さないと落ちる事があるので注意
0059デフォルトの名無しさん
2009/02/19(木) 17:28:260060デフォルトの名無しさん
2009/02/19(木) 18:30:000061デフォルトの名無しさん
2009/02/19(木) 19:56:36bug fixが中心だろう
Memory allocation error in for_array_copy_in()
とか
11.0 message catalog produces garbled messages with Japanese compiler
とかあったから、とりあえず入れといた
0062デフォルトの名無しさん
2009/02/19(木) 23:33:4111.0でコンパイルしなおしただけで計算時間が1/3になってワロタ
別に最適化オプション何も使ってないデフォルトなんだけどな
これでOpenMPで最適化したら一体どうなるんだ
0063デフォルトの名無しさん
2009/02/20(金) 04:02:32黒カエラーはまだ直ってないなｗ
誰もIntelに通報してないのか？
0064デフォルトの名無しさん
2009/02/20(金) 13:43:56７になった時にですごく良くなった気がする。
あとはあまり変化は大きくない気がする。
・・・といって今使ってるのが９のままなので10,11でまた進歩があったのかな？
SSE絡みでは順次良くなったそうだけど
明示的にそれを使うことは少ないからなあ。
0065デフォルトの名無しさん
2009/02/20(金) 14:23:19見違えるように良くなったって話だが
0066デフォルトの名無しさん
2009/02/21(土) 00:20:220067デフォルトの名無しさん
2009/02/21(土) 10:32:37みなさんはエラー系の出力はwrite(0,*)を使っていますか？
それとも、標準出力もエラー出力も全部print *を使っていますか？
0068デフォルトの名無しさん
2009/02/21(土) 19:42:57標準出力とエラー出力があるっていうのはUNIXの話だから、
メインフレーム上がりでUNIXをおもちゃと馬鹿にしているおっさんたちは、
そんなもの存在も知らず、相変わらず６番に出力していると思う。
0069デフォルトの名無しさん
2009/02/22(日) 00:56:09integer, allocatable :: imem(:)
msize = 1000*1000 * 200
call alloc4i(imem, msize)
c :
stop
end
subroutine alloc4i(imem, msize)
integer, allocatable :: imem(:)
if ( .not. allocated(imem) ) then
allocate( imem(msize), stat=ier )
if (ier .ne. 0) stop
else
stop
endif
return
end
こんな感じでメモリ割付をするサブルーチンを作りたいのですが、サブルーチン側の
配列宣言でコンパイルエラーになってしまいます。このように親側でallocatableにした
配列は、サブルーチン側でどのようにallocateするのでしょうか？
0071デフォルトの名無しさん
2009/02/22(日) 01:47:54外部サブルーチンだからエラーが出んじゃない？
内部サブルーチンやモジュールで使ってるけど
エラーは出ないよ
007269
2009/02/22(日) 02:58:23program memtest
integer, allocatable :: imem(:)
msize = 1000*1000 * 200
call alloc4i
stop
contains
subroutine alloc4i
if ( .not. allocated(imem) ) then
allocate ( imem(msize), stat=ier )
if (ier .ne. 0) stop
else
stop
endif
return
end subroutine
end program
レスどうもです。こんな感じでしょうか？
たしかにコンパイルは通りました。ただ、ライブラリ化
したいため、分割コンパイルしたいのです。何とか
外部プログラムで実現できないでしょうか・・・
0073デフォルトの名無しさん
2009/02/22(日) 04:39:47MODULEに入れるか、本体側でINTERFACEを書けばいいんでね？
ALLOCATABLE属性をつけても、７７式に裸にしておくと、引数がマッチしているのか分からんから
怒られるんだと思うが。
007469
2009/02/22(日) 09:30:19module malloc
integer, allocatable :: imem(:)
contains
subroutine alloc4i(imem, msize)
integer::imem(:)
integer::msize
allocate ( imem(msize), stat=ier )
if (ier .ne. 0) stop
return
end subroutine
end module
program memtest
use malloc
integer::msize
msize = 1000*1000 * 200
call alloc4i(imem, msize)
stop
end program
う～ん・・・ ｆ90に疎くて、こんなコードしか書けませんでしたorz
ホントは、下みたいな感じで、program側で複数の配列確保に使える
サブルーチンが欲しいのですが・・・
program memtest
：
call alloc4i(imem1,msize)
call alloc4i(imem2,msize)
：
■ このスレッドは過去ログ倉庫に格納されています