FORTRAN III

.SUFFIXES:
.SUFFIXES:.f .c .o
を頭のほうにいれてみるとかどう？

>>952
どうもありがとう。やってみましたが、うまくいきませんでした。

main.o:main.c
main.c:main.f

subr.o:subr.c
subr.c:subr.f

を加えればうまくいくのですが、実際に全部加えるのは大変。。。

>>950のだと*.oと*.f間に依存関係が記述されてないからなんじゃない?
>>953みたいに依存関係を記述するのが正解だと思うけど。

そうしないとすると*.fから*.oをつくるようにするとか。
確認していないけど、こんなイメージ。
.f.o:
$(F2C) $(F2CFLAGS) $< -o hoge.c
$(CC) $(INCLUDES) $(CFLAGS) hoge.c -o $@

多分もっといいやり方があると思う。

newmann条件の熱方程式を差分法を使って解きたいんですけど、うまくできません。
(du/dt) = (d^2u)/ (dx^2) (0≦x≦1)
du/dx=0 in x=0,1
u=cos(pi*x) in t=0
どなたかご教授お願いします。
program sabunhou
parameter (ndim=101)
dimension u(ndim)
write(*,*) 'INPUT : nfunc='
read(*,*) nfunc
write(*,*) 'INPUT : n='
read(*,*) n
write(*,*) 'INPUT : dt='
read(*,*) dt
write(*,*) 'INPUT : tmax='
read(*,*) tmax
h=1.0/float(n)
flam=dt/h**2
do 1 i=1,n+1
fi=i-1
u(i)=f(fi*h,nfunc)
1 continue
write(*,100) 0.0,(i,u(i),i=1,n+1)

100 format(1h ,'t=',1pe12.41h ,5(2x,'i',5x,'U(i)',4x)(1h ,5(i3,1pe12.4,1x)))
nstep=0
2 continue
uleft=u(-1)
do 3 i=2,n
unew=(1.0-2.0*flam)*u(i)+flam*(u(i-1)+u(i+1))
u(i-1)=uleft
uleft=unew
3 continue
u(n)=uleft
u(n+1)=u(n-1)
nstep=nstep+1
t=dt*float(nstep)
write(*,100) t,(i,u(i),i=1,n+1)
if(t+dt.lt.tmax+0.1*dt) go to 2
end
function f(x,nfunc)
pi=3.1415927
f=cos(pi*x)
end

>>952

%.o:%.f

をいれると.fと.oの間のデフォルトルールが消されて
うまくいくっぽい。

>>955
何がうまく行かないのか、もう少し具体的に書いてくれないと答えられない。

世間的には、いきなりソース見せて、『おかしいの教えて！』と言ったら。
『おかしいのはお前の態度だ！』といわれちゃうｗ

>>949
いまいち理解できないのだが、データそのものを乱数で生成したいということか？
年号に関してはポアソン分布をなすように、それに対応する６個の値は別の分布ということ？

それならば、それぞれの分布をなす関数が用意してあればDO LOOPをまわすだけだ。
それらをサブルーチンかファンクションで作ればいい。

>>955
境界の位置のuの値が固定されてしまっていました。
なんとか解決することができました。
ありがとうございました。
煮詰まった余りおかしな態度だったことをお詫びします。

>>959
OKｗｗｗｗｗｗ
頼りになる人は気難しいし、気さくな人は頼りにならん。
奴らの操縦法は下手に出ておだてる事と、適切な質問を与えることだ。

g95　と　ｇｆｏｒｔｒａｎ　の関係について教えてください　（＞＜

できれば３行で！

無
関
係

>>926
いや、先祖はｇ９５で同じだろｗ　
ｇ９５の中の人がＧＣＣと一緒になるのを嫌ったのかな？
詳しいことしらん。

gfortran はGCC(GNU Compiler Collection)の中のFortran95 コンパイラ。
g95 はGCCを元にしたFortran95コンパイラ。
g95の進捗があまりにも遅かったのでgfortranが別プロジェクトとして分岐した(が、今のところ完成度はg95の方が上)。

ttp://gcc.gnu.org/wiki/TheOtherGCCBasedFortranCompiler

Mersenne Twisterで乱数を発生させたいのですが（fortran77）、
http://www.math.sci.hiroshima-u.ac.jp/~m-mat/MT/VERSIONS/FORTRAN/TAKANO/mt19937.f
（http://www.math.sci.hiroshima-u.ac.jp/~m-mat/MT/VERSIONS/FORTRAN/fortran.html）
上記のプログラムだと[0,1]の一様乱数発生になります。
これを(0,1)の一様乱数発生に変えるには、プログラムをどのように変えればいいのでしょうか。


簡単に変えられるようでしたら、(0,1]と[0,1)も教えてください。

>>965
メルセンヌ・ツイスターの本家ページにＣ版の半区間を見た記憶があるが、幻かもしれない。

C版にはあるみたいだね。
http://www.math.sci.hiroshima-u.ac.jp/~m-mat/MT/MT2002/CODES/MTARCOK/readme-mt.txt
C版の関数を呼ぶなり、移植するなりするということで。

mt19937ar.fは全部入り

>>966-968
ありがとうございます。
mt19937ar.fを使ったらうまくいきました。

そろそろ次スレが要るな。

次スレ

FORTRAN �W
http://pc8.2ch.net/test/read.cgi/tech/1163319215/

　

　

　

　

　

ちゃんと最後まできっちり使えよ！

FORTRANIVはすなわちFORTRAN66だなｗ

　

　

　

　

　

CONTINUE

　

　

　

　

　

　

　

　

　

　

　

　

　

　

　

do 1000 i =1 ,n

STOP

このスレッドは１０００を超えました。
もう書けないので、新しいスレッドを立ててくださいです。。。