fortran66のブログ

fortran について書きます。

2014-01-01から1年間の記事一覧

最初のBug リレー中の蛾

Stroustrupの講義スライドから #22 Programming - Lecture slides http://www.stroustrup.com/Programming/lecture-slides.html First actual case of bug being found.

FortranでAA白黒 改

入力フォント データ例 はてなのアップローダの仕様でjpgに変換されていますが、本来は24bit bmpです。 font4x8.bmp font8x8.bmp module m_bmp implicit none type :: t_bmp_file_header sequence integer(2) :: bfType = transfer('BM', 0_2, 1) ! BitMap i…

Fortran で ASCII ART カラー版

少し Fortran2003/08 の機能を試してみることにして、やっつけ仕事でカラーの BMP も出力するようにしてみました。1バイト整数問題は、変数をキャラクターにして、iachar で4バイト整数に直すようにしてみました。これで unsigned int(1) として扱えます。…

Fortran で ASCII ART

Fortran で BMP -> ASCII 変換のプログラムを作ります。基本的な原理は、フォントのビットマップを(非直行、非完全、非独立)基底ベクトルとして、画像をフォントのサイズごとに細切れにして、その小さなビットマップをベクトルとみなして、各フォント基底…

SX-ACE 続報

2.NECが東北大などからSX-ACEスパコンを受注 http://www.geocities.jp/andosprocinfo/wadai14/20140531.htm NEC、ベクトル型スーパーコンピュータ「SX-ACE」を東北大学、大阪大学、国立環境研究所から受注 http://jpn.nec.com/press/201405/20140529_03.ht…

OpenCL kernel in Fortran

Intel CPU 内蔵グラフィックスアクセレレータで用 OpenCL が Intel Fortran でも使えるようです。詳細は不明。https://software.intel.com/en-us/forums/topic/515753

Modern Fortran Reference Card

Fortran90 reference card の作者ゲルツ氏が Fortran2003/08 も含んだ改訂版を出してくれました。http://michaelgoerz.net/refcards/直リン http://michaelgoerz.net/refcards/fortran_refcard_a4.pdf

昔のイタリアのセンス良すぎw

Mina _ Bossa nova medley _ live 1968 - YouTube ELIS REGINA (MEDLEY) ITALIA 1972 - RAI ...

乞食速報 Mac 用 PGI Fortran Compiler 無料!

Free PGI — HPC Compilers for Apple OS Xhttp://www.pgroup.com/products/freepgi/index.htmFortran2003完全対応 & c99 コンパイラ付き! 乞食マカーの皆さん急げ!

SX-ACE

東北大と阪大の計算センターの次期機種が NEC SX-ACE に正式に決まったようです。https://www.hpci-office.jp/pages/isc_tohoku_2014-2 https://www.hpci-office.jp/pages/cmc_osaka-u_2014-5東北大は最大構成クラスターを二台、阪大はレーザー研と共同で最…

GTK3でGUI

GTK3を用いて Fortran で GUI をやる方法やツール開発を紹介しております。 http://fortrandev.wordpress.com/

Interactive IPython

IPython が 2.0 になって新機能が増えたようです。スライダーがほぼ何もせずに付きます。これは楽だわ。 subroutine cmandel(c, m) integer, parameter :: maxiter = 255 integer, parameter :: kd = kind(1.0d0) complex(kd), intent(in ) :: c integer , i…

Damian Rouson によるビデオ講義

麻生元首相も通ったスタンフォード大学が Fortran ビデオ講座を開講したようです。講師は、Fortran でのオブジェクト指向本を書いている方で、そちら方向の講義のようです。http://jolts.stanford.edu/72/modern_fortran_for_scientists_and_engineers目次 S…

メモ帳 parameterized derived types 例

Modern Fortran Explained 13章 を参考に試してみました。ソース・プログラム program test implicit none type :: t_mat(ikind, n, m) integer, kind :: ikind integer, len :: n, m real(ikind) :: x(n, m) end type t_mat ! type, extends(t_mat) :: t_la…

The new features of Fortran 2008 改定版

The new features of Fortran 2008 改定版が出ていました。The new features of Fortran 2008 (Reid) - supersedes N1828 ftp://ftp.nag.co.uk/sc22wg5/N1851-N1900/N1891.pdfちょっと見には、どこが改定されたのかよくわかりません。

ACM SIGPLAN Fortran Forum 4月号

イギリス Fortran 同人誌最新号出ました。ACM SIGPLAN Fortran Forum Volume 33 Issue 1, April 2014http://dl.acm.org/citation.cfm?id=2594488&CFID=304828140&CFTOKEN=91798490今月は Arjen Markus による書評と記事の2本立てでした。

急げネットイナゴ! Beta for Intel Developer Tools 2015 始まる。

Intel Fortran 次期バージョンのベータテストが始まりました。 登録すると今年の9月後半まで気前よく全部入りのベータ版が試せるようです。Fortran に関していえば Parameterized Derived Types と block..end block 構造の導入と が目玉でしょうか?Paramet…

中央二項係数

Pascal の三角形によって二項係数が求まります。その中で *1 にあたる要素だけを取り出します。 これを中央二項係数と呼びます。Pascal の三角形の中央に並んだ列に当たります。次にこの中央二項係数を用いて、級数 http://oeis.org/A086463 を計算してみま…

あいまいみー二期はマジ!あびゃー!!

コミックナタリー - 「あいまいみー」2期決定!エイプリルフールの予告はガチ あいまいみー / ちょぼらうにょぽみ / まんがライフWIN

無駄話

最近 amazon の洋書の値段が高いような気がします。以前は米 amazon + 送料と同程度だったのが、最近は送料払っても米 amazon が安いような。でも米 amazon で重い本頼むと箱が壊れて、本の角がつぶれて背中が割れて送られてきたりするから困ります。 あとこ…

組み合わせ

1~m迄の数字をn個重複しないように並べる。 実行結果 m=6、n=3 1 2 3 1 2 4 1 2 5 1 2 6 1 3 4 1 3 5 1 3 6 1 4 5 1 4 6 1 5 6 2 3 4 2 3 5 2 3 6 2 4 5 2 4 6 2 5 6 3 4 5 3 4 6 3 5 6 4 5 6 Press any key to continue . . . Fortran2003 版 mod…

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 を作る…