fortran66のブログ

fortran について書きます。

2014-03-01から1ヶ月間の記事一覧

Modern Fortran Explained の Errata の新版が出てました。

PDFに直リン ftp://ftp.numerical.rl.ac.uk/pub/MRandC/edits.pdf 親ディレクトリ ftp://ftp.numerical.rl.ac.uk/pub/MRandC/

FORTRANステッカー

http://www.amazon.com/CafePress-FORTRAN-Bumper-Sticker-Standard/dp/B00FOMUXMO

Fortran information & resources @ Fortranplus

Fortranplus のサイトにある情報まとめが小さいながらツボを押さえている気がします。 http://www.fortranplus.co.uk/fortran_info.html この文書も最近の Fortran 規格の差分情報の胆が要領よくまとめられていると思います。 http://www.fortranplus.co.uk/…

Fibonacci number 

fibonacci 数は漸化式の形で与えられているので、逐次的にしか計算できないように見えますが、漸化式を行列で書くと、行列のべき乗を求めることで非逐次的に解くことが可能です。 実行結果 行列法と漸化式を逐次求めたものの結果を比較してみます。当然同じ…

JAMSTEC の Fortran90 入門講座

テストに間違いが多すぎw http://www.jamstec.go.jp/esc/simschool/f90learning/index.html

Python で有理数 6j symbol

メモ帳 とりあえず定義式通りに。計算量が少ない変形もあるんですが・・・ import sympy import numpy def fa(n): if n < 0: return 0 else: return math.factorial(n) def dicsort(d): return sorted(d.items(), key=lambda x: x[0]) def pr(i): a = sympy.…

Python で 3j-symbol

Python は確かに楽に有理数・多倍長計算ができる。 フォーマット文にあたるものがよくわからなくて整形出力できない・・・IPython では TeX 表記もできる様だから、数学記号的に出力させられるのではないかとも思うのだが、よく分からんw In [9]: chk3(4,4,…

Python で有理数計算

よくわからなくて苦しいw ライブラリの関係も謎すぎる、どの程度信じていいのか?w ThreeJ symbol の m=0 の時でチェック。 In [190]: chk(10,6,6) -0.0884757071422 -6*sqrt(2028117)/96577 ({2: 2, 3: 3, 7: 1}, {17: 1, 19: 1, 13: 1, 23: 1}) Python …

Ruby から Fortran

メモ帳 Ruby全く知らないw インストール Ruby 64bit は何故かエラーが出て動かないので、32bit版で。DevKit もいれると GFortran が使える。 https://www.ruby-lang.org/ja/installation/#rubyinstaller gem install ffi で何とかなる。 Fortran DLL を作る…

Fortran の変数に Python からアクセスする。

メモ帳 参考サイト http://docs.python.jp/contrib/ctypes/tutorial_jp.html http://slashdot.jp/journal/445869/python-ctypes%E4%BD%BF%E3%81%A3%E3%81%A6%E3%81%BF%E3%81%9F%E3%80%82 http://docs.python.jp/2/library/ctypes.html Fortran ソース・プロ…