トップページ⇒tech
991コメント434KB

FORTRAN Ⅳ

■ このスレッドは過去ログ倉庫に格納されています
0001デフォルトの名無しさん2006/11/12(日) 17:13:35
FORTRAN総合スレッドです。
FORTRAN77/90/95/2003, HPF等、FORTRAN全般に関する話題を書き込んで下さい。

・FORTRAN関連情報
お約束だが、まずは自分で調べよう
http://www.google.co.jp/

x86で動作するFortranコンパイラ
http://www.nminoru.jp/~nminoru/programming/x86-fortran.html

The Fortran Company
http://www.fortran.com/
Polyhedron Software
http://www.polyhedron.co.uk/
Google Directory - Fortran
http://directory.google.com/Top/Computers/Programming/Languages/Fortran/

・前スレ
FORTRAN III
http://pc8.2ch.net/test/read.cgi/tech/1104724162/
FORTRAN II
http://pc5.2ch.net/test/read.cgi/tech/1068351911/
Fortran
http://pc2.2ch.net/test/read.cgi/tech/1003214738/
0478デフォルトの名無しさん2008/11/15(土) 00:55:13
>>476
lは１と見間違いやすいからあまり使わないが、貴重な整数型なのでもったいない。
0479デフォルトの名無しさん2008/11/15(土) 05:22:43
>>478
Lは大文字で使う
0480デフォルトの名無しさん2008/11/16(日) 00:36:10
>>479
ボーランド記法？
昔PASCALに散々いじめられたから嫌ｗ
0481デフォルトの名無しさん2008/11/16(日) 12:49:36
>>477
c を複素数に振るときには character 型は s で始めるようにしてるな。
>>478
l は先頭でだけ使うようにしてる。 先頭が数字の変数ははじかれるから
混同しなくてすむ。
0482デフォルトの名無しさん2008/11/26(水) 01:24:43
東京で先週開かれたＩＳＯのＦｏｒｔｒａｎ部会でＣｏＡｒｒａｙＦｏｒｔｒａｎに関する扱いで
何か進展はあったのでしょうか？

つまりＣＡＦを排除できたのでしょうか？
0483デフォルトの名無しさん2008/11/29(土) 16:57:56
>>482
東京会議の結論がうｐされたようだ。
http://www.nag.co.uk/sc22wg5/


0484デフォルトの名無しさん2008/12/06(土) 18:55:27
The shapes of the array expressions do not conform.
って出てくるのですが、どういうことですか？
0485デフォルトの名無しさん2008/12/11(木) 16:23:19
配列間違ってないか？
0486デフォルトの名無しさん2008/12/13(土) 02:11:30
C Pad for Salford FTN77
でCPU TIME を出力できるようにする方法はありますか？
0487デフォルトの名無しさん2008/12/13(土) 18:02:53
>>486
HIGH_RES_CLOCK@

使い方は Library Reference に載ってる。
0488デフォルトの名無しさん2008/12/23(火) 09:11:38
>>487
亀ですが
ありがとうございます。
出来ました！
0489デフォルトの名無しさん2009/01/24(土) 19:01:34
関連スレが新しくなりました。

くだすれFORTRAN（超初心者用）その4
http://pc11.2ch.net/test/read.cgi/tech/1232789521/l50
0490デフォルトの名無しさん2009/02/12(木) 00:23:26
うちの研究室の人たちがFortranでプログラミングしているけれど
お前らgoto使いすぎ・・・
というか自分の提出したプログラムがgoto文に書き換えられている・・・

subroutine、functionを使いすぎるのはスパコン計算に向かないと言っていたけれど
勉強してきた構造化プログラミングの手法自体がfortran向きでないということなの？
0491デフォルトの名無しさん2009/02/12(木) 02:32:22
>>490
インラインされるようにしろｗ

コンパイラが安心してグローバル最適化できるように、付加的な情報を与えてやれ。
グローバルな構造が人間が見て分かりやすいならコンパイラも見て分かりやすい。

ＩＮＴＥＮＴ指定はつける。
ＦＵＮＣＴＩＯＮは、ＰＵＲＥ　（可能ならＥＬＥＭＥＮＴＡＬ）にする
とかすればDO..LOOP内で関数呼び出しとかしても並列化やベクトル化してくれるはず。
コンパイラによるが。

実際にベンチマークを取って、自分が勝利すれば、後は文句いわれないべｗ
0492デフォルトの名無しさん2009/02/12(木) 12:52:18
>>491
intent属性つけると微妙に速くなるよね。
特にloop内で頻繁に呼ばれるやつだと全体で5から20%くらい違ってくる。
最近のコンパイラはF77で書いててもintentをわかってくれるから楽だ。

