fortran66のブログ

fortran について書きます。

2012-06-01から1ヶ月間の記事一覧

面白いバグ

http://www.fortran.bcs.org/2007/jubileeprog.php Bugs I Have Known and Loved - Ron Bell, AWE Aldermaston http://www.fortran.bcs.org/2007/jubilee/bugs.pdf 常数も変数のようにメモリー上に置かれ、かつ最適化すると共有されてしまいます。ところが、…

forall 命令のダメっぷりについて

Implementing the Standards - including Fortran 2003 Malcolm Cohen, Numerical Algorithms Group Ltd http://www.fortran.bcs.org/2007/jubileeprog.php http://www.fortran.bcs.org/2007/jubilee/f50.pdf こんなダメダメな forallちゃんですが、取り柄も…

Intel® Visual Fortran Composer XE 2011 for Windows* Update 11, Version 2011.11.344 のお知らせが来ました。

Fortran2003 機能の拡張は無いようです。update では対応せず、次の Version で一気に残りを片付けるのではないかと思います。最近では、CRAY、IBM、PGI が Fortran2003 完全対応をうたっているので、そろそろ Intel や NAG 先生にも頑張っていただきたいと…

Features deleted in Fortran

Features deleted in Fortran95 Non-integer do indices Assigned go to and assigned formats Branching to an end if statement The pause statement H edit descriptor Features deleted in Fortran2003 Carriage control 消すには惜しい人材たちw

ソース・コード

module m_mandel implicit none integer, parameter :: kd = kind(0.0d0) integer, parameter :: maxiter = 50 contains pure elemental integer function imandel(z) complex(kd), intent(in) :: z complex(kd) :: c c = (0.0_kd, 0.0_kd) do imandel = max…

実行結果

80x24 136x45

キャラクター版のマンデルブロ図形

ネットでよく見るキャラクター版マンデルブロ図形を Fortran で描いてみました。参考:ASCII Haskell Mandelbrot http://yannesposito.com/Scratch/en/blog/Haskell-Mandelbrot/intel visual fortran の visual studio のアウトライン表示は、forall 1行の…

N1921 Special mathematical functions in Fortran (supersedes N1688)

将来 Fortran 規格に加わるかもしれない特殊関数のリストが出ています。 ftp://ftp.nag.co.uk/sc22wg5/N1901-N1950/N1921.pdf ftp://ftp.nag.co.uk/sc22wg5/N1651-N1700/N1688.pdf オイラーの常数、ベッセル関数、ルジャンドル関数等、物理数学の教科書に出…

Fortran 2003+2008 概略

規格に直接関わっている高田正之氏によるまとめ。 http://www.itscj.ipsj.or.jp/forum/forum2008SC22fortran/Fortran20032008.pdf

添え字変数を使わない方向で素因数分解

DO LOOP のループ変数をできるだけ使わない、あるいは配列の添え字変数を使わない方向で、有限個の素数のテーブルにしたがって素因数分解し、テーブルにない素数は余りとして別変数に積を取ってためるプログラムを書いてみます。 試しなので、さほど一貫した…

open 文に於ける newunit

Fortran では外部ファイルの unit 番号の決定が問題になります。そもそも unit 番号が何番まで使えるのか、何番がシステムに予約されているのかすら、規格に定められていないためによくわかりません。 Fortran2008 では、その問題を解決するために、open 文…

ワロス

天安門事件の日付が上海株指数に…検索制限騒ぎ【香港=竹内誠一郎】上海株式市場で4日、市場全体の値動きを示す総合指数の下落幅が、中国の民主化運動が武力弾圧された天安門事件の日付1989年6月4日と同じ64・89ポイントとなった。 事件からちょ…