FORTRAN III

FORTRAN総合スレッドです。
FORTRAN77/90/95/2003, HPF等、FORTRAN全般に関する話題を書き込んで下さい。

・前スレ
FORTRAN II
http://pc5.2ch.net/test/read.cgi/tech/1068351911/
Fortran
http://pc2.2ch.net/test/read.cgi/tech/1003214738/

・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/

・フリー／無料評価版のコンパイラ
The G95 project （Fortran95）
http://www.g95.org/
GCC (g77) （今のところ FORTRAN77）
http://gcc.gnu.org/
Salford FTN77/95 Personal Edition （Fortran77/95, Windows, 要登録）
http://www.salfordsoftware.co.uk/software/downloads/compilers.html
FTN77/95のインストール方法
http://www.coastal-env.k.u-tokyo.ac.jp/koibuchi/fortran.htm
the F Programming Language （Fortran90のサブセット, Win/Linux/Solaris）
http://www.fortran.com/F/compilers.html
Open Watcom Fortran Compiler （FORTRAN77）
http://www.openwatcom.org/
Intel Fortran Compiler （Fortran95, 無料評価版(Win)／非商用バージョン(Linux), 要登録）
http://www.xlsoft.com/jp/products/download/download.html
Lahey/Fujitsu Fortran （期限付きトライアル版, Win/Linux, 要登録）
http://www.lahey.com/downloads.htm

・ドキュメント
Compaq Visual Fortran オンライン・ドキュメント
http://www.xlsoft.com/jp/products/intel/cvf/docs/index.html
Intel Fortran Compiler 日本語版マニュアルのダウンロード （要登録）
http://www.xlsoft.com/jp/products/intel/download.html
富士通Fortranオンラインマニュアル
http://www.cc.nao.ac.jp/fsunman/japanese/Fortran/index.html

FORTRAN77 standard
http://www.swcp.com/~walt/F77_std/f77_std.html
http://www.fortran.com/fortran/F77_std/rjcnf0001.html

JISデータベース
http://www.jisc.go.jp/

・チュートリアル等
Fortran90プログラミング
http://www.ip.media.kyoto-u.ac.jp/htomita/
Fortran 90 言語マニュアル
http://www.rs.kagu.tus.ac.jp/yama/
やさしいフォートラン講座
http://www.g.dendai.ac.jp/kouza/top.htm
Fortran ResQ - Fortran初心者に捧ぐ
http://homepage.mac.com/galois21/fortran/index.html
FORTRANプログラミング初級編
http://svr01.damp.tottori-u.ac.jp/~akimoto/Fortran/
FORTRAN 77 入門
http://www-aos.eps.s.u-tokyo.ac.jp/~takagi/f77-enshu/index.html

・その他
Fortran Source Check （オンライン文法チェック）
http://www.lahey.com/check.htm
ftnchek
http://www.dsm.fordham.edu/~ftnchek/
c2f
http://www.unics.uni-hannover.de/rrzn/gehrke/c2f.html
F2C
http://www.vector.co.jp/soft/win95/prog/se056813.html

Fortranの歴史
http://members.at.infoseek.co.jp/kitaurawa/fortran.html

The Real Programmer Stories （日本語訳）
http://www.genpaku.org/realprogrammerj.html

前スレが落ちたようなので、立ててみました。

というわけで、後は任せた。
オレはもう寝る。

とりあえず乙>>1

>>1

print *,"乙"

配列変数の宣言なんですが、引数の範囲を実数で指定すればうまくいく。
引数の範囲を同じ値の変数で指定すれば範囲が大きいときにstackoverflowがでる。
という現象が起こっています。
なぜでしょうか？

>>10
>配列変数の宣言なんですが、引数の範囲を実数で指定すればうまくいく。
>引数の範囲を同じ値の変数で指定すれば範囲が大きいときにstackoverflowがでる。

いまいちよく解らないんだが。

SUBROUTINE FOO(A, N)
REAL A(N)
・・・

とかやってるってこと？

どうも失礼します・・・ fortran77使ってます。
temp = v(im,jm,km)*(x(number)-r(1))*(z(number)-r(2))*(y(number)-r(3))
& +v(im,jm,kp)*(x(number)-r(1))*(z(number)-r(2))*(r(3)-y(number-1))
& +v(ip,jm,km)*(r(1)-x(number-1))*(z(number)-r(2))*(y(number)-r(3))
& +v(ip,jm,kp)*(r(1)-x(number-1))*(z(number)-r(2))*(r(3)-y(number-1))
& +v(ip,jp,km)*(r(1)-x(number-1))*(r(2)-z(number-1))*(y(number)-r(3))
& +v(ip,jp,kp)*(r(1)-x(number-1))*(r(2)-z(number-1))*(r(3)-y(number-1))
& +v(im,jp,km)*(x(number)-r(1))*(r(2)-z(number-1))*(y(number)-r(3))
& +v(im,jp,kp)*(x(number)-r(1))*(r(2)-z(number-1))*(r(3)-y(number-1))
この括弧の数はあってると思うんですが、最後の行で、
Syntax error, found END-OF-STATEMENT when expecting one of: , )
& +v(im,jp,kp)*(x(number)-r(1))*(r(2)-z(number-1))*(r(3)-y(number-1))
-----------------------------------------------------------------------------------^
というエラーが。。。括弧以外に原因があるのでしょうか？？

>>12
72桁超えてる。

>>11
全くそのとおりです。説明が不足でした。
サブルーチン呼び出し時に値を渡し、
サブルーチン内でその変数を使って配列を宣言しているのです。

>>10
>引数の範囲を実数で指定すればうまくいく。

整数だよね

>引数の範囲を同じ値の変数で指定すれば範囲が大きいときにstackoverflowがでる。

日本語がよく分からんが、配列の次元を馬鹿でかくするという意味なら、関数またぐ時だし、仕様かな。
ただ、ソースも具体的な数字もないので、何とも言えんが。
http://www.xlsoft.com/jp/products/intel/cvf/docs/vf-html/pg/pg04_02_43.htm

想像するに 10 は

「配列変数の宣言なんですが、添え字の範囲を定数で指定すればうまくいく。
添え字の範囲を、引数で渡される変数で指定すれば範囲が大きいときにstackoverflowがでる。
という現象が起こっています。
なぜでしょうか？ 」
といいたいのでは？

であれば、
・前半は、配列は、ふつーの配列であり、ヒープからメモリが割り当てられる
・後半は、配列は、自動割付配列となり、スタックからメモリが割り当てられる


対処法は
・>>15にあるように、スタックサイズをあげる
・F90で書き、配列をALLOCATEによって確保する
あたりか？

ここの過去ログはどこかに置いていませんか？

>>17
Fortran
http://pc5.2ch.net/tech/kako/1003/10032/1003214738.html

「FORTRAN II」（２スレ目）の方はもうしばらくしてhtml化されるのを待つか、●を買え。

>>18
どうもありがとうございます。

３次関数の最小自乗法による係数決定のやり方はどこかにありませんか？

http://pc5.2ch.net/test/read.cgi/tech/1103561159/

>>20
>>2の一番上にあります。

unformattedなファイル読み込みルーチン作製中です。
unformattedで書き込んだファイルのreadは全く問題ないのですが、
openの際に、うっかりそうでないファイルを指定した時の振舞で悩んでます。

やりたいことは「ファイル先頭をチラっと見る。意図してないファイル形式なら
その旨のエラー出して終了」です。ですが、ウチの環境で不適切なファイルを
unformattedでopenすると、readの段で強制的にプログラムが停止してしまいます。
これではエラー出すヒマが無く、一方で(私の知識の範囲内では)ファイルの判断には
readをするしかないため、どうにもジレンマに陥ってしまいました。

何かイイ方法御存知の方いらっしゃいますでしょうか？
当方 RHEL3(amd64) + Intel Fortran 8.1 であります。

>>23
iostatでチェック

データファイルから読み込むのに、一行に一つのデータでファイルを作ると
読み込むのに、一行に ","　区切りで複数のデータを与えたら一つ目のデータ
しか正しく読めず、あとのデータが????に成ります。
どうしたらいいのでしょうか。
Salford FTN77です。

Fortranで
　・Cの仕様書工房
　・Cppunit
みたいなツールってないでしょうか？

変数 / 関数の型を表示したり、変数の共有のされ具合が分かる
emacs lisp を探しています。
fortran はハンパじゃないので etags だけではちょっと無理です。。。

>>25
データの数かFORMATの指定が間違ってない？
うちの環境では","で区切っても正しく読み込めたけど。

基礎的なことなんですがｘ＾ｎを求める関数は何でしょうか？

x**n

Fortranの規格名称のことでちょっと教えてください。
『Fortran 95』というふうに95の前にスペースを入れるのと，
『Fortran95』というふうにスペースを入れないのと，
どちらが正しいのでしょうか？

どちらも正しくない。それは通称だ。

特に決まりは無いと思うよ。ぐぐる君の意見だと『Fortran 95』の勝ちだね。

"Fortran 95" の検索結果 約 107,000 件
"Fortran95" の検索結果 約 12,800 件

>>32-33
THANX!

>>31
ちなみに Fortran 95 の規格は
　ANSI X3J3/96-007
　JIS X 3001:1998
　ISO/IEC 1539-1:1997
内容は全部同じ。

>>35
わ，そうだったんですか！
詳しいフォロー，ありがとうございました。

WinXP のパソコンで Fortran 90 のコーディングをしたいのですが，
Salford FTN95 を使うか，それとも，うちの職場にあった唯一の
PC用Fortranコンパイラである Fortran PowerStation 4.0 を使うか，
迷っています。
（他の商品版を買うのは高くて，今のところは躊躇しています）

Win95の頃に発売された Fortran PowerStation 4.0 は Win XP に
インストールしてもコンパイルなどは正常に行われるのでしょうか？

>>37
WindowsXP上で、FORTRAN77のソースをコンパイル、正常実行できた記憶があるよ

> 37
cygwin上でgcc-4.0(CVS版ですが)をmakeしてみる. gfortranができるが
これがf90. 簡単なやつでテストしてみましたが, よっぽどf90の機能を
使ってない限り, 大丈夫そうでした.

g95 の linux バイナリがあるということは，
cygwin か coLinux を入れればフォートランランランってことなんじゃないかな．

g95なら >>3 にWindows用のバイナリも置いてあるわけだが。

>>38-41
レスありがとうございました m(_ _)m"

Windows XP をサポートしている商用版の Fortran 90 或いは Fortran 95 コンパイラで
お薦めのものがあったら教えてください。

↑
できれば日本語マニュアルがあって，便利なライブラリも付属しているものが良いですが‥‥。
ご存じのかた，よろしくお願いします m(_ _)m"

>>43
Intel の Visual Fortran

何をするのに「便利な」ライブラリが欲しいのか知らないが
Windows上でグラフィックが使える QuickWin や
Pro版にはIMSLも付いている。

また付属ではないが、Windows上の事実上の標準なので
MKL や Lapack(95) 、 f90gl 等、有名どころのライブラリは
大抵使える。

>>45さん，ありがとうございました。
ライブラリの件ですが，目指しているのは地球のダイナモの成因となる
外核のMHDシミュレーションなので，その数値計算に便利なライブラリが
あればと考えています（地球流体電脳ライブラリが既に存在してはいますが）。

g95って，ホームページには There are still bugs buried here and there. （あちらこちらに
バグがまだ埋まっている）と書いてありますが，どういうバグなのでしょうか？
コンパイルや計算が正しく行われなかったりするのでしょうか？

>> 43
ttp://www.polyhedron.co.uk
にベンチマークも含めてのっていますよ

>>48
情報ありがとうございました。
早速見てみます。

あるsubroutineの中からその属するsubroutineと同じsubroutineを呼び出したのですが
どうすればいいでしょう？

下記のようなプログラムを書きたいのですが
コンパイルできません。

program main
call hoge
end
subroutine hoge
write(6,*) "hello world"
call hoge
return
end

>>50
f90の場合だと再帰的呼び出しにすれば大丈夫でしょう。サブルーチン宣言の頭に recursive をつけるだけでいいんじゃないでしょうか？

recursive subroutine hoge

>>50
ちゃんと再帰から抜けるようにしとけよ。
無限ループやるなよ。

ちょっとお聞きしたいのですが，write文でcsv形式を書くにはどうしたらいいでしょう?

write (10,'(50d15.7)') TT,(U(j)/K(j), j=1,NMax)

上のようなwrite文を書くとします．その際，csv形式で出力したいのですが，'(50d15.7)'というような書き方をするときに「カンマ」をどうやって紛れこませるか?--ということをお聞きしたいのです．

write(10,1000) TT,(U(j)/K(j), j=1,NMax)
1000 format(d15.7,',',d15.7,','d15.7)
というような形式で書くのならわかるのですが...何十回も繰りかえすのはメンドイもので．御存じの方よろしくお願いします．

>>53
自己レスです．

write (10,"(49(d15.7,','),d15.7)") TT,(U(j)/K(j),j=1,NMax)

と書けばうまく逝きました．
失礼しました．

すみませんご教授下さい。
Fortranで作られたプログラムWindowsで（DOSで）動かしたいんですが、一応データファールとかも揃えて走らせると、
最初の１行だけ読んで「cannot open file "F90.EER" #__」みたいなメッセージが出て計算されません。
これはいったいどういうことなのでしょうか？

ちなみにファイルはここの「TWINSPAN」です。
↓
http://cc.oulu.fi/~jarioksa/softhelp/softalist.html

>>55
もっと詳しく

>>56
すんません。

昔、同じプログラムをNEC98（５インチ時代）に使ったことあってそん時は問題なかったんですよね。
DOSプロンプトのあとにexeファイル名（TWINSPAN）を入れて実行させると、inputデータファイルを尋ねてくるのでそのファイル名を入れてリターンすると
「Cannot Open Error file "F90.EER", error #151」
てのが出て終わってしまいます。

中身はそのまま触ってないので、これ以上のことはわからないのですが・・・

>>57
> 「Cannot Open Error file "F90.EER", error #151」
"F90.EER" ではなく "LF90.EER" では？

Run time error messages from LF90
ttp://agbu.une.edu.au/~kmeyer/errmsg.html
> 151 File specified STATUS= "OLD" doesn't exist (see "Input/Output" in the Essential Lahey Fortran 90 Reference).

ようするに、入力ファイルが無い or 名前が違う。

>>58
ありがとうございます。

そうです「LF90.EER」でした。済みません。
入力ファイルがないとは、inputデータファイルが読めていないということでしょうか。
データ行の１文目は読んでいるんですが。
それか、もはや元のファイルがないのかな。
済みませんあまりに素人で。

gfortran(GCC4.0)でnamelist読みこませた人います？
&abc
de=5
fg=3.1
&end
だと読み込まず
&abc de=5 fg=3.1/
でもだめでした。
方言があるんですかね。

↑です。
使っているプログラムでうまくいかなかったのでテストプログラムを作ってみました。
その結果、結果が異なりましたので訂正します。
以下のプログラムコンパイルして実行するとエラーになります。
FC3　gcc4-gfortran-4.0.0.29.x86_64.rpm


programe main
integer i,j,k
dimension j(10)
character*20 c
namelist /n1/i,j,k,c
read(5,nml=n1)
write(6,n1)
end


./hoge
&n1
i=20
k=3
j=1,2,3,4,5,6,7,8,9,0
c='test'
&end
cannot match a namelist variable

g77だとコンパイルして実行すると、正常に動きます。
いろいろ探したらg95に行き当たりました。これだと正常にコンパイル,実行できました。
g95(gcc 4.0.0 20050124 (g95) Mar 12 2005)

gcc4.1コレクションも試しましたが、このバグは、直っていませんでした。

ずいぶん前のことですみませんが，ご存じのかた教えて下さい。

MS-FORTRAN Ver.5.1 は Fortran 90 に似て非なる物だったという記憶があるのですが，
ふと思い立って，8年くらい前に買った森北出版の 「FORTRAN 77 と Fortran 90」 の
「２部：Fortran　90編　Fortran　90とMS・FORTRAN」を見て
自分で確かめようとしたら，本棚から行方不明になってて
確かめることができませんでした。

MS-FORTRAN Ver.5.1 は Fortran 90 と同一のものだったのか，
あるいは Fortran 90 に独自拡張を加えただったのか，
それとも，似ているけど違うものだったのか
ご存じのかた，教えて下さい。

>>62
実際にMS-FORTRANを使ったことはないですのが
手元に「FORTRAN 77 と Fortran 90」があったので
確認してみたところ、MS-FORTRANは
"Fortran 90 の一部機能を含む独自拡張仕様"
のようです。

F90との主な違いは
・継続行の指示記号は "&" ではなく "-"
・構造体の宣言は TYPE ではなく SURUCTURE
・ポインタ、WHERE文、モジュール等がない
などだそうです。

>>63
解説情報をアップして下さって，本当にありがとうございました。助かりました。
日曜の朝早くから妙な話題につきあわせてしまい，すみませんでした。

再度すみません。
もうひとつ教えていただきたいことがあるのですが，「FORTRAN 77 と Fortran 90」 の本で
扱われている MS-FORTRAN について，本の中でバージョン番号は書かれていますでしょうか？

>>65
書籍では「Microsoft FORTRAN V5.1」を用いたということで、巻末の参考図書に
　マイクロソフトの「Microsoft FORTRAN - Reference -」(1991)
が挙げられています。
ちなみに Fortran 90 の方は「NAGWare f90 Compiler」でテストしたということです。
（2003年5月20日発行の第１版第７刷で確認）

>>66
詳細な情報をいただき，感謝申し上げます。
二度もお答え下さり，本当にありがとうございました。

GUI　WindowsアプリにFORTRANのソースを組み込みたいのだけれど、どうしたらいい？

>>68
�@VFで、GUI WIN32アプリをFortranで組む。
�AC/C++でGUI組んで、C/C++からFortranの関数を呼ぶ。

>>68
VFならDLLをつくるという手もある。

>>68
CVFならCOM DLLをつくるという手もある。

>>68
ほとんどのコンパイラが DLL を作れます．
あとはシンボルを調べれば好きな言語から呼べます．

Intel fortran compiler 8.0 for linuxを使っています。
プログラムの中でサイズの大きいファイルを開こうとすると、
forrtl: File too large
forrtl: severe (30); open failure , unit 10, file 〜
というエラーが出て止まってしまいます。
このような場合はどう対処すればいいのでしょうか。
どなたかご存知の方がおられましたら、どうかよろしくお願いします。

>>73
100Gのファイルでも開こうとしたのか？

約4Gです。

2G以下のファイルに分ける

例えば、「Input No.:」と出力させた直後に改行せず入力したい場合には
どのようなコードにすればいいのでしょうか？
いろんなテキストを見たけど載ってない…。

あと、もう一つあります。
「write(*,*) 'へろー,わーるど？'」と
「print*, 'へろー,わーるど？'」とでは本質的にどう違うんですか？

>>77
advanc='NO'を出力時に指定する。

>>78
77か90か

write(*,*)　顔に見える
print*,　顔に見えない

>>79
ヒント、アリガトウございました。使っているのは 90 の方です。
「character NO」を宣言しなければならないんですね。
あと、「advance」ですよね。

「write(*,'(A)',advance='NO')」のように使っていますが、「'(A)'」は
文字列の指定を意味しているんですか？
例えば「'(I)'」の時なんかは整数かな、と推測できるのですが「'(A)'」に
ついては馴染みがないんで教えて欲しいです。

２GB以上のファイルを読むためのオプションみたいなものはないものでしょうか。

>>82
一般的な32 ビットプラットフォームのLinuxでは，最大ファイルサイズが2Gバイトに制限されているので、

�@64 ビット環境でコンパイルする。
�A32 ビットアーキテクチャなら、LinuxのカーネルをLarge File Supportが有効な新しいものに変える。

アドバイスありがとうございます。
こちらでも調べてみましたが、Large File Supportはカーネル2.4から組み込まれているみたいです。
今私が使っているlinuxのカーネルは2.6なのですが、Large File Supportは有効にはなっていないのでしょうか。
それとも他に何か原因があるのでしょうか。

http://vs6.f-t-s.com/~otonavip/src/vip0658.mp3

休憩しようよ

FORTRAN やらないか？

http://members.at.infoseek.co.jp/kitaurawa/cgi-bin/wiki.cgi

fortran用のeclipseプラグインだって

http://www.photran.org/

>>87
ダウンロード画面にベトナム料理屋が表示されるね。
フォーのレストランてことか。

スレ汚しですみませんが…… Fortranプログラミングに煮詰まったときの
現実逃避で『 諸君 私はFortranが好きだ 』というネタ演説を作ってみました。
先月，私がガイドライン＠2ch板に投稿したものに少々手直し。
（ >>63 & >>66 さん，先月は失礼いたしました。実はあの朝はこれを打っていました。）

ttp://www.geocities.co.jp/Technopolis/4404/i_love_fortran.html

スレ汚し，大変失礼いたしました m(_ _)m"

↑をギコナビじゃなく IE で見て，いま初めて気がついた‥‥
何故 h が抜けてたんだか‥‥ orz

演説 『 諸君 私はFortranが好きだ 』
http://www.geocities.co.jp/Technopolis/4404/i_love_fortran.html

再度のスレ汚し，本当に！ 本当に失礼いたしました m(_ _；)m"

課題で、適当な数字をキーボードから入力して、
文字列に変換して出力するサブルーチンを作れ
といわれたんですが、だれか教えてくらさい。

Fortress

http://www.sun.com/presents/minds/2005-0302/

>>91
ヘンな課題ですね。

>>91
VAX FORTRAN という Fortranコンパイラの文法を解説した Webページの
http://www3.jpn.hp.com/upassist/ktree/cr/vms/v030/v030-04-08.html
で，「メモリ内部入出力」という方法で行う 数値型 ⇔ 文字型 の型変換が
紹介されていました。

要するに，WRITE文を使って「メモリ内部入出力」というのを行うことで
数値型から文字型へ変換しています。すなわち，

WRITE(文字型変数名,FORMAT文番号) 数値型変数名

というようにコーディングすると，数値型変数の値が文字型に変換できるというものです。
以下にプログラム例を示します。
（WindowsXP上の Cygwin の g77 で動作確認しました。
ただし，他のどんな Fortranコンパイラでも動く普遍的な Fortran文法なのか
どうかは知りません。ご存知の方，フォローして下さると助かります。）

● 整数の場合
例えばマイナス符号を含めて５文字以内の整数なら
      WRITE(moji,100) seisuu
  100 FORMAT(I5)
とコーディングして型変換します。

【プログラム例１】

      PROGRAM main
      INTEGER seisuu
      CHARACTER*5 moji
      CALL int2char(seisuu,moji)
      STOP
      END

      SUBROUTINE int2char(seisuu,moji)
      INTEGER seisuu
      CHARACTER*5 moji
      WRITE(*,*) 'マイナス符号を含めて５文字以内の整数をキー入力せよ'
      READ(*,*) seisuu
      WRITE(moji,100) seisuu
  100 FORMAT(I5)
      WRITE(*,*) 'あなたが入力した整数は',moji,'です'
      RETURN
      END

● 実数の場合
例えばマイナス記号を含めて整数部が４文字以内，小数部が３桁以内の実数なら
      WRITE(moji,100) jissuu
  100 FORMAT(F8.3)
とコーディングして型変換します。

