差分

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

--- fortran:tips [2014/01/15 19:40]
saito [コンパイラオプション]
+++ fortran:tips [2017/10/03 12:55] (現在)
@@ ライン 2: / ライン 2: @@
   * intelのコンパイラは[[fortran:ifort]]参照のこと．
   * [[http://www.owlnet.rice.edu/~ceng303/manuals/fortran/]] が役に立ちそう．
+  *
+===== よくある間違い =====
+コーディングにおいてよくある間違いをまとめる．
+==== べき乗の順序 ====
+<nowiki>計算の優先順位は
++- < */ < **
+である．</nowiki>
+<code fortran>
+print*,-1**0
+</code>
+<nowiki>
+これは-(1**0)と解釈されるので，出力結果は−1．
+-1に対してべき乗する場合は，(-1)**0．
+</nowiki>
+==== 論理演算の順序 ====
+<nowiki>
+論理演算子の優先順位は
+.eqv., .neqv. < .or. < .and. < .not.
+である．
+</nowiki>
+<code fortran>
+logical :: a=.false.
+num :: i=1
+if(a.eqv..false..or.i==1) then
+print*,'if'
+else
+print*,'else'
+end if
+</code><nowiki>
+これは，a.eqv.(.false..or.i==1) と解釈される．
+(.false..or.i==1) はtrueなので，上のコードはelseを吐き出す．
+前後の条件をorでつなぐ場合は，
+if((a.eqv..false.).or.(i==1)) としてカッコで順序を明示する．
+</nowiki>
+ちなみに
+  * (.ture..or..false.)は.true.（どちらかが真なのでtrue)
+  * (.ture..and..false.)は.false.（どちらも真でないのでfalse)
+になる．
+===== allocated =====
+すでにallocateしたかどうかわからない時は，allocatedで確かめてからallocateする．
+https://gcc.gnu.org/onlinedocs/gcc-4.1.0/gfortran/ALLOCATED.html
+<code fortran>
+  if (allocated(x).eqv..false.) allocate(x(n))
+</code>
+===== construct name =====
+入れ子のループで，cycleやexitを使って外側に飛びたいときは，do ループにラベルをつける．
+これで文番号を付けずに好きなところへ飛べる．
+  - do の前にラベルとコロンをつける．
+  - cycleやexitの後にラベルをつける．
+  - end do の後ろにラベルをつける．
+例
+<code fortran>
+outer : do i=1,n
+  do j=1,n
+    if(...) cycle outer
+  end do
+end do outer
+</code>
 ===== 読み込み専用，書き込み専用で開く =====
 読み込み専用で開く場合．
@@ ライン 46: / ライン 110: @@
-===== makefileでコンパイル =====
+===== makefileでコンパイル その２=====
+以下の記事は古い．[[fortran:makefile|別のページ]]を参照せよ．
+moduleファイルを含む場合，依存関係をはっきりさせる必要がある．
+依存関係を変えるたびにmakeファイルを書き変えるのは面倒なので，依存関係を見てくれる
+ f90_mod_deps.py を使う．
+http://lagrange.mechse.illinois.edu/mwest/f90_mod_deps/
+コンパイル不要のpardisoやompのモジュールを読み込んだり，同一ファイル内に複数のモジュールが存在し循環参照が起きる場合は，f90_mod_deps.py のwrite_depsに以下のように修正．
+<code python>
+def write_deps(outf, filename, deps, mods):
+    filebase, fileext = os.path.splitext(filename)
+#>>>
+    loadmod="MKL_PARDISO.mod"
+    if loadmod in deps:
+        deps.remove(loadmod)
+    loadmod="omp_lib.mod"
+    if loadmod in deps:
+        deps.remove(loadmod)
+#<<<
+    outf.write("%s.o: %s\n" % (filebase, " ".join(deps)))
+    for mod in mods:
+#>>>
+        if mod in deps:
+            deps.remove(mod)
+#<<<
+	outf.write("%s: %s\n"
+		   % (mod, " ".join(deps)))
+</code>
+以下，Makefile のサンプル．
+<code bash>
+TARGET = a.out
+DEPS=hoge.deps
+OBJECTS = piyo.o foo.o bar.o main.o
+F90 = gfortran
+FFLAGS=-Wall -ffpe-trap=invalid,zero,overflow,underflow -fbounds-check -fimplicit-none -fbacktrace -fopenmp
+#FFLAGS=-O3 -fopenmp
+#F90 = ifort
+#FFLAGS=-check all -warn all -std -gen_interfaces -fpe0 -ftrapuv -traceback
+#FFLAGS = -mkl -O3
+SOURCES=$(OBJECTS:.o=.f90)
+.SUFFIXES :
+.SUFFIXES : .o .f90 .f
+.f.o:
+	${F90} -c $< ${FFLAGS}
+.f90.o:
+	${F90} -c $< ${FFLAGS}
+.PHONY:all
+all: ${DEPS} ${TARGET}
+${DEPS}: ${SOURCES}
+	python f90_mod_deps.py --output $@ $^
+${TARGET} : ${OBJECTS}
+	${F90} -o $@ $^ ${FFLAGS}
+-include ${DEPS}
+%.o %.mod: %.f90
+	${F90} -c $< ${FFLAGS}
+	touch $(patsubst %.f90,%.mod,$<)
+clean :
+	rm ${TARGET} *.mod *.f90~ ${OBJECTS}
+</code>
+===== makefileでコンパイル =====
 複数のファイル（サブルーチンetc.）で構成されるプログラムを動かすのに便利．
 以下の文を「makefile」という名前で作り，コンパイルしたいファイルと同じディレクトリに置く．
@@ ライン 89: / ライン 224: @@
 </code>
-===== データファイルの読み込み =====
+===== 何行あるか分からないときのデータファイルの読み込み =====
 データファイルの中に何行あるか分からないときにはiostatかendを使う．
@@ ライン 211: / ライン 346: @@
 $ sh hoge.sh
 </code>

Applied Mechanics wiki

ユーザ用ツール

サイト用ツール

差分

ページ用ツール