くだすれFORTRAN（超初心者用）

くだすれFORTRAN（超� ..

830:デフォルトの名無しさん
06/08/08 16:22:18
質問があります。
ある文献のFORTRAN77で書かれたプログラムを、練習がてら
そのまま書き写してbuildすると

libc.lib(crt0.obj) : error LNK2001: unresolved external symbol _main
test.exe : fatal error LNK1120: 1 unresolved externals
Error executing link.exe.
test.exe - 2 error(s), 0 warning(s)

とりあえず、build errorの中を調べてみたんですが、全く解決できません。
コンパイルエラーはゼロです。

どのような解決方法がありますか？教えてください。

831:デフォルトの名無しさん
06/08/08 18:02:10
>>830
program文がないんだよ。

FORTRANのプログラムは program 文で始まり end program 文で終わる。
program文をちゃんと書いたか確かめてみな。

832:デフォルトの名無しさん
06/08/08 19:02:10
program 文は一応あってもなくてもコンパイル可能。普通は書くけど。
end は大事だけど、end がないだけだとそういうエラーにはならない。

このエラーが出てくるのは、
関数とかサブルーチンとかしか定義してないファイルの場合。
まあ、それを指して「program 文がない」と表現してるのかもしれないけど。

833:デフォルトの名無しさん
06/08/09 13:06:31
age

834:デフォルトの名無しさん
06/08/09 22:51:51
>>831-832
うーん。なるほど。
なんとかがんばってみます。ありがとうございます。

835:デフォルトの名無しさん
06/08/09 23:39:19
>>831-832
おまいら、ここは初心者スレなんだからもっと優しく教えてやれよ！

さぁ次はしょんべんだ！ってかんじでさー

836:デフォルトの名無しさん
06/08/10 09:54:51
しょんべんワロタ

