差分

この文書の現在のバージョンと選択したバージョンの差分を表示します。

--- fortran [2009/07/10 16:13]
10.249.202.217
+++ fortran [2009/12/01 13:59] (現在)
saito
@@ ライン 1: / ライン 1: @@
 ====== Fortranのtips ======
+  * intelのコンパイラは[[ifort]]参照のこと．
 ===== コンパイラオプション =====
@@ ライン 31: / ライン 31: @@
 </code>
-このファイルをperlに食わせて処理すると便利かも．
+処理は[[movie]]参照のこと．
-===== ifortで改行を抑制する =====
-ifortでは出力時に勝手に改行する仕様になっている。改行を抑制するためには、Format文を使用すればよい。
+===== makefileでコンパイル =====
+複数のファイル（サブルーチンetc.）で構成されるプログラムを動かすのに便利．
+以下の文を「makefile」という名前で作り，コンパイルしたいファイルと同じディレクトリに置く．
+（以下，main.f90, modules.f90[モジュールファイル], sub1.f90, sub2.f90をコンパイルする場合）
+<code>
+TARGET = a.out
+OBJECTS = modules.o main.o sub1.o sub2.o
+F90 = ifort
+FLAGS = -check all -warn declarations -CB -fpe0 -traceback
+COMMON_MOD = modules.f90
+.SUFFIXES :
+.SUFFIXES : .o .f90
+.f90.o:
+	${F90}	-c $< ${FFLAGS}
+${TARGET} : ${OBJECTS}
+	${F90}	-o $@ ${OBJECTS}
+clean :
+	rm  *.o ${TARGET} *.mod
+${OBJECTS} : ${COMMON_MOD}
+</code>
+あとは,makeコマンドを打てばコンパイル完了．そして実行（上の例では./a.outを打てば良い）．
+<code>
+$ make
+$./a.out
+</code>
+また，以下の様にmake cleanで実行ファイル(a.out)及びオブジェクトファイル(*.o)を消去できる．
+<code>
+$ make clean
+</code>
+以下，makefileの変数について簡単に説明する．
+<code>
+TARGET = 実行ファイル名
+OBJECTS = モジュールファイル プログラムファイル (すべて拡張子は.oとする．)
+F90 = コンパイラコマンド (gfortranを使うなら「gfortran」)
+FLAGS = コンパイラオプション
+COMMON_MOD = モジュールファイル (拡張子は.f90)
+</code>
+===== データファイルの読み込み =====
+データファイルの中に何行あるか分からないときにはiostatかendを使う．
+iostatの方がスマートなやり方かも．
+==== やり方その１ ====
+read文にiostatのオプションを付ける．
+iostatが正の数を返す場合はエラー，負の値を返す場合は行の終わりを意味するらしい．
 <code fortran>
-integer n
+program hoge
-real(8),allocatable::a(:)
+  implicit none
-n=9
+  real a,b
-allocate(a(9))
+  integer i,inputstatus
-do i=1,n
+  do i=1,10
-   a(i)=i
+     read(*,*,iostat=inputstatus) a,b
+     if(inputstatus>0) stop "error"
+     if(inputstatus<0) exit
+     write(*,*) a,b,a**2-b
+  end do
+end program hoge
+</code>
+==== やり方その２ ====
+endを使う．
+以下サンプル．プログラムの中では10回ループを回すが，読み込むのは5行分だけ．
+<code bash>
+gfortran hoge.f90
+./a.out < hoge.dat
+</code>
+hoge.f90
+<code fortran>
+program hoge
+  do i=1,10
+     read(*,*,end=999) a,b
+     write(*,*) a**2-b
+  end do
+continue
+end program hoge
+</code>
+hoge.dat
+<code>
+. 1.
+. 4.
+. 9.
+. 16.
+. 25.
+</code>
+===== NaNの検出 =====
+浮動小数点演算例外を検出するコンパイラオプションがあるので、それを使用する。
+===== Makefile =====
+保存するためのmakefile
+<code bash>
+DATE=  `date +%y%m%d%H%M`
+TAR = backup/backup${DATE}.tgz
+backup :
+     tar czvf ${TAR} Makefile *.f90
+</code>
+===== implicit none =====
+宣言文で宣言した変数以外を利用しないようにする。
+長いプログラムを書く場合はミスを避けるためにも使用するべき。
+自由形式の場合はuse文の後に使用。
+<code fortran>
+program hoge
+  use hoge_lib
+  implicit none
+  real(8) a,b,c,
+  ...
+</code>
+===== 装置番号 =====
+write(?,*)やread(?,*)の?に入る番号。
+通常、５が標準入力（キーボード）、6が標準出力（ディスプレイ上の端末）に使用されるので、
+と6を使ってopenするのは避けるべき。
+slatecのI1MACHを使えば、標準入力、標準出力の番号が調べられる。
+装置番号にはcharacterも入る。
+===== 連番のファイル名をつける。 =====
+ファイル名をcharacter変数に代入して使う。ファイル名の変更にはwrite文を使う。
+逐次計算に使える。
+<code fortran>
+character(len=12) filename
+...
+do i=1,100
+  write(filename,'(a,i4.4,a)') 'hoge',i,'.dat'
+  open(10,file=filename)
+  ...
+    write(10,*) ...
+  ...
+  close(10)
 end do
-write(*,'(<n>f)') a
-end program
 </code>
+===== 簡単なシェルスクリプト =====
+連続してプログラムを実行したいとき、shファイルが便利。
+端末で入力する内容をあらかじめファイルに書いておき、それを実行すれば、上から順に実行してくれる。
+windowsで言うバッチファイル。
+<code bash>
+gfortran hoge.f90
+cd case10
+../a.out
+cd ../case20
+../a.out
+cd ../case30
+../a.out
+</code>
+<code>
+$ sh hoge.sh
+</code>