で・・・・pure,elementalの事はよくわからないので詳しく！
0493デフォルトの名無しさん2009/02/13(金) 00:04:13
>>492
PUREもELEMENTALもFortran95で導入されたキーワードで、関数やサブルーチンの
副プログラムの性質をコンパイラに教える修飾語。

PURE FUNCTION unko(x) などのように宣言の頭につけて使う。

PUREというのは、大雑把には副作用がないこと。それゆえ、並列に動作させても
値が変わらないので、コンパイラは並列化など最適化をやれる。
（PUREを満たさねばならない条件はいくつかあるので、ｋｗｓｋはマニュアル見てくれ）

ELEMENTALはPUREの中でさらにいい性質を持ったもの。
本来は１変数で定義してあっても、配列を引数にとって、配列を返すように使ってもいいもの。

たとえば、
REAL :: a(10), b(10)
b = SIN(a)
こういう使い方が組み込み関数では出来るわけだが、自作の関数でもこうやって使っても
いい場合にELEMENTAL　（要素毎の意味）　をつける。　

これもコンパイラが利口なら、最適化をかけてくれるはず。

普通、DO..LOOP内に関数やサブルーチンの呼び出しがあると、呼び出し先で何が
あるかわからないので、最適化が抑止されるが、PUREやELEMENTALがあれば、
安心してコンパイラが最適化できる。（原理的にはｗ）

0494デフォルトの名無しさん2009/02/13(金) 08:29:24
>>493
ありがとう。pure 属性は使い易そうだ。
04954942009/02/13(金) 11:01:13
loop内で頻繁に呼ばれる４個つくらいの function や subroutine に pure 属性つけただけで
並列(OpenMP)の時にちょびっと速くなった。わーい。
もともとdo-loopで一つずつ処理しているので elemental を使える状況は少ないけど、
もう少し理解してから使ってみよう。
inlineを少し深くするようになったからかコンパイル時間がすこし長くなったけどまあ気にしない。
0496デフォルトの名無しさん2009/02/17(火) 22:12:50
>>491
なるほど、ためになったわ。サンキュー。PUREを付けるだけで簡単に速くなったわ。

担当のは基本的に古い人間だから、あまりコンパイラ任せにするなとか言うんよな…。
Fortran90っぽい77形式の文も、「完成した後で直します」とか言ってそのままだったり。
implicit realもcommon文もどうにかしたいけれど、きっと文句言うんだろうな。
0497デフォルトの名無しさん2009/02/18(水) 02:22:54
>>496
PURE 属性はもともと HPF で導入されたものなので、それなりのメーカならコンパイラ
にも最適化のノウハウの蓄積があるはず？
ELEMENTAL　も一変数用にサブルーチン書いて、配列を渡せばいいので、DO..LOOPを
あからさまに書かずにすむ。時々ソースがすごくすっきりして楽になることがある。

でも、ついでだから古い人間を擁護してみるかｗ

FORTRAN77 にも、むきだしの良さがあるｗ
Fortran90 以降だとソース上は簡潔になっても、影で動的にテンポラリ配列とか取るので
気持ち悪いところもある。長時間並んだJOB がStackOverflow とかで止まっていたり
すると殺してやりたくなる。古い時代の人はその辺が嫌なのかもしれないｗ

COMMON もグローバル変数だと思うとむかつくだろうが、 自分でガベージコレクション
しながら、変数領域を自由自在に再利用して操っていると思うと、結構気持ち良い。
メモリーが貴重だった時代には、利用価値が高かった。

IMPLICIT に関していえば IMPLICIT NONE をつけた上で、頭文字を守ればいい。
C と違ってローカルスコープの捨て変数が使えないから、頭文字を守らないと
いちいち宣言部に戻って型を確かめないといけないのでめんどい。
Fortran では数値の精度が致命的に重要で、倍精度計算中に１個単精度が混じるだけでも
すべての結果が単精度のゴミになるのだから、数値精度に無頓着の他のおしゃれ言語の
連中の言説を信じると痛い目に会わされるｗ