837:デフォルトの名無しさん
06/08/10 16:55:46
　　　　　　　＿＿＿＿＿＿　　　　　　 |
　　　 ,,..-‐";;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;` ､　　　i　　　　　　あ
　　／;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;､　ヽ　つ　　あ
　 /;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;|　　〉　ぎ　　　：
　 |;;;;i "'`~　 "`~ `i||i"　''　ﾞ`　"　|;;;;;;|　 /　　は　　　：
　 |;;;;| 　　　　　　ヽ`　　　　　u　　|;;;;;|　| 　　シ
　 .|;;|　,-;;;;;;;;;;"ﾌﾉ　ヾ`;;;;;;;;;;;;;;;ヽ　|;;;;|　| 　　　ョ
　,,ト;| ',,_＝＝-､く　　>ﾞ-＝＝､　 |/ i　|　　　ン
　|i ､|　　　'￣"彡|　　　　　　　　　||　|　| 　　べ
　|'. (| 　　　　彡|　　　　　　　　 |)) |　.| 　　ン
　 ! ､|　　　　　i,"(_ ,, ､,　　　　　　|" i　　| 　　だ
　ヽ_|　　　　　　　 `　　　　　　.|_/　　ノ_
　　　.|ﾞ　　　　､,.－-‐ ､,,　　　　|　　　　　ヽ､＿,,,､＿ノ
　　　.i ﾞ､　　　　' ￣ニ￣　　　　 /|
　　　| 　､　　　　￣￣　　　　, '　|
　　　|　 i　` ､　　　　(　　　 , "　　 |
　　 |　　　　 ` ー---― "|　　　　|
　　　|　　|　　　　　　　　　 i　　　　 |

838:デフォルトの名無しさん
06/08/10 21:33:56
>>835
ソースがなければ、これ以上どうにもならん。

839:デフォルトの名無しさん
06/08/11 00:10:30
ソースがなければ、しょう油をかければ良いじゃない！
文句があるならベルサイユへいらっしゃい！！

840:デフォルトの名無しさん
06/08/15 16:48:28
自宅のWindowsPCでfortranを使いたいのですが
インストール用のディスクを入れると

This kit does not support the currently
running Operating system or Hardware

と表示されます。「あきらめろ」てことでしょうか。
動かす方法があれば教えてくださいっ

841:デフォルトの名無しさん
06/08/15 17:06:54
そのFORTRAN処理系に未練があるのなら、
処理系の動作環境を確認して、
それに応じたOSとハードを用意すればいい。

842:デフォルトの名無しさん
06/08/15 21:08:04
>>840
基本的にはあきらめろってことだｗ

しかし、せめてOSの種類とインストールを試みているコンパイラの名前くらい
書かないとまともな返事は期待できないｗ

843:デフォルトの名無しさん
06/08/17 01:11:50
>>840
つ URLﾘﾝｸ(www.google.co.jp)

844:デフォルトの名無しさん
06/08/20 21:02:50
さあ、楽しい夏休みも残り二週間だ。
みんな、宿題はもう終わったかな？

845:デフォルトの名無しさん
06/08/22 06:28:00
Absoft のコンパイラーをWinXP上で使ってもらったソースのコンパイルをしようとしています。
現在、ソースのコンパイルは出来るのですが、そのプログラムが壊れているようです。
調べてみると以下のような　getarg の部分でつまずいているみたいです。

．．．
character*100 extflag, newbin, oldbin, schnam
integer*2 status
．．．
ii = 1
call getarg(ii, extflag, status)
ii = ii + 1
call getarg(ii, iname, status)
．．．

getargについてすこし調べたのですが説明はすべて二つのパラメターを用いることを前提としているようなので
この部分が何をしているのかわかりません。何をしようとしているのでしょうか？
また、この部分のコンパイルの問題を解決する方法はありますか？
お願いします。

846:デフォルトの名無しさん
06/08/22 10:38:20
GETARG はコマンドライン引数を受け取るサブルーチン。
Fortran2003で正式の規格に入ったが、それ以前のバージョンではメーカーごとの
独自拡張仕様である。が、たいていは存在している。

解決法
１．コマンドライン引数を使わないならコメントにして飛ばす。
２．マニュアルの後ろのほうにある、非標準メーカー拡張の章をしらべて書き直す。
（たいていはPOSIX規格互換になっているので、微妙な修正で解決できると思う。）

847:デフォルトの名無しさん
06/08/30 18:23:24
保守党

848:デフォルトの名無しさん
06/08/31 08:34:44
コンパイル時にエラーがないのに実行できないのはプログラムに問題があるという事なのでしょうか？

849:デフォルトの名無しさん
06/08/31 10:29:29
はい、そうです。

850:デフォルトの名無しさん
06/09/01 12:06:38
実行環境の方に問題がある場合もある。

851:初心者
06/09/02 06:26:17
くだらない質問ですみません・・・
Compaq Visual FortranとFortran77を用いる際、プログラムの違いはあるのでしょうか？？

あと、前者を用いて、分割数を上げていき、円周率を求めるプログラム例がわかりません。
どなたかご教授願えないでしょうか？

852:デフォルトの名無しさん
06/09/02 08:13:32
module oppai
contains
real function pi()
pi = 3.1415926535897932384626433832795028841971693993751058209749445923078164
return
end function pi
end module oppai
program main
use oppai
implicit none
print *, pi()
stop
end program main

853:デフォルトの名無しさん
06/09/02 10:11:39
>>848
chmod a+x

854:デフォルトの名無しさん
06/09/02 10:24:05
program hoge
implicit none
integer, parameter :: WP = SELECTED_REAL_KIND(15)
integer :: i
do i = 8, 160, 8
print *, "分割数:", 4*i, ", π =", oppai(i)*4
end do
stop
contains
pure function f(x)
real(kind=WP), intent(in) :: x
real(kind=WP) :: f
f = 1 / (1 + x**2)
end function f
pure function oppai(n)
integer, intent(in) :: n
integer :: i
real(kind=WP) :: oppai, dx, s
dx = 1.0_WP / n
s = 0
do i = 1, n/2
s = s + f(2*(i-1)*dx) + 4*f((2*i-1)*dx) + f(2*i*dx)
end do
oppai = s * dx / 3
end function oppai
end program hoge

855:デフォルトの名無しさん
06/09/07 23:09:34
すみません。フリーのFortran90/95のプログラム作成ツールありませんか？

それで質問なんですけど。定積分をするときにf(x)=3*x**2などのような簡単な関数は

処理の手前に文関数として定義することで認識されて計算できるのですが、

f(x)=erfc(a*x)+(4*a*c)*exp(2*a*k)+・・・・・・・・・・・・・・・・のような長い関数を

文関数として定義してやると、エラーが出ます。エラーの内容は配列とみなされどうのこう

のといったないようです。どうしたらよいでしょうか

856:デフォルトの名無しさん
06/09/08 03:27:36
>>855
>フリーのFortran90/95のプログラム作成ツールありませんか？
URLﾘﾝｸ(www.g95.org)

>それで質問なんですけど。
URLﾘﾝｸ(ace.phys.h.kyoto-u.ac.jp)

857:デフォルトの名無しさん
06/09/08 20:50:48
ファイルをオープンして入力処理をしているプログラムを安全に中止する方法
はありますか？ctrl+cで強制終了するとファイルが壊れる恐れがあることを聞きました。
linux環境で実行しています。

どうか、よろしくお願いします。

858:デフォルトの名無しさん
06/09/08 23:13:53
>>857
読み込みだけならCtrl-C(SIGINT)を受け取って終了する時にファイルが壊れることは無いと思うが
心配ならSIGINTを無効にするなり、自分でシグナルハンドラを書いて適切な終了処理をさせれば良い。

処理系依存でSIGNALとか言う感じの関数なりサブルーチンが用意されている（ことが多い）

859:デフォルトの名無しさん
06/09/09 02:24:31
858さん

教えて頂いてありがとうございます。
知識不足なので、これを手がかりにやってみます

860:デフォルトの名無しさん
06/09/12 21:24:00
TBasicなのですが･･･
キーボードから西暦年を入力し、その年が閏年かどうかを判定するプログラムを作れ。
但し、if文は一つだけ使うこと。

1.西暦が4で割り切れると閏年である。
2.ただし西暦が4で割り切れても100で割り切れる年は閏年としない。
3.ただし西暦が100で割り切れても400で割り切れる年は閏年とする。

if ((a mod 4 = 0) or (a mod 400 = 0)) and (a mod 100 <> 0) then

こうすると2000年が閏年ではなくなってしまうのですがどうすれば良いのでしょうか？
どなたか教えてください。お願いします。

861:デフォルトの名無しさん
06/09/12 21:56:29
スレ違い

すれ立てるまでもない質問はここで第78刷
ｽﾚﾘﾝｸ(tech板)

862:デフォルトの名無しさん
06/09/12 22:52:28
>>860
結論から先に言うと２０００年はうるう年ではない。カレンダーを見て確かめろ。

つぎにこの課題はすでにこのスレで2回ぐらい出てきている。過去レス嫁。

最後にここはFortranスレなのでBasic野郎は失せやがれ！ｗ

863:デフォルトの名無しさん
06/09/13 02:55:02

>　結論から先に言うと２０００年はうるう年ではない。カレンダーを見て確かめろ。　

え?
ホントに2000年のカレンダー見たの?

864:デフォルトの名無しさん
06/09/13 03:48:58
>>863
グレゴリオ暦じゃないんだろｗ

865:デフォルトの名無しさん
06/09/13 08:01:29
2000 年は閏年だよー。
普通は 100 年ごとに閏年じゃないけど 400 年ごとに閏年なんだよ。

866:デフォルトの名無しさん
06/09/13 13:08:27
そろそろ本題に入ろうよ＞＜

867:860
06/09/13 15:56:18
ご迷惑をお掛けしてすみませんでした。今後気をつけます。

868:デフォルトの名無しさん
06/09/21 13:10:11
はじめまして。
授業でfortranの授業をやっているのですが、さっぱりわかりません。。
ソフト？は７７でやってます。
この問題を誰か解いてください涙。。

行列の計算

　　１　２　３　　　　　－１　－２
A= 4　５　６ B=　－３　－４
　　７　８　９　　　　　－５　－６
A*B
この行列のプログラム、教えていただけないでしょうか。
お願いいたします。

869:デフォルトの名無しさん
06/09/21 14:15:44
>>868
まず、手でA*Bの計算はできるんだろうな？？

870:デフォルトの名無しさん
06/09/21 14:16:44
できます！！！！！！

871:デフォルトの名無しさん
06/09/21 14:48:20
AとBを2次元配列に入れる
ex. A(1,1)=1, B(3,2)=-6

A*B=Cとするなら、
C(1,1)=A(1,1)*B(1,1)+A(1,2)*B(2,1)+A(1,3)*B(3,1)
とうい具合にC(1,1)～C(3,2)まで6個式を書けばいい。

872:デフォルトの名無しさん
06/09/21 14:52:13
do ループをつかわなければならないんです；。

873:デフォルトの名無しさん
06/09/21 14:58:55
すみません。どうにかできそです＾＾

874:デフォルトの名無しさん
06/09/21 21:18:24
>>868
mixiでも質問してだだろ。

875:デフォルトの名無しさん
06/09/21 22:06:02
>>873
がんばれ！

うまくいかなかったら、また聞きに来い

そのかわりちんちんうｐ

876:デフォルトの名無しさん
06/09/29 13:40:50
ルンゲクッタつかおうとおもったら
functionの形じゃないとだめですかね？
というかfunction以外でできる方法あるの？

877:デフォルトの名無しさん
06/09/29 20:54:07
>>876
つ SUBROUTINE

878:デフォルトの名無しさん
06/09/30 00:22:52
簡単すぎて申し訳ないとは思うんですがどうしていいのか...

real dI,dV,dR3,R1,R2,R3,Rx
write(*,*) 'dR3 ?'
read(*,*) dR3
write(*,*) 'R1 ?'
read(*,*) R1
write(*,*) 'R2 ?'
read(*,*) R2
write(*,*) 'R3 ?'
read(*,*) R3
write(*,*) 'Rx ?'
read(*,*) Rx

dI=R2*2*dR3/(R1*R2(R3+Rx)+R3*Rx(R1+R2)+9090(R1+R2)(R3+Rx))

write(*,*) 'dI=',dI

dV=dI*9090

write(*,*) 'dV=',dV
stop
end

879:デフォルトの名無しさん
06/09/30 00:24:32
をコンパイルするとこんなエラーが出ました。
どうすればいいのですか？

0013) dI=R2*2*dR3/(R1*R2(R3+Rx)+R3*Rx(R1+R2)+9090(R1+R2)(R3+Rx))
0014)
*** R2 is already in use as a local name
1 ERRORS [<MAIN@>FTN77 Ver 4.03]

*** Compilation failed

880:デフォルトの名無しさん
06/09/30 00:32:21
13行目途中から掛け算の*が消えているのを直せばいいんじゃない？

始めに消えたR2(R3+Rx)のとこで、R2が実数ではなくて関数とみなされたんで、
もうR2ってのは実数として定義してますよっていうエラーだと思う。

881:デフォルトの名無しさん
06/09/30 00:35:20
いけました！ありがとうございます！

882:初心者
06/10/04 21:11:02
フォートランでサブルーチンを使ってフーリエ変換をするプログラムをつくりましたが、横軸の値のとり方がわかりません！！わかる方いらっしゃいましたらどうか教えてください！

883:デフォルトの名無しさん
06/10/04 22:18:19
>>882
その質問では、そなたがフーリエ変換をよく理解していないという以外分からない。

もう少し、詳しく書き込め　３行で

884:sage
06/10/05 05:43:12
>>882
横軸は、周波数か周期では？

885:初心者
06/10/05 18:46:37
作ったプログラムを書きます。
C
USE NUMERICAL_LIBRARIES
PARAMETER(NN=3000000)
C
REAL XN(0:NN),XX(0:NN),N(10000),COEF(0:NN),FFTN(0:NN)
INTEGER NNF
CHARACTER DATAFILE*35
C
DATAFILE='experience.csv'
OPEN(100,FILE='data.csv')
DO 110,I=1,1024
READ(100,*) XN(I)
110 CONTINUE
CLOSE(100)
C
N=I
NNF=I
CALL FFTRF(NNF, XN, COEF)
DO 200 K=0,(NNF)-1
FFTN(K)=ABS(COEF(K))/DBLE(NNF/2)
200CONTINUE
C
OPEN(200,FILE='OUTPUT.dat')
DO 210,K=0, NNF-1
WRITE(200,*) 　FFTN(K)
210CONTINUE
CLOSE(200)
END
縦軸はFFTN(K)ですが、横軸の周波数をどのようにすればいいのでしょうか？？

886:デフォルトの名無しさん
06/10/06 03:08:16
>>885
プログラムがどうこう言う前にフーリエ変換の基礎から勉強し直せ。

887:デフォルトの名無しさん
06/10/07 00:34:18
関数COTH(X)を使うために、実行文の前に
COTH(X) = 1.0 / TANH(X)
と定義したのですが，コンパイルしてみると
COTH is not an array
とエラーが出ます。どんな原因が考えられますか？

888:デフォルトの名無しさん
06/10/07 02:17:17
>>887
コンパイラは何ですか？　
文関数は最新の規格では廃止になっているので厳格なコンパイラならエラーになります。

その配列と間違えるエラーは、実行文より後に使うと出るものです。

889:887
06/10/07 10:57:16
レスありがとうございます。
プログラムは本当に初心者でコンパイルがどれをさすのか
分かってないのですが、silverfrostのFTN95の個人用を使っています。
関数の定義は、変数の型宣言の直後に書いています。
他に簡単なプログラムを書いてみて、試してみたところ、
うまくいきました。エラーはここ以外に出ていませんが、
warningでIn a previous call to COTH, the first argument was of type
REAL (KIND=1), it is now INTEGER (KIND=3)
と出ています。

890:887
06/10/07 12:36:34
自己解決しました。
他の部分に不備があったようです。

891:デフォルトの名無しさん
06/10/08 15:07:12
f(x)の区間[0,b]の積分を台形公式で近似した
∫f(x)dx=(h/2)*(f(x_0)+2(x_1)+2f(x_2)+…+2f(x_n-1)+f(x_n))
(ただしh=b/nはn分割したときの刻み幅でx_i=h*i)
という式を数値積分をしようと思うのですが、
f(x)=((x^4)*exp(x))/((exp(x)-1)^2)
なのでf(0)が発散してしまいます。
　　こういう場合、どう処理したらよいのでしょうか？
思いついたのは
(f(x_0)+2(x_1)+2f(x_2)+…
という部分を
f(x_1)+2(x_2)+2f(x_3)+…
とする方法（端のほうを無視する）なのですが、他に方法はないのでしょうか？
　　みなさんでしたら、どうするのか、あるいは他に方法はない(思いつかない)のか、
どうか助言お願いします。

892:デフォルトの名無しさん
06/10/08 20:31:02
>>891
変数変換するべし。

893:891
06/10/08 22:58:27
なるほど！
どうもありがとうございます。
少しまた考えて見ます。

894:デフォルトの名無しさん
06/10/10 00:20:42
>>893
というか、そもそもf(+0)は0なのでは？
確かに、f(0)そのものは分母がゼロになるので、計算できないが、
f(x) x->0　はゼロなので、普通に端っこだけ除けば計算できる。
（もちろん、十分結果が収束するまで分割数は多くしないとダメだけど）

それか、計算区間を[0,b]ではなく、[ε,b](εは十分小さい数)で計算するとか。
同じことだけど。

895:デフォルトの名無しさん
06/10/11 22:29:36
ちょっとお聞きしたいのですが
atanとatan2とdatanの違いはなんなのでしょうか？
検索してみたのですがいまいちわかりません（バカですいません）
プログラム書いてるのですが、どう考えてもここしか間違いがないのですが・・・

896:デフォルトの名無しさん
06/10/11 23:10:24
　　　　　引数の型(個数) 関数の型
atan 実数型(1) 実数型
datan 倍精度実数型(1) 倍精度実数型
atan2 実数型(2) 実数型
datan2 倍精度実数型(2) 倍精度実数型

897:デフォルトの名無しさん
06/10/12 00:38:18
>>895
ATANではよく、ATAN(Y/X)　みたいな計算をするが、この場合X＝０の時、
０割の例外が起きて困る。だが実際はATAN(∞）＝π／２だ。

こういうのを避けるのに使うのがATAN2（Y,X)だ。

898:デフォルトの名無しさん
06/10/12 01:50:11
>>896,897
ありがとうございます。
もしかしてatan(-1/-1)=atan(1/1)じゃないんでしょうか？

899:デフォルトの名無しさん
06/10/12 01:58:35
>>898
それはFortranの話ではなく高校数学の話。
atan(-1/-1)=atan(1/1)=atan(1)=π/4

900:デフォルトの名無しさん
06/10/12 02:46:10
>>899
ありがとうございました！
なぞがすべて解決しました。深夜にサンクスでした

901:デフォルトの名無しさん
06/10/12 18:43:21
高速化したい場合は、どのような場合でも-O5をつけてOKですか？

902:デフォルトの名無しさん
06/10/12 23:16:29
>>898
それはいい質問だ。
ＡＴＡＮ（ｙ／ｘ）の場合、象限に関する情報が失われてしまうのも問題で、
ｘ＝０．０のゼロ割りの時と並んで厄介なものである。

ＡＴＡＮ２の場合、その問題にも少し対応している。

ＡＴＡＮの場合は値域が－π/２～π/２だが、ＡＴＡＮ２の場合象限に関する付加的な
情報があるので、返す値域が－π～πに拡大している。

>>901
その質問は難しい。一概には答えられない。
まぁ最適化の手引きのようなものを読んで勉強するしかない。

その上でマニュアルを読んで、使用しているコンパイラの最適化がオプションごとに
どうなっているかを知ったうえで、自分のプログラムと合わせて考えねばならない。

単純には、人間が見ても自明なソースはコンパイラも素直にコンパイル＆最適化
できるので、正しく最適化されることが多い。

903:デフォルトの名無しさん
06/10/13 02:46:34
なんだ、昨日はいろいろなスレにatan厨が湧いたのか

904:デフォルトの名無しさん
06/10/15 15:04:15
FortranのDLLをC#もしくはBasicから呼び出そうとしています。
↓このあたりを参考にコードを書いているのですが、FatalExecutionEngineErrorというエラーが出てしまいます。

URLﾘﾝｸ(www.atmarkit.co.jp)
URLﾘﾝｸ(www.xlsoft.com)
URLﾘﾝｸ(dotnet247.com)

文字列をあつかわない関数に関しては成功しているので、文字列の引渡しがうまくいっていないようです。
DllImport文には以下のようにCharSetの記述をしています。

[DllImport("*****.DLL", CharSet=CharSet.Unicode)]

文字列の渡し方は、stringで渡す、stringBuilderで渡す、char配列で渡す、の3通りを試しましたがいずれもエラーとなりました。

もとのFortranコードおよび
自分が具体的に書いたコードを以下に記します。
なにかお気づきの点などあれば教えていただけないでしょうか？

905:904
06/10/15 15:05:00
オリジナルのFortranコード

subroutine SETUPdll(i,hfld,hfm,hrf,ierr,herr)
implicit double precision (a-h,o-z)
implicit integer (i-k,m,n)
dll_export SETUPdll
parameter (ncmax=20)
character hfld*10000,hfm*255,hrf*3,herr*255,hfmix*255
character*255 hf(ncmax)
do 10 j=1, ncmax
10 hf(j)=' '
if (i.eq.1) then
hf(1)=hfld
else
m=1
do 20 k=1,i
j=index(hfld(m:10000),'|')
if (j.ne.0) then
j=j+m-1
hf(k)=hfld(m:j-1)
m=j+1
endif
20 continue
endif
hfmix=hfm
call SETUP (i,hf,hfmix,hrf,ierr,herr)
end

906:904
06/10/15 15:06:52
自分で書いたコード（char配列を利用）

[DllImport("*****.DLL", CharSet=CharSet.Unicode)]
private extern static void SETUPdll(ref int nc, ref char[] hfld, ref char[] hfm,
char[] hrf, ref int ierr, ref char[] herr);

static void Main(string[] args){

int nc = 1;
char[] hfld = new char[10000];
char[] hfmix = new char[255];
char[] hrf = new char[3];
char[] herr = new char[255];
string strBuff = "fluids/nitrogen.fld";
strBuff.CopyTo(0, hfld, 0, strBuff.Length);
strBuff = "fluids/hmx.bnc";
strBuff.CopyTo(0, hfmix, 0, strBuff.Length);
strBuff = "DEF";
strBuff.CopyTo(0, hrf, 0, strBuff.Length);
strBuff = " ";
strBuff.CopyTo(0, herr, 0, strBuff.Length);
int ierr = 0;

SETUPdll(ref nc, ref hfld, ref hfmix, hrf, ref ierr, ref herr);
}

907:904
06/10/15 15:09:49
自分で書いたコード(stringBuilderを利用)

[DllImport("*****.DLL", CharSet=CharSet.Unicode)]
private extern static void SETUPdll(ref int nc, ref char[] hfld, ref char[] hfm,
char[] hrf, ref int ierr, ref char[] herr);

static void Main(string[] args){

int nc = 1;
StringBuilder hfld = new StringBuilder(10000);
hfld.Append("fluids\\nitrogen.fld");
StringBuilder hfmix = new StringBuilder(255);
hfmix.Append("fluids\\hmx.bnc");
StringBuilder hrf = new StringBuilder(3);
hrf.Append("DEF");
StringBuilder herr = new StringBuilder(255) ;
herr.Append("");
int ierr = 0;
SETUPdll(ref nc, ref hfld, ref hfmix, hrf, ref ierr, ref herr);
}

908:904
06/10/15 15:12:50
自分で書いたコード(stringを利用)

[DllImport("*****.DLL", CharSet=CharSet.Unicode)]
private extern static void SETUPdll(ref int nc, ref string hfld, ref string hfm,
ref string hrf, ref int ierr, ref string herr);

static void Main(string[] args){

int nc = 1;
string hfld = "fluids\\nitrogen.fld";
string hfmix = "fluids\\hmx.bnc";
string hrf = "DEF";
string herr = "";
int ierr = 0;

SETUPdll(ref nc, ref hfld, ref hfmix, ref hrf, ref ierr, ref herr);

}

909:デフォルトの名無しさん
06/10/15 20:03:53
>>904
難しい質問だなｗ
Fortranコンパイラはなんなの？

一般的にC言語などは尻にｘ００のNULL文字が入る内部構造なので、
Fortan側で文字列をいじったあとは尻にNULLをつけたりしないと駄目だったりするが・・・
DECの拡張では尻にCと書くことでNULL文字が付け加わる。

910:904
06/10/17 08:48:25
出張で返信遅れました。
すみません。

コンパイラはSalfordのftn77を使っています。
Salfordのftn95を使って.NET用にコンパイルしたところ、
一応は使えるようになりました。
ただ、一定時間広告が出てしまうのでできるならばftn77でなんとかしたいところです。

!DEC$　の末尾にCも付け加えましたが、やはり同一のエラーが出ます。
もう少し検討してみます。

911:デフォルトの名無しさん
06/10/17 09:49:08
で、公開関数はいつはいるの？

912:デフォルトの名無しさん
06/10/17 19:01:46
intelfortran ver.9.0を使っています。
教官からもらったプログラムをコンパイルしたところ
npionkai.f(2803) : Error: This character is not valid in a format list. [R]
6000 FORMAT(/42H TOO MANY VARIABLE PARAMETERS. YOU REQUEST,I5/47H THIS

-^
compilation aborted for npionkai.f (code 1)
といわれました。
これはformat文の中に[R]を使うなボケェってことなんでしょうか
それならなんか理不尽です

913:デフォルトの名無しさん
06/10/17 21:35:29
>>912
プログラムの"6000 FORMAT ..."の行を晒しては如何。
その書式は教官以上にエッチにならないと使えないらしいよ。

914:デフォルトの名無しさん
06/10/17 21:41:49
ホレリスの文字数がずれてたり、１行当たりの文字数制限でいつの間にかお尻が切れてたりする可能性がある。

915:デフォルトの名無しさん
06/10/17 23:00:53
>>912
たぶん継続行のところでホレリスの数えがずれたんじゃないかな。
その行だけじゃなく、文全部をだせや

916:912
06/10/18 00:11:18
format文全文は
WRITE(2,6000) NINT,MAXINT
6000 FORMAT(/42H TOO MANY VARIABLE PARAMETERS. YOU REQUEST,I5/47H THIS
1 VERSION OF MINUIT IS ONLY DIMENSIONED FOR,I4//)
です。
このformat文があるサブルーチン内で IMPLICIT REAL*8(A-H,O-Z) を宣言しています。

よろしくお願いします。
やはりホレリスずれで最後のFORの[R}が文句いわれてるんですか

917:912
06/10/18 00:15:09
連投すみません
47Hを48Hにかえたらいいんですね

918:デフォルトの名無しさん
06/10/18 00:31:10
>>916
6000 FORMAT(/42H TOO MANY VARIABLE PARAMETERS. YOU REQUEST,I5/47H THIS

ホレリスの途中で行を分けてるのが不味いな。
行末の「THIS」の後ろにスペースが入ってるとそれがカウントされて文字数がずれる。

919:デフォルトの名無しさん
06/10/18 08:04:08
Fortran77を使っているんですが，
乱数を発生させようと思い

real a

do i= 1,10

a = rand()

write(*,*)a

enddo

stop
end

このようなプログラムを作ったところ，
乱数は書き出されるのですが毎回同じものになってしまいます．

実行するたびに違う乱数を得るにはどうしたらいいのでしょうか．

920:デフォルトの名無しさん
06/10/18 09:34:30
>>919
乱数を使用する前に SRAND で乱数の種を設定する。

PROGRAM HOGE
REAL A
CALL SRAND(TIME())
DO I=1, 10
A = RAND()
WRITE(*,*) I, A
END DO
STOP
END

但し、RAND/SRAND/TIME はFORTRAN77標準ではないので
もし使っている処理系にない場合は同等のサブルーチン/関数を探すこと。

921:デフォルトの名無しさん
06/10/18 11:19:26
>>920さんのプログラムを入力してみたのですが
やはり何回やっても同じ値しか返してきません。

>もし使っている処理系にない場合は同等のサブルーチン/関数を探すこと

なんでしょうか。もう少し調べてみます。

922:デフォルトの名無しさん
06/10/18 12:42:59
>>921
F90にはRANDOMIZEというのが標準であったはずだが・・・

しかし、この手の問題はBASIC全盛時代のゲームの乱数の頃からのFAQだｗ
時計で乱数の種を作る。

宿題程度ならともかく、まじめなシミュレーションをやりたいなら、乱数の専門サイトへ行って、
乱数マニアの御託を我慢して読んで、目的にあった適切な擬似乱数のルーチンを貰ってくるのが吉。
大体言語内蔵の乱数はウンコなことが多いという経験的な事実がある。

金が余っているならハードウェア的に乱数を吐き出すボードもあるようだ。
熱雑音か放射性物質の崩壊を利用して天然物の乱数を出してくれる。

923:デフォルトの名無しさん
06/10/18 21:00:08
fortran77で計算結果をfile.csvに記録するとき

WRITE(3,*)a(1),',',a(2),',',a(3)・・・
（a(n)は数値）

とカンマ区切りで出力しているのですが、70文字目くらいで折り返されて
記録されているため、Excelで開いたときに1行で結果が見られません。
これを1行で記録するにはどうすればいいのでしょうか？

924:デフォルトの名無しさん
06/10/18 21:45:57
>>923
横着せずに書式を指定する

write(3, '(100(F8.3,:,'',''))') (a(i), i=1,100)

925:923
06/10/18 22:11:40
>>924
出来ました！
ありがとうございます。

926:ふつふつ
06/10/19 23:25:37
fortran77を使っているんですけど、得られた６つの値を１行１列、２行１列、３行１列
１行２列、２行２列、２行３列の順番にformat文を使って表示するにはどのようにしたら
いいでしょうか？初心者でぜんぜんわからなくて。どなたかお願いします。

927:デフォルトの名無しさん
06/10/19 23:33:15
>>926
>得られた６つの値
どんな形式で与えられるのか、詳しく説明せい。

928:ふつふつ
06/10/19 23:41:23
>>927
説明不足ですみません。E25.16E3です

929:デフォルトの名無しさん
06/10/19 23:44:46
掲示板のHTTML教えてください

930:ふつふつ
06/10/20 00:45:03
ほんとごめんなさい。どこに書いてあるのですか？htmlって？

931:デフォルトの名無しさん
06/10/20 21:41:25
>>926
データ数とフォーマットによる表示数が合ってない。
>>927の書式もおかしい。

問題文を忠実に写せ。このままでは答えようが無い。

932:ふつふつ
06/10/22 02:30:06
Nx=256,Ny=256で、以下のプログラムでやっています。これを一列目に２５６個の結果を出し、２５７個目から
２列目に結果を表示していき、行２５６、列２５６にしたいんです。E25.16E3は変えられません
　　　　DO 250 I=1,Nx
　 DO 250 J=1,Ny
　write(30,260) U(I,J)
250 continue
260 Format (E25.16E3)

933:デフォルトの名無しさん
06/10/22 03:04:50
>>932
なんか >>926 と言ってることが違うんだが…
>>932のプログラムから察するとこういうことか？

WRITE(30, 260) ((U(I,J), J=1,256), I=1,256)
260 FORMAT(256(256E25.16E3,/))

934:ふつふつ
06/10/22 17:45:27
実はこのプログラムが実際のものだったんです。逆にややこしくしてしまったようです。
それにもかかわらず答えていただきありがとうございました。
明日試してみます。

935:デフォルトの名無しさん
06/10/23 20:21:09
学校でfortranやらされる事になった・・・

936:デフォルトの名無しさん
06/10/23 20:29:55
>>935
おめでとう！

937:デフォルトの名無しさん
06/10/23 23:10:38
>>935
FORTRANはいいぞぅ！

1950年代のコードが未だに生き残っている。

ソース再利用では最高水準の言語だｗ

938:デフォルトの名無しさん
06/10/23 23:21:30
7カラムにぴったり揃えたあの美しさ。インデント糞喰らえ。
これに匹敵できるのはアセンブリしかない。

939:デフォルトの名無しさん
06/10/23 23:27:21
おいおい、そんなこと言ってるとCOBOLerの皆さんが…

940:デフォルトの名無しさん
06/10/24 10:45:52
× 匹敵できる
○ 匹敵する
○ 比肩できる

941:stack overflow
06/10/24 14:17:33
severe(170): Program Exception - stack overflow というエラーが出るのですが、解決方法が分かりません。
過去の書き込みに同様のものがありますが、もう1度スタックのサイズの増やし方を教えてください。
あと、配列はすべてallocatableで宣言していて、deallocateを繰り返しているのですが、stack overflowになってしまうものなんでしょうか？
エラーの出る行がcall文だし、そもそも、どうするとメモリが消費されるのか分かりません。
ど素人にご教授願います。

942:デフォルトの名無しさん
06/10/24 20:07:34
fortran4で書かれたコードを使いたいのですがgccではコンパイルできませんでした。
fortran4でもコンパイルできる今手に入るコンパイラはないでしょうか?

943:デフォルトの名無しさん
06/10/24 22:24:00
>>941
>スタックのサイズの増やし方
Windows上のCompaq Visual Fortranなら
df hoge.f90 /link /stack:0x12ab000

Intel Visual Fortranなら
ifort /F0x12ab000 hoge.f90

ほかは知らん。

>エラーの出る行がcall文
(FORTRAN77形式の)次元と範囲が明示された配列を引数に取るサブルーチンに
(Fortran90以降の)部分配列やポインタを渡して呼び出そうとすると
暗黙のうちに一時的な配列が作られ、これがサブルーチンに渡される。
このとき要求される配列のサイズが大きいとスタックがあふれる可能性がある。

944:stack overflow
06/10/24 23:06:34
エラーの出るサブルーチンの引数に配列は含まれていません。
ループを繰り返していると、見えないところで、メモリが切られて、スタックがあふれてしまう？？？っぽいんですが・・・

945:デフォルトの名無しさん
06/10/24 23:44:35
動的メモリー割付をとる領域には、ヒープとスタックの２種類がある。
普通はヒープはOSの管理下に、スタックはプログラムの管理下にある。
ゆえにスタックサイズはリンカーで指定することになる。

allocateはヒープ領域からメモリーを取られるのが一般的だ。

スタックの方は、オートマチック変数などに取られる。
呼び出し側のサブルーチンで、オートマチック変数を取っていないか？
これがでかすぎるとスタックオーバーフローが起こる。

コンパイラの実装にもよるが、スタックオーバーフローはａｌｌｏｃａｔｅではないんじゃないか？

946:デフォルトの名無しさん
06/10/24 23:55:32
>>942
ＣＰ／Ｍのエミュレーターを入れて、その上でＦＯＲＴＲＡＮ　ＩＶのコンパイラを動かす
ことは出来る。エミュレータとソフトはネットに流れているので拾える。

というのは冗談だ。

ＦＯＲＴＲＡＮＩＶはＦＯＲＴＲＡＮ６６と同じものだと考えてよい。
ＦＯＲＴＲＡＮ６６は、ごく些細な部分を除けばＦＯＲＴＲＡＮ７７に含まれている。
したがって、（色々落とし穴はあるが）若干の修正で動くはずである。

ＦＯＲＴＲＡＮ９０→９５→２００３　と進むにつれて徐々に古臭い機能を切り落としているが、
現実の実装では機能は生き残らせてあるのが普通だ。
新しいコンパイラでも動くだろう。

ＤＥＣ／ＩＮＴＥＬ　ＦＯＲＴＲＡＮなどではＦＯＲＴＲＡＮ６６形式のＤＯ廻り
（条件にかかわらず必ず１回は回る）とかをオプションで指定できる。
あと拡張ＤＯ　ＬＯＯＰもＤＥＣ時代には対応していた。

ＧＣＣでどうなってるかは知らん。

947:942
06/10/25 00:16:11
>>946
回答ありがとうございます。intelで試してみます

948:デフォルトの名無しさん
06/10/25 01:33:30
大学でFortranを使うことになったのですが、詰まってしまいまして。
以下のようなプログラムを作るのですが……。

１．すでにAと言うデータファイルがあり、2列で構成されている。
　　それぞれの列が粒子のx座標、y座標に対応で、1行1粒子に当たります(座標情報は実数です)

２．ここから新たに別のファイルBを作成して、
　　Bに原点からの距離(sqrt(x**2+y**2))を計算して入力していく。

というものです。
試しにプログラムを作り入力データが１行のAを用いたところ、できたBの中に距離を表す数値が
一つだけ表示されていて成功したのですが、行数が２行以上になると、どう試してもできず。
Aの各行の値がばらばらでも、Bで表示される数値が全部の行で同じに……orz

行数が少ないなら自力で計算するところなのですが、Aが数千万行クラスのファイルなので
手計算ではどう考えても時間がかかりすぎます。

内容的に初歩の物だとは思うのですが、幾つかのサイトを見てみても複数行のデータを
うまく読み込み、計算して書き出す方法がわからず、先に進めません。

どなたか、御力をお貸しくださいませ。よろしくお願いしますm(_ _)m

949:デフォルトの名無しさん
06/10/25 02:08:57
>>948
program main
implicit none
real(kind=SELECTED_REAL_KIND(10)) :: x, y, d
integer :: ios
open(unit=11, file="A", status="OLD")
open(unit=12, file="B", status="REPLACE")
do
read(unit=11, fmt=*, iostat=ios) x, y
if (ios /= 0) exit
d = sqrt(x**2 + y**2)
write(unit=12, fmt=*) d
end do
close(11)
close(12)
end program main

950:948
06/10/25 15:46:23
>>949
すぐに答えて頂いたのに、返答が遅くなってしまい申し訳ありませんm(_ _)m
さっそく上記のプログラムを用いたところ、無事計算されたファイルが
出力されました。
これでようやく先に進めそうです。本当にありがとうございました!!

ただ、一つ疑問点が……。
3行目のreal以下の（）文を書いたままコンパイルすると、

Invalid declaration of or reference to symbol
`selected_real_kind' at (^) [initially seen at (^)]

という警告文がでて、コンパイルが止まってしまいました。
中身を含めて（）を消してみたところうまくコンパイルできたのですが。

もしよろしければ、（）内の意味を教えていただけないでしょうか。
重ね重ね申し訳ありませんが、よろしくお願いします。

951:stack overflow
06/10/25 23:16:42
返信遅くなりました。昨日(945)の続きお願いします。
今日は寝ないつもりです。(もちろんデフォルトの名無しさんには強制しません)
オートマチック変数ってどのような変数のことをいうのでしょうか？
いろいろ調べてみたのですが、いまいちです。
結局、オートマチック変数が原因の場合どうすればよいのでしょう？

952:デフォルトの名無しさん
06/10/25 23:28:30
>>950
あれはＦＯＲＴＲＡＮ９０以降で導入された精度保障のための関数だ。
kind=SELECTED_REAL_KIND(10)
これは有効数字１０桁を保障するような型に対応する種類（ＫＩＮＤ）を返す。
普通のＦＯＲＴＲＡＮは短精度と倍精度を持つが、短精度の有効桁は７～８、
倍精度は１４～１５なので、倍精度に相当する値を返す。
ＦＯＲＴＲＡＮ７７で言えばＲＥＡＬ（８）に相当する。

これでコンパイラのエラーが出るとすると、ＦＯＲＴＲＡＮ９０完全準拠で無いと思われる。

953:デフォルトの名無しさん
06/10/25 23:39:38
>>951
オートマチック変数または自動変数とは、サブルーチンが呼び出されたときに
一時的に割り付けられる変数だ。

SUBROUTINE unko(a, n)
IMPLICITE NONE
REAL, INTENT(IN) :: a(n)
REAL :: work(SIZE(a))
云々

とあったとすると、配列ｗｏｒｋが自動変数にあたる。
これはあまり大きなものを取ることを想定しておらず、スタック領域に割り付ける。

こういう使い方で大きな配列などを取っていると、スタックオーバーフローがおきる。

サブルーチンを呼び出すと、戻る場所などの情報をスタックに保存することになる。
またサブルーチンで一時的に使用するメモリー等もスタックに取る。
サブルーチンの呼び出しが過剰に多いか（再帰関数が無限ループとか）、
一時変数が領域を食いすぎるとオーバーフローが出る可能性がある。

エラーメッセージからするとＤＥＣ系のコンパイラの感じだが、トレースバック情報などの
デバッグ情報を吐き出させるオプションをつけてコンパイルし実行してみるのも手だ。

まぁ、もう少し詳しい情報が無いと助言のしようが無い。

954:stack overflow
06/10/26 00:19:00
programを省略して書くと、

(宣言文)
open(21,file='filelist')
do
allocate
read(21,*,end=1000)unko
open(51,file=unko)

　(計算)

close(51)
deallocate
end do
1000 close(21)
stop
end
データのファイル(具体的には地震動の波形)がたくさんあって
1つ１つに同じ処理をして出力ファイルを作成したかったので、
cmdで入力ファイルのパスのリストを作って、リストのファイルが読み終わるまで、
doループで繰り返すように動かしたいのです。
大体2000回くらいループを繰り返すとstack overflowします。(動かして20分ぐらい)
mainで使う配列はすべて、doループのなかでallocateとdeallocateを繰り返しています。
subroutine内でも引数でない配列は、allocateとdeallocateをしています。
こんな風な動かし方ってできないんですかね?

OSはWindowsXP、programはvisual studioで動かしています。

955:デフォルトの名無しさん
06/10/26 00:43:27
>>954
その範囲内では問題ないと思われる。
２０分はちゃんと動いているとすれば、結構厄介なBugかもしれない。
Ｖｉｓｕａｌ Studioを使っているという事は、コンパイラはＩｎｔｅｌのものでいいのか？

昔のＭＳPowerStationだとＩ／ＯにメモリーリークがあるのでＩ／Ｏを繰り返すと
死にメモリーが増えて行きどこかで死亡することになっていた。
これはタスクモニターでメモリー使用量を観察していると分かる。

Ｉｎｔｅｌではそういう事は無いと思うのだが、もしかしたらそういう可能性が無いともいえない。

ALLOCATEとDEALLOCATEの対応が狂って、きちんとメモリーが解放されていない場合も
タスクモニターでメモリー使用量を見ていると、どんどん使用メモリーが増えてゆくので
見て取れる。

ところで実行はDEBUG MODEでためしているだろうか？

956:stack overflow
06/10/26 00:47:04
プログラム概要を記します。

(宣言文)
open(21,file='filelist')
do
allocate
read(21,end=1000)unko
open(51,file=unko)
read(51)
(計算)
close(51)
deallocate
enddo
1000　close(21)
deallocate
stop
end

データ(具体的には地震動波形(6000stepから10000stepぐらい))のファイルがたくさんあって、
１つ１つに同じ処理をしたいので、cmdでファイルのパスのリストを作って、このリストを読み終わるまでdoループを繰り返しているのですが、
2000回くらい繰り返すと(大体20分ぐらい)stack overflowになります。
mainで使う配列はすべてdoループのなかでallocateとdeallocateを繰り返しています。
subroutineでも引数でない配列はallocateとdeallocateを繰り返しています。

こんな使い方しませんか？

OSはwindowsXPで、visual studioでプログラムを作ってます。

957:stack overflow
06/10/26 00:47:56
プログラム概要を記します。

(宣言文)
open(21,file='filelist')
do
allocate
read(21,end=1000)unko
open(51,file=unko)
read(51)
(計算)
close(51)
deallocate
enddo
1000　close(21)
deallocate
stop
end

958:stack overflow
06/10/26 00:48:58
957の続き

データ(具体的には地震動波形(6000stepから10000stepぐらい))のファイルがたくさんあって、
１つ１つに同じ処理をしたいので、cmdでファイルのパスのリストを作って、このリストを読み終わるまでdoループを繰り返しているのですが、
2000回くらい繰り返すと(大体20分ぐらい)stack overflowになります。
mainで使う配列はすべてdoループのなかでallocateとdeallocateを繰り返しています。
subroutineでも引数でない配列はallocateとdeallocateを繰り返しています。

こんな使い方しませんか？

OSはwindowsXPで、visual studioでプログラムを作ってます。

959:stack overflow
06/10/26 00:50:11
957の続き
データ(具体的には地震動波形(6000stepから10000stepぐらい))のファイルがたくさんあって、
１つ１つに同じ処理をしたいので、cmdでファイルのパスのリストを作って、このリストを読み終わるまでdoループを繰り返しているのですが、
2000回くらい繰り返すと(大体20分ぐらい)stack overflowになります。
mainで使う配列はすべてdoループのなかでallocateとdeallocateを繰り返しています。
subroutineでも引数でない配列はallocateとdeallocateを繰り返しています。
こんな使い方しませんか？
OSはwindowsXPで、visual studioでプログラムを作ってます。

960:stack overflow
06/10/26 00:51:33
すいません。
同じ内容を何度も書き込んでしまいました。

961:stack overflow
06/10/26 01:17:54
ファイルリストがいくつもあるのでバッチファイルで動かしていました。
filelistをcall getargで読み込ませています。
DEBUG MODEでやった方がいいですか？
タスクモニターを見ながらもう一度動かしてみようと思います。

962:stack overflow
06/10/26 01:43:26
つかぬ事をお聞きしますが、タスクモニタのどこをみればよいのでしょうか？
見ながら動かして見ましたが、あまり変わりません。

コンパイラはIntel FORTRANです。

963:948
06/10/26 01:58:24
>>952
なるほど、精度保障のための関数でしたか。勉強になりました。
一度、自分の使っているFORTRANのバージョンを確認してみます。
お答えいただき、どうもありがとうございました！

964:stack overflow
06/10/26 02:15:31
!!!!なぜか、急に動くようになりました。????
もう少し長い時間動くかどうか、やってみます。結果は明日報告します。

ご丁寧に説明していただき、ありがとうございました。お時間取らせてしまってすいません。
and同じ書き込み何度もして、すいません。

965:デフォルトの名無しさん
06/10/26 02:39:52
>>961
一応DEBUG　MODEで1回は動かしておくのがいい。
時間はかかるので小さめのデータを用意して確かめておくのがマナー。
いきなりリリースモードに行くのは、体を洗わずに湯船につかるようなものｗ

>>962
概略ならパフォーマンスのところを見るとグラフが書かれているのでそれで傾向を知れる。
プロセスのところでMEMUSAGEをみれば、数値でどのくらい使われているかが分かる。

966:デフォルトの名無しさん
06/10/26 05:22:04
f90glを使いたいのですが、導入の方法からつまづいています。

URLﾘﾝｸ(math.nist.gov)
を取ってきたあと、どうすれば導入＆コンパイルできるのですか？
環境はWindowsXP、cygwinでインストールしたOpenGLを使っています。
コンパイラは何を使うのでしょう？gfortran？g95？

967:stack overflow
06/10/26 09:07:18
>>965
お騒がせしました。program無事に終了しています。
stack overflowが出る前は、debug modeで試すというのはやっていましたが、
そのあとの修正が反映されてなかったってことなのかな・・・
コンパイルはしてたはずだったんですが・・・
初心者なもんでお許しを
何がともあれできてよかったっす。ありがとうございました。

968:stack overflow
06/10/26 10:53:04
できたと思ったのですが、やっぱりできていませんでした。
昨日か一昨日、どこでメモリが食われているのか、調べるために
出力の部分をコメントアウトしたりしていたのですが、
コメントアウトしたままで実行したので動いていたようです。
やはり、出力時にメモリが使われているようです。
出力はバイナリ形式なんですが、なにか関係があるのかなあ・・・

969:デフォルトの名無しさん
06/10/26 22:03:07
>>966
まず、f90glから使える GLUTが必要。
これは（少々バージョンが古いようだが）ソースとWindows用にコンパイルした物が
f90gl Softwareのページから拾える。

次に f90gl-*/INSTALL, mf_keyをよく読んで自分の使っているOS/C/Fortranコンパイラに
合ったMakefile(mf*)またはバッチファイル(mf8n*.bat)を探す。
適切な物が見つかれば、それを自分の環境に合わせて修正し
makeまたはバッチファイルを実行し、ライブラリとモジュールファイルを作る。

適切な物がなければ、似たような物から自分でMakefile/バッチファイルを作るか
Makefile/バッチファイルの処理内容を見て手動でコンパイルする。

f90gl-1.2.12.zipにはcygwin gcc + g95/gfortran用のmf*は含まれていないようなので
自分でMakefileを作るしかないようだな。

次ページ