【プログラム例２】

      PROGRAM main
      REAL jissuu
      CHARACTER*8 moji
      CALL real2char(jisuu,moji)
      STOP
      END

      SUBROUTINE real2char(jisuu,moji)
      REAL jissuu
      CHARACTER*8 moji
      WRITE(*,*) 'マイナス記号を含めて整数部が４文字以内で'
      WRITE(*,*) '小数部が３桁以内の実数をキー入力せよ'
      READ(*,*) jissuu
      WRITE(moji,100) jissuu
  100 FORMAT(F8.3)
      WRITE(*,*) 'あなたが入力した実数は',moji,'です'
      RETURN
      END

文字に変換してなんか意味あるのかな？

訂正というか補足というか蛇足というか。
>>94 での

WRITE(文字型変数名,FORMAT文番号) 数値型変数名

の「FORMAT文番号」の部分は正しくは，

WRITE(文字型変数名,FMT=書式識別子) 数値型変数名

ですね。FMT= を省略してコンマのあとにすぐ書式識別子を置いても可。
なので，>>95 のプログラム例１での

      WRITE(moji,100) seisuu
  100 FORMAT(I5)

の二行は

      WRITE(moji,'(I5)') seisuu

また，>>96 のプログラム例２での

      WRITE(moji,100) jissuu
  100 FORMAT(F8.3)

の二行は

      WRITE(moji,'(F8.3)') jissuu

というように一行にまとめてもOKです。
むしろ 文番号を使わないで済むので，特に Fortran 90 以降なら
こちらのほうがすっきりして良いかもしれません。

>>97
僕も最初そう思ったんですが，ついさっき，
「commonで定義した変数を文字型に変換してファイル名に使いたい」という例を見つけました。
http://okweb.jp/kotaeru.php3?q=797367

>>91さんには，その課題の出された背景を教えてほしいです。
どんな使われ方をするのか興味があります (^ ^)

>>94
> ただし，他のどんな Fortranコンパイラでも動く普遍的な Fortran文法なのか
> どうかは知りません。ご存知の方，フォローして下さると助かります。

‥‥に対する自己レス （うざくてすみません）。

WRITE文とREAD文の装置番号のところに文字変数を置くやり方は，
培風館から出版されている「Ｆｏｒｔｒａｎ９０プログラミング」の
著者・冨田博之先生の同名サイトの中の
http://www.ip.media.kyoto-u.ac.jp/htomita/sec7.html#7.2.
で，「内部ファイルの方法」として紹介されていました。
↓
（装置番号が）文字変数：プログラム中で使われている文字変数に読み書きできる。
　　　　　　　　　　　　　　　 「内部ファイルの方法」といい，必ず書式を指定すること。

＞＞８９さん
サンクスです。助かりました。課題は、fortranを覚えるために出された数問のなかにありました。
何に使うかはわからないのですが、私がc言語を知っているという背景を配慮して、
cとfortranを関連付けるためなのかとも思いました。私自身もまだよくプログラムに関しては
よく理解してないのでわかりませんが、研究とかにはまず使わないと思います。

本当になんでこんな変な問題がでたんでしょうね。
とにかくありがとうございました。