0498デフォルトの名無しさん2009/02/18(水) 13:19:26
pure属性付けたサブルーチン内で
pure属性付けた他のサブルーチン呼ぶ時には
interfaceで引数部分の性質を教えて置かないといけないんだね。
f77書式で書いてるのに interface intent が
現れるのでムズムズするけどムズムズするだけなのでそのまま。
0499デフォルトの名無しさん2009/02/26(木) 16:08:37
OpenMPを試してみてるんだけど、いまいちうまくいかないんだが・・・
ifort -QParallel でコンパイル＋実行するとうまくいくんだけど、
同じソースをifort -Qopenmp でコンパイル＋実行すると
Program Exception - stack overflow
っていうエラーが出る
ソース中にOpenMP構文は use omp_lib　だけ
これってどうしたらいいのかな？
05004992009/02/26(木) 16:22:10
おｋ事故解決
-auto_scalar　オプションでいけた
0501デフォルトの名無しさん2009/02/26(木) 16:49:38
>>500
だたの
-auto
や
-static
でうまく行くこともあるよ。あとは実行時の環境変数の
setenv KMP_STACKSIZE 16M
とかデフォより大きめに取ると幸せなことが多い気がする。
0502デフォルトの名無しさん2009/02/26(木) 16:50:50
>>500
ていうか、
-auto_scalar
ってデフォルトになってない？
0503デフォルトの名無しさん2009/02/26(木) 17:23:54
-openmpだとデフォで動的（-auto）になるんだっけか
>>499は単にスタック不足なだけでは？

ところで、KMP_STACKSIZEってunlimitedに設定できないの？
0504デフォルトの名無しさん2009/02/26(木) 17:58:56
動的っていってもヒープとスタックがあるよね。
Fotranの場合ってどういう風に割り当てられるんだろう・・・そういやあまり意識したことなかった。
0505デフォルトの名無しさん2009/02/27(金) 02:43:33
>>504
元々はAUTOMATIC変数の類はスタックに、ALLOCATE命令のときはヒープだった気がする。

ただ歴史的経緯のせいでサブルーチン側で巨大な配列をとっているプログラムが多くて
それがSTACKOVERFLOWを引き起こすから、IntelFortranなんかはデフォールトでは
スカラー変数はスタックで、配列の類はヒープで確保するようになった気がする。

インテルの掲示板に、デフォールト・オプション切り替えの話が出ていた旗もするのだが
昔のことなので記憶がはっきりしない。

OpenMPを指定すると、オプションが切り替わってしまう話はこのスレに出ている。
http://software.intel.com/en-us/forums/intel-visual-fortran-compiler-for-windows/topic/59627/


FortranはFORTRAN66時代がSTATICに変数をとっていた関係で、
FORTRAN77で規格の上ではサブルーチン類の変数はAUTOMATICでありうるように
なったにもかかわらず、互換性の観点からデフォールトでSTATICに変数を
とっていることが多くて、それはFortran90でもしばらく続いていたという歴史的経緯がある。
0506デフォルトの名無しさん2009/02/27(金) 11:54:19
>>505
そのスレのSteveさんが言ってるlocal scalar variablesって、配列じゃない（次元１の）フツーの変数のことだよね？
だとしたら、
If you say /Qopenmp (or /Qopenmp-stubs), this changes to /automatic - ALL local scalar variables are allocated on the stack.
のALL local scalar variablesはALL local variables（つまりローカルスカラ変数だけじゃなく、ローカル配列変数も含む）の間違いじゃないかな？


勘違いしてたらすいません
0507デフォルトの名無しさん2009/02/27(金) 13:12:01
>>506
そうなんじゃないかな？
配列までスタックにとるからoverflowしていると考えるとつじつまが合う？
05085062009/02/27(金) 13:52:41
>>507
そうだよね。SteveさんIntelのスタッフなんだし書き間違わないでくれよぉ・・・

OpenMPの場合にデフォルトで全てのローカル変数（スカラも配列も）がスタックになるのは
スレッドセーフのためってあるけど、逆に言えばそれらをstaticに扱ったらスレッドセーフに
ならないってことだよね。それってやっちゃまずいことなんだろうか？

OpenMPまともに使った経験がないもんで・・・識者の方コメントおながいします
0509デフォルトの名無しさん2009/02/27(金) 14:48:49
-openmp
つけてコンパイルすると、仮にautomaticは嫌いだっと
-noauto
つけても消されるはず。なのでopenMPを読んだ状態だと
有無をいわさずにぜんぶautomaticになる・・・と理解しているけど。
shared属性をつけてかつcommonブロックに貼り付けるとかして
やや無理なことをしないとstatic領域には配置されないと思ふ。
-staticはその名のイメージとは違って共有ライブラリを使わない
(というか実行ファイルに埋め込んじゃう)オプションだし。

