FORTRAN専用スレッド

>>128
どうもです。ありがとうございました。

>>128
廃止"予定"事項なのでFortran95ではまだ使えます。

Fortran90で廃止予定事項として上がっていた物のうち
実数型制御変数を使ったDOループやPAUSE文などは
Fortran95で廃止されています。

ただ、規格上は廃止されていても実際のコンパイラでは
これらの機能を使える物がほとんどのようですが。
（オプション等でエラーや警告が出るようにはなります）

至急、質問があります。
ＦＦＴ（高速フーリエ変換）のプログラムを作ってます。
今まで処理データ数が１６０００個くらいで、順調に計算できていました。
しかし、データ数を７００００個くらいに増やすと、出力ファイルに
『*』印がたくさん出てきます。
ちなみに言語はフォートランで出力はＣＳＶ形式で出力してます。
これはフォーマットの設定が悪いんでしょうか？

>>131
マルチポストするような奴には教えてやらん。

申し訳ありませんが質問なのですが。

character*10 cc
real array(10)

array(5)=5.0

write(cc,'(A)') 'array'

として、character の　cc　を使ってarrayの中身を表示することは可能でしょうか？
イメージ的には

write(*,*) cc(5)

のような感じです。（上のままではコンパイル通りませんが）

>>133
最終的には，文字列が加工できるから使う変数を動的に変えられるだろうな，って話?

>134
その通りです。
でも、やり方がわかりません。何か別にうまい方法があるのかもしれませんが。

>>135
こういう方法を使わないアルゴリズムを作れませんか?
高級言語なら，文字列からシンボルを引くことができます．
C 言語でも似たようなことができ，関数へのポインタと文字列をメンバにした
構造体を用意します．

fortran はよく知りませんが，関数へのポインタを変数に入れられましたよね?
文字列の配列を作って，さらに関数へのポインタを入れる配列をインデックスを揃えて
作るとかは? ．．．できるのかな ;-D

>>136
すでに存在する巨大なfortranプログラム＋巨大なinclude fileをベースに
一部変更する必要があり、その中で上のような事ができれば便利だと
感じたので質問してみました。

何か別の方法を考えてみます。

>>136-137
EXTERNALを使えば関数へのポインタを引数に取れるけど
変数に代入したり、配列にするのは無理。

intrinsic cos
external bar2, bar3
write(*,*) foo(cos, 0.)
a = 2.
write(*,*) foo(bar2, a)
write(*,*) foo(bar3, a)
end

real function foo(f,x)
external f
foo = f(x)
end

real function bar2(x)
bar2 = x**2
end

real function bar3(x)
bar3 = x**3
end

あとはこのページが参考になるかな？
ttp://user.ecc.u-tokyo.ac.jp/~tt37103/tips/f90/function_pointer.html

>>138
できないとなると悲しいな．

DIMENSION A(100)というのの100のところは最大で、どのくらいまでできるんでしょうか？

DIMENSION A(131072)までは、プログラムが実行されるのに、
DIMENSION A(262144)にすると、
「Stack Overflow Re-ling program with bigger stack value.
stack:reserve,commit.Will attempt to trace back」
「Stack Overflow: at address=541718,esp=547f6
00541718 routine at address 541718[+0000]」
って出てプログラムが実行できないんです。

>>140
処理系に依る。
FTN95だと
DIMENSION A(300000)
とかでもできた。

http://www.hulinks.co.jp/software/pf_mac/

>>140
コンパイル・オプションで、スタックのサイズを変えられないの？
あとは、動的にメモリを取得するとか。

>>140
コンパイル・オプションで、スタックのサイズを変えられないの？
あとは、動的にメモリを取得するとか。

>>139
Fortran 2003ならprocedure pointer （関数・サブルーチンへのポインタ）が使えるんだけどね。
いかんせん、Fortran 2003準拠の処理系がまだないのが…

>>145
お、いつのまにかFortran2003なんて出てたんですね。っていうか2000は？

>>146
ttp://www.itscj.ipsj.or.jp/committees/sc22/sc22fortran.html
>従来，通称Fortran 2000と呼んできたが，2003年度に開催されたWG 5会議において，
>通称Fortran 2003と呼ぶことが決議された．

規格制定が遅れたんで名前が変わっちゃいました。

社会保険庁のシステムがFORTRANで組まれているという話をよく聞きますが、なんででしょう？
この手の業種はかつてはCOBOLが一般的だったような気がするのですが。。。

