差分

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

--- fortran [2009/09/18 01:30]
122.103.100.45
+++ fortran [2009/12/01 13:59] (現在)
saito
@@ ライン 1: / ライン 1: @@
 ====== Fortranのtips ======
+  * intelのコンパイラは[[ifort]]参照のこと．
 ===== コンパイラオプション =====
@@ ライン 31: / ライン 31: @@
 </code>
-このファイルをperlに食わせて処理すると便利かも．
+処理は[[movie]]参照のこと．
-===== ifortで改行を抑制する =====
-ifortでは出力時に勝手に改行する仕様になっている。改行を抑制するためには、Format文を使用すればよい。
-<code fortran>
-integer n
-real(8),allocatable::a(:)
-n=9
-allocate(a(9))
-do i=1,n
-   a(i)=i
-end do
-write(*,'(<n>f)') a
-end program
-</code>
@@ ライン 92: / ライン 78: @@
 COMMON_MOD = モジュールファイル (拡張子は.f90)
 </code>
+===== データファイルの読み込み =====
+データファイルの中に何行あるか分からないときにはiostatかendを使う．
+iostatの方がスマートなやり方かも．
+==== やり方その１ ====
+read文にiostatのオプションを付ける．
+iostatが正の数を返す場合はエラー，負の値を返す場合は行の終わりを意味するらしい．
+<code fortran>
+program hoge
+  implicit none
+  real a,b
+  integer i,inputstatus
+  do i=1,10
+     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
+</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>