spec.orgにあるOMP2001の情報のところに試した人たちが
実際に使ったオプションが載ってるので、自分のところのシステムに
近い(or 同じ)マシンがあったらそれと同じモノをだいたい使っているけど、
まあ大間違いではないと思う。お薦めできるのかは分らないけど・・・。
0510デフォルトの名無しさん2009/02/27(金) 15:30:30
え、OpenMPだとサブルーチン内ではstatic配置ってできないわけ？
数値計算では馬鹿でかい配列は可能な限りstaticにするのが基本じゃない？
サブルーチン渡る度に配列をコピーとかやってられんでしょ？

その代わりメモリのアロケーションはコーディングする側がきちんと考えてあげないといけないけど。
この辺がリソースをある意味富豪的に使うことが前提とされるJavaやC#などのオブジェクト系言語との違いだね。
0511デフォルトの名無しさん2009/02/27(金) 16:56:15
510はちょっと尋ね方が変だな。
副プログラムが並列のloop内にあるのか、
副プログラム内で並列のloopがあるのかで答は変わる。
05125102009/02/27(金) 19:38:52
的外れな質問してたのかな、ごめん・・・。
>>511の両方のケースとも知りたいので、よかったら教えてほしい。
流体とかの数値計算では行列の反復法とかはサブルーチン内でやるけど、
そういうケースでメモリのアロケーションはどうすべきなのか、という視点で興味があります。
0513デフォルトの名無しさん2009/02/28(土) 00:58:18
>>512
引数に載ってれば、Fortranは参照呼出しなんだから配列のコピーが生じるはずもない。
CやPASCALとは違う。（まぁ最近はVALUE属性も出来るようになったが）

ここでは副プログラム内で新たに割り付けられる配列の話だろう。
それが動的に取られるとして、スタックに取られるか、ヒープに取られるかの違い。

05145102009/02/28(土) 01:28:49
>>513
＞副プログラム内で新たに割り付けられる配列
あ、そういう話か・・・ありがとう。

自分は一時作業用の配列も全部メイン側で用意してやって引数渡ししてるんだけど、
そういうのはサブルーチン内でスタックやヒープで取ってやる方がいいのかな？
皆さんどうしてますか？
0515デフォルトの名無しさん2009/03/03(火) 06:16:03
allocatableを副プログラム内で作ってそれ(先頭アドレスやらサイズその他)を
メイン階層や他の副プログラムに行き渡らせるのが割りと
面倒臭いのでメイン階層でやってる。
メイン階層が長くなるのは正直あまり好みではない、けど
まあ好みの問題なだけなので気にしないようにしている。

allocatable使うと場合によってはプログラム内で明示するサイズより大きいメモリを扱えたり
(コンパイラはハンドルできるメモリサイズ ＜ OSがハンドル出来るメモリサイズ、なので)、
格子数(データ数)の違う処理を同じ実行イメージで出来る(コンパイルしなおさなくて良い)のが
いいけど、最近ではサイズをプログラム内で固定変数で明示していちいちコンパイルしている。
その方がミスは少ない気がする。気のせいかもしれないけど。
0516デフォルトの名無しさん2009/03/03(火) 12:11:34
何がスタックやヒープに割り当てられるのか分からんかったのでググった。
ブログだけどｗ

ttp://monologuemidnight.blogspot.com/2009/01/blog-post_03.html

module中の大域的な変数、save属性付きの変数→.bss
allocatable属性な配列→ヒープ
普通にサイズが宣言されている自動配列→スタック
手続き中の変数（mainでも同様）→スタック

ってことらしいが…。
この辺は規格ではないっぽいので、コンパイラの仕様をちゃんと確認したほうがいいな。
0517デフォルトの名無しさん2009/03/11(水) 18:46:48
ifortで複数の数字を出力すると勝手に改行されるんだけど、改行を抑制するオプションってある？
gforntranでは改行されないんだけど。
real(8) a(5)
write(*,*) a(:)
とするとa(3)とa(4)の間で改行されちゃう。
0518デフォルトの名無しさん2009/03/12(木) 04:57:59
>>517
気持ち悪いけど、下記のようにしてもだめかな？
クラシックな書き方だけど。

　　write(*, '(100f10.5)') a(:)
0519デフォルトの名無しさん2009/03/12(木) 09:37:42
>518
ありがとん。
write(*,'(100f)') a(:)
でいけた。
0520デフォルトの名無しさん2009/03/15(日) 01:01:35
改行抑制なら、write文ならadvance='no'を付ければできるはず
出力するものが決まってるなら、改行してしまう書式の箇所に\（バックスラッシュ）を入れてもいい
0521デフォルトの名無しさん2009/03/15(日) 11:34:42
>520
advanceはフォーマット文がいるから、結局519と同じでは？
0522デフォルトの名無しさん2009/03/17(火) 03:12:33
>>521
'(100f)' なんて書くよりはスマートな気がするけど？
ま、どうでもいいことだな
0523デフォルトの名無しさん2009/03/18(水) 04:24:12
do i = 1, 5
  write(*, '(f10.5,$)') a(i)