おいおい、ところどころ宿題スレになってんな。
もっとマ板らしいネタスレにしようぜ。

俺はFORTRAN使ったことない

おれはある

俺は今でも使ってる

俺はFORTRANで書かれている

>>152
名前は６文字以下ね

FORTRAN77(の方言)で、メモリを動的確保するサブルーチン/関数
はありますか。

Fortranでいいフリーのデバッガってないですかね？

gdbも一応Fortranに対応しているようですが，世の中の情報が少なくて困ります．

あと，静的なツールもフリーでご存知の方いらっしゃったらお願いします．当方はctagsを使っているのですが...

open watcom

わからない問題がいくつかあるんですが暇つぶしに解いてくれる猛者はいませんか？？

いないか(´Д⊂ｸﾞｽﾝ

sdf

助けてくれるんですか？？

分からない問題があるのですが助けてくれる方いませんか？？

>>161
タダじゃいや。

以下のようなプログラでは、配列Aの部分配列（スライス）をポインタPと結合しています。

このとき、ポインタの下限値 lbound(P)が1、上限値 ubound(P) が 4 になってしまうんだけど、なんとかポインタの下限値を結合先と一致させる方法はありませんか？ポインタの下限値と上限値をそれぞれ 0 と 3にしたいのです。

program main
integer,dimension(0:3,0:3,0:3,0:10),target :: A
integer,dimension(:,:,:),pointer :: P
p => A(:,:,:,2)
print *, lbound(P), ubound(P)
print *, lbound(A), ubound(A)
end program main

ちなみに、以下のように、配列A全体を結合させると、結合先と添字の下限値が一致しますが、目的はこのようなコードではありません。

program main
integer,dimension(0:3,0:3,0:3),target :: A
integer,dimension(:,:,:),pointer :: P
p => A
print *, lbound(P), ubound(P)
print *, lbound(A), ubound(A)
end program main

フリーのfortranコンパイラを探しています。
Windows用の。

フリーのfortranコンパイラを探しています。
Windows用の。

フリーのfortranコンパイラを探しています。
Windows用の。
















全裸で

ここはプログラ「マ」板ですから，>>164 は釣りだと思います．
いや，「探しているプログラマ」というのを演出しているのかも．

フリーのfortranプログラマを探しています。
Windows用の。

フリーの100万円探しています。
Windows用の

まだVisualFortranってあるの？

Power Stationって何

>>171
1. かつてMicrosoftで出してたWindows 上のFORTRAN

Microsoft Fortran
↓
Microsoft Fortran Powerstation
(↓ 実は別物)
DEC Visual Fortran
↓
Compaq Visual Fortran
↓
(HPに会社ごと身売り HPブランドではVFは放置）
(↓ 実は別物?)
Intel Visual Fortran




2. 昔のIBMのEWS(RS/6000)も、POWERSTATION という名前だったような気がする。

>>172
おお、サンクス

流れが良くわかったよ

>>172
CVFとIVFの関係は、フロントエンド部がCVFを元にしており、
バックエンドがIntel独自という関係ではないかな？
つまり半分別物。

>>163
一般的には、部分配列を取ると、コンパイラ側でいったん非明示的ワーク配列に展開してそれを改めて引き渡したりしている。
今の状況では、本来メモリー的にはすでに連続しているので必要ないのだが、多分ワークに展開されていると思う。

非明示的ワーク配列の添え字が１から始まるという事情で、ポインター側の添え字も１から始まることになってしまうのだろう。
１個目と２個目のプログラムの挙動の違いはそういうことだと推測される。

したがって、手動で明示的に３次元の０から始まるワーク配列に展開して、それにポインタを割り付けるしかないと言うことになる。
多分オーバーヘッド的にはあまり変わらないと思う。

ただし、以上はいくつか他の場合の挙動に基づく推測なので、正しくないかも。

>>175 の補足
program main
integer,target :: A(0:3,0:3,0:3)
integer,pointer :: P(:,:,:), q(:,:,:)
p => A
q => A(:,:,:)
print *
print *, lbound(P), ubound(P)
print *
print *, lbound(q), ubound(q)
print *
print *, lbound(A), ubound(A)
end program main


出力結果　ｂｙ　DVF
0 0 0 3 3 3

1 1 1 4 4 4

0 0 0 3 3 3

♪おおブレネリ　あなたのおうちはどこ
わたしのおうちは　スイッツァランドよ
きれいな湖水の　ほとりなのよ
ヤッ〜ホ〜

FORTRAN IV YOU !