>>101 egaさん
いえいえ，お役に立てたようで何よりです。
こちらこそ，今まで知らなかった WRITE 文と READ 文の使い方が判って勉強になりました。
僕は Fortran と VB がメインなので C のことはよく判らないのですが，
VB なら Str 関数で一発だから，Fortran にも似たのがあるんだろうと思って
軽い気持ちで調べ始めたら，これが見つからない見つからない (^ ^；
やっぱり Fortran は科学技術計算用の言語なんだなと改めて思いました（笑）

でも、write文で文字列変数をファイルに見立ててカキコするなんてヘンなん仕様だな。
全角で数字入れたらどうなるんだろう？

>>103
WRITE(moji,100) seisuu
100 FORMAT(I5)
に行く前にすぐ前の
READ(*,*) seisuu
のところでエラーになるんじゃないの？

>>90 の Fortranが好きだ演説、笑ったｗ
そういや意外と今まで無かったんだな。Ｃ言語が好きだバージョンなら見たことあるが。

ところでＦ言語って初めて知ったが、拡張子は何なん？

ちなみに、内部ファイルはFORTRAN77でもある。
それ以前でもあるかどうかは知らない。

文字から数値への変換はたまに使った。
CSV読む時とか。

ちょっと質問、

複数の実数型多次元配列(a, b, c,...)を１つの１次元配列にパックしたり、
アンパックしたいのだけど、できれば、元の多次元配列の形も一緒にパックして

{ aの次元数, size(a,1), size(a,2), ... , a(1,1,...), a(2,1,...), ...
bの次元数, size(b,1), size(b,2), ... , b(1,1,...), b(2,1,...),...
cの次元数, size(c,1), size(c,2), ... , c(1,1,...), c(2,1,...), ....}

みたいにしたい。 Fortran90では無理？

win上で動く早いコンパイラで無料のってなんですか？
cygwinは遅いと聞いたので不安になって・・・・

>>108

>>2 の 「x86で動作するFortranコンパイラ」 に評価版などの紹介がある。
http://www.nminoru.jp/~nminoru/programming/x86-fortran.html

あと，Salford の FTN95 Personal Edition （英語版） なら 個人用途限定・非商用用途・サポート無し で 無料・無期限
http://www.salfordsoftware.co.uk/software/downloads/compilers.html

>>108

>>3 に 「フリー／無料評価版のコンパイラ」 一覧があったや。

もちろんそれらのサイトは見ました。
質問がわかりにくくてごめんなさい。

現在数値計算をしているんですがcygwinは
FORTRANコードをF2Cみたいなもので一度C言語に書き換えて
それをコンパイルするみたいなこと（よって遅い）を聞いたので、
無料の中でも早いもの。。はなんですか？という質問です。

ごめんなさい。

>>111
PolyhedronのベンチマークではFTN77/95ともにg77より良好な結果が出ている。
OpenWatcomもほとんどの場合、g77よりも速いようだ。

実際にどの位速くなるかは自分の環境で試してみないとわからない。
#物によってはほとんど変わらない／かえって遅くなるということもなきにしもあらず。

とにかく速い物が欲しいのであれば財布と相談して商用の物を購入するか
Linuxを入れてIntel Fortranの非商用版を使う方がいい。

>>112
ありがとうございます。
FTN95入れてみます。

>>112
すみません，Linux版 Intel Fortran の非商用版（無期限フリー？）はどこで入手できるのでしょうか？
Intel（日本）とエクセルソフトのサイトを見ても情報が見つからなかったので‥‥。

>>114
Free Non-Commercial Software Development Products for Linux
ttp://www.intel.com/software/products/noncom/index.htm

ここから登録とダウンロードが出来る
…はずなんだけど、今つながらない？

>>115
アクセスできました。情報ありがとうございました m(_ _)m"

ftn77やftn95を使ってMS-DOS上のコマンドプロンプトでプログラムを
コンパイル、実行したいのですが「ftn77(ftn95) 〜〜.f」
とうつだけでは〜.exeの実行ファイルができません。〜.obj
というファイルはできるようなのですが。どのようにしたら
コマンドプロンプト上でコンパイル・実行できるのかどなたか
教えていただけますか？

>>117
ググッてみたところ，FTN77 についてなら，
ttp://www28.tok2.com/home/kasayan86/fortran/fortran.htm
の ページ中ほどの 『２　Salford FTN77』 で簡潔に説明されていました。

>>118
この方法でftn77のみでなくftn95に関しても実行ファイルを
作成してプログラムをまわすことができました。
問題も解決して気持ち良く寝れそうです。大変ありがとうございました。

>>119
FTN95 については先ほど http://www.silverfrost.com/23/ftn95/ftn95_documentation.asp から
Salford FTN95 User Guide (Including Clearwin+ and Plato3 documentation) をダウンロードして
ざっと見てみたんですが，
・The /LINK compiler option
・Compiler options on a command line
という２つのオンラインヘルプによると，

FTN95 <pathname> /LINK

としてやれば，コンパイラが <pathname>.F95 (or .F90) という名前のファイルを探してコンパイルした後に
リンカが自動的に起動して，<pathname>.EXE を生成するらしいです （この場合は .OBJファイルは保存されないらしい）。

例えば，MYPROG.F95 から MYPROG.EXE を生成したいなら

FTN95 MYPROG /LINK

但し，こちらの手元では試していないので間違ってたらすみません。
また，僕の訳も間違っているかもしれないので，念のためにご自身で上述のユーザ・ガイドをダウンロードして
・The /LINK compiler option
・Compiler options on a command line
その他，コンパイルに関係しそうなヘルプを検索して，一読されることを強くお薦めします。

>>120
確かにftn95 〜〜.f /link とすると直接実行ファイルが生成されました。
まだ、プログラムを始めたばかりでまだまだ分からないことだらけで。。。
ご丁寧に説明して頂きありがとうございました。
ヘルプを一読してみたいと思います。

FORTRAN90で、

INTEGER :: ARRAY(5)
ARRAY(1) = 1
ARRAY(2) = 2
ARRAY(4) = 4
DO I = 1,5
WRITE(*,*) ARRAY(I)
END DO
END

のように、ところどころ値のない配列をＤＯループで出力
（値の無いときは空行）したいのですが、値のないところで
Reference to undefined character
のエラーを出さないように対処するにはどうしたらよいかどなたかご教授下さい。
コンパイラはSalfordのFTN95です。

-999とかですべてを初期化して、-999なら空行を書けばいいだけでは
-NaNが入ってるともかぎらんしね。

>>122
write文の引数に「err=行番号orラベル」を追加して、
そこでエラー処理(今の場合は空行出力)。

>>122
INTEGER,PARAMETER :: IUNDEF = HUGE(1)
INTEGER :: ARRAY(5) = IUNDEF
ARRAY(1) = 1
ARRAY(2) = 2
ARRAY(4) = 4
DO I = 1,5
IF (ARRAY(I) /= IUNDEF) WRITE(*,*) ARRAY(I)
END DO
END

コンパイラの未定義チェックの多くは、へんてこな値を未定義に予約して
それを参照したりすると、警告を出す仕組み。上記のように、自前でやってもいい。

普通、不定値があるような配列を操作しないだろ。データ構造から勉強したほうがいいかと。

初心者質問に答えてくださりありがとうございました。たすかりました。

Opteronつかってるのですがインテルのfortranコンパイラでも最適化は適切な物になるのでしょうか？

>>128
あまり速くならない。
最新のAMD64ならよいかも。
それよか、gfortranかg95をつかったら。
gfortran(gcc 4.1)は、g77の1.5倍位速い。

gfortranは、バグがある。
ひとつは、namelistで配列を使うと正しく読めない。
もうひとつは、清水氏の作った、column-row読み込み速度プログラムが
g77 col 254MB/s row 641MB/s
g95 col 227MB/s row 612MB/s
なのに
gfotran col 37MB/s row 35MB/s
col 65MB/s row 58MB/s
などばらつくし、桁が違う。
どっかおかしい。

手許のプログラムはgfortran-4.0ではコンパイル出来なくて，
g95だと出来る．g95のサイトから，毎日のようにコンパイル
済みバイナリをダウンロードしている．

自作PGの精度をH/W、コンパイラ両面でいろいろ評価しています。

CompaqコンパイラやIntelコンパイラの場合、
実数や整数を0で(勝手に)初期化してくれてるんですが
PGIやXL Fortranだと初期化してくれない場合があって
精度があいません。

もちろん初期化処理を入れてないのが悪いっちゃあ悪いんですが、
このようなコンパイラの細かい仕様についての情報がある
サイトご存じの方おりましたら教えてください。

ヨロシクお願いします。

130だけど
5月1日版のgcc-4.1のgfortranでは、namelistのバグは、なおっていた。
2つめの、配列のアクセス速度のほうは、まだなおっていない。

>>132
初期化しないおまえが悪い。
Fortran90なら宣言文で初期化できるので、修正は簡単だよ。

integer :: i=0, j=0

とかね。

＞PGIやXL Fortranだと初期化してくれない場合があって
＞精度があいません。

それは変数初期化ではなく、
使用する浮動小数点フォーマットとか、ライブラリの違い
に起因するものだと思った。

FORTRANのよい入門書を教えてください

>>136
入門ということで、

ttp://www.amazon.co.jp/exec/obidos/ASIN/4894713217/qid=1115760491/sr=8-1/ref=sr_8_xs_ap_i1_xgl/250-9355971-7357068
読みやすいが、付録の組み込み手続きの説明が、親切じゃない。

ttp://www.amazon.co.jp/exec/obidos/ASIN/4890524932/qid=1115760948/sr=1-9/ref=sr_1_10_9/250-9355971-7357068
非常に良いが、廃刊なのか？

はじめまして、教えていただきたいことがあり書き込みします。
計算が途中でおかしくなりデータが NaN になった場合，プログラムが終了するようにしたいのですが、良い方法はないでしょうか？
データの上限下限を設定し, その値からはずれた場合に止まるようにプログラムを書き直すことも考えましたが、
NaNが現れた場合に強制終了するための、もっと簡単な方法がありましたら教えていただけると助かります。
よろしくおねがいします(ﾉД`)　

>>138
コンパイルオプションで一発。
どういうオプションかは、コンパイラのマニュアルに書いてある。

＞139さん
-fpe0 で良いのでしょうか
最適化の後ろにオプションをつけてみたのですが、何かうまくいかなくて
すみません、よく分かってませんで･ﾟ･(ﾉД`)･ﾟ･｡　

すみません、また138です。もう少し勉強してみました･ﾟ･(ﾉД`)･ﾟ･｡　
-qflttrap で終了させることができることまでわかったのですが、
最適化と同時に使った場合、0/0では終了するのですが、1/0で素通りしてしまいます。
最適化と同時に使う場合、どうすればよいのでしょうか。
なんども申し訳ありませんがよろしくお願いします(⊃д⊂)

>>141
1/0はNaNじゃなくて、Infじゃないのかな？

っていうか、もまえがどこのコンパイラ使ってるかわからんし、
俺の家に、お前が使っているコンパイラのマニュアルがあるわけじゃない。

がんばって、マニュアル読んでくれ。

>>138
つ IF (ISNAN(x)) STOP

>>142
がんばって読んでみてます･ﾟ･(ﾉД`)･ﾟ･｡　
コンパイルオプションでできることがわかっただけで大収穫です。
ありがとーー

>>143
それをやると、すべての変数に対して全部書いていかないといけなくなるので･ﾟ･(ﾉД`)･ﾟ･｡　

質問です。
intel Visual Fortran ver.8をインストール後
ハロープロジェクトで試そうとしたのですが、
以下のようなメッセージが出て終了してしまいます。
誰か詳しい人お願いします。

-out:hello.exe
-subsystem:console
-entry:mainCRTStartup
hello.obj
LINK : final error LNK1104: コンパイラは、ファイル'ifconsol.lib'を開くことができません。

LNK1104というのはどうもVisual C++のiostream関連のエラーであるらしいのですが、よく判りません…

ソースコードは以下のものです。

PROGRAM HELLO
INTEGER I
DO I=1,3
WRITE (*,*) 'HELLO WORLD!'
END DO
END PROGRAM

ifconsol.libはIntel付属のコンソールライブラリ。
LNK1104とは単にライブラリファイルが開けないと
VCが言っているだけ。ライブラリのパスが通っていない
とかVCとFortranのバージョンが違うとか。

>>145
つ 「インストールと入門」 ガイド 3-8〜3-9

簡単にマルチスレッド化する方法を教えてください

研究者、大学院生、プロのプログラマのいずれかの方にお聞きします。
興味本位の質問ですが、貴方の職業（ご身分）もあわせて教えていただ
けるとうれしいです。

１）いま仮に何かの物理現象を扱っている1000行ぐらいのフォートラン
プログラムを「ホイ」と渡されたとき、このプログラムの中で実際
何がどういう風に処理されているか、読み解き、必要に応じて
中身をいじれるくらいになるまで、どのくらいの時間がかかりますか？
数時間？数日？数週間？数ヶ月？永遠に無理？

２）上記の作業を効率よくするのにどのようなツールを使っておられ
ますか？例えば、印刷したものをじっくり見るとか、テキストエディタ
で画面を見ながらとか、デバッガーみたいなのとか、その他いい方法が
あったら教えてください。

３）こういったプログラミングの業界は世界と比べて日本は進んでいる
と感じますか？

ちなみに自分は学生なのですが、ある程度フォートランの文法の基礎は
分かってるつもりです。でも、1000行の見ただけで、もうついてけない、
と思ってしまう根性なしで最近落ち込み気味です。

hpfの良書を教えてください

ファイルの読み柿がよくわかりません。
どこかによい解説はないでしょうか？

自由形式でかけるFORTRANはありませんか？

千行の道も一行から。

FTN77 Personal Edition Compilerをダウンロードしたいんですが
download FTN77 Personal Editionをクリックしても
名前やメールアドレス，住所などを入力するフォームが表示されないんですが
どうすればいいんでしょうか？

>>155
http://www.polyhedron.co.uk/salford/products/ftn77_personal/

>>156
ありがとうございます！

教えて君ウザイ

intel fortran (Linux版)の質問です。
ランタイムオプション -check をつけてコンパイルして、実行すると、

forrtl: warning (402): fort: (1): In call to [サブルーチン名], an array temporary was created for argument #1

という警告がでます。マニュアルによると、
『ルーチンが呼び出される前に、実引数が一時的な記憶域にコピーされる場合、ランタイム情報を含むメッセージを発行します。』
ということらしいのですが、何故引数がコピーされるのでしょうか？Fortranって call by reference (call by value の対義語)じゃないの？また、コピーされることの弊害として、どのような事が考えられるでしょうか？

>150
まあ似たような立場のモンです。

あっしもいくつかプログラム渡されて数週間解読しましたでつ。
やったこととしては印刷してじっくり読むのは基本として

文法が分かるのなら
・似たような処理をしてるところがあったら吐き出す結果が違わないように
注意しつつ出来る限りサブルーチンや関数にまとめる。
まあ、これはあっしが渡されたプログラムは結構スパゲッティだったから通用したのかも
知れませんが…

あと、
・研究で使ったプログラムならその結果を元に論文が書かれてると思うのでその論文や
そこに上げられた参考文献も読んでみる。
というのも重要かも。あっしはこれで転送行列法なるものを知ったです。

ついでに質問です。

……
real(8),dimension(1:3)::bar
……
call foo(…,bar)
……
end program

subroutine foo(…,bar)
real(8),dimension(:),intent(INOUT)::bar
……
end subroutine

という感じにプログラムを組んだんですが、実行してみるとbarの値や形状が
正常に受け渡しされてないことが判明しました。原因がさっぱり分かりません。
何か分かる方がいらっしゃいましたらご教授お願いします。

使用コンパイラはIntel Fortran Compiler8.1 for Windowsです。

>>161
interface文書いてやらんと、整合配列は通らないんじゃないの？

>>149,153
つ [ifort]

>>151,152
つ [google]

>>159
>何故引数がコピーされるのでしょうか？
サブルーチン側が連続した配列を想定しているのに
連続していない（可能性がある）配列が渡されたから。

program test
　integer, dimension(5,4) :: a = 1
　integer, dimension(3,2) :: b = 2
　call foo(a(2:4,2:3))
　call foo(b)
　print *, a
　print *, b
contains
　subroutine foo(x)
　　integer, dimension(3,2), intent(inout) :: x
　　x = 3
　end subroutine foo
end program test

>Fortranって call by reference (call by value の対義語)じゃないの？
この場合、ちゃんと「一時的な配列の先頭アドレス」が渡されます。

>また、コピーされることの弊害として、どのような事が考えられるでしょうか？
配列のコピー・書き戻しによるオーバーヘッドの増加、使用メモリ量の増加等。

>>164
わかりやすい例示をありがとうございます。
具体的には、以下のコードで困っています。
ちなみに、character型配列ではなく、integer型にすると、一時配列は作られません。
また、以下の例では、b は2次元配列ですが、1次元配列に変更しても一時配列は作られます。
--------------------------------------------------------
program test
character,dimension(:,:),pointer :: b
allocate(b(3,2:3))
call foo(b, size(b))
print *, b
end program test

subroutine foo(x,nmax)
integer,intent(in) :: nmax
character, dimension(nmax),intent(inout) :: x
x = '3'
end subroutine foo
--------------------------------------------------------
実行すると以下のようになります。

$ ifort -check test2.F90
fortcom: Info: Feature not yet implemented: check pointer
fortcom: Info: Feature not yet implemented: check shape
$ ./a.out
forrtl: warning (402): fort: (1): In call to FOO, an array temporary was created for argument #1

333333

x = '3'

のとこが

x(:) = '3'

と解釈され実行されたんじゃないの？
つまり、２＊３の配列すべてに’３’が書き込まれたと。

character,dimension(:,:),pointer :: b

って明示的に書くと

character(len=1), dimension(:,:),pointer :: b

なんだけど、これでいいんですか？

>>165
character型は特殊なんだよ。とても変な仕様になってる。

gcc4.0.0で同じようなコードをコンパイルしてもfortranの方がC++より高速に感じるのだが

厳しいフォートラン講座はどこかにありませんか？

>>169
つ ttp://garota.fismat.umich.mx/mn1/tutor_fort/index.html

>>170
たしかに厳しい。何かいてあるのかまったくわからん。

>162
この件についてはご教授の通りで解決いたしました。
どうもありがとうございます。

>>168
そりゃそうでしょ。

alphaプロセッサを使っている場合、速度の最適化を最高にするオプションは何ですか

>>174
-fast -arch host -O5
ただ、これが最速とは限らない。

fortranはO5まであるのね。C++はO3までしかないね

C++は-O5まであるよ。
ただマニュアルには
-O4 and -O5 have the same effect.
ってあるので、実質-O4まで。

LINUXでFORTRAN７７のコンパイル〜実行ができません。

[k@localhost k]$ gedit test.f
[k@localhost k]$ g77 test.f

となって何もなりません。
a.outは作成されるのですが、クリックしてもなにもなりません。
よろしくお願いしします。

$ ./a.out

>>179
出来ました！
学校のだと　./　必要ないので分かりませんでした。
ありがとうございました。

>>180
おやおや。その学校のマシンって結構危ないのではなかろうか。

なんで./が必要なのかを知ってれば別にいいんじゃね。

しかし、、「クリック」なのかぁ、、。

>>182
すいません！
知りませんｗ

カレントディレクトリはコマンドパスに入れてはいけない。

変数を NaN で初期化するには、どうしたら良い？

>>185
NaNでそんなｱﾌｫなことを？

>>186
いや、つい出来心でつ。
NaNをフラッグにして、条件分岐できるかなっと。

フラッグを -1 とか 999999 とかにすると、
いわゆる2000年問題的な問題がおきるしね。
フラッグの変数(boolean)をわざわざ作るのも面倒だし。

>>167
character型は腐ってるということで、良いですかね。

>>188
Cだとchar型の配列と変数は整合性があるけど、fortranの場合はそれが無い。
つまり、

character(len=10) choge

と

character(len=1) choge(10)

は全く違う内部構造になってるんだよね。

ファイルの読み書きがよくわかりません。誰か教えてくだすれ

f77で動的にメモリ領域を割り当てるにはどうしたらよいのでしょうか？

>>190
open, read, write を勉強しる。

>>191
できない。f90からできる。

HP FORTRANだと構造体を使ったプログラムがf77でもコンパイル
出来たっけな。あくまでHPの独自拡張だけど。

>>189
FortranのCHARACTERは長さの属性を持っている。
「文字型」と言っても実質的には文字列型になる。

高級言語という観点からすれば整数型 charの配列を文字列とする
Cの方が手抜きである、とも言える。

Fortranのsubroutineをc言語から呼び出すプログラムをWINDOWS XP用のフリーのコンパイラを
使用してコンパイルしたいです。
どのコンパイラを使ったらよいですか？

>>195
cygwin

OPEN文でファイル名を指定すれば、新しいファイルを作ることはできますが、
ディレクトリを作成するにはどうしたらよいのでしょうか。
どなたかお教えくだされたまえ。

>>197
call system('mkdir AAA')

めちゃくちゃ処理系に依存しそうだな

http://www.imslab.co.jp/products/products.html#F*FLOWER

こんなツールでf90対応のソフトを知りませんか？

組み込み関数としてRAND()が認識されないのですがなんでなんでしょうか？
ヘルプを参照したりググッてみたけどわかりませんでした。

使っているのは
Microsoft Developer Studio 97
Visual Fortran Professional Edition5.0A
というやつです。

http://web.archive.org/web/20040223051246/http://www1.jpn.hp.com/products/software/developme

198さん。できました。どうもありがとうございました。

>>198
call systemを使わないやり方はありますか？

初心者の質問です．
LAPACK95を入れてテストしたらこんなエラーが出るんですけど
詳しい人おしえて

>g95 sample01.f90 -llapack95
In file sample01.f90:2

　use la_precision, only: WP => SP
　　　　　　　　　　　　　　　　　　 　 1
Fatal Error: Can't open module file 'la_precision.mod' at (1)
for reading: No such file or directory

ファイル la_precision.mod は c:\g95\lib にちゃんと存在するんですけど

>>204
おそらくないと思うよ。

でも、たいがいのコンパイラで
call system('mkdir AAA')
または
ii=system('mkdir AAA')
のどちらかはサポートしてると思う（mkdirがない場合は知らんが）

>>204
処理系依存のサブルーチンを使う。
Visual Fortran なら MAKEDIRQQ とか。

>>206,207
どうもです。以前、中身の不明なディレクトリ内のファイルを全て読む
プログラムを作った時も、命令を見付けられなかったので仕方なく

call system("ls hoge/ > list")
open(io,file="list")

とかやっていました。lsが無いMS-DOSプロンプトだったら実行出来ない
ですね。

処理系依存性を吸収するためには、プリプロセッサを使う。

おれなら、ディレクトリ操作なんかはPythonスクリプトなんかでやって、本質的なところでFortranを呼び出す形で解決するけどな

Python つかうと、逆にポータビリティーが下がらないか？
基本は sh だろう。

>>205
モジュールファイルのパスはライブラリとは別に指定します。

>g95 -Ic:\g95\lib sample01.f90 -llapack95

または
set G95_INCLUDE_PATH=c:\g95\lib
としてから
>g95 sample01.f90 -llapack95

>>211
じゃあ、ディレクトリツリーを全て巡回しながらデータファイルを探し出し、それを読み込んでfortranで処理させるのshで作れよ。

#!/bin/sh
~/bin/getfilelist ~/dat > filelist.txt
~/bin/hoge filelist.txt

shの機能以外を使ってるから×

Pythonっての知らないんだけどさ、
例えばどんなスクリプトになるか
見せてちょ。

210=213ですか？
211が言っているのはshのみを使うというよりは
sh+一般的にどこでも使えそうなコマンド(ls等)で作るという意味では。

外部プログラム呼ぶのがｵｹｰなら、Pythonスクリプト呼ぶのと大差ないでしょ？
つまり、sh自体にディレクトリ巡回機能が付いてるのかどうかなんですが

シェルスクリプトは複数のプログラムを組み合わせてバッチ処理を行うための
ものだと思うんだが。

ようは目的に合わせて自分の環境で使える道具を選択すればそれで良いんだよ。

結論を急ぐと、

Fortran単体では柔軟な処理は不可能

ということです

一昔前ならともかくFortran単体で全てやる必要性もないし。
Fortranに限らず何でもかんでも単一の言語でやろうと余計な策を弄するのは本末転倒だし。

それがFortranクオリティ

フォートラン　(・∀・)ｶｯｸｲｲ

Fortranは真の高級言語だな。

Fortranは真の高級言語だな。

Fortranは真の高級言語だな。

Fortranは真の高級言語だな。

Fortranは真の高級言語だな。

Fortranは真の高級言語だな。

そうですが何か？

Fortranは真の硬旧言語だな。

ＳＵＢＲＯＵＴＩＮＥに二次元配列が渡せないよ。なぜ。

>>221
だな。
ポータビリティーのことも考えると、
なるべく使いそうな計算環境にありふれた言語を使うことかな？

以前、Cで書いたコードをFortranからリンクしたら、
共同研究者から怒られた。Cが読めないんだって（笑）
そいうのは、例外的事例ってことで。

複数の人間でつくる時はどんな言語を使うかも含めて
事前にルールを決めとかないと…

これからＦｏｒｔｒａｎをやろうとするのは無謀ですか？
Ｃ言語はコンソールアプリケーション程度、しかもあまりできない。
Ｂａｓｉｃはアプリケーション開発をしばらくやってたけど理解できてるのかどうかわからなくてやめた。

何故やりたいのか
ちょいと興味があった。さわり程度に使ってみたい。
後はＭｓ−Ｄｏｓのシンプルな画面に惚れている。

>>235
無謀じゃないよ。古臭い手続き型言語だから、本やWebでいくつか例題やればすぐ習得できる。
でも、あんまり面白い言語じゃないし、エレガントでもない。

formatなどを使って、２×２の行列を表示させるにはどうしたらいい
のでしょうか？ほぼ初心者です。

FORMAT文の勉強をするといいと思います

勉強してもわからんとです。

宿題は自分でやろうね

Intel Fortran Compiler 9.0 出たようです。
ttp://www.intel.com/software/products/compilers/flin/flinux.htm
ttp://www.intel.com/software/products/compilers/fwin/fwindows.htm

何かよくなったんやろか、、。

Linux超初心者です
大学の卒研でFortran(今のところ77で十分)が必要になったんですが、インストールの段階で早くもつまづき、どうにも前へ進めません
gccが何のことかすらよくわかりません
Linuxユーザの方HELP ME

Linuxのインストールがうまくいかないという話ならLinux板で質問したほうが
いいと思うよ。

>>243
つ http://www.google.co.jp/search?hl=ja&q=g77+linux&lr=lang_ja

>>241
ヒロシです。
9.0 のリリースの件，ダイレクト電子メールで勤務中に知ったとです。
3月に 8.1 を買ったばかりだったので，ちかっぱショックだったとです。
ヒロシです‥‥

安値でアップグレードできるんだろうか？
今まではどうだったか，ご存じのかた教えて下さい。

>>246
インテルにユーザー登録すれば、購入後１年間は無償で最新版をダウンロード出来ます。
２年目以降は保守ライセンス契約の更新（有償、新規価格の約半分）という形で
このサポートサービスを継続していくことが出来ます。

ttp://www.xlsoft.com/jp/products/intel/register.html
ttp://www.xlsoft.com/jp/products/intel/support/index.html

>>247
あっ Σ（￣□￣ ；
そういえばそうでした！ そのことをすっかり忘れていました。
情報ありがとうございました。大感謝 m(_ _)m"

Fortranで作ったプログラムの中に外部ソフトを組み込むことって可能ですか？
exeでしか存在していない外部ソフトです。
自分で作ったプログラムならソースがあるから、include出来ました。
ソースが無くても可能なのでしょうか？

>>249
dllをリンクとかで出来るでしょ

>>250
返信ありがとうございます。でも、その方法がわかりません。
状況としては、
Ａ（自作Fortranプログラム）とＢ（外部）があって、
Ａで計算した結果を用いて、Ｂで計算をする。そしてＢでの結果をまたＡに返して計算。
Ａは起動させたまま、ということを考えています。

dllはどうやればリンク出来るのか教えていただけますか？

>>251
>Ａで計算した結果を用いて、Ｂで計算をする。そしてＢでの結果をまたＡに返して計算。
>Ａは起動させたまま、ということを考えています。

大抵の処理系では外部コマンドを起動するためのサブルーチンが用意されているので
それを使えば良い。

>>252
ごめんなさい。ＶＦ６．６でオンラインヘルプとか見たのですが、わかりません。
普通のサブルーチンみたいに、
call (B.exe)
は無理ですよね？
どう記述するのですか？またコンパイル時に何かやるのですか？

VFならこんな感じ↓

program test
　use DFPORT

　integer :: ios
　integer(kind=4) :: r
　character(len=80) :: str

　r = system("dir > flist.txt")
　print *, "終了コード = ", r

　open(unit=1, file="flist.txt", status="OLD")
　do
　　read(unit=1, fmt="(A)", iostat=ios) str
　　if (ios /= 0) exit
　　if (str(1:1) /= " ") print *, trim(str)
　end do
　close(1)
end program test

詳細はオンラインヘルプ等で「SYSTEM」または「SYSTEMQQ」を参照のこと。

g77で乱数の発生のさせ方を教えてください！

　　　　　　real x
integer l
do=1,5
call rand(x)
write(*,*) x
enddo
end

とすると、サブルーチンの使い方がおかしい！みたいなことを言われてしまいます。
よろしくお願いします！

自己解決しました！
　　　　　　real x
　　　　　　integer l
　　　　　　call srand(1234)
　　　　　　do=1,5
　　　　　　　call rand(x)
　　　　　　　write(*,*) x
　　　　　　enddo
　　　　　　end
こんな感じであってますよね？
失礼しました

あれ？コンパイルできない・・・。
　　　　　　real x
　　　　　　integer l
　　　　　　call srand(1234)
　　　　　　do=1,5
　　　　　　　x= rand()
　　　　　　　write(*,*) x
　　　　　　enddo
　　　　　　end
とするとできるのですが、普通のサブルーチンと扱いが違うのですかね？
あたふたしてすみません。

今現在ＣＰａｄで一人で勉強中。
このコンパイルの瞬間がＣ言語でも味わったけど、好きだ。
無機質にＤｏｓ画面で表示されている瞬間なぜか技術者気分だ。
コード自体は入門サイトからのだからそんな難しいものではない。
だけどいつかは大学の言語解析ソフトくらいはつくれるようになりたい。
・・・というかプログラムか。
なんか英語辞書が組み込まれててテキストの中に入っている文書から英単語を読み出して、それの個数を読み出すソフト。
そしたらメモ帳に出力させてた。
一定のスペースをつけて読み出してるからエクセルでも処理可能。
それくらいは作ってみたい。
ただそのソフト、日本語はサポートされてないからひとつの文書にある文字ベースのエントロピーが求められん・・。
日本語入れた瞬間バグ起こしたもんで教授から大目玉とはこれいかに・・。

>>257
　　　　　　do=1,5

のとこじゃない？

古いフォートランのコードで、Ｃに直したいのですがコードが読めません。
もし分かる方いらっしゃったら確認お願いします。
以下コード

FUNCTION RAND(K)

INTEGER K, M, CONST1
REAL RAND, CONST2
PARAMETER (CONST1 = 2147483647, CONST2 = .4656613E-9)
SAVE
DATA M /0/

IF(M .EQ.0) M = K
M = M * 65539
IF(M .LT. 0) M = (M + 1) + CONST1
RAND = M*CONST2

END

以下私の理解している事： 参考にしたページ： http://www.g.dendai.ac.jp/lecture/fortran/fort12.html

任意のＫが関数に渡される。
関数は Ｋ＊６５５３９＊ＣＯＮＳＴ２をいつも返す。

でもいくつか疑問点があります。
どうして型宣言が行われていないのでしょうか？
どうして関数は特定の値をいつも返すのに、ＩＦの分岐が必要なのでしょうか？
困ってます。アドバイスをよろしくお願いしますｍ（＿＿；）ｍ

「暗黙の型宣言」てやつで、RAND自体は単精度実数になってる。
関数が持ってるMという変数の値はSAVE宣言されてるから、関数から抜けても値が保存される。
呼び出されるごとに関数が返す値は違う。IF文は乱数発生アルゴリズムを理解すれば意味がわかる。

>>260
・DATA 文による初期化は、プログラムの実行前(?)に１度しか実行されない

・SAVE 文があるので、関数内の変数はすべて static

（例外はあるが、FORTRANのローカル変数はSAVEがあろうが無かろうが
すべて static と考えたほうがよい）

> （例外はあるが、FORTRANのローカル変数はSAVEがあろうが無かろうが
> すべて static と考えたほうがよい）
>

え−っ！　すごい昔の話じゃないの？

address errorがでるときの原因を教えてください！
お願いします。

>>264
OSやコンパイラの種類を教えてくれないと分らん。

それでもエスパーなら、エスパーなら何とかしてくれる

配列の範囲外をアクセスしようとしたとか

>>265
ウィンドウズです。
コンパイラは７７だと思います。

それはコンパイラとは言わない

>>264,268
これらの情報から考えられる原因：

お前のプログラムがカスだから

address errorというのはコンパイル時に出るものなのか、
実行時に出るものなのか。address errorというだけで自明の
ことなのだろうか。

【科学技術計算】C言語とFortran　どっちがいい？
http://science3.2ch.net/test/read.cgi/sim/1042715803/476
476 名前：名無しさん＠３周年[] 投稿日：2005/06/22(水) 17:26:14
fortran使って作ったプログラムがコンパイルできるんだが
実行するとaddress errorって出るあるよ
どうしてあるか？

>>263
たしかに、保証はされてないですね。
そのための save 属性だし。

>>261
>>262
ありがとうございました！！周りに知ってる人がいなくて困ってました。
m（＿　＿）m　重ね重ねお礼申し上げます。。

質問なんですが、Windows で動くフリーの Fortran9x コンパイラ(FTN95など) で
Lapack をライブラリにして使用することってできるんでしょうか？

現在はcygwinでg77+Lapackっていう組み合わせで使ってますが、
もともとはC/C++使いだったんで、Fortran77が使いにくく感じて仕方ないです…。

Intel Fortran Compiler買えばいいのに…

>>276
Intelのコンパイラは速いらしいですね。
Math Kernel Library なんていうありがたいものまで販売されてるし。
しかし、私の所属する研究室（私は大学院生です）では
行列演算をするのは自分だけだし、他の人はFortran77で満足してるみたいなので微妙です。
一人でFortran77派の先生を説得できそうにはないですね･･･。

>>275
MinGW版のg95でBLAS/Lapack/Lapack95のコンパイルは一応エラーなしで通った。
（ライブラリのみ）
ただ、テストはしてないので正常に機能するかどうかまではわからん。

>>277
Linux使えば？compilerもMKLも学生ならただで使えるよ。

レスありがとうございます。
>>278
g95プロジェクトって遅々として進んでないような印象があるんですけど、
いつになったら完成するんでしょうね。
完成or完成間近で、バグがほとんどないようなら使いたいと思ってるんですが･･･。

>>279
MKLもタダでしたか！
これを機にLinuxとWinのデュアルブートにしちゃおうかな…。

ttp://www.g95.org/
を見ると毎日のように更新されているね。
Linux上で動かしている手元のプログラムは、gfortran-4.0では
コンパイル出来なくて、g95では大丈夫。実行するとAlphaの
Compaq Fortranとほぼ同じ結果を返してくれる。Intel Fortran
でもコンパイル出来たけど、実行すると異常終了して、デバッグを
諦めた。

Intel Fortran9 買った。Visual Studio .NET上で使うと訳分らん。今までCVF6.5使ってたから。

俺もVF6.6使ってたけど、まだ慣れないよ。

IntelというかVC.NETで
6列目と73列目以降の色分けがしたいんだけど、
できるのかな・・

>>283
固定フォーマット(*.for, *.f)ならVS.NETで桁の色分けしてくれてますが。

　　WindowsでプログラミングをしてみようとSalford FTN95 Personal Editionをインストールしたのですが、
コマンドラインでコンパイルしようとして「ftn95 test.f」のように入力してもEXEなどの実行ファイルが作成できないんですがこのやり方ってなんかまずいのでしょうか？
なんか画面みると問題ないようなのですが肝心の実行ファイルがないんです。
しかも「/Lgo」というコンパイル成功時に一度実行させるコマンドラインオプションをつけるときちんとプログラムは実行されているんです。
　　なんかものすごく初歩的な質問のような気がするのですが、調べても逆に初歩過ぎるのか解決方法が見つかりません。
すみませんが解決方法を教えてください。お願いします。

>>285
/LINK オプションは試した？

>>286
　解決しました。ありがとうございます。
　意外と簡単なところに限って躓くものですね…。
いままでunixで使っていてそのコンパイラーがオプション無しで実行ファイルa.outを作ってくれてたので盲点でした。

コンパイルして./a.outで出力すると「セグメンテーション違反です。」とか出るんですけど、解決策は無いですかね？

>>288
再現出来る最小限のソース

逆にＣをＦｏｒｔｒａｎへ変換するソフトってあるんですかね？

>>290

>>6 の C2F

Salford FTN95で実行ファイルを作ろうとしてCommand Promptで、
「ftn95 "main.f90" "sub.f90" /link」
と打ったら、
more than one source file specified
　　1 error [main.f90] compilation failed
とかいうエラーメッセージが出たんですが、
複数のファイルに分かれた場合のコンパイル法ってこれじゃまずいんですか？
とりあえず、Salford以外じゃ問題ないんで、ファイル自体には問題なくて
単純にSalfordの操作ミスのようなのですが。

>>292
個別にコンパイルして最後にリンク

ftn95 main.f90
ftn95 sub.f90
slink main.obj sub.obj -FILE:prog.exe

面倒くさかったら mk32やバッチファイルで処理。

OPEN文を使って次のプログラムを実行してみたのですが、うまくいきません。
(「中断しました。」とでます)
どなたかよろしくおねがいします

OPEN(1,FILE='701')
READ(1,*)X,Y
S=0.0
DO 100 I=1,15
READ(1,*)U,V
S=S+(U-X)*(Y+V)/2.0
X=U
Y=V
100 CONTINUE
WRITE(*,*)S
CLOSE(1)
END

ちなみに計算したいデータは↓のものです。
ヒントなどでもうれしいです。

ttp://lcofao.hp.infoseek.co.jp/Output.pdf

>>294
見た目間違ってなさそう。

それだけ短いプログラムだったら
計算過程を全部writeで出力したら？
変数の計算式やreadの式の直後に全部writeで
write(*,*)"S=",S
みたいに出力すればいいんじゃん？

>>296さん
返信ありがとうございます。
大学でしかフォートランできない状態なのでまた明日確認してみます

ひょっとしたらOPEN対称ﾌｧｲﾙがダメなのかもしれないですよね。
別テキストで数字だけを羅列していったのですが。。。
↓みたく。
0,10.5
125,8.2
307,7.5
420
以下同様

>>294

READ(1,*,END=200)U,V

.
.
.

100 CONTINUE
200 CONTINUE

>>293
ありがとうございました。

>>298さん
返信頂きありがとうございます。
READ(1,END=200)ですよね?
明日さっそく試してみたいと思いますが、見てみた感じ
READ文の中にこのような命令が入っているのは初めてで、
本やウェブでいろいろ見てみたのですが、どういう処理をしているのかよくわかりません。

あるところで、READ文とCLOSE文の間にWRITE文があるとうまくいかないことがある。とあったのですが、
それと関係あるのでしょうか?

>>300
ttp://globe.nagaokaut.ac.jp/~kumakura/cecp/b3_1/fortran/io.shtml
ttp://globe.nagaokaut.ac.jp/~kumakura/cecp/b3_1/fortran/

どうもありがとうございます。
今日学校でやってみたところ、解決しました。

>>302
>>298の方法で実行時エラーが出なくなったと言うことは
データの個数が足りないか、正しく読み込まれていない
可能性があるのでは？

>>303
実はOPEN(1,FILE='701')
の部分ですが、保存していたのが701.TXTだったので、
それがアボートの原因だったのかも知れません。

そこを変更して、教えていただいたやり方でやっても、元の形でやっても
うまくいきました。

すいません。
自分でなんとか書いてみたんですが、合ってるかどうか見てもらえますか？
本当ならこれを学校のFortranでやってみて、確かめて、gnuplotでグラフ。
とまでやりたかったんですが、学校のが使えませんでした。（家も使えません。）
どなたか数値計算していただけない・・ですよね。。
もし計算結果がでたらエクセルでグラフを書けばなんとかできるかなと思っているんですが、
↓こんなんです

program prog036
implicit real*8
write(*,*) 'a, b, n, dt = ?'
read(*,*) a, b, n, dt
write(*,*) 'a, b, n, dt = ',a,b,n,dt
c
x1 = a
x2 = b
y1 = 0
y2 = 0
c
do k = 0,n
call subroutine bibunkoshiki(x,y,z,dt)
end do
end
c ***********************************************************
subroutine bibunkoshiki(x,y,z,dt)
c ***********************************************************
x1 = x1 + y1*dt
x2 = x2 + y2*dt
y1 = y1 + z1*dt
y2 = y2 + z2*dt
c
write(*,*) x1,x2
end
c ***********************************************************
function z(x)
c ***********************************************************
implicit real*8(a-h,o-z)
z1 = (-11*x1 + x2)/10
z2 = ( x1 - 11*x2)/10
end

>>305
[Salford FTN77/Win32 v4.03, Copyright (c) Salford Software Ltd. 1988-1998]
Licensed to: FTN77 Personal Edition
Department: Non-commercial use only
0002) implicit real*8
*** Invalid list
1 ERRORS [<PROG036>FTN77 Ver 4.03]
0025) end
WARNING - Z2 has not been given an initial value (see -undef and
-only_undef options)
WARNING - Y2 has not been given an initial value (see -undef and
-only_undef options)
WARNING - X2 has not been given an initial value (see -undef and
-only_undef options)
WARNING - Z1 has not been given an initial value (see -undef and
-only_undef options)
WARNING - Y1 has not been given an initial value (see -undef and
-only_undef options)
WARNING - X1 has not been given an initial value (see -undef and
-only_undef options)
NO ERRORS,6 WARNINGS [<BIBUNKOSHIKI>FTN77 Ver 4.03]
0032) end
0033)
WARNING - X1 has not been given an initial value (see -undef and
-only_undef options)
WARNING - X2 has not been given an initial value (see -undef and
-only_undef options)
WARNING - Z has not been given an initial value (see -undef and -only_undef
options)
NO ERRORS,3 WARNINGS [<Z>FTN77 Ver 4.03]

*** Compilation failed

うはｗひどいなｗ
>>307さんありがとうございますｗ
なおそなおそ（ｿｿｸｻ..

（´・ω・`）ｻｯﾊﾟﾘ

>>305

つ >>3

積分を行うときに区分数を大きくしすぎると、逆に誤差がでてきてしまうのは
どうしてですか?

サブルーチンの並びはこんな感じでいいですか?
IMPLICIT REAL*8(X,Y)
REAL*8 FUNCTION FUNC(X,Y)
FUNC=SIN(X)+COS(X)
RETURN
END

>>311
つ ttp://maya.phys.kyushu-u.ac.jp/~knomura/education/numerical/text4/node1.html

>>312
つ ttp://www.g.dendai.ac.jp/lecture/fortran/

>>313さん
ありがとうございましたが、
REAL*8 FUNCTION FUNC(X,Y)
IMPLICIT REAL*8(X,Y)
FUNC=SIN(X)+COS(X)
RETURN
END
となるんですね?

てか関数の場合サブルーチンいいませんね。

ケタ落ちは整数同士の割り算のとき、あまりが切り捨てられること。→×割り算だけでない
ケタ溢れはおもに掛け算で発生するが、加減や割り算で発生することもある。→○
整数の加減ではケタ溢れが起こらない限りは正しい結果である。→○


ですよね？

ふつう、実数として浮動小数点演算を使うとき、気をつけるものだと思うが・・・

独学で初心者ですが，「外部関数を用いてeuler法により微分方程式を解く」
といったプログラムを作ってみましたが，エラーがでてしまいます．
エラーを見たのですが，どうしてエラーなのかいまいちわかりません．
使用しているコンパイラはvisualfortranです．どうしてエラーが出るのかわかる方，回答よろしくお願いします．

program euler
implicit none
integer :: i
real :: h

interface
function f(x,y)
implicit none
real :: f
real,dimension(1:1000),intent (in) :: x,y
end function f
end interface

x(1)=0.0
y(1)=1.0
h=0.1

do i=1,1000
y(i+1)=y(i)+h*f(x(i),y(i))
x(i+1)=x(i)+h
write(*,*) 'i=',i
write(*,*) 'x=',x(i)
write(*,*) 'y=',y(i)
end do

end program euler

function f(x,y)
implicit none
real :: f
integer,dimension(1:1000),intent(in) :: x,y
f=0.5*(1.0+x)*y*y
return
end function f

長くてごめんなさい．エラーはこんな感じです．
C:\Program Files\Microsoft Visual Studio\MyProjects\7_11_2\7_11_2\euler.f90(14) : Error: This name does not have a type, and must have an explicit type. [X]
x(1)=0.0
--------^
C:\Program Files\Microsoft Visual Studio\MyProjects\7_11_2\7_11_2\euler.f90(15) : Error: This name does not have a type, and must have an explicit type. [Y]
y(1)=1.0
--------^
C:\Program Files\Microsoft Visual Studio\MyProjects\7_11_2\7_11_2\euler.f90(19) : Error: This name has not been declared as an array. [Y]
y(i+1)=y(i)+h*f(x(i),y(i))
--------^
C:\Program Files\Microsoft Visual Studio\MyProjects\7_11_2\7_11_2\euler.f90(20) : Error: This name has not been declared as an array. [X]
x(i+1)=x(i)+h
--------^
C:\Program Files\Microsoft Visual Studio\MyProjects\7_11_2\7_11_2\euler.f90(32) : Error: The shapes of the array expressions do not conform.
f=0.5*(1.0+x)*y*y
--------^
Error executing df.exe.

7_11_2.exe - 5 error(s), 0 warning(s)

>>317
･･･エラー見たら英語分からなくても、変数XとYがおかしいことくらいわかるでしょ？
参考までにexciteでエラーを翻訳しといたよ。

This name does not have a type, and must have an explicit type
→ この名前には、タイプを持たないで、明白なタイプがなければなりません。
This name has not been declared as an array.
→ この名前は配列として宣言されていません。
The shapes of the array expressions do not conform.
→ 配列式の形は従いません。

2番目の訳がexciteにしてはうまく訳せてるよね。
これを参考にXとYの何がおかしいか考えてみて。
頑張れ！

>>321
ありがとうございます．
確かに変数xとyに問題があることはエラーからわかるのですが，
配列や型は宣言しているのに，なぜエラーなんだろうと思ってしまいまして・・・
もう少し粘って考えてみます．

>>322
そういう疑問だったのか。
じゃあ function f(x,y) の2つの引数x,yを、 function f(a,b) のように別の文字a,bを使って
書き換えてみればエラーの原因が分かるかもね。

ちなみにこれ、iを1000まで回してるけど、i=20程度でy(i)が発散するよ。

>>323
ありがとうございます．
色々チャレンジしましたがだめでした・・・
このような簡単なプログラムが組めないようでは
まだまだ勉強不足ですね．
また明日から頑張ります．

>>315
あってますか?
過去問にあったので教えてください。

>>315,325
○×はあってるが、1問目の理由が違う。
1問目は「丸め誤差」、桁落ちではない。

tp://e-words.jp/w/E6A181E890BDE381A1.html

>>326
ありがとうございました。
ややこしいですが、大分とわかってきた気がします。

最初にキーボードからintegerを読み込ませて(例えば100とします)、
保存対象のファイルをこのintegerによって変えたいのですが(この場合'100.dat'に保存したい)、
下のようにするとうまくいきません。どこがおかしいのでしょうか？
教えてくださいお願いします。
implicit none
integer :: i
character*16 :: name
read(*,*) i
name = (char(i))//'.dat'
open(11,file = name)
write(11,*) 'ここに書き込む内容をかく'
close(11)
stop
end

非対称行列の固有値と固有ベクトルを求めるprogramを
作りたいのですが、どこかにプログラム例などありませんか？
ありましたら、教えてください。

>>329
LAPACK ( http://www.netlib.org/ ) は？
固有値問題の数値計算について研究してるんじゃなかったら、
自分で作ると言わずそのまま使った方がいいよ。
作るの面倒だし、作っても確実にLAPACKの方が速いし。

>>328
Fortranには整数型や実数型の変数を文字型変数に変換する関数はありません（多分）。
代わりに「内部ファイルの方法」というのを使います。
こういうのです
↓
WRITE(文字型変数名,書式識別子) 数値型変数名

例えば整数型変数 i ( = 100 とする) を，３文字の文字型変数 kazu に変換するとしたら，

implicit none
integer :: i
character*7 :: name
character*3 :: kazu
read(*,*) i
write(kazu,'(i3)') = i
name=kazu//'.dat'
open(11,file = name)
write(11,*) 'ここに書き込む内容をかく'
close(11)
stop
end

ただ，こちらの手元に今，Fortranコンパイラがないので，試せていません。
ダメだったらごめんなさい。

【参考】
http://www.ip.media.kyoto-u.ac.jp/htomita/sec7.html#7.2.

>>329
ARPACKとか使えば？

>>328
char()って、「文字コード→文字(≠文字列)」な変換関数じゃない？
ASCIIコードなら、100→「d」。

>>331
内部ファイルを使うなら >>328のプログラムで
> name = (char(i))//'.dat'
の部分を

WRITE(name,'(I3,A)') i, '.dat'
または
WRITE(name,'(I3,''.dat'')') i

とすればＯＫ

>>329 です。
LAPACKをインストールするまでに時間を使ってしまいそうですが、やってみます。

>>331 >>333 >>334
みなさんありがとうございます。

メモリの動的割り当てと解放について、初心者でもわかるようにわかりやすく説明されたサイト
はないか教えて下さい。
テストがあるんで急ぎです

>>337
宿題は自分でやろう。
これを見てゆっくり考えよう。
http://h18009.www1.hp.com/fortran/docs/lrm/dflrm.htm

サブルーチンの並びを考える問題なのですが、これで問題ないですよね?

REAL*8 FUNCTION FUNC(X,Y)
IMPLICIT REAL*8(X,Y)
FUNC=SIN(X)+COS(X)
RETURN
END

よろしくおねがいします

試してから書き込めばいいじゃん

IMPLICIT は NONE 以外使うな。
変数宣言使え。

implicit none 以外の implicit の使い方初めて見た希ガス

age

プログラムを動かすと、REAL*8 FUNC(A,A)のAの値を変えるごとにプログラムの速度が違う様なんです。
配列数を変えることによってプログラム実行速度はかわってくるのでしょうか？
6000＞10秒
5800＞5秒
5600＞10秒
5400＞4.8秒
もう何が何やらorz

環境はWindows pen４　メモリ１Gの環境です。

>>344
配列各要素に操作をしているなら実行速度に正の相関関係が生じるのは当然ですが、
秒数を見るとそれに起因したものではなさそうですね
気になるのは「その秒数はどうやって／何を計ったか」なのですが。
マルチプロセスのOSで応答時間を測定する事はあまり意味がない　ということをご存知か？

秒数は目安としてdtimeで計測しています。

繰り返し計算で一回の計算が終わるたびに値をDos窓に表示させているのですが
それを見ていても明らかに遅いのが判るので・・・・・。

同じ配列数場合、いつも実行しても、時間は変わらないのでしょうか？

はい、変わってないですね。6000で10秒くらい、5400で5秒付近です。
dtimeなので信用はできませんが
明らかに時間差があるため・・・謎ですね。

サブルーチンの時間計測してるのですがやはり一から書き直したほうが良さそうですねorz
解はまちがってないんだけどなぁー。

解決しました！

サブルーチンのプログラムの中にある配列数を曖昧に今まで取ってたんですが
一番余分に配列組んでそうなところを少なくしてあげたら6000にしようが5100にしようが
4秒代で回るようになりましたorz
きちんとした原因がそこかは判らないのですがコンパイルの種類、f77,f90にしたりといろいろ
変化させた結果、それで安定したので・・・。

>>345、347　ありがとう。

キャッシュのせい？

夏休みも残りわずか。
皆さん、宿題は終わりましたか？

>351
宿題っつうか学会発表のためのプログラム(fortran製)の
ぶん回し中ですが何か？

>>352
何学会？ 一寸興味あり （・∀・）

すみませんがちょっと教えて下さい。
学生向けだと無料になる Linux 用 Fortran コンパイラがあった気が
するんですが‥‥ありませんでしたっけ？
Intel だったかなと思ったのですが，9.0 のアカデミック価格は \46,725 （税込）でしたし‥‥

>353
応物

>354
>3

研究で、Fortran IVで書かれた古いプログラムを使う必要が有ります。
プログラムといえば、関数電卓のBASICくらいしかやったこと無いので、
自分でそのプログラムを書き換えることは無理だと思います。
Fortran IVのプログラムが動くソフトって有りますか？

Windows XPか、Linux上で

コンパイルできればどんなOS上でも動くんじゃないの？

>>356
>Fortran IVのプログラムが動くソフト
それは、通常コンパイラと呼ばれるソフトです。
私があなたなら、該当プラットフォーム（OS）上で動くFortranコンパイラを販売している会社を調べ、
その会社に問い合せる。

ファイルの中身を全部コピーしたいんですけど、
read文を使う場合、読み込ませる変数を細かく設定しなければいけませんよね？

例えば、
"do i=1,10"
を読ませたかったら、
read(1,*)a1,a2となりますよね？

毎回ファイルの中身が変わるので、それは面倒なんですよ。
とにかく、そのままコピーしたいんですけど、いい手ありませんか？

>>360
call system("cp hoge fuga")

>>360
ファイルの中身をコピーしてから「何をしたい」のかを
示さないと、>>361のような恐らく望んでいない回答を
されて、再度質問をすることになりはしないか？

＞361,362
申し訳ありませんでした。
ファイルに、新たな内容を追加したいということです。
A.datというファイルがあり、そのファイルの最下部にフォートランで計算
した値を追加したいのです。

ですから、新たなファイルをopen文で作成して、A.datファイル全体をコピーして
、そして追加内容をwrite文で書き込もうかと考えています。そのコピーの方法
を教えていただきたかったということです。

もし他にわかりやすい方法があれば、是非教えていただきたいです。

ちなみに、コンパイラはCompaq Visual Fortran 6.6.0使用してます。

>>363
単純に A.datに追加書き込みするだけで良いのなら

open(UNIT=1, FILE='A.dat', POSITION='APPEND')
　　：
write(1, *) hoge
　　：
close(UNIT=1)

>>360
ttp://www.uopmu.ees.osakafu-u.ac.jp/~yabu/soft/fortran.html
12. 文字型変数に対する read
15. read 文実行時のエラー処理

でどうですか？

>>360
別のファイル名で作成して、call systemを使わないなら
FORTRAN90で

character(len=100) :: buf
open(io1,file='hoge',status='old')
open(io2,file='fuga',status='new')
do
read(io1,'(a)',iostat=iostat) buf
if(iostat/=0) exit
write(io2,'(a)') trim(buf)
enddo
close(io1)
write(io2,'("tekito")')
....
close(io2)

配列の大きさを求めた後で配列宣言したいのですが、
プログラムの途中で　real::配列　としてもエラーになります。
どうすればよいですか？　

自己解決しました

>>368
エスパー募集をやめる。
エラーと言っても様々なエラーがある訳で。

work_dir に a.out と 01,02…,099,100というディレクトリが既にある
として下さい。その時に、

do j=1,100
i = system("./a.out > memo")
i = system("mv memo ./j")
end do

のような事をしたいのですが、何か上手い方法を御存知の方はいないでしょうか？
system=("")の中で変数を使えないかという事なのですが、どうでしょうか？
ちなみに訳合ってshell scriptで、というのは無理なんです。

>371

内部WRITEで文字列を生成してsystemにわたせば？

372さん、書き込みありがとう御座います。
仰ってるのってこんな感じですか？

character ccc*3
do j=1,9
write(ccc,10)j
i = system("./a.out > memo")
i = system("mv memo ./ccc")
end do
10 format('00',i1)

でもこれでは普通に ccc ってファイルができてしまいます・・・。
変数であることをどう認識させるのですか？？
$ccc でも @ccc でもうまくいってません。
基本的な質問ですいません。。

system("mv memo ./"//ccc)
か
character*100 cbuf
cbuf = "mv memo ./"//ccc
call system(cbuf)
とかでいけるかな

自信はないけど、こんな感じでいけるんじゃね

374様
Thanks a lot !!
これで納期に間に合いそうです！ｗ
ｍ(__)ｍ

VBからFortranのdllを呼び出すサンプルプログラムをうｐしていただけませんでしょうか？
変数（配列になっているやつも）の引渡しの部分も書かれているやつで。

ネットでサンプルプログラムを探しているのですが、なかなか見つからなくて困っております。

>>376
ttp://www.xlsoft.com/jp/products/intel/cvf/docs/vf-html/pg/pg18_05.htm
ttp://www.xlsoft.com/jp/products/intel/cvf/docs/vf-html/pg/pg08.htm

fortran で表形式のデータに対して　SQLの sort やjoin みたいなことが出来るライブラリ等は無いでしょうか？
行列演算と、データベース的な処理の両方が必要なのですが、fortranは行列に偏っているように思われます。
（科学技術分野では当然か）

>>378
使用したことはないので目的に合うかどうか判りませんが
Windows用ならODBCデータベースにアクセスするライブラリがあるようです。

f90SQL
ttp://www.canaimasoft.com/f90sql/index.htm

Windows　Ｆｏｒｔｒａｎによるデータベースの利用 （3月2日 5-10）
ttp://platypus.lns.tohoku.ac.jp/Program.htm#5th

>>379 ありがとうございます。しかしDB操作で無くて、プログラム中でSQL的なデータの加工がしたいので、
それは使えません。自作するしかないかもしれません。

Intel FORTRAN9を購入した。
自動で並列化してくれつという宣伝につられた。
確かにタスクマネージャーをみると２個とも計算しているっぽい。
しかし、COMPAQのCVF8でコンパイルしたものと比べて（これは１つのCPUが動いている）
全然変わらない。
試したのはCAVITY流れの非定常計算（自作）で、同じマシンで両者を比較。
PentiumD 820を残った少ない予算で購入したのに・・・。
これまで使っていたPentium4 560のほうが速い。意味がなかった。
FORTRAN77形式で書いたのがまずいのかな。
（でも並列化してくれているし）

自動並列化がうまくいく確率なんて相当低そう。
PenDを無駄にしたくないならOpenMPで並列化してみ。
たいして難しくないはず。

>>381
自動並列化は全然役に立たない。多分、非常にシンプルなプログラムに限ってその効果が現れると思われ。
382の言うようにOpenMPで自分で並列化していくしかない。

>>383
多分、並列化の技術の使いかたをイメージできてないんじゃない？
自動並列化は役に立たないんじゃなくて、上手く並列化されるようにプログラムを変える
数式を解きやすいように変形するのと一緒

>>384
それでうまくいけば誰もOpenMPやpthread使わない。

数百行のソースなら簡単だけど。

並列化計算は概要を読んだだけでほとんど勉強していない。
インテルのホームページを読むとこれまでのプログラムのままでも
多少の効果は期待できそうなことが書いてあったので。

自動並列化: インテルVisual Fortran コンパイラには、
ループを自動的にスレッド化する自動並列化機能が搭載されています。
このため、並列化によるマルチプロセッサ上でのアプリケーション・パフォーマンスの改善を
手軽に実現できます。自動並列化機能は安全に並列実行できるループを検出し、
これらのループに関してマルチスレッド・コードを自動的に生成するため、
イタレーションのパーティショニング、データの共有、スレッドのスケジューリング、
同期化といった低レベルの細かい作業をユーザが行わなくても並列化が可能となり、
マルチプロセッサ・システムやハイパー・スレッディング・テクノロジ対応システムの強力な
パフォーマンスを活かすことができます。並列実行ループのベクトル化機能や、
ソフトウェア・ベースのスペキュレーティブ・プリコンピュテーション新機能など、
ハイパー・スレッディング・テクノロジ対応インテルプロセッサでのパフォーマンスの
向上につながる最適化の改良が行われました。


それにしても一応デュアルで動いているのに効果が全くないのは
どういうこと？思ったほど高速化しなかったな、というのは
よくあることだけど、計算が振り分けられれば多少は高速化すると
思うのだが。
タスクマネージャーをみる限りでは２個ともフルに動いていそうなのだが。

次のような単純ループのプログラムを自動並列化オプション付きでコンパイルしても、
CPU占有率は75パーセントがいいとこだったよ orz...

program Console1
implicit none
integer i, j, k
integer, parameter :: NLOOP=1000
real(8) :: dum

dum=0.0d0
do i=1, NLOOP
do j=1, NLOOP
do k=1, NLOOP
dum=dum+sin(k*1.0d0)
enddo
enddo
enddo
write(*,*) dum

end program Console1

sinなんて関数呼んだら、ループが並列化していいかどうか
判断できないんじゃ？

行列の積だと１００％になったよ。
ただし、行列×ベクトルだとコンパイルした時点で並列化してくれない。
もっとも多い計算が行列×ベクトルなのに・・・。
たとえば、これだと１００％になった。


PROGRAM REI
IMPLICIT NONE
INTEGER ::I,J
REAL(8),DIMENSION (5000,5000) :: A,B,C
DO I=1,5000
DO J=1,5000
A(I,J) = 1.1d0
B(I,J) = 2.1d0
END DO
END DO
C = MATMUL(A,B)
STOP
END PROGRAM REI

>>389
MKLに入ってるblasだったら並列化されてるから、それ使えばいいんじゃない？

コンパイラに金を出したみたいだし、タスクマネージャとかいってるから
Winユーザでしょ?
てことはMKLも有償じゃぁ…

>>389
って事はベクトルは行列にしてから計算しろってことじゃね？

>>391
お察しの通りWindows版です。
Linuxは扱い方がよくわからず、いろいろ不便なので
（インストールしてネットワークの基本的な設定ができる程度）
軟派ですがWindowsを使っています。
>>392
B(5000,1)にしろということでしょうか？（単純すぎ？）
この方法では並列化しませんでした。
うまいアルゴリズムを考えなくては。
並列化とOpenMP勉強します。
皆様、アドバイスありがとうございました。

>>393
なんとなくだがB(5000,5000)のB(x,x)だけ使えばいいんじゃね？

>>388
自動並列化がそんなにバカなら使い物にならないな。

ご教授おねがいします．

module variables
　　　　real(8),dimension(:,:),allocatable :: x
　　　　integer :: i,j
end module
としてモジュールをつくり，
program main
　　　　use variables
　　　　implicit none
　　　　read(*,*) i,j
　　　　allocate( x(i,j) )
　　　　!ここで配列xに値を代入する操作を行う
　　　　call test()
end program
subroutine test()
　　　　use variables
　　　　implicit none
　　　　!配列xに対する操作・演算を行う
end subroutine

というプログラムを実行した場合，
配列xに代入された値は，サブルーチン内から確実に参照できるのでしょうか？
（正確な値が参照できる保障はあるのでしょうか）
また，モジュール内でallocatable属性の変数を宣言し，
プログラム内で動的確保，サブルーチン等で共有することは
プログラミングの作法として許されているものなのでしょうか？
なお，コンパイラはCompaq Visual Fortran 6.1を使用しています．
よろしくお願いします．

>>396
普通はmodule内に、allocateするためのルーチンを作っておいて、それを呼びます。例えば

module variables
implicit none
　　　　real(8),dimension(:,:),allocatable :: x

contains
subroutine alloc_x(n1,n2)
implicit none
integer, intent(IN) :: n1, n2
allocate(x(n1,n2))
end subroutine

end module

みたいな感じで作っておいて、呼び出し側では

use variables
read(i,j)
call alloc_x(i,j)

で配列確保します。メリットはC++ライクなカプセル化ですかね。

>>397
ありがとうございました．

Fortran90で，とりあえず動くプログラムを作ることばかりを考えてやってきてたので，
一度プログラミング作法やFortran90の仕様などをきっちり勉強したいと思います．

fortran90でOOPぽいプログラミングするなら

http://www.cs.rpi.edu/~szymansk/OOF90/F90_Objects.html
http://www.amath.washington.edu/~lf/software/CompCPP_F90SciOOP.html

も読んどいたらいいかもね。

>>399
そういえば、Fortran2000だったか2005だったかで、OOP支援機能がなんか入るんじゃなかったっけ？

>>396
module variables
　　　　real(8),dimension(:,:),allocatable :: x
　　　　integer :: i,j

CONTAINS
subroutine test()
　　　　use variables
　　　　implicit none
　　　　!配列xに対する操作・演算を行う
end subroutine

end module
とするのがこの場合の本来的な型。

>>396
module variables
　　　　real(8),dimension(:,:),allocatable :: x
　　　　integer :: i,j
!
CONTAINS
subroutine test()
　　　　use variables
　　　　implicit none
　　　　!配列xに対する操作・演算を行う
end subroutine
!
end module
とするべし。
ａｌｌｏｃａｔｅに関しては問題なし。

>>402
subroutine test() 内の use文は不要じゃない？

後、配列の大きさは size() で調べられるから
モジュール内の i, j は不要な気もするけど
そこら辺は好みの問題か。

variable::test() みたいに呼べればいいのにな。

質問

学生時代ずっとf77ばっかりで
会社に入ってようやくf90をやりました。

配列式、メモリの動的確保、do-enddoとか
便利なところは感じますが疑問もあります。

moduleってどんなメリットがあるのでしょう？

使い方も難しくて利点がよくわかりません。

commonとincludeを使わないで良くなるからmoduleのほうが好きだ。
コードもきれいに書けるし。

>>405
メリットとしてはcommonに比べて安全に変数を共有できる、とか
モジュール内のサブルーチンを使う時にコンパイラが引数の型チェックをしてくれる
などかな。

使い方はここらへんを参考にして
ttp://ace.phys.h.kyoto-u.ac.jp/~tomita/education/fortran90/sec9.html#9.3.
ttp://www.mri-jma.go.jp/Project/mrinpd/coderule.html

研究でVAX FORTRANという言語を使わないといけないんですが
開発環境はどこでダウンロードできますか？
探したんですが見当たらなくて　誰か教えてください。。。

>>406
>>407
ありがとう。

理解したとは言えないけど、なんとなくわかった気がします。

学生時代、研究でドイツ人の書いたf90ソースを改良しなくてはならないのに
f77しか知らずに、moduleとかすごい嫌いになったんですよね、、。

勉強してみようかな、、。

>>407
すいません、追加で質問です。
commonはなぜ危険なのでしょうか？

どこで変更しているかが明示されないから

commonとequivalenceはgoto文と同じくらい凶悪です。

>>408
はい、これ↓。Hobbyistなら非商用に限り無料で使えるOpenVMS
ttp://www.openvmshobbyist.com/

AlphaかVAX用だけどね。

VAX互換で良ければIntel Fortranあたりでも使えば？
（一応、DEC→Compaq→HP→Intelの流れでVAX FORTRANがご先祖様だし）

Microsoftがfortranに手を出さないのは何故でしょうか？

>>414
つ 【Microsoft Fortran PowerStation】

>>413
ありがとうございます
助かりますっ

>>415
なんで廃れたのかね？

売れないからでしょ

姫野ベンチってさ、FortranのよりCの方が30パーセントくらい速いんだけど、そんなもんなの？

どのFortranやねん

intel visual Fortranの購入を考えているんですが、これのマニュアルにはネット接続環境でないと
アップデートができないとありますが、インストールしたあとでアップデートを受けたい場合には
ネットに接続し、手動でオンライン認証を受ける必要があるんでしょうか？

Visual FortranオンリーでGUIなプログラムつくっている人っているのかな？

>>421
アップデート方法は基本的に
インテルのプレミアサポートから新しいバージョンを丸ごとダウンロードして
古いバージョンをアンインストール→新しいバージョンをインストール

プレミアサポートにログインするためには予め登録しておいたIDとパスワードが必要。

c2fはどこにあるんでしょうか?
>6のリンク先のリンク先は切れているし
c2f.zipでググってもリンク切ればっかりで見当たりませんでした

>>424
ttp://home.earthlink.net/~dave_gemini/c2f.zip

ttp://www.talkaboutprogramming.com/group/comp.lang.fortran/messages/161252.html
> I lost my copy of C2F.ZIP file when I changed computers/ISP in last year,
> and now have a request to make it available again.

・・・（´Д｀；）

>425
ありがとうございます

簡単な質問ですいませんが、正の実数Xを入力して
その小数部を出力するにはどうすればいいんですか？

簡単な質問ですいませんが、正の実数Xを入力して
その小数部を出力するにはどうすればいいんですか？

(Xの小数部) = (実数X) - (Xの小数部を切り捨てた値)

x = x - real(int(x))

文字処理でcharacter型変数に格納されている文字数を取得したい場合
どうすればよいでしょうか？
C言語ではヌル文字があって比較的簡単だったのですが，Fortranでは無いようなので…

WRITE(*,*) LEN_TRIM(' Hello, World! ')

>>432
即レスありがとうございます．
できました！

>>429さん>>430さん
ありがとうございます。
お礼を言うのが遅くなってすいません…。

Fortran90/95
の入門書としてお勧めの書籍があれば教えてください．

http://www.amazon.co.jp/exec/obidos/ASIN/4894713217/
一応これがいいかなって思っているんですが．

IntelのMKLって試用期間過ぎてもスタティックライブラリをリンクして使えてるんだが、こんなもんなの？

>>435
俺もその本使ってる。
分かりやすくはないんだけど、本屋で立ち読みした他の本は
もっと分かりにくかったのでその本を買ったよ。

でも、ちょっとFortranに慣れてきたらその本より
ttp://www.rs.kagu.tus.ac.jp/yama/f90/f90-lang.html
の方がよっぽど役に立つよ。

>>437
おおありがとう．
リンク先も参考にさせて頂きます．

エラーが出てしまい、自力での解決が難解なので、皆さんの知恵を貸してください。

プログラムの一部に

E(I)=1.0D00/DSQRT(DREAL(X(I,I)))

という行があるのですが、新型のPC&コンパイラにしたところ、今までは出てこなかったエラーが出てきてしまいました。
エラー表示は

fortcom: Warning: ファイル名.f line 2401: This argument's data type is incompatible with this intrinsic procedure; procedure assumed EXTERNAL. [DREAL]
E(I)=1.0D00/DSQRT(DREAL(X(I,I)))
------------------------------^
/tmp/ifortIC65Ir.o(.text+0x18f5d): In function `secula_':
: undefined reference to `dreal_'
/tmp/ifortIC65Ir.o(.text+0x19264): In function `secula_':
: undefined reference to `dreal_'

といった感じです。dreal、dsqrtの両方のdを外せば解決はしますが、精度が低くなりすぎて使えません。
dをつけた状態でコンパイル時にエラーが出なくなるにはどうしたらよろしいでしょうか。

PC構成は、OSがsuse9.1、コンパイラが無料版のintel fortran 9.0です。
古いPCですと、OSがredhat6.2、コンパイラがabsoftのpro fortranです。

よろしくお願いします。

dをつけた状態と書きましたが、精度が失われなければ問題ありません。
あるいはより高い精度になる分にも一向に構いません。

ただ、これよりも精度が劣ることがないようにしていただきたいのです。

>>439
種別型パラメタがあってないのかも。
・E(I)=1.0D00/DSQRT(DREAL(X(I,I),KIND=8))
・E(I)=1.0D00/DSQRT(DBLE(X(I,I)))
のどっちかでコンパイル通ったりしない？

種別型パラメタって、
REAL*8の"*8"とかREAL(8)の"(8)"ってやつね。

>>441
ありがとうございます。
下の方法にてコンパイルが成功いたしました。

現在実際に計算しています。
精度については結果待ちなので明日にならなければ分かりませんが、
いい結果が出ることを祈ってます。

精度に関しても特に問題ありませんでした。
本当にありがとうございました。

E・∇・ﾖﾉｼ　＜444ｹﾞｯﾄ♫

form='binary' で open して読み書きしたいんだけど g95 だとコンパイルエラーになる・・・
何か対策無い？

>>445
form='binary' のかわりに access='stream' を使うとか。

普通、 unformatted使わないか？

>>446
ありがとー。
access='stream' で望んだ結果が得られました。

これいいかも

http://www.nag-j.co.jp/FortranBuilder.htm

WinXPユーザーなんですが、
Salford FTN95でコンパイルしたものと
g95でコンパイルしたexeのサイズがぜんぜん違うんです。（１０倍ぐらい）
（30行くらいのものです）
Salford FTN95はいちいち最初にメッセージが出るんでコンパイラーを変えたいんですが
g95をこのまま使い続けてもいいんでしょうか。

>>450
どっちかがスタティックリンクをしているから大きくなる？

>>451
そうなんですか。気にしないでおきます。たいしたプログラムは組みませんし。
ありがとうございました。

ところで、gfortranとg95はどちらが優れているのでしょうか。
とにかく早く計算したいんです。
別にこの二つに限らないのでWinXPで動くフリーのFortran90コンパイラーを紹介してください。

「早く」ということは、実行速度はともかく短期間で
実行可能なプログラムを作りたいということと考えれば
g95でないかな。gfortranではコンパイル出来なかった
プログラムがg95では出来たことがある。
フリーのコンパイラの件は知らん。

>>452
貧乏なら無償のIntel FortranをLinuxにインストール。
金があるなら、Windows版Intel Fortranを買え。アカデミック版なら大した金額じゃない。
お遊び程度なら、Linux上でgfortranかな。

>>454
Linux上で、っていう要件を満たせるなら、遊びでも素直にIntel Fortran使ったほうが吉かと。

>>452
coLinux上でIFC

>>456
IFCってIntel Fortran のバージョン7までだよね？
8からはifort。7系列をダウンロード出来る場所を
見付けられなかったよ。premierサポートに登録して
いて8と9はダウンロード出来た。

いやバージョン関係なくIntel Fortran Compilerの略でIFCって使ったんだけど。
８以降はifortって言わないといけないのか。

みなさん、ありがとうございます。
LinuxのIntel製コンパイラーが優秀なのは知っているのですが、
なぜかインストールに失敗してしまいます。
Linuxに関してはかなり初心者なのでインストール方法があっているかも
自信はありません。
Vine Linux3.2はとっつきやすそうなので挑戦したいと考えてはいます。
LinuxになれたらIntelコンパイラーを入れてみようと思います。

すいません、わからないので教えてください。
使ってるのは「fortran90」です。

「0度から90度まで5度刻みで、小数点以下5桁の精度の
三角関数sinとcosとtanの表を作れ。ただし、tan90°の
値には・・・・(中点4つ)を表示させること。」

途中まで作ったんですが…、最後の「tan90°の
値には・・・・(中点4つ)を表示させる」がわかりません…。

real, parameter :: rad=3.1415926/180.0
real :: r,s,c,t
integer :: k
print*, "角度 sin cos tan"
do k=0,90,5
r=k*rad
s=sin(r)
c=cos(r)
t=tan(r)
print'(1x,i3,2x,3f10.5)',k,s,c,t
end do
end

ＩＦ文使えば？

FORTRANというより数学の問題じゃないか。
なんで90度のtanだけ特別扱いになるか分かってるか？

こんなの宿題か何かだろ？いちいち助言してたらキリがないよ。
>>460も自分で考えな。
プログラムを書くことからじゃなくて、アルゴリズムを先に考えるんだぞ。

そういえばFortran 2003 には Floating-point exception handling が有るそうだが
どんな感じになるの？
調べてもようわからんかった・・・

intelのサイトしか見てないけど
intelのwindows版コンパイラはvisual fortran版しかないの?
linux経験がないので非商用版を使う前に評価版で性能とかを
win上で確かめたかったけど

>>465
どういう意味？
Windows版でもコマンドライン環境で使えるよ。

>466
http://www.xlsoft.com/jp/products/intel/compilers/index.html
みたらwin版はvisual fortran ってかいてて
linux版はfortranとしか書いていないから
ちがうものかと思ってました。
なんかlinux版にもvisual fortranがあると思いこんでました

質問です。
今までCompaq Visual Fortran 6.6（以降CVF）を使っていましたが
最近Intel Visual Fortran（以降IVF）に乗り換えました。

ところがデバッグの仕方がわからず四苦八苦しています。

例えば
real*4 x(100)
に対して
do i = 1, 101
x(i) = 0.0
enddo
と書きます。
CVFを使って、この状態でデバッグ実行すると、配列オーバーで止まり
その時点での変数iの値などを見れました。

ですが、IVFだと止まることは止まるのですが
変数の中身が見れません。
Visual Studio関連のソースだと思うのですが
Cのソースコードが表示された状態からどうしたらよいかわかりません。

IVFでこの問題解決する方法ってありますか？

マニュアル読むにしてもキーワードもわからないので
何かアドバイス頂けたら、、と思います。

数値解析や3Dグラフの描画にフォートランを使おうかと考えています。
いろいろと調べてみたところ、ヒューリンクスのプロフォートランが
良いかと思いました。
これを使っている人はいますか
http://www.hulinks.co.jp/software/pf/plplot.html

>>468
Intel Visual Fortranでよくありますね。
エラーが出てプログラムが止まった時にDOS窓に戻ったら、
どこでエラーが出たかメッセージが書いてませんか？
サブルーチン名みたいなのと行番号みたいなのがあわせて書いてあると思うのですが。

それがないようなら、PAUSE文をたくさん挿入するなり、
「カーソルの行まで実行」で順に見ていくなりして、
様子を見るしかないと思います。

cvfも残しておいたほうがデバッグがしやすいよな

>>470
>>471
レスありがとうございます。

機能的に無理なのですかね〜。
上位互換ということでCVFでできたことは全部できるかと予想していたのですが。


行番号などは確かにプロンプトに出ていますね。
それがループの中とかじゃなければ、その情報でも作業できるんですけどね〜。

g77のときには無事コンパイルできていたのですが、gfortranになってからコンパイルエラーになります。

INTEGER MASKA
DATA MASKA/'377'0/
END
某プロフラムのソースで見かけたのですが、勉強不足なのか、このDATA行の意図することが分かりません。
MASKA
はINTEGERで、DATAにある
'377'
は「377というCHARACTER」を入れようとしているように読めるのですが…。
「377*0」でもなさそうですし、「3770」なら、そのまま書けばいいように思うのですが、どういった意味になるのでしょうか。

>>473
自己レス

「'377'0」ではなく「'377'O」でした。
で、これは8進定数 (octal constant) というものでした。
gfortranでは「O"377"」と書いてやらないとコンパイルエラーを起こします。

ちょっとお聞きしたいです。fortran77です。
format文でデータを５個横に羅列し、
改行してまた５個、以後繰り返しというプログラムのはずなんですが、４個目で改行、
５個目がひとつ書いてあって改行、以後繰り返し、という感じになってしまいます。
format文なしでやっても一緒でした。
一行あたりの文字量の問題かとも思ったのですが・・
なにかヒントいただけたら幸いです。

>>475
WRITE(*,*) 0.1, 0.2, 0.3, 0.4, 0.5
と、書式を指定しない場合は、（コンパイラ依存かもしれないけど）
1行72文字 or 80文字で自動的に改行が入る。
(ファイルへの出力、標準出力いずれでも結果は同じ。)

WRITE(*,'(4F12.6)') 0.1, 0.2, 0.3, 0.4, 0.5
と、Format文で4つ分の書式しか指定していなければ
5つ目のデータ(ここでは0.5)が改行されてから出力される。

WRITE(*,'(5F12.6)') 0.1, 0.2, 0.3, 0.4, 0.5
と、ちゃんと5つ分の書式をしていれば、必ず改行なしで出力される。
80桁を超えてももちろんOK。

なんで、Format文でデータ5つ分の書式を指定していることをチェックしてみれば？

>>476
ありがとうございました！５つ分の式が書いていませんでした。。
というか、出力６個あったのに５個と勘違いしていました。。
非常に助かりました。ご丁寧にありがとうございました。

Sun、C/C++およびFortranの開発環境「Sun Studio 11」も無償化
http://www.itmedia.co.jp/enterprise/articles/0511/17/news025.html

太っ腹ですね。

>>478
お、すげー！x86でSolaris試したかったんだけど、CコンパイラはともかくFortranコンパイラがないのがネックだったんだよね。

rand()　を使わずに乱数を作りたいのですが、
誰かわかるかた教えていただけませんか？お願いします。

>>480
「乱数」伏見正則，東京大学出版会がいいと思う。

乱数発生するPCIカードとかあるよ。

>>481
ありがとうございます！調べてみます
>>482
できればプログラムの書き換えだけで済ませたいのです、、、。
でも書き込みありがとうございます！

>>481
Numerical Recipe に解説+プログラムが
数パターン載ってる。

なんとか解決できそうです。
みなさんありがとうございました！！

>>485
どう解決するの？

学校でFortranを習っているのですが、レポートがわかりません。

..以下のプログラムを作成しメールで提出
..「級数展開によりsin(xπ) を求め、表示する」
.. x とn は実行時にキーボードから入力
.. 実数の計算はすべて倍精度で行う

という内容です。級数展開の内容は次の通りです。
ｓｉｎ（ｘπ）　≒　ｘπ　-　（（ｘπ）^3）/3!　+　（（ｘπ）^5）/5!　-　（（ｘπ）^７）/7!　+　（（ｘπ）^9）/9!　-　・・・・　（（ｘπ）^（2n-1））/（2n-1）!

授業で習った文法（？）で使えそうなものは、doとifぐらいです。
作ってみたのですが、これで合っているかどうかわかりません。
御指導お願いします。


program taylor_expansion
implicit none
interger :: i, j, n, kai
real(8) :: pi, x, result
intrinsic atan, sin

pi = 4.0D0 * atan(1.0D0)

write(*,*) "input x"
read(*,*) x
write(*,*) "input n"
read(*,*) n

if (n<=0.0D0) then
write(*,*) "error! n must be natural number!"

else

kai = 1
result = 0
do i = 1, n
do j = 1, 2*i-1
kai = kai*j
end do
result = result + (-1.0D0)**(i+1.0D0) * (x*pi)**(2.0D0*i-1.0D0) / kai
end do

write(*,*) "sin(x*pi) is about",result
write(*,*) "sin(x*pi) is",sin(x*pi)

end if

stop
end program

>>487-489
kaiの初期化の位置がおかしい。
このままでは i のループの２回目以降で前回の値が持ち越されてしまう。

それから n を大きくすると階乗の値がオーバーフローする恐れがあるので
kai は倍精度実数の方がいいかも。

>>490
考えてくださって本当にありがとうございます！
「kai」の初期化の位置については全く気がついていませんでした。

result = 0
do i = 1, n
kai = 1
do j= 1, 2*i-1

と書くべきですよね。
明日学校でコンパイル・実行して試してみます。
ありがとうございました。

級数展開でsin関数を計算する意図が良く分からない。手計算しかできないときは役立つけどさ。
その先生、頭悪いんだろうな。
ご愁傷様。

>>492
ただの練習問題なんだろう。
そんなのの内容にまで一々ツッコミ入れてたらきりがないよ。

>>487
それよりも、ここに書き込み出来る環境があるなら
机上だけで悩んでないでフリーのコンパイラでも入れて
自分でも試すようにした方がいいよ。

>>492
287を積んでない286マシンで極座標とか円筒座標の座標計算するときには役に立つと思う。
今はありがたみを感じないかも知れないけど。

たぶんセンセはその時代の人なんでしょう。

>>494
コプロ世代か。

つか、その世代だとテーブルに格納してたと思うけど。

>>490
今日学校でコンパイルして実行してみると、うまく動いてくれました。
まだプログラミングを始めて間もなく、自分で作ったプログラムの中ではいちばん難しかったのでとてもうれしかったです。
アドバイスしていただきありがとうございました。

>>492-493
ただの練習問題です。
フリーのコンパイラを入れて自分のパソコンでもコンパイルできるようにしようと思って試したのですが、
コンパイラの設定がうまくいってないのか、挫折してしまいました。
また今度やってみます。

sage

質問いたします。
fortranからmatlabを呼び出して関数を使えるらしいのですが、うまくいきません。
コンパイル時に undefined symbols として engopen などが認識されないようです。
当方はfortran77 , matlab7.0.1 を使っていますが、そのせいなのでしょうか？
↓参考にしました
ttp://www.cybernet.co.jp/matlab/support/helpdesk/r14/toolbox/matlab/matlab_external/ch_engi5.shtml#25618

real A( 5, 2 )

call s( A )

subroutine s( A )
real A( 4, 2)

・・・

end

上のようなプログラムを作った場合、例えばメインプログラムのA(1,2)と
サブルーチンのA(1,2)は違うものになってしまいます。C言語と配列の
メモリ上配置が違うためというのは分かったのですが、これに対処するために
どうすればよいでしょうか？

・配列の添え字を逆にする。
・最後の添え字の次元をサブルーチンに渡す。

この2つしか思いつかないのですがどちらもわずらわしいです。
なにかいい方法は無いですか？

>>500
最初の次元の長さが固定なら、
PARAMETER (N=5)
とか書いたヘッダファイルparm.hを用意して、
PROGRAM MAIN
INCLUDE 'param.h'
REAL A(N,2)
...
END

SUBROUTINE S(A)
INCLUDE 'parm.h'
REAL A(N,*)
...
END
とか書けば?

Ｃから書き直してるなら添え字を逆にしないと性能が大きく変わる
可能性があることは知っておくべし。

>>499
ライブラリがつながっていないのではないかな？
HELPのUNIXでのコンパイルとリンクないしWindowsでのコンパイルとリンクのところをみて
ライブラリのリンクするための環境変数の設定を確かめてみれば？

>>500
FORTRAN90以降ならサブルーチン側での宣言を
real A(;,;)
で次元だけ合わせる形にしておけばとりあえず呼び出し側と同じサイズに自動的になる。

FORTRAN77なら整合配列を使う。呼び出し側で配列の大きさも渡す。
call(n1, n2, a)
...
sub(a)
real a(n1, n2)
...

ただし、FORTRAN77はサイズのチェックをしないから、サブルーチン側で擬寸法仮配列として
宣言しておけば、添え字の最後の大きさを明示的に渡さなくてもいい。
real(5,4)
...
call(5, a)
...
sub(n1, a)
real a(n1, *)
...

さらにいえば、FORTRAN６６時代に書かれている既存のプログラムの場合、最後の次元を
擬寸法で与えることも無く、任意の整数値を入れてあることがある。（１とか１０のことが多い。）

>502

解決しました。どうもありがとうございました。

プログラム全体を通してのパラメータの定義をする方法というのは
ありますか？
サブルーチン、関数ごとに
parameter ( a=10, b=20 )
という風に何回も書いているのですが、バグの元になっています。

use か include を使う。

>>502
遅くなりましたが、ありがとうございます。
調べてみます･･･。

　999 CONTINUE

条件分岐で
例えば
変数 i,j,k,l
が1から10までそれぞれとりうるとき
その各々について
分岐させたいんだけど
なにかいい命令文ない?

それぞれで内容がおおきく異なったり、おなじになったりするから

if,select case,do,where

どれも適当なものにならない
(書く量が莫大になる)
なにかいい命令文ない?

ちなみにfortran90
です

>>508
いまいち状況がよく分からないけど、
num=(i-1)*1000+(j-1)*100+(k-1)*10+(l-1)
ってな感じでひとつの変数にまとめてから、
select caseを使って分岐するってのはだめ？

if文だらけの状態でここにサンプル書いてみれば？
効率的にできるアドバイスもらえるんじゃね？

Microsoft Visual .NET でFortranコンパイラを使っているのですが，
実行すると
　ifort: Command line warning: ignoring option '-extfor'; argument required

　ifort: Command line warning: unrecognized source type
　　'C:\Documents and Settings\aoyama\My Documents\Visual Studio
　　Projects\yabai\shallowgtemp3d8_new'; object file assumed
　ifort: Command line warning: no action performed for specified object file(s)

というエラーメッセージが表示されてしまいます。
これはどういう意味なのでしょうか？
またその解決法があれば教えていただきたいです。
お願いします。

>>512
ソースのファイル名が「-extfor」になってない？
ファイル名と拡張子（固定形式なら .for か .f、自由形式なら .f90）を確認してみて下さい。

>>513
ファイル名を.f90にしたところ上のようなエラーメッセージがなくなりました。
ありがとうございます！

10個の数字を大きい順に並べる方法を教えてくれい

>>515
「ソートアルゴリズム」とGoogle先生に聞けば即答してくれるよ。
しかも何通りも出てくる。

f2cを使ってfortranのソースをｃ言語のソースに変換したいんだが、redumeが英語で
使い方が良く分からん！ｆ2ｃってどうやって使うの？

>>517
f2c -A -c hoge.f

>>518
？　よくわかんないッス

>>519
せめて自分の環境と何をどこまでやったのかくらい書いてくれや。
>>517だけじゃ回答する方もよくわかんないッス

プログラム終了時にsendmail等を使ってメールを送信することは出来ますか？
(fortran77)

シェルスクリプト組めば何でもできるでしょ。

>>522
具体的にソースにどのように書けばいいのかを教えて欲しいのですが。

#! /bin/sh

（fortranの実行ファイル名）
cat finish.txt | mail -s 'pugaaa' someone@somewhere


とでも書いておけ。つか、Unix勉強汁

>>524
出来ました！
どうもありがとうございます！

intel visual fortran 9.0と　visual c++ .net 2003の連携って
どうやったらできるのでしょうか？
いきなりの質問ですいません
教えてくださいm(__)m

IVFに連携用のインストーラあるから、それ入れると自動的にできるよ。

Fortranがタダで手に入りますよ

サン、Solaris 10対応開発環境最新版「Sun Studio 11」の無償提供を開始
http://japan.zdnet.com/news/devsys/story/0,2000052522,20093324,00.htm?ref=rss

＞Sun Studio 11は、C/C++およびFortranを使用するアプリケーション開発ツール。
＞最新のマルチコアUltraSPARC、x64およびx86プラットフォームを含むサンの
＞最新ハードウェア上で稼働する32ビットおよび64ビットの高性能アプリケーションを開発できる。

×　Fortranがタダで手に入りますよ
○　商用Fortranがタダで手に入りますよ

>>528-529
Solaris版ね。
Linux版の方にはコンパイラは含まれていない。

電気工学の研究をしている者です。
非常に基本的な質問で恐縮なのですが、
誰かご存知の方がいましたら教えてください。

ある数値データファイルおいて、
数値データファイルの１行目から２０行目まで、
そのデータの測定情報が記載されていて（つまり文字データ）、
２２行目から数値データの羅列が始まっている数値データファイルがあります。

この数値データファイルを読み込ませて処理をする場合、
２１行目までを無視して、２２行目から読み込むようにする方法はあるでしょうか？

自分でも色々やってみましたが、結局同じファイルをコピーして、
文字データ部分を削除したファイルを作成し、使用しています。
これは面倒だし、ファイル数も多くなるので、
直接元のデータファイルを処理できれば便利なのです。

>>531
行頭から、ASCII文字が21文字詰まっているものとして読み込み＆無視

>>532
ごめん、行じゃなくて列だと思ってた。

＞ASCII文字がつまっているものとして処理

すいません。具体的にどうやればよいのでしょうか。
これを応用すれば、1行目から21行目まで、
IF文で「数値データでないとき、次の行にうつる」
といったことができそうですね。

>>531
WRITE文で後ろに変数を書かなければ
自動的に1行スキップしてくれる。
> OPEN(10,FILE='hoge.dat')
> DO i=1,21
> WRITE(10,*)
> END DO
> (…このあとデータ読み込み)
> CLOSE(10)

あと、>>532 が言っているように行頭の20文字を無視したい場合は
CHARACTER型の変数を用意して
> CHARACTER*20 chr
> WRITE(10,*) chr, (データを格納する変数)
って感じに書けばOK。

>>535
writeでやると、ファイルに空白行書き込んじゃうことにならない？

>>535
>> WRITE(10,*) chr, (データを格納する変数)
WRITE→READの間違いだとして
書式指定なしでやると途中でスペースやカンマが出てきた場合
そこで文字列の読み込みを止めて変な値を読んじゃうよ。

READ(10, '(A20,3G8.0)') chr, x, y, z
または
READ(10, '(20X,3G8.0)') x, y, z

間違いだらけだし…。すまん。

>>536
WRITE→READの間違いです。

>>537
それは知らなかった。指摘ありがと。

みなさん。ありがとうございました。
みなさんの意見を参考に、以下のようなプラグラムをつくり、解決できました。


OPEN(10.FILE='file.csv',STATUS='OLD')
DO 100 K=1,LEN+21,1
IF(K.LE.21) THEN
READ(10,*)
ELSE
READ(10,*) T,C
HAIRETSU(K-21)=C
END IF
100 CONTINUE
　　　 …
CLOSE(10)

LENというのは、数値データの個数です。
読み込むデータファイル「file.csv」は、
一列目が空白で、２列目に数値データファイルが入ってるので、
一列目にTというダミーの変数を用い、２列目のCの値を配列に格納するプログラムです。

ここでまた質問なのですが、このプログラムだと、1列目の空白をTという変数に読み込むという
無駄な操作をしていることになり、効率が悪いと思われます。
この操作をなくして、直接２列目の数値データだけを読み込ませる方法というのはあるのでしょうか？

>>539
データが馬鹿でかくて読み込みにものすごい時間がかかってるのでないなら気にすんな。
データが馬鹿でかくて読み込みにものすごい時間がかかってるとしても、
そのプログラムを何回も実行しないなら気にするな。
そのプログラムを何回も実行するなら素直にデータ変換しる。

>>539
書式指定して読み込め。効率はほとんど変わらないけど。

この場合の書式指定のやり方がよくわからないです。
でもたしかにこの部分は一瞬で終わるので、気にしないことにします。
ありがとうございました。

　100 CONTINUE

GOTO 543

/*******************
ここへは絶対来ないw
********************/

失敗

! ここへは絶対来ないとか
c ここへは絶対来ないと書くべきでは？

STOP 'ERROR'

>>539 単にこれで良いのではないか？
OPEN(10.FILE='file.csv',STATUS='OLD')
DO 100 K=1,LEN+21,1
　　READ(10,*) 　HAIRETSU(K-21)
100 CONTINUE
CLOSE(10)

HALT

マックを買おうと思うんですが、マックで使えるフリーの95コンパイラってありますか？

g95とか gfortranとか

>549

駄目です。１列目を考慮する部分が抜けているからだと思います。
そのプログラムだと、１列目（空白）が配列に読み込まれるから、
配列にはすべて「0 もしくは0.000…0E+00」が読み込まれます。

ところで、２１行目までを無視する部分は無くても
正しい結果は出力するみたいですが、
私の使用している富士通のコンパイラだと実行時に、
英語の文章がずらずらと表示されて不快なので、
（よく読んでないですが、どうも入力装置10や、文番号100に関係するようです）
やはり２１行目までを無視する部分は必要だと思います。

すいません、どなたか教えていただきたいのですが

Fortranでｆ（ｘ）の積分グラフの面積を台形公式により求め、
さらにその誤差を求めるプログラムなのですが・・・
どのようにしたらよいのでしょうか。
わかる人どなたか、お願いします。
m(._.)m

そんな典型的な例題も探せないのか？人生やりなおせ

どなたか教えてください

例えば

INTEGER A,B,S,P

DO 4 A=2,7
DO 4 B=2,7
DO 3 S=-2,2
DO 3 P=-2,2
WRITE(6,*) A,B
WRITE(6,*) S,P

3 CONTINUE
4 CONTINUE
STOP
END

を実行してA,B,S,Pを表示させてみると、Aが6から始まったり、Pが−1から始まったりしてしまいます。
何がまずいのでしょうか？

>>554
ttp://atlas.shinshu-u.ac.jp/class/Kclass/F-2.html

>>556
まずはファイルに出力させて結果を確認しよう。
後、ループの終端を共有させるのはやめれ。

>>549
配列はみだしてるｗ

>>556
多分問題は全然違うところにあるｗ

ところでFORTRAN77で行くなら、暗黙の型宣言は死んでも守ったほうがいい。バグの元。
悪いことは言わないから　A,B,S,P　はそれぞれ　IA,IB,IS,IP　にしておけ。

>>556
宣言してあってもAやBが整数なのは気持ちが悪い。

…問題の本質とは関係ないが。

ところでコンパイラは何？

>>558
REAL*8 HAIRETSU(-20:LEN)
とかしとけば大丈夫！

…文字列読み込んだとたんにruntime errorになりそうだが

>>560
10 READ(10, *, err = 10) HAIRETSU(K-21)

これで乗り切る手もあるがｗ　ホンモノのエラーが出ると無限ループへｗｗｗ

素人ばっかだな

>>562
玄人さまｗｗｗｗ　なんか芸を見せてみろｗ

print *, 'Pugyaaaaa!'

どなたかintelのコンパイラでintelのmklのLAPACKについておしえてください．
自前のテスト用プログラムや，exampleのファイルをコンパイルしようとすると，'Error in opening the Library module file'というエラーが出ます．
対処方法がわかりません．それどころかライブラリのインストールや設定に問題がありそうなのか，プログラムのコードに問題がありそうなのか，どのあたりが怪しいのかも判断できません．
ちなみにライブラリを利用しないコンパイルは問題ないです．どのあたりを確認していけばよさそうでしょうか．

環境：
Linux fedora 2.6.11-1.27_FC3 i686 athlon i386 GNU/Linux

http://www.xlsoft.com/jp/products/intel/perflib/mkl/lin/Getting_Started.htm#Linking_Your_Application_with_Intel_MKL
を参考に、
$ifort gesvd.f90 -L/opt/intel/mkl/8.0/lib/32/ -lmkl_lapack95 -lmkl_lapack -lmkl_ia32 -lguide -lpthread -o ~/f/a.out
とやってみました．

fortcom: Error: gesvd.f90, line 27: Error in opening the Library module file. [MKL95_PRECISION]
USE mkl95_PRECISION, ONLY: WP => SP
----------^
fortcom: Error: gesvd.f90, line 28: Error in opening the Library module file. [MKL95_LAPACK]
USE mkl95_LAPACK, ONLY: GESVD
…(略)…
compilation aborted for gesvd.f90 (code 1)

intel_mklのマニュアルに書いてあった設定は済ませています．といっても環境変数だけです．
$echo $PATH
/opt/intel/idb/9.0/bin:/opt/intel/fc/9.0/bin:/usr/kerberos/bin:/usr/local/bin:/bin:/usr/bin:/usr/X11R6/bin:/home/user/bin:/usr/java/jdk1.5.0_01/bin
$echo $LD_LIBRARY_PATH
/opt/intel/mkl/8.0/lib/32:/opt/intel/fc/9.0/lib
$echo $INCLUDE
/opt/intel/mkl/8.0/include

>>565
F95用のインターフェース定義Moduleがみつかんねーと言っているので、
F77のサンプルを試してみれば、ライブラリの設定に問題があるのかどうか分かるのでは？

と無責任に言ってみるｗ

>>565
Fortran95インターフェイスを使うには、まず自分でインターフェイス部分をコンパイルする必要があるみたいだぞ。

http://www.intel.com/software/products/mkl/docs/mklgs_lnx.htm

貧乏人はマニュアルとにらめっこしてろ

金持ちのオレもマニュアルとにらめっこするぞｗ

古くからある他人の作ったプログラムを改良する作業をしているのですが
common使いまくりでmoduleに変更しようとしています。

使い方そのものはわかりますが、
commonのメモリの使い方とmoduleにした場合でどう違うか？など
内部的な話を知りたいのです。
それらを知ることができるようなサイトご存じであれば教えてください。

>>570
ＣＯＭＭＯＮはメモリー領域の共有割付なので、ＥＱＵＩＡＶＥＮＣＥをかけているのと
同じようになっている。

たちのいいプログラムなら、COMMON文が共通の名前になっていてModule化は楽だが、
変数名がサブルーチン毎に違っていると、それ相応の手当てが必要になる。

また昔のプログラムは、メモリー領域の節約の為にCOMMONを違ったサイズに切って
再利用したりしているので、機械的な修正ができないこともある。

COMMONの割り当て法に関してマニュアルで調べれば必要なことは分かると思う。

>>572
● .● ● ●● ●● . ..●
● ● ● ● ● ● ● ●
● ● ● ● ● ● ● ●
● 　　●●●●●●●● ●
● ● ● ●
● ● ● ●
● ● ● ●
● ● ● ● ● ●
● ● ● ● ● .●
● ● ● ● ● ..● ●
● .● ● ● ●● ● ●
● 　●●●●●●●● ● ● ● ● ●
● ..● ● ● ● ● ● ●
● .● ● ● ● . ●

>>573
● .● ● ●● ●● . ..●
● ● ● ● ● ● ● ●
● ● ● ● ● ● ● ●
● 　　　　●●●●●●●● ●
● ● ● ●
● ● ● ●
● ● ● ●
● ● ● ● ● ●
● ● ● ● ● .●
● ● ● ● ● ..● ●
● .● ● ● ●● ● ●
● 　　　●●●●●●●● ● ● ● ● ●
● ..● ● ● ● ● ● ●
● .● ● ● ● . ●

>>574
● .● ● ●● ●● . ..●
● ● ● ● ● ● ● ●
● ● ● ● ● ● ● ●
● 　　●●●●●●●●●● ●
● ● ● ●
● ● ● ●
● ● ● ●
● ● ● ● ● ●
● ● ● ● ● .●
● ● ● ● ● ..● ●
● .● ● ● ●● ● ●
● 　　●●●●●●●●● ● ● ● ● ●
● ..● ● ● ● ● ● ●
● .● ● ● ● . ●

>>575
● .● ● ●● ●● ●
● ● ● ● ● ● ● ●
● ● ● ● ● ● ● ●
● ●●●●●●●●●● ●
● ● ● ●
● ● ● ●
● ● ● ●
● ● ● ● ● ●
● ● ● ● ● ●
● ● ● ● ● .● ●
● .● ● ● ●● ● ●
● ●●●●●●●●● ● ● ● ● ●
● .● ● ● ● ● ● ●
● .● ● ● ● ●

>>576
● .● ● ●● ●● ●
● ● ● ● ● ● ● ●
● ● ● ● ● ● ● ●
● 　　●●●●●●●●●● ●
● ● ● ●
● ● ● ●
● ● ● ●
● ● ● ● ● ●
● ● ● ● ● ●
● ● ● ● ● .● ●
● .● ● ● ●● ● ●
● 　　●●●●●●●●● ● ● ● ● ●
● .● ● ● ● ● ● ●
● .● ● ● ● ●

プログラミングに関して全くの初心者なのですが、
大学でfortranを使った計算をすることになりました。

Compac Visual Fortran Pro.Edit V6.6
で計算していたプログラミングファイルを、そのまま
absoft Pro Fortran 8
で計算することは可能なのでしょうか？

>>577
基本的には可能。
だが実際にはなんともいえない。

元のプログラムがどれぐらいベンダー依存性を持つかによる。

ベンダー依存性・・・ですか。
すみません、初めて聞く言葉でよくわからないのです。

コンパイルするまでの手順も違うので、戸惑っています。
とりあえず、マニュアルに従ってコンパイルしたら、
local variable mx never referenced.
local variable my never referenced.
というエラーメッセージが出ました。

>>579
それは宣言された変数がプログラムの中で
使われていないことに対する「警告」では？

>>579
ベンダー依存とは、コンパイラー製造メーカーごとの独自仕様のこと。
元のプログラムを書いた人が、標準の範囲内で書いていれば基本的には問題なし。


しかし、最近は浮動小数点もIEEE形式に統一されたからいいが、
ちょっと前まではハードごとメーカーごとに仮数部と指数部のBit数が違うので
収束やオーバーフローの範囲が異なり色々面倒が起きた。

MacとWinだとEndianの違いがある。バイナリーのファイルをやり取りしようとすると
問題になる。でもまあコンパイラーのオプションで対応できるだろう。

秒単位で表された時間ｔを読み込み、何日何時間何分何秒になるかを表示するプログラムって
どうやってやればいいですか？本当初心者です。。。

>>582
考え方が解らないのか、コーディングが解らないのか？
ともかく単なるわり算の問題だぞ。

>>582
１日＝２４時間
１時間＝６０分
１分＝６０秒
です。

>>582
時間や日付を取得する関数はプラットフォームに依存するからなあ…。
コンパイラ何使ってるかくらい書かないと

本当初心者なんでコーディング？プラットフォーム？って感じです・・・
学校の課題で出てるんですが・・・s=int(x/y)とr=amod(x,y)を組み合わせて
考えろって書いてあるんですがどうやってもうまくいかなくて。

書いたところまで貼ってみな

>>582
あー、なんか分かってきた。秒の単位の変数ｔを日にちと時間に分解するやつだろ？
紀元０年からの秒数かな？うるう年とか絡んでるだろな。
amod()使うのはうるう年の検出のためかな？

>>588
ちがうだろｗ
分と秒　時間と分　日と時間にわけるためだろ

Fairfield公式 でググると幸せになれるかもしれない。

>>588,590
ユリウス日とかそんな高級な問題ではないと思います。
単に切り捨ての割り算（整数の割り算）と余りの問題だと思います。

PROGRAM test
IMPLICIT NONE
INTEGER :: ntime, isec, imin, ihour, iday
READ(*, *) ntime
isec = MOD(ntime, 60)
ntime = INT(ntime / 60)
imin = MOD(ntime, 60)
........
STOP
END PROGRAM test

残りは繰り返しだから自分で完成させろｗ

秒が整数で与えられると仮定したので、amodでなくてmodにした。
その辺は自分でハッテンさせてくれ。

ま、宿題は自分でやれよｗ

>>583　〜　>>592

みんな凄く親切だな・・・。

卒論生を指導したことがあるが，
「次の作業を指示してください。」という依頼に対して，
「自分で試行錯誤してみるのも，研究だよ。
　いろいろ試してみて，わからなければ聞いてね」
と門前払いしたことがある。

作業内容が明確にわかる指示に対してしか反応できない学生が
多すぎるのかな・・・。


>>582
教科書や補助プリントが全くないとは思えないから，
参考になりそうなプログラムを調べてみるということは，
思いつかないか？

「教科書例題13（入出力）と例題25（関数）を組み合わせれば，
　できそうだが，組合せ方がよくわからない。」
とか，
「こんなエラーがでて，バグが発見できない。」
という質問なら，可愛げもあるのだが。

>>580,581
ありがとうございます。やはり素人には難しいです・・。
重ねて質問してもいいでしょうか。

mx と my が使われていないということは、
それらをプログラムファイル上で書き換えないといけないということでしょうか。
mx と my はなぜ使われないのでしょう。

プログラムファイルは先生が書いたものなのですが、
ここに来る以前にエラー文を見せた際『分からない』と言われました。
先生がfortranを学んだのは20年近く前になると思います、
『標準の範囲内』を超えている可能性もあるのでしょうか。
Winを使っています。

自然数 N が与えられたとき、
1 から N までの数字を N 個並べる組み合わせをすべて
列挙するプログラムは FORTRAN で書くとどうなりますか？

例えば N = 3 のとき
1 2 3
1 3 2
2 1 3
2 3 1
3 1 2
3 2 1
となるようにです。
よろしくおながいします。

>>595の質問は Pascal スレでもみたんだけど、最近流行ってるの？

>>594
> mx と my はなぜ使われないのでしょう？
　変数名のタイプミスの可能性もある。
または，（動作確認用などで）とりあえず変数宣言したが，
動くようになってからは，コメントアウトしてしまった。
も考えられる。
「警告」は，エラーではなく，「動くけど，間違いてないか？」
という，（時にはお節介だが）ありがたい注意だ。

impicit none
をプログラムの最初に書くと，吉かもしれない。
　# 未宣言エラーが大量にでるかもしれないけど。

>>597
正　implicit none
誤　impicit none
俺がタイプミスをしている。

>>597
Compac Visual Fortran v6.6では計算ができるファイルなのですが、
タイプミス等と考えられるのでしょうか。
エラーではないのですね。

本当に基本的なことばかりですみません。

implicit none
書いて試してみます！未宣言エラーって予想もできないので怖いですけど。

>>599
２０年前に書いたとするとFORTRAN７７だから、暗黙の型宣言を使っているはず。
IMPLICIT NONEはかえって大混乱を及ぼしかねないｗ

昔のプログラムの書き方では変数を使わないことは結構あるので気にするな。
たぶんmx,myはサブルーチンの仮引数だろう。

>>600
大正解です。

implicit none を書いたら、
シンタックスエラーチェックをした際にエラーが100個出てきました。
軽くパニックです。

なので、mxとmyを何とかしようと、意味も重要性も分からないまま
思い切ってmx→xm，my→ymに全部書き換えてみたところ
同じ箇所で警告がでました。

real mx(mm,2),my(mm,2),ox(mo,2),oy(mo,2),oz(mo,2),refc(4)
^
このmxとmyが駄目みたいです。

でも、気にしないで良い程度のことだと分かり
それだけでとても安心しました。ありがとうございます。

>>595　余りエレガントでないｗ　
判定条件を簡素化出来ればいいのだが、２乗和を比較するとか。
でも一意だと証明できないｗ

PROGRAM vip
IMPLICIT NONE
INTEGER, PARAMETER :: n = 3
INTEGER :: i, k(n)
LOGICAL :: qpr
k = 1
DO
IF ( ALL(k == n) ) EXIT
DO i = 1, n
IF ( k(i) > n ) THEN
k(i) = 1
k(i + 1) = k(i + 1) + 1
END IF
END DO
!
qpr = .TRUE.
DO i = 1, n
IF ( ANY(k(i) == k(i + 1:)) ) qpr = .FALSE.
END DO
IF (qpr) PRINT *, k(n:1:-1)
k(1) = k(1) + 1
END DO
STOP
END PROGRAM vip

すいませんフォートランとは直接関係ないかもしれませんがアイディアをいただければ。。

三角形が無秩序に並んでる（各点の座標は分かってます）なかに任意の点（これも座標は分かってる）
をおいたらその点はどの三角形に入っているかということなんですが何かアイディアはないでしょうか？

三角形の最小x最大x最小ｙ最大ｙ出してその範囲内にあるやつってやってたんですが途中でダメなことに
気づきまして、ちょっとしたヒントないでしょうか。

>>603
丸め誤差や例外判定などをまったく無視すれば，
ベクトルの三角形内部判定を使えばいいと思う。

実際にコーディングしようと思うと，
境界線問題や，0割など例外処理が頻発するので，

FORTRAN 計算幾何プログラミング / 杉原厚吉
岩波コンピュータサイエンス
http://www.amazon.co.jp/exec/obidos/ASIN/4000077082/
などに目を通した方が吉かもしれない。

>>603
各頂点から、極座標表示で残りの2頂点と判定すべき点をあらわして、
判定すべき点の位相角がすべての頂点から見て、２頂点の角の間に
入ってればいいんでね？

>>603

検証してないので間違ってたらすまんが

既知の三角形ＡＢＣと点Ｐがあるとして
∠ＡＢＣ＜∠ＡＰＣ
∠ＡＣＢ＜∠ＡＰＢ
∠ＢＡＣ＜∠ＢＰＣ
が全部成り立てばいいのでは？

>>604-606さん
ありがとうございます。すぐに３つも候補をいただけるとは感謝です。
全然思いつきませんでした。
みなさんの考えを参考にさせて頂きましたが結構簡単で計算負荷も軽そうな>>605
さんのやり方で作ってみることにします。できなかったら604さんや606さんのもやってみます。
微妙にスレ違いでしたがみなさんご親切にありがとうございました。

すみません。
エロ動画をダウンロードするための連番リストを作成しようとしています

do i=1,10
write(6,*) "http://www4.xtokyo.com/jpmo/j171/jpmo_1",i,".mpg"
end do

これでループを回すとリダイレクトで連番ファイルが作成されるのですが、
jpmo_ 1 .mpg
のように、切れ目で空白がひと文字分はいってしまいます。

jpmo_1.mpg
のように連続的に文字列をつなげるにはどうすればよいのでしょうか？
よろしくお願いします。

>>608
PROGRAM ero_list
IMPLICIT NONE
CHARACTER (LEN=50) :: fmt
INTEGER :: i
DO i = 1, 10
WRITE(fmt, '(a, i1, a)') '(a, i', INT( LOG10(REAL(i)) ) + 1, ', a)'
WRITE(6, fmt) "http://www4.xtokyo.com/jpmo/j171/jpmo_",i,".mpg"
END DO
STOP
END PROGRAM ero_list

出力結果
http://www4.xtokyo.com/jpmo/j171/jpmo_1.mpg
http://www4.xtokyo.com/jpmo/j171/jpmo_2.mpg
http://www4.xtokyo.com/jpmo/j171/jpmo_3.mpg
http://www4.xtokyo.com/jpmo/j171/jpmo_4.mpg
http://www4.xtokyo.com/jpmo/j171/jpmo_5.mpg
http://www4.xtokyo.com/jpmo/j171/jpmo_6.mpg
http://www4.xtokyo.com/jpmo/j171/jpmo_7.mpg
http://www4.xtokyo.com/jpmo/j171/jpmo_8.mpg
http://www4.xtokyo.com/jpmo/j171/jpmo_9.mpg
http://www4.xtokyo.com/jpmo/j171/jpmo_10.mpg

解説
　スペースが入ってしまうのは、フォーマットを指定せずに自由形式を使っているから。
フォーマットは、隙間が開かないように文字列、数字の長さ、文字列を指定すればいい。
文字列は長さを指定しなければ自動的に与えた文字列長になる。すなわち　a　。
問題は数字で、今回の場合１桁の場合と２桁の場合が混じっているのであらかじめ
決まらない。したがって動的にFormat文を生成しなければならない。
（IF　文で場合分けしてもいいが）
生成したいFormatは、数字が１桁のとき '(a, i1, a)' ２桁の時　'(a, i2, a)'

動的にFormatを作りたい場合は普通、内部ファイルという物を使って作り上げる。
それが最初のWRITE文。１０の対数をとって１を足しているのは変数iの桁数を求めている。
fmtという文字列の中に、'(a, i1, a)' ないし　'(a, i2, a)'　を書き出している。

Format文はWRITE文の中に文字列として直接書いてもいいし、文字変数でもいい。
FORTRANの教科書などでは行番号を用いる形式を良く使うが、情報を局所に展開
させるほうが好ましいいので避けるべし。

>>609
あ、すごい。早速ためしてみます。
ありがとうございましたm(_ _)m

WRITE(fmt, '(a, i1, a)') '(a, i', INT( LOG10(REAL(i)) ) + 1, ', a)'
なんてものがあるんですね…

今まで出力形式に付いてはあまり拘らずに、データをログファイルに記述できれば
それでことたりてたので、こういうのは初めて見ました…

出力形式についても>>610の解説を参考に勉強してみようと思います。

内部ファイルについてはレスの>>89-106あたりでも触れられている。

FORTRAN６６時代は、ENCODE、DECODEというのを使っていたので、
FORTRAN７７以降の機能だと思う。

いろいろいじくって遊んでみました。

WRITE(fmt, '(f7.5)') REAL( LOG10(REAL(i_count)) ) + 1
WRITE(6,*) fmt

こんな感じで実数形式でまわして書き出してみると
1.00000
1.30103
1.47712
1.60206
1.69897
1.77815
1.84510
1.90309
1.95424
2.00000

ってなりました。
対数と頭のINTをつけて形式を変換することによる
小数点以下切り捨てによってi1とかi2とか区別していたんですね。
勉強になりました。

あと、ifのベンチもとってみました。
Iを1000000でまわして
if (i < 10) then

WRITE(fmt, '(a, i1, a)') '(a, i', INT( LOG10(REAL(i)) ) + 1, ', a)'
WRITE(6, fmt) "http://www4.xtokyo.com/jpmo/j171/jpmo_00",i,".mpg"

else

WRITE(fmt, '(a, i1, a)') '(a, i', INT( LOG10(REAL(i)) ) + 1, ', a)'
WRITE(6, fmt) "http://www4.xtokyo.com/jpmo/j171/jpmo_",i,".mpg"

end if

のときと

WRITE(fmt, '(a, i1, a)') '(a, i', INT( LOG10(REAL(i)) ) + 1, ', a)'
if (i < 10) then

WRITE(6, fmt) "http://www4.xtokyo.com/jpmo/j171/jpmo_00",i,".mpg"

else

WRITE(6, fmt) "http://www4.xtokyo.com/jpmo/j171/jpmo_",i,".mpg"

end if

で比較すると
前者4.169u 0.225s 0:04.52 96.9% 0+0k 0+0io 0pf+0w
後者6.366u 0.239s 0:07.01 94.0% 0+0k 0+0io 0pf+0w

となりました。ifの前で読み込むよりも、それぞれのif文のなかで読み込ませる方が速いようです。

DO i = 1, 10
WRITE(6, '(a, i3.3, a)') "http://www4.xtokyo.com/jpmo/j171/jpmo_",i,".mpg"
END DO

エロリストの数字を007の形にしたいなら、FORMATの I3.3 をつかえばいいお。
これなら ｉ の桁数で場合分けしたりしなくてもいいお。

http://www4.xtokyo.com/jpmo/j171/jpmo_008.mpg
http://www4.xtokyo.com/jpmo/j171/jpmo_009.mpg
http://www4.xtokyo.com/jpmo/j171/jpmo_010.mpg

でもエロ画像はＦｌａｓｈＧｅｔやＩｒｉａ（の派生）で落としたほうがいいお。

iを001〜075とかでまわしても、頭の0が無視されてたのでifでやるべきなのかなと思ったのですが
ifなくてもできるんですね…

上と同じような条件のベンチは
3.176u 0.333s 0:03.62 96.6% 0+0k 0+0io 0pf+0w
if文はできるだけ使わない方が良いみたいですね。

エロ動画も落とせてfortran力もついて一石二鳥。
すごいおいしいです。

>>593
大学教育としては問題があるが、より広い意味での問題解決力としては、
頭を下げて人に教えてもらって解決するというのも一つの重要なテクニック。

分からないなりに理解しようと努力しているのであれば
手助けしてやろうという気も起こるだろう。
何も考えずに丸投げするような奴は論外。

erolist インリン　受けたｗ
今日はおもろいねた無いのか？

この流れで抽象的な質問をするのはもうしわけないのですが、
Salford FTN95で、よそからもってきた固定形式(.for)のプログラムを
コンパイルしたところ、

Attempt to call a routine with argument number twenty containing
too few array elements

というエラーが出て困っています。
このエラーが出る場合、どういう問題が起きているのでしょうか。
問題となっているサブルーチンは、引数を23個もっているので、
その内の20個が問題あり？どういったことを試せばよいかも
わからず途方にくれています。どうかアドバイスをお願いします。

×コンパイルしたところ、
○コンパイルして実行したところ、

でした。

>>621
２０番目の引数の配列の大きさが小さすぎると言っている。

ソースが古いプログラムなら、意図的にやっている可能性もある。
ＦＯＲＴＲＡＮは参照渡しだから、配列を引数としても先頭要素のアドレスを
渡しているだけなので、昔風のプログラミングではサブルーチン側の配列の引数を
Ａ（１０）とかと適当な数で宣言して、ルーチン内でははみ出して使っていた。

>>621
２０番目の引数の配列の大きさが小さすぎると言っている。

今の情報ではこのくらいしかいえない。

ちなみにソースが古いプログラムなら、意図的にやっている可能性もある。
ＦＯＲＴＲＡＮは参照渡しだから、配列を引数としても先頭要素のアドレスを
渡しているだけなので、昔風のプログラミングではサブルーチン側の配列の引数を
Ａ（１０）とかと適当な数で宣言して、ルーチン内でははみ出して使っていた。



もっと詳しく知りたければ、ＣＡＬＬしている文と、ＣＡＬＬされているサブルーチンの
頭から引数の宣言文までを見せてくれないとなんともいえない。

>>624さま
レス本当にありがとうございます。試しに20番目の配列の大きさを小さめに
して実行してみたら、その部分はエラー無く動きました。
しかし、同じルーチンの他のいくつかの引数についても同様のエラーが
出たので、次々に配列を小さくしていったのですが、今度はそれらの
引数をwriteするところで、format/data mismatchが出てしまいました。

これらの次数とformatの指定を丹念に調節する以外になにか方法は
ありますでしょうか。古いルーチンに合うコンパイラが必要？
以下は、call文とサブルーチンの頭です。サブルーチンは、
インプットファイルからデータを読む働きをするものです。

CALL READT(TVA,RAR1,RAR2,ASA,INVECT,TVB,IDEG,NL,V,VL,JJS,TST,
& TSTS,THETA,FI,LMMAX,NPSI,ES,PHSS,PHSS2,L1,IPR,NEL)

SUBROUTINE READT(TVA,RAR1,RAR2,ASA,INVECT,TVB,IDEG,NL,V,VL,JJS,
& TST,TSTS,THETA,FI,LMMAX,NPSI,ES,PHSS,PHSS2,L1,IPR,NEL)
COMPLEX VL(NL,2)
DIMENSION ARA1(2),ARA2(2),RAR1(2),RAR2(2),ASA(10,3),ARB1(2)
DIMENSION ARB2(2),RBR1(2),RBR2(2)
DIMENSION CPARB1(2),CPARB2(2),PHSS2(NPSI,80)
DIMENSION V(NL,2),JJS(NL,IDEG),ES(NPSI),PHSS(NPSI,80)
DIMENSION IT1(5),THDB(5),AM1(5),FPER1(5),FPAR1(5),DR01(5)
DIMENSION DRPER1(5),DRPAR1(5)

>>625
ＦＴＮ９５コンパイラの仕様と、ソースプログラムの信用度にもよるのだが、
動くんだったら下手にソースをいじらず動かしたほうがいい。

元のソースが正しいと仮定しての話だが・・・・
単にＦＴＮ９５のチェックが厳しすぎるだけで、オプションで実行時チェックを切れば動くはず。


配列の最後の要素の大きさはＦＯＲＴＲＡＮでは気にしないのが伝統だった。
理由は参照呼出しと多次元配列のメモリー上での取り方の性質による。
詳しい説明はメンドイので誰かに聞いてくれｗ

オプションで実行時のチェックが切れるのですか。
恥ずかしながらそんなことも知りませんでした。
今日は用事があるので明日また調べてやってみます。
本当に有り難う御座いました。

>>627
どういう環境で実行しているのか知らないが、統合環境を使っていていると
デフォールトではＤｅｂｕｇモードでＢｕｉｌｄされるとおもう。その場合はチェックが
効いているだろう。

コマンドラインでコンパイルしているなら、自分でオプションをつけない限りは
チェックは入らない。したがって、この場合はコンパイラがきつい仕様になっていると
いうことで諦めて、他のコンパイラで試してみるべし。ｆ７７なら昔のＷＡＴＣＯＭとか
ｇ７７とか色々ただの奴がある。ＦＴＮ７７も個人使用はＦｒｅｅになっているようだ。


・記憶領域上の多次元配列に関して。〜特に２次元の場合〜
ＦＯＲＴＲＡＮでの多次元配列は記憶領域上では連続して１次元に展開されている。
このとき、第１列、第２列、・・・、第Ｎ列、・・・、最終列というように列ごとに展開される。
つまり配列Ａ（１０，８０）の場合、Ａ（１，１）の次はＡ（２，１）であり、Ａ（１０，１）の両隣は
Ａ（９，１）−Ａ（１０，１）−Ａ（１，２）とＡ（９，１）とＡ（１，２）である。

サブルーチンに配列を含む変数を渡すとき、ＦＯＲＴＲＡＮは参照渡し（アドレス渡し）
なので、実際は先頭要素のＡ（１，１）の記憶領域上のアドレスを渡しているだけである。
それがどこまで続いているかは、実行時にはＦＯＲＴＲＡＮは気にしていない。したがって、
はみ出しさえしなければ、プログラマーの責任でサブルーチン側の宣言サイズに必ずしも
囚われる必要は無い。（が、普通はそれなりに敬意を払う。）

ＦＯＲＴＲＡＮの数値計算ライブラリなどを見れば分かるが、行列の列の大きさは
引数として要求されるが、行の大きさは要求されないのはこの事情による。

vine linux 3.2 でgcc-3.3.2 を使ってます．
fortranプログラムのコンパイル時に
-r8 のオプションを使いたいのですが，
そんなオプションはないと怒られます．
何かかわりのオプションはあるでしょうか?

>>629　マニュアル嫁
GNUとかじゃなく、まともなコンパイラ使え

>>630
すんまそん
お金無くて…

Intel fortranのLinux版タダじゃん。

>632
サンクス

ちょっとためしたみました
エラー: 依存性の欠如:
libsvml.so は matfor_f-3.2-0 に必要とされています
libscml.soなるものが必要だそうだ．
探してみます．

>>629　ｇ７７にはなさげ。

Increasing Precision/Range

Some compilers, such as f2c, have an option (-r8 or similar) that provides
automatic treatment of REAL entities such that they have twice the storage
size, （中略）

g77 does not yet offer any option that performs these kinds of transformations.
Part of the problem is the lack of detailed specifications regarding exactly how
these options affect the interpretation of constants, intrinsics, and so on.

ttp://programming.ccp14.ac.uk/g77-fortran/~michael/os2/g77/g77.htm#hd_670

ｆ２ｃ＋ｇｃｃという昔風にやるか、IMPLICIT REAL*8 (A-H, O-Z)を各サブルーチンごとに
いれるか・・・でも定数がｒｅalのままなだと解決にならんかも。

■
ｇ９５はいつまでたっても完成しねな。もうFORTRAN２００３の規格も決まったんだが。

GNUのfortranもどっかの商用メーカーが倒産してそのフロントエンドでも
貰わないとまともに動かないんでね？

Mac,OpenSource系は信者様の調子のいい話が多くて判断つかねぇ。
連中はすぐ頭に血が上るからうっかりした冗談もいえねえしｷﾓｽ。

>634
親切にどうも
ｆ２ｃ＋ｇｃｃでなんとかなりそうです。

FORTRAN素人で困っています。

intel visual fortran を初めて使ってみているのですが、
libm.lib を利用しようとして座礁しています。
use libm
とすると、ファイルが見つからない！と怒られます。
多分、libm.modを探しに行っているのかなと思うのですが･･･。
もちろん libm.modはありません。
さて、どうしたらいいのでしょう？

助けてください。

>>636
いまいち状況が分からんが、多分勘違いをしていると思う。
本来の付属ライブラリならその行は要らないはず。

もしくは・・・

USE　で用いる .MODファイルというのは、C言語でいう .H ファイルのようなもので
サブルーチン等のインターフェース情報が入っている。
それはソースから自動生成されるのもなので、ソース（インターフェースのソース）があって
自分でコンパイルするか、はじめから提供されているかしなければどうしようもない。

それでも、やる必要があるなら　USEのかわりに　各個のサブルーチンを　external 文で
並べておけば、コンパイル時のインターフェースのチェックは無しで、LINKERがLINKしてくれる。
まぁ７７時代のようなもの。

>637
早速のアドバイス、ありがとうございます。
明日、再トライしてみて、結果をご報告します。

以前質問した603です。
いろいろやってみたんですが604さんの三角形内部判定というやつでできました。
正直こういう数学屋っぽい証明がいるような物は実用性ないだろうと思ってたんですがやってみて
あまりのシンプルさにびっくりしました。教えていただいてありがとうございました。

90が77用のコードを実行できたように、95も90や77で作ったコードを実行できますか？

>>640
９５は９０を含んでいるので出来る。
わずかだが７７の一部の機能は廃止になっている。
なので規格に完全に忠実だと動かない物も出るはず。

しかしながら大抵のコンパイラーメーカーは過去との互換性の為に規格から削除された機能を
サポートし続けているのでまず大丈夫だろうｗ　ウォーニングぐらいは出るかもしれんが。

>>640
Metcalf, Reid の　Fortran Explained 90/95 に、規格から削られた部分、
次に削る予定になっている部分などが記述してある。

http://www.amazon.co.jp/exec/obidos/ASIN/0198526938/249-3009705-2437946

http://www.amazon.co.jp/exec/obidos/ASIN/0198505582/249-3009705-2437946

>>642
この本には，「Other Features　使わない方がよい機能」という，
お節介な（でも有益だと思う）章がある。
「bit 別冊　詳解 Fortran 90」の改訂復刊を希望する。

Solaris(UNIX)でFortran(コンパイラ f77:SC3.0.1)プログラムを作成しています。
このFortranからC言語(コンパイラ cc:SC3.0.1)の関数を呼び出しています。

これらをWindows XPに移植しようと思っているのですが，UNIX/PCの
Fortranの互換性や違いを教えていただけないでしょうか。

また，移植が簡単に行えるFortran/C言語のコンパイラの組み合わせ
があれば教えていただけないでしょうか。移行ツールでも良い。

例えば，インテル Visual Fortran + Microsoft VC++ .NET とか。

以上，よろしくお願いします。

全然問題なし。マニュアルに書いてあるし。

>>644
Lahey-FujitsuやABSOFTのFORTRANはprofessional versionを買えば
Cコンパイラが付いてくる。

Compaq(DEC)のVisualFortranは、MS VisualC++V６とVisualStudio環境で統合される。
日本語版と英語版が入り混じったりするが・・・・

IntelVisualFortranは、IntelC++とVisualStudio2003/.NETで連携するが、
MSのC++との関係はワカンネｗ

g77はgccと組めたはず。

現段階ではこの手の機能はベンダー依存だが、FORTRAN２００３ではCとの相互呼び出し
機能が規格として規定されているので気長に待つという手もあるｗ

FortranにはCで言うところの終端文字というのが無いから、文字列渡すときだけめんどくさいな

FORTRAN2003のオブジェクト指向用の拡張がうんこ過ぎる。
あれを考えた奴は出てきて、どういう使用状況を考えて用意したのか説明しなさい。

ISOの偉い人達が集まって
「どーせFortranでOOPなんてロクに使われないんだから、適当でいーんじゃねーの？（鼻ﾎｼﾞﾎｼﾞ」
って決めますた。

いろいろありがとうございます。

647さんの情報のようにUNIX/PCのFortranの互換性や違いが書いてある
サイトはあるのでしょうか。

>>650
FORTRAN Explained 2003にFORTRANとCでの違いのまとめがある。
データ型の互換性や文字列の終端記号、構造体、call by value, call by reference
などの問題。

それ以外はコンパイラーのベンダー毎に独自の解決をしているので、個別の問題点を
まとめていることはあっても、一般的なまとめは無いのではないかな。マニュアルの相応の
ページを読むしかないと思う。

大文字・小文字を区別するCとしないFORTRANでは、LINK時のサブルーチンの内部
表現名をどうするかなど、細かい落とし穴がいっぱいある。

>>639
役に立てたようで，とてもうれしい。
ただ，誤差対策をちゃんとしているかな？
行列式の0割防止，隣接三角形との所属判定など，
面倒な例外処理が必要だと思うのだが。

>>652
ご指摘ありがとうございます。
調べた感じ正直そこが一番大事そうでしたがとばしてしまいました。。
ただ僕がやろうとしていることの性質では一応線分上にはほぼ来ないはずと踏んでいたし、
使う目的がその点の三角形の各点の変数からその点の参照値を出すということだったので、
近いところが分かればある程度いいやと考えていたので
大丈夫かなと思っていました。
実際に全部でやるべき１１２５点をやるべき三角形全部に判定させてループさせましたが
線上判定かぶりの点はなくどこかの三角形１つに入ってくれました（ちゃんと入るべきところに入ってるかはやってませんので
意味ないかも。。）
線上でかぶってたら
ただ言われてみれば判定する点が三角形の各点に近いのでそのあたりがとても不安です。
（ちょっと勉強不足なのですがやっぱり点に近いと誤差どうこうでおかしくなるのでしょうか？）

なんか書いてるうちに不安になってきました。。（笑）教えてくださった参考書をちょっと明日でも借りてきます。

>>653
簡潔に説明できないが，
EPSL（微少量），s,t（ベクトルの係数）として，
すべての点が
　EPSL < s < 1-EPSL
　EPSL < t < 1-EPSL
　EPSL < s+t < 1-EPSL
に入っているということかな？
だったら，問題ない。

嫌らしいのは，-EPSL < t < EPSL などの境界上にあり，
隣接三角形との取り合いがある場合。
みかんのへたのような，複数三角形が密集する頂点近傍に，
判定すべき点があると，とても神経質になりそう。

領域判定の理論は、ゲーム関係のサイトに載ってるよ。

>>653
すまないが，例にあげた参考図書
「FORTRAN　計算幾何プログラミング」
ttp://www.amazon.co.jp/exec/obidos/ASIN/4000077082/
には，そのものずばりの解答は載っていない。

無誤差の整数にしてしまう，というとても遠回りな手法をあげている。
「誤差を真剣に検討するなら，このくらいやるべきかも」
というつもりで参考図書にした。

とても遠回りには違いないが，誤差無しという最強兵器を持っているので，
誤差検討に神経を使うなら，意外に近道だと思う。

>>655 氏のあげている，ゲームグラフィックスは非常に参考になるはず。

「ゲームプログラミングのための3Dグラフィックス数学」
Eric Lengyel，狩野智英
http://www.amazon.co.jp/exec/obidos/ASIN/4939007375/
を資料費で購入してもらおうと思ったら，タイトルで怒られた。

LU分解のプログラムを誰か…
ピボッティングはなしで，LとUに分解するだけでいいんで…

>>657
LAPACKやMKLに入ってるやつ使えば？

>>657
精度比較用に「無しでなければならない」ということか？
ピボッティングありのプログラムは，
宮崎大などですぐに見つかるから，
ピボットの意味を考えて，がんばれ。

余談だが，f77 の3重do文をf90 化すると，
配列操作の練習になるかもしれない。

open(10,File='test.txt')
write(10,*) 'Hello,world!'
close(10)
stop
end

これでカレントディレクトリにtest.txtが作成され，中にHello,world!が書き込まれるのは分かるのですが，
プログラムからカレントディレクトリに新しくディレクトリを作成するにはどうすればよいでしょうか？
言語はFortran90です．．

>>660
ＦＯＲＴＲＡＮ９０／９５ではそれは規格に無い。
ベンダー依存の機能なのでマニュアルの当該部分をみるべし。

ＦＯＲＴＲＡＮ２００３ではＰＯＳＩＸ互換のルーチンを標準化するようになっていたと
思うので、共通のやり方ができるようになるだろう。

>>657
すまない。嘘をついた。
宮崎大は，ピボッティングなしだった。
「学外の方には有用ではない」という注釈を忘れていた。

>>657
ＮＵＭＥＲＩＣＡＬ　ＲＥＣＩＰＥＳのＨＰに行けばソースが見られる。
ＬＵ分解はバグ無かったはずｗ
ピボッティングしていたが。

http://www.nr.com/

暗黙の型宣言は使うなとか，common文は使うなとか，
変数名はこう命名しろとかいうような，より良いコードを書くために必須な知識が纏められている書籍
って無いですかね？
C++で言えば，Effective C++に相当するような・・・

>>661
ﾚｽありがとうございます．
マニュアルを参照してみます．

>>664
書籍じゃないけど、こんなのはどうかな？

Fortran 標準コーディングルール
ttp://www.mri-jma.go.jp/Project/mrinpd/coderule.html

必死さから、一瞬卒研やってる学生かと思ったが、
さすがに今の時期にLU分解を組み込んでるようなら違うかな。

>>666
�d
参考にさせて頂きます

>>667
今月から本腰入れるっていうところも多いから強ちはずれでもないかもよ

>>664
FORTRAN77時代のものなら
http://www.star.le.ac.uk/~cgp/prof77.html
これがまとまっていたと思う。

お聞きしたいのですが
サブルーチン文を２回別の変数を持って呼び出すときに、１つ目の変数をAと
したとき二つ目の変数を-Aにしたいのですがこのとき
call ****(-A,B,.....)
としていいのでしょうか？
やっぱりコールする前にA=-Aとしないといけないのでしょうか

>>670
FORTRAN９０の場合、INTENT(IN）属性をつけたものならマイナスをつけてもOK
INTENT(OUT)、INTENT(INOUT)の場合は駄目だとおもう。

それはFORTRANはアドレス渡しで変数を渡しているので、マイナスが付いたときの
変数の戻し方に任意性が生じるから。

FORTRAN７７時代はどうだったか・・・・
あるいはINTENT属性をつけない場合は・・・・
まぁ昔のFORTRANでは何事も気にしないことが多いので、多分そのまま通るだろうｗ
そして結果はコンパイラーに依存するだろうｗ

Intel Fortran Compiler for Mac OS ベータテスト開始
ttp://www.xlsoft.com/jp/products/intel/mac_beta.html

すいません質問です

Windows用 Intel Fortran EM64T対応のを使ってます。
諸事情でメモリ2G制限を越えるプログラムを作らなくてはならなくなったのですが、
これって静的領域に確保することは可能なのでしょうか？
どうも2G以上の領域を確保しようとするとリンカでハネられてしまうのですが…

>>674
単一の配列で２G超えですか？複数で合計２G超えですか？
OSはWindows64ですか？32bitWindowsでターゲットがwin64用ですか？

でも、ここで聞くよりマニュアルを調べた方が早いと思う希ガスｗ

フォートランをマスターするとどんな幸せな事にめぐり会えますか？

>>672
ありがとうございました。つけて走らせてみたらできました。おかげでめんどくさいことせずにすみました。

>>676
FORTRANはいいぞ〜
最古の高級言語であって、コンピューター界のヒエログリフといってもよい。
地層のように積みあがった言語構造は、神聖文字、神官文字、民衆文字、コプト文字
を思わせて美しい。

　　　　 　 '　「　　　　　　´　{ハi′　　 　　 　 　 }　 l
　　　　　 |　 |　　　　　　　　　　　　　　　　　　　 |　 |
　　 　 　 |　 !　　　　　　　　　 　 　 　 　 　 　 　 |　 |
　　　　　 | │　　　　　　　　　　　　　　　　　　　〈 　 !
　　　　　 |　|/ノ二__‐──ｧ　　 ヽﾆニ二二二ヾ } ,'⌒ヽ
　　　　 /⌒!|　 =彳ｏ｡ﾄ￣ヽ　　　　　'´ !o_ｼ｀ヾ　| i/ ヽ !　　　　や　ら　な　い　か　？
　　　　 ! ハ!|　　ー─　'　　i　　!　 　 `' 　 '' " 　 ||ヽ　l |
　　　　| | /ヽ!　　　　　　　　|　　　　　　　　　　　 |ヽ i　!
　　　　ヽ { 　|　　 　 　 　 　 !　　　　　　　　　　　|ﾉ　 /
　　　　　ヽ　 |　　　　　　　 _ 　 ,､　 　 　 　 　 　 !　, ′
　　　　　　＼ !　　 　 　 　 '-ﾞ　‐ ﾞ　　　　　　　　レ'
　　　　　　　 ｀!　　　　　　　　　　　　　 　 　 　 /
　　　　　　　　ヽ　　　　　ﾞ　￣　 ￣ ｀　　　　 / |
　 　 　 　 　 　 |＼　　　 　 ー ─‐ 　 　 　 , ′ !
　　　　 　 　 　 |　 ＼　　　　　　 　 　 　 ／　　 |

ウホッ！ いいスパゲティ！ orz

>>676
彼女ができる・・・（と，いいな）。

よく言われるが，計算ライブラリーなど過去遺産が豊富。

また，1980年以前に隆盛があった学術分野（気象，機械，土木など）や，
電気回路など複素数変数が便利な分野では，
Fortran が未だに主流ではないかな？

話題はずれるが，>> 90 の演説も面白い。

>>675
えっと、結論から言うと動的に確保することで上手くいってしまいました。

で、実は気になるのは速度の問題なんです。
以前動的に確保することで速度が遅くなったんで動的に確保するのを躊躇ってたのですが
今回はそれほど遅くなってませんでした。

恐らく前回はポインタを使ったために
最適化が上手くいかなかったんだと思うのですけど
もしかしてpointer属性とallocatable属性って違うんですかね？

任意の生年月日と今日の日付を入力すれば、
産まれてから今日まで生きてきた日数を出力するプログラムを作れ。
さらに、任意の整数ｎを入力して、産まれてからｎ日目が、
いつかを出力させるプログラムを作れ。
ちなみに、閏年は以下のようにして決められる。

１）基本的に西暦が4で割り切れたら閏年。
２）しかし、西暦が100で割り切れたら閏年ではない。
３）しかし、西暦が400で割り切れたら閏年とする。

レポート作成
予備的考察、プログラム、プログラムの説明、計算結果、考察

フォートラン使用･･・

>>681
ALLOCATABLEの方が最適化が効きやすくていいと思う。
POINTERの場合、実際のメモリーが割り当てられているのか、TARGETを指しているのか
コンパイラが判断するのが難しいが、ALLOCATABLEは実際にメモリーが割り当てられて
いると仮定できるから最適化をかけやすい、と何かに書いてあったが思い出せないｗ

まぁ、話半分情報ということで・・・

>>682
まず正月からある日までの日数を数えるプログラムを組め。
うるう年は無視していいぞ。

これが出来たら、残りを教えてやるｗ　

>>684
すみません　
まったくわかりません・・・わかろうとしてないかも・・・
意味不明なんです･･･
UNIXが、フォートランが･･・　

教えないほうが本人の為になりそう。自分で考えない馬鹿は学校行かなければいいのに

頼む！　金が無いんだよ･･･　留年するための・・・

>>678
>>680
ありがとうございました。
暇をみて手を出してみます。

>>687
なんでそのガッコにいったの？

とりあえず大学へ･･･　の流れですorz
ごめんなさい、卒業したいです；；

こんな問題も解けない奴に学士号を与えて世間に出されてはこちらが迷惑です。
お願いですからもう一年がんばって下さい。

板汚してすみませんでした。

とりあえず自分で考えられるところまで作って↓で聞いてみれば？

くだすれFORTRAN（超初心者用）
http://pc8.2ch.net/test/read.cgi/tech/1138063703/

いきなり閑古鳥が鳴いてるけど…

>>690
マジレスすると、今すぐ大学やめて違う専門に入り直した方がいい
進路の修正は早い方がいいからな

>>690
大学では壁にぶつかった後、立ち止まって進みも退きもせずに
引きこもりになるより２ｃｈで答えを得るほうがよっぽどいい。

カンコー学生服は丸洗いOK！　VIPPERなら丸写しOKｗｗｗｗ

Fortran Source Code Converter ってどうよ？

>>696
それなーに？ｋｗｓｋ

program Calculation of Response
implicit none
integer i,j,n1,n2,neq,hnda,nda,odd,st,p1,p2,p3
real*8 h
parameter(n1=4096,n2=8,neq=2,h=4.d-3)
real*8 x(n2),y(n2),yd(n2),a(n1,n2),b(n1,n2),k1(n2),k2(n2),
&k3(n2),k4(n2),f(n2),m,c,k
character*20 infile,outfile
c
write(6,*)'input >> input file name'
read(5,*) infile
write(6,*)'input >> output file name'
read(5,*) outfile
write(6,*)'input m,c,k'
read(5,*) m,c,k
open(1,file=infile,status='unknown')
read(1,*) nda
read(1,*) p1,p2,p3
do 10 i=1,nda
read(1,*) (a(i,j),j=1,2)
10 continue
close(1)
で始まるプログラムをコンパイルしようとすると下のようなエラーが出てしまいます。
Undefined first referenced
symbol in file
rk4_ /var/tmp//ccMMIE6P.o
ld: fatal: Symbol referencing errors. No output written to shuuron2.exe
collect2: ld returned 1 exit status
infileの名前はちゃんとcharacterで定義しているので何が悪いのか分かりません。
どなたか宜しくお願いいたします。

>>698
ローダーのエラーなので、ファイル名の問題ではないだろう。
＞symbol in file
これが正しいならinfileｆはない。エラーメッセージがこれ全部なのだろうか？
未定義変数でも使っているんだろう。

こまぎれでわかんねが　rk4　か　first　の問題でね？

shuuron2 って・・・ 今頃プログラムつくってて大丈夫なのか？

>>700
なかなかチェックが厳しい。
中間報告（春の学会発表とか）でがんばっているM1 だろう。

>>698
F77形式の上に、暗黙の型宣言を守っていないので、bugは中々とれず
周囲の経験者やプログラム相談員からは相手にされず、苦労するだろう。

要約すれば・・・・・・
『地獄に落ちる！』

　　　　　　　　　　　 ,,.　 -‐―-　､
　　　　　　　　　 ／ _..-‐―- 、　　.ヽ
　　　　　　　　 / ／-‐―-..__..ヽ　_..、ヽ
　　　　　　　　.!三三三´　　　 ｀´　　:　i
　　　　　　　　|.!三三　　 __　　￣ _,. .;. i　暗黙の型宣言を守らない
　　　　　　　　|ii　i!"　　´ ｪｪ｀　l´rｪ､ i..!　　　　　　>>698は地獄へ逝く
　　　　　　　　 !ｉ､ i!:::..... -￣r　　ヽ　　.!
　　　　　　　　　!ｉｉ!.|:::::::　　/ ｀−´ヽ　/
　 _,,..　===＝"/%;;;::::.　　　'∠ﾆﾆ」　 /、　　　　　　　　　　　　_　-‐‐-､
/´　　　　　　 | .%　:::;;;;:　,,,.,,,:. ⌒: ./´　￣＝=-- 　_,,.. -‐‐'"´　_,..-‐ ´
　　　　　　　 |　 %　　::::::;;;;;;;;;:.""%.／　　　　　　ｒ"　(<二￣￣>
　　　　　　　.|　　%｡ﾟ+。　　。+% ／　　　　　　　｀ヽiﾉ<二.￣￣>
　　　　　　　|　　　%+*ﾟ+。*ﾟ%／　　　　　　　　　（　<二￣￣>

>>702
デジャブかな？どっかで同じような光景が…

>702
implicit none をつけたとしても、暗黙の型宣言は守ったほうがいいな。
FORTRANの場合、型の宣言がルーチンの頭部に限定されているので、
いちいちそこを参照するのが面倒くさい。つまりローカルに解決しない。

C言語みたいに捨て変数を｛｝の単位内でローカルに宣言できるなら、
まあいいかもしれないが・・・

本物のプログラマはフォートランを使うと聞いてとんできました。

リスパー＞コボラー＞Fortranコーダー

"Real Programmers Don't Use PASCAL"
http://spike.scu.edu.au/~barry/RealProgrammers1.html

リスプとコボルが上にいるの？
フォートランはブロンズセイントなの？

歴史と伝統ではFORTRANが一番に決まってる。
が、その使用者は向上心がないからな。

ももしきや　古き軒端にしのぶれど　なほあまりある　昔なりけりｗ

エスパー様。お助けください。
Salford FTN77のコンパイラでInvalid character in field
というエラーが出ます。原因はなんでしょうか。

少し状況をいうと、READ文でINTEGERの数値をテキストから読むときに、
入力値を1〜7の数値にするとエラーになり、8以上だと普通に動きます。

>>711
ファイルからのデータ読み込みで、実行時エラーが起きている状況でおｋ？
たぶんデータの並びと、FORMAT文がずれている。

もしくは１〜７の数字が全角になっているかｗｗｗ　または空白が実はTabとか００ｈとかｗｗｗ


そのread文のformatとエラーのでるファイルの行を見せてくれ。

>ファイルからのデータ読み込みで、実行時エラーが起きている状況でおｋ？

そのとおりです。問題となるREAD文は

READ (4,112) NSYM,NSYMS,ASTEP,VSTEP

で、この文のフォーマットは、

112 FORMAT (2I3,2F7.4)

で、読みこみファイルは

1 5 0.02 0.20 NSYM NSYMS STEP VSTEP

です。２番目のNSYMSを７以下の数値にすると
実行時にRUN−TIME ERRORとなります。
はじめ、ウインドウズのワードパッドで入力をつくっていたのですが、
空白がタブの可能性を疑ってメモ帳で保存しなおしたりもしましたが。
うまくいきませんでしたorz。

>>713
２ｃｈの特性で空白が飛ばされているせいかもしれないが、入力のフォ-マットおかしくね？

３桁整数　３桁整数　７桁実数　７桁実数になってないとﾔｳﾞｧｽ

..1..5...0.02...0.02
（ . はスペースのつもり）
のようにしないとﾏｽﾞｽ

>>713
ごめｗｗカラムはだいじょぶそうかなｗｗｗ

READ文の次に　
PRINT *, NSYM,NSYMS,ASTEP,VSTEP
を入れて何が読み込まれているか見て味噌

>READ文の次に　
>PRINT *, NSYM,NSYMS,ASTEP,VSTEP
やってみましたところ、入力した通りの数値が
そのとおりDOS窓に表示されました。
1 5 　2.000000E-02 0.200000 等
しかし、数字が７より小さい場合に限り、そのあとやはり同じエラーがでます。
エラーのダイアログには、

00424588 READPL [+0799]
00401000 （プログラム名） [+1244]

と表示されます。READPLは今問題になっているのとは別のサブルーチンです。

NSYMSはたぶん、場合分けを指定する数字で、READPLもファイルを
読むルーチンなので、そこでINVALID　CHARACTERがあるのかも
と思いつきました。READ文を教わった方法で一つづつ調べてみようと思います。
おかげさまで解決のとっかかりができました。ありがとうございます。

>>716
それはエラーの場所がそのREAD文では無いということだｗ

トレースバックの行番号799行をミロ！

799行はCで始まるコメント行でした。1244行もです。orz.
READPLは9000行くらいに本体があり、呼び出しは517行めでした。

全角で数字打ち込んでるんだろ

>>717

00424588 READPL [+0799]
00401000 （プログラム名） [+1244]

（これが普通のトレースバック情報だとしたらの話だが）
これの意味は１２４４行にあるメインルーチンのCALL文で呼ばれた
サブルーチンREADPLの中の799行にある文でエラーが起きているということだ。

多分そこのREAD文のデータとフォーマットが合っていないのだろう。

>>719　うは、ごめｗ　機械語のアドレスだったかｗ
ただ、コンパイラーのオプションでソース行を表示させることが出来るかもしれない。

たとえばＣｏｍｐａｑのVisualFORTRANではオプションつけないと
トレースバック情報でソース行を吐き出さない。

少し進展したので報告します。
問題のREADPLをcallしている前と後で
PRINT*,'行番号'
をしたら、前は出力されましたが後ろは出ませんでした。
やはりREADPLの中のようです。もう少し調べてみます。ありがとうございます。

>>720
見なおしましたが数字はすべて半角のはずです。
しかし、ワードパッドで打ったので空白などが
へんなのかもと思っています。

解決しました。
NSYMSによって、READPLの中でREADする行数が変化する部分が
あったため、FORMATが一致したりしなかったりするのが原因でした。
ありがとうございました。m(_ _)m

>>725
乙！
中々バグが取れないときは、発生箇所や原因を勘違いして決め付けているときだからなｗ
犬や猫と心が通じるように、エラーメッセージの読み方を学ぶとFORTRANとも心が通じる。
相手（コンパイラー）を理解する気持ちが大切だｗ

エスパー様。また来てしまいました。orz
バグとりの質問ばかりですがどうか許してください。

上で尋ねたのとは別のプログラムを実行するときに、
まったく同じエラー「Invalid character in field」
が出て困っています。コンパイラはFTN77で、エラーダイアログは、
00436b78 READE [+0552]
0042d818 EXPAN [+0273]
00401000 （プログラム名） [+1345]
です。READEはファイル読みこみのサブルーチン。
EXPANは計算のサブルーチンです。
教わったように、PRINT*，'(行番号）'などをDOSに
出力すると、プログラムのごく最初の方にあるOPEN文の前まで
しかPRINTの出力ができません。そのOPEN文は
OPEN (UNIT=2,FILE='search.s')
というものです。このserch.sというファイルは入力値を入れるものではなく、
プログラム自身が生成し、値を書き込むためのものです。
この文の前にはREAD文はないため、原因個所がわかりませんでした。
そこで、OPTIONのことを教えてもらいましたので、マニュアルをみて、
/check　というオプションをつけてみたところ、こんどは、
Integer arithmetic overflow
というエラー表示に変わり、
LOOKUP2 - in file tleed2.for at line 2307 [+0203]
（プログラム名） - in file tleed2.for at line 376 [+13f7]
と表示されました。LOOKUP2は計算サブルーチンです。
2307行は、そのLOOKUP2の前に記述してある別のサブルーチンのEND行
（そのENDのあと、コメントがあり、そのあとLOOKUP2の記述が始まる）
です。行376は、そのLOOKUP2のCALL行で
CALL LOOKUP2(CPVPOS,ILOOK,ILKBD,NLAY2,NST1,NSTEF,NSYM,
& NLTOT2,NSHIFT,NCL,NLMAX,IROT,IMIR,1,FPOS,POSS,POSS2)
というものです。
原因個所を探がすのに、どのようなことを試みるとよいでしょうか。

あ、プログラム名伏せてたのにばれてしもたOTL。。

>>727
難しいな。コンパイラは奇妙な反応を返しているし。
FTN95に慣れていないと、コンパイラの息遣いを感じることも出来ないｗ

風の息づかいを感じていればのガイドライン
http://ex13.2ch.net/test/read.cgi/gline/1135748314/


配列のはみ出しか、引数の不整合が、起きている。もしくは配列のはみ出しで
引数の不整合が起きているｗ　
check で argment check や subscript check も見ていると思うが、option調べて味噌。

inputが間違っているので配列のはみだしが起きているのだと思うがｗ

CALLの直前直後、および呼ばれたサブルーチンの先頭で、大きな配列でない引数の値を
printしてみるくらいしかないかな。このパターンではCALLの直前しか書かないと思うが。

>>729

周りに頼れる人がいないため、レスが大変ありがたいです。
未解決ですが過程を報告します。
あれから、マニュアルでオプションをもう少し調べて
/UNDEFと/BREAKというのをためしてみました。
/UNDEFでは初期値が与えられていない部分を指摘されました。
しかしこの部分は、オプションをつけなければそのままでも初期値０として
動いてくれるようでした。（言い忘れてましたが、サンプルのインプットがあって、
それでは問題なく動く）
/BREAKでは、デバッガが起動してびっくりしました。無手勝流ですが、
今、これを使って変数に入ってる数値を調べているところです。

>CALLの直前直後、および呼ばれたサブルーチンの先頭で、大きな配列でない引数の値を
print

なぜか始めの方のOPEN文より後の部分に、print文を書いても出力が出ないんですorz.

>inputが間違っているので配列のはみだしが起きているのだと思うがｗ

まったくその通りです。ルーチンは動く保証は一応あって、
INPUTで四苦八苦してるんです。

バグとりの質問などは面白くないと思いますが、また進展があれば報告させて
いただきます。

報告。
デバッガで入力変数をまえから順に調べていったところ、
前述のREADEというルーチンの中のあるところから先が値が入ってなかった
ので、そのFORMATをしらべたところ、入力とのズレを修正することができました。
それで、現在、プログラム動差中（たぶん）なのですか、前述の通り、
なぜか、プログラム中のあちこちに書いたはずのPRINT文が、
途中から出力されません。（そのため、本当に動いてるのか不安です）
数十分動かしてプログラムも終了する気配が無いし、なんだか悪い予感がします。
PRINTが出力されないのはいったいなぜなんだろう。。

>>731
反応が返ってこないのは暴走。
最後にprintされたところと、次のprint文の間でどっかへ行っている。
print文の範囲を縮めていって、問題箇所を特定する。


OPEN文に関しては、status=unknownを入れておくといいかも。
あと、暴走を途中で止めた場合ゴミファイルが残っているかもしれないので
消して置くといいかもしれない。

>>732
レスありがとうございます。

>最後にprintされたところと、次のprint文の間でどっかへ行っている。
>print文の範囲を縮めていって、問題箇所を特定する。

ご指摘で、どうやら不具合がわかったっぽいです。PRINTが無効になる
境目だけを別プログラムとして実行してみました（以下）

PRINT*,'263'
OPEN (UNIT=2,FILE='search.s')
OPEN (UNIT=15,FILE='prdir.d')
PRINT*,'265'
OPEN (UNIT=10,FILE='restart.d')
PRINT*,'267'

END

すると、263だけがDOSに出力され、なぜかsearch.sの中に265 267が書かれていましたOTL....
もしかして、UNIT「2」と「*」って重複なんでしょうか。
ともかく、PRINTは使いたいので、UNIT番号を書きかえてやってみます。
今日は午後は予定があるのでしばらく2ちゃんを離れます。
重ね重ねご助言ありがとうございます。

>>733
思いがけない展開だが、たまにあるパターンだｗ

>>もしかして、UNIT「2」と「*」って重複なんでしょうか。
*は普通は６番だと思うのだが、STDOUTとSTDERRで分けていたりするのかもしれない。
１桁や９０番台のUNITはコンパイラーによって予約されていたりする。
これは、コンパイラーの息遣いを感じる以外ないｗｗｗ


安直には２番を２２番とかに書き換えることだろう。

FTN77 User’s Guide (UserGuide.pdf) 102ページ
> by preconnection; the standard input and output are preconnected to units 1 and 2
> respectively (units 5 and 6 are respective alternatives).

息遣いを感じるのも悪くはないが、処理系の取説もよく読んだ方がいい。

だれかこのプログラミングを教えてください。お願いします。
ｘ=0.1 , 0.2 , 0.3 , 0.4 , 0.5
ｙ=1.228 , 1.005 , 0.823 , 0.674 , 0.552
近似式　ｙ=a*exp(bx)
上の関数とデータで最小自乗法で近似した式を求めよ。
この問題のプログラムがわかりません・・・。お願いします。

>>736
こっちで頼んでみろ

くだすれFORTRAN（超初心者用）
ttp://pc8.2ch.net/test/read.cgi/tech/1138063703/l50

Fortranで一度動的に確保した配列のサイズを途中で大きくしたい場合
どのような方法を使うのが一般的？

自分は、ベタに
1. 配列A(N)を確保
2. 同じ大きさの配列B(N)を確保
3. B(1:N) =A(1:N)
4. A解放
5. Aを新規に大きめに確保

くらいしか思いつかないのですが。

できればポインタ使わない方法が嬉しいのですが。

>>736
オマエが考えなきゃいけない問題なんでは？
ここで宿題しても身にならんだろ。

>>736
どうすればaとbを求められると考えているの？プログラムを作るのはそれから

>>738
基本的にはそれしかないだろう。
ただし、ちょっと増えるたびにその操作をすると時間がかかるから、あらかじめ増やす単位を大目に取っておく。
その単位は問題によるけど

そうですか〜、残念。

確保＆解放繰り返すと、断片化が起きそうで怖いんですよね・・。