enddo
write(*,*) ''

はダメかな？
改行するなよ、の＄は非標準だけどDEC以来ほぼ標準と思っていいのではないだろうか
0524デフォルトの名無しさん2009/03/18(水) 21:51:40
ちょっとテストしてみた＠Intel Fortran 10.1

integer::i
integer,parameter::imax=5
real(8)::f(imax)
do i=1,imax ; f(i)=i ; enddo

print *, "print *,f"
print *, f
print *, '---'

print *, "write(*,'(f)') f"
write(*, '(f)') f
print *, '---'

print *, "write(*,'(f)') (f(i),i=1,imax)"
write(*, '(f)') (f(i),i=1,imax)
print *, '---'

print *, "write(*,'(f\)') f"
write(*, '(f\)') f
print *, '---'

print *, "write(*,'(f)',advance='no') f"
write(*, '(f)',advance='no') f
print *, '---'

print *, "write(*,'(f)',advance='no') (f(i),i=1,imax)"
write(*, '(f)',advance='no') (f(i),i=1,imax)
print *, '---'

end
0525>>524の結果2009/03/18(水) 21:52:34
print *,f
1.00000000000000 2.00000000000000 3.00000000000000
4.00000000000000 5.00000000000000
---
write(*,'(f)') f
1.0000000000000000
2.0000000000000000
3.0000000000000000
4.0000000000000000
5.0000000000000000
---
write(*,'(f)') (f(i),i=1,imax)
1.0000000000000000
2.0000000000000000
3.0000000000000000
4.0000000000000000
5.0000000000000000
---
write(*,'(f\)') f
1.0000000000000000 2.0000000000000000 3.0000000000000000 4.0000000000000000 5.0000000000000000 ---
write(*,'(f)',advance='no') f
1.0000000000000000
2.0000000000000000
3.0000000000000000
4.0000000000000000
5.0000000000000000 ---
write(*,'(f)',advance='no') (f(i),i=1,imax)
1.0000000000000000
2.0000000000000000
3.0000000000000000
4.0000000000000000
5.0000000000000000 ---
0526デフォルトの名無しさん2009/03/18(水) 22:07:03
ごめん。
print文かwrite文か、また出力並びをfとするか(f(i),i=1,imax)とするかで、
出力結果が違うのかと思ってやってみたが、違いはないっぽい。
というわけで、>>524をもっとすっきりさせた。

integer::i
real(8)::f(3)
do i=1,3 ; f(i)=i ; enddo

print *, "write(*,*) f"
write(*,*) f
print *, '---'

print *, "write(*,'(f)') f"
write(*,'(f)') f
print *, '---'

print *, "write(*,'(f\)') f"
write(*,'(f\)') f
print *, '---'

print *, "write(*,'(f)',advance='no') f"
write(*,'(f)',advance='no') f
print *, '---'

end
0527デフォルトの名無しさん2009/03/18(水) 22:19:08
↑の結果

write(*,*) f
1.00000000000000 2.00000000000000 3.00000000000000
4.00000000000000 5.00000000000000
---
write(*,'(f)') f
1.0000000000000000
2.0000000000000000
3.0000000000000000
4.0000000000000000
5.0000000000000000
---
write(*,'(f\)') f
1.0000000000000000 2.0000000000000000 3.0000000000000000 4.0000000000000000 5.0000000000000000 ---
write(*,'(f)',advance='no') f
1.0000000000000000
2.0000000000000000
3.0000000000000000
4.0000000000000000
5.0000000000000000 ---


結論
①書式省略すると、適当な位置で改行してくれる（規格で決まってるんかな？）
②書式指定すると、書式の中身が終わる度に改行
③書式に\を付けると改行抑制
④書式にadvance='no'（デフォルトはyes）を付けるとwrite文が終了したときに行う改行を抑制（②との違いに注意）

って感じかな。
改行抑制の\と改行指定の/が混ざるとどんな出力になるんだろう。
0528デフォルトの名無しさん2009/03/18(水) 22:20:24
ごめん、>>526は
real(8)::f(3)
do i=1,3 ; f(i)=i ; enddo

は↓に修正。

real(8)::f(5)
do i=1,5 ; f(i)=i ; enddo
0529デフォルトの名無しさん2009/03/22(日) 01:30:43
＞①書式省略すると、適当な位置で改行してくれる（規格で決まってるんかな？）
コマンドプロンプトの右端で改行だと思う
違うかなー
0530デフォルトの名無しさん2009/03/22(日) 10:23:20
>>529
他のコンパイラでは改行しないので、コンパイラのせいだと思います。
>>524-528
ありがとうございます。
参考になりました。
write(*,'(f\)') f
が良さそうですね。
0531デフォルトの名無しさん2009/03/23(月) 15:14:13
>>529
Fortranに限らず、プロンプトのサイズは関係ないよ
0532デフォルトの名無しさん2009/03/27(金) 05:59:41
亀だけど、Ifortなら<>がお勧め。
多次元配列の最初の数を入れることが多いです
例
program main
implicit none
integer,parameter :: num = 9
integer :: ii,jj
real :: arry(num,num)
do ii=1,num
do jj = 1,num
arry(ii,jj) = ii*jj
enddo
enddo

write(6,'(<num>F)') arry
end program
0533デフォルトの名無しさん2009/03/27(金) 09:10:33
>532
勉強になった。
ありがとん
0534デフォルトの名無しさん2009/03/28(土) 00:52:37
>>531
WindowsはDOS窓のサイズで改行位置が変わるぞ。

0535デフォルトの名無しさん2009/04/03(金) 21:43:51
>>532
<>で整数くくるとその数が使えるんだね・・・知らなかった。昔からあった機能？

>>534
そうなんだ。>>531はLinuxの場合ルだけかな
0536デフォルトの名無しさん2009/04/03(金) 22:09:05
>>535追加

<実数変数>でも勝手に整数にキャストしてくれるね。
あと他のところにも、例えば
num1=5
num2=1
write(6,'(<num>F<num1>.<num2>)') arry
みたいにして使えるね。
0537デフォルトの名無しさん2009/04/04(土) 02:01:36
>>536
Fortran2008　では、　*( i3, '-' ) 見たいな感じで　＊　が無限反復の指定子になるらしい。
0538デフォルトの名無しさん2009/04/08(水) 09:50:52
*( i3, '-' )
を、スキー帽がぶったバカボンのパパっぽい大人と
すました子供に見えるようになると、やばいね。
0539デフォルトの名無しさん2009/04/08(水) 12:40:58
やばいな
0540デフォルトの名無しさん2009/04/17(金) 23:47:22
Fotranはじめようと思ってインストールして、
簡単なプログラムを実行しようとしたんだけど、

ld: crt2.o: No such file: No such file or directory

ってエラーが・・・。ググってもよくわからず・・・。
助けてください。
05415402009/04/18(土) 00:12:28
アンインストールして、はじめからやり直したらでけた。
0542デフォルトの名無しさん2009/04/18(土) 20:15:44
>>542
おまえ、プログラムのセンスあるから、今後も缶ばれ！
0543デフォルトの名無しさん2009/04/19(日) 12:16:56
>>542
自画自賛?
0544デフォルトの名無しさん2009/04/23(木) 16:31:22
windowsXPでフォートラン90をフリーで使いたいのですけど、
cygwin以外でどうにか出来ますか？
0545デフォルトの名無しさん2009/04/23(木) 22:10:01
free Fortran でぐぐれ
0546デフォルトの名無しさん2009/04/24(金) 00:32:48
>>544
g95だな・・・
0547デフォルトの名無しさん2009/04/24(金) 23:53:11
windowsでVMWareを立ち上げてUbuntuいれてgfortran
0548デフォルトの名無しさん2009/05/13(水) 00:24:34
n次元配列をなめるルーチンを作成できる?
reshepe使って一次元に落とし込んで処理しようとしたら、segmentation faultになってしまった。
コードはこんな感じ。
-----------------------------
program main
integer :: arr(3,4,5)
integer :: i,j,k
do i=1,3
do j=1,4
do k=1,5
arr(i,j,k) = i+j+k
enddo
enddo
enddo
call check(arr)
end program main

subroutine check(arr)
integer,intent(in) :: arr(:,:,:)
integer,allocatable :: tmp(:)
integer :: i,elem=1
do i=1,size(ubound(arr))
elem = elem * (ubound(arr,i)-lbound(arr,i)+1)
enddo
allocate(tmp(elem))
tmp = reshape(arr,shape(tmp))
do i=lbound(tmp,1),ubound(tmp,1)
print *,tmp(i)
enddo
end subroutine check
-----------------------------
0549デフォルトの名無しさん2009/05/13(水) 01:22:27
>>548
PRINT *, arr
でだめなのか？
0550デフォルトの名無しさん2009/05/13(水) 20:37:56
>>548
ほらよ。つーかこれテンプレにしてちょ


・形状引継ぎ配列を引数に使う
・割り付け配列を引数に使う
・ポインタを引数に使う
・引数キーワードやoptional属性を使う
・関数の返値を配列にする

↑をしたいなら、副プログラムは明示的引用仕様にしなきゃならない。

■明示的引用仕様
・副プログラムを内部副プログラムにする（呼び出し側でcontains）
・副プログラムをモジュール副プログラムにする（モジュールでcontains）
・副プログラムが外部副プログラムならば、それを呼び出す側でinterfaceを用意する
■暗黙的引用仕様（FORTRAN77的なやり方）
・副プログラムが外部副プログラムで、それを呼び出す側でinterfaceを用意しない

なお、副プログラムとは関数やサブルーチンを指す。
0551デフォルトの名無しさん2009/05/13(水) 23:25:40
>>548
多次元配列を一次元にするには PACK UNPACK を使うのが一般的だと思う。
Fortran90 で導入されている。

例題はマニュアルをみてちょｗ
0552デフォルトの名無しさん2009/05/14(木) 23:40:08
FORTRAN90で引数の整合性をチェックする方法ってないかな？
よくそこでバグがおきるんだ。
05535482009/05/14(木) 23:50:44
>>549
print の部分は、本番仕様ではガチャガチャ変更したいところなのでございます。
0554デフォルトの名無しさん2009/05/14(木) 23:51:21
>>552
>>550
05555482009/05/14(木) 23:51:43
>>550>>551
って～ことは、呼び出し側でinterfaceを用意して、packで一次元化するということか。
------------------------------------------
program main
integer :: arr(3,4,5)
integer :: i,j,k
interface
subroutine check(arr)
integer,intent(in) :: arr(:,:,:)
end subroutine check
end interface
省略
call check(arr)
end program main

subroutine check(arr)
integer,intent(in) :: arr(:,:,:)
integer,allocatable :: tmp(:)
integer :: i
tmp=pack(arr,mask=(arr .ne. 0))
do i=lbound(tmp,1),ubound(tmp,1)
print *,tmp(i)
enddo
end subroutine check
------------------------------------------
これでできました。ありがとうございます。
余談だけど、gfortranでは期待通りに処理されるのに、intel fortranだとダメみたい。
pack()の返り値が壊れてしまうっぽい。バグか?
0556デフォルトの名無しさん2009/05/15(金) 00:45:08
>>555
tmpをALLOCATEしてないからかな？

ALLOCATE(tmp(SIZE([arr .ne. 0])))

これを入れておくと動く。
配列を返すものを、未割り当ての割りつけ配列に入れたときの挙動が、
F90の範囲ではプログラム側に問題があって悪いとおもうが、
新しい規格ではかなり自動で調整してくれるようになったので、
これがBugなのかどうか文法書を読み直さないと分からんｗ

[ ] はF2003で、f90なら(/ /) に対応。
05575562009/05/15(金) 00:49:27
ああ、Fortran95/2003　Explained 17.5.2 に そのものずばり、PACKの結果を
未割り当てALLOCATABLEに突っ込む例題が載っていた。

Fortran95までは出来なくて、F2003で出来るようになるものだそうだ。

あと、>>556で [ ] は無くても良かったｗ　
0558デフォルトの名無しさん2009/05/15(金) 01:22:11
>>556.
＞新しい規格ではかなり自動で調整してくれる
へーそうなんだ
いちいちallocateしなくてすむのは嬉しい
05595482009/05/17(日) 22:24:17
>>556>>557
フォロー、ありがとうございます。
allocateでintel fortranでも無事動くようになりました。
てか、2003いいね。期待大。
0560デフォルトの名無しさん2009/05/21(木) 00:38:15
スレ違いだったら失礼します

FORTRANで他からデータを読み込むのに、そのデータの最初にデータ概要が書いてあるので
読み飛ばす方法があったら教えてください

初心者なので分かりにくくてすみません
0561デフォルトの名無しさん2009/05/21(木) 01:30:02
>>560
適当な変数を定義して、読み込みだけさせて、その中身を
使わなければいいぢゃない。嘘データで良いから実例だし
てくれ。

# リアルなのを晒すとマズいかもしれんし。
0562デフォルトの名無しさん2009/05/21(木) 01:54:07
>>561数字を読み込んで並べるだけなのですが


読み込むデータが

○月×日
気象データ
東京都日黒区

0.0001
0.0025
0.0025
0.0105
0.1568
1.5385
1.4657
12.368

こんな感じで数字の部分のみ読み込みをしたいです
0563デフォルトの名無しさん2009/05/21(木) 02:05:29
データ構造が既知ならば読みたいところまで回していけばいいし、
レコード長が分かってれば指定のレコード部を直接読みに行ける、はず

もっと良い方法あったらﾌｫﾛｰよろ
05645632009/05/21(木) 02:08:22
あ、sequentialで書いたファイルをdirectで読む（あるいはその逆）って、できないんだっけ？
自信なくなってきた＞＜；
0565デフォルトの名無しさん2009/05/21(木) 02:15:57
>>563すみません、もう少し詳しくお願いできますか？
0566デフォルトの名無しさん2009/05/21(木) 09:45:54
>>562
いくつか方法はある。
最初のデータ概要の行数が決まっていて一定ならば簡単で、
CHARACTER TEXT*256
DO i = 1, 4
READ(7, *) TEXT
END DO
的に読み飛ばせばいい。

手抜きな方法としては、数値変数に文字列を読み込むとエラーが出ることを利用して、
10 CONTINUE
READ(7, *, ERR= 10) dat
で無理に読み飛ばすことも出来る。
本物のエラーのときや、データ概要が数値として読めてしまうと困るｗ

やや高度な方法としては、いったん文字列に読み込んで、その文字列の中身をしらべて
文字列を内部ファイルとして入力に使って数値を読むという方法もある。

データ開始地点まで適当にすすんで、ＢＡＣＫＳＰＡＣＥで適切なところまで戻る手もある。
0567デフォルトの名無しさん2009/05/21(木) 12:20:07
>>566
２つ目の手抜き法はオイラもよく使う。
念のため、

icount = 1
20 continue
read(7,fmt=*) string
read(string,fmt=*,err=10,end=999) dat(icount)
icount = icount + 1
goto 20
10 write(*,*) string
goto 20
continue 999

みたいに文字列を一度かまして
読み込みは内部ファイルで処理して、
飛ばされたケースは表示させて目でチェックするとか。
あまりほめられないスパゲッチ～だけど。
0568デフォルトの名無しさん2009/05/23(土) 20:59:00
>>562
本気でスピード求めるならdirectアクセスがベストだが
多分あなたの使い方じゃそうじゃなさそうだから、普通に次々読んでいって
使わないデータは破棄しちゃえばいいじゃん

>>562がASCIIデータだとすると、5行目からの数字が欲しいなら、
do i=1,5
read(unit) tmp
enddo
とやって捨てちゃえばいい
0569デフォルトの名無しさん2009/05/23(土) 21:16:21
読み捨てるなら

DO I = 1, 5
READ(NUNIT, *)
END DO

でおｋ
0570デフォルトの名無しさん2009/05/23(土) 23:01:20
そっか、捨てるんだから何もいらないか
0571デフォルトの名無しさん2009/06/09(火) 12:40:21
FORTRANがCOLUMN-MAJORなことを調べてみたいのですが、
変数のポインタ値を知る方法はありますか？
0572デフォルトの名無しさん2009/06/09(火) 12:49:09
>>571
無いが、サブルーチンに送ればわかる。
つまり、２次元配列を、サブルーチンの１次元形式引数の配列に渡せばよい。
0573デフォルトの名無しさん2009/06/09(火) 13:18:41
ありがとうございます。確かにWRITEなんかにわたすことで
順番を知ることはできるのですが、メモリ上の位置を知ることは
できないのですね
0574デフォルトの名無しさん2009/06/09(火) 13:32:54
>>573
標準文法ではないが、LOC 関数などがベンダー拡張として用意されていることが多い。

またリンカーでマップを出力させると静的な変数ならばアドレスが分からなくも無い。
配列のならびに関しては、そのままでは分からないが、さらに EQUIVALENCE でも
かければ分かるのかもしれない。
0575デフォルトの名無しさん2009/06/20(土) 14:16:58
IFort10.1+マルチコア環境で使ってる人いるかな
OPENMPで並列化してるんだけど、
デュアルコア環境でトータルCPU使用率が50%から上がらない理由わかる？
１コアしか使ってなくて2コア目がアイドルになっちゃう
0576デフォルトの名無しさん2009/06/20(土) 14:21:55
OSによるだろうけど、
setenv なんたら、
をちゃんとしてる？
05775752009/06/21(日) 16:49:20
setenvはよくわからないんだけど、OPENMP+自動並列化で5スレッド走ってる
環境はOSがXPProSP3で、Core2DuoE8500
■ このスレッドは過去ログ倉庫に格納されています