fortran66のブログ

fortran について書きます。

2011-12-01から1ヶ月間の記事一覧

DAVID GOLDBERG 『浮動小数点について知っておきたいこと。』

元々Sunのサイトにあったんですが、Oracleに買収されてリンクが移ったのでメモ。What Every Computer Scientist Should Know About Floating-Point Arithmetic. PDF直リン1 PDF直リン2浮動小数点演算について PDF直リン 付録D

CMPLX の KIND 指定

組み込み関数 CMPLX() は REAL(4) がデフォールト値なので、引数が倍精度でも第二引数で KIND を倍精度に指定しないと、単精度で変換されます。 実行結果 ソース・プログラム program Complex implicit none complex (kind = 8) :: c real (kind = 4) :: a r…

DO LOOP 終了後の DO LOOP 変数の値

FORTRAN66 では DO LOOP 終了後の DO LOOP 変数の値は不定ですが、FORTRAN77 以降は終了時の値です。

Householder法による三重対角化

昔Householder法とQR法のサブルーチンを外人からもらいましたが、EISPACKの tred と tql で F66 で書かれていました。のちに Numerical Recipes の F77 版に置き換えたけれども、基本的に LAPACK のものだと思います。微妙に変えてあって、tql は旧版で収束…

クリントン国務長官の論文

クリントン(雌の方)がForeign Policy誌10月11日号に発表した論文“America's Pacific Century”がネットで話題だというので暇つぶしにななめ読みしたが、日米安保破棄も何も書いてないw このおっさんの記事『日米安保破棄を真剣に検討し始めた米国』はなん…

講談社学術文庫版 宇野哲人『論語新釈』432頁 憲問第十四の二十

子言衛霊公之無道也、康子曰、夫如是、奚而不喪。孔子曰、仲叔圉治賓客、祝鴕治宗廟、王孫賈治軍旅。夫如是、奚其喪。[通釈]孔子が衛の霊公の無道なことを語られたので、季康子が問うて曰うには、「お話のように無道なのに、どうして位を失わないのですか。…

LU分解 Crout法 部分Pivot付き

Pivot選択のためにちょっと余分な計算をしなければなりません。 実行結果 ソース・プログラム module m_mat implicit none integer, parameter :: kd = kind(0.0d0) contains subroutine lu2(a, b, x) real(kd), intent(in) :: a(:, :), b(:) real(kd), inte…

LU分解 Crout法

実行結果 ソース・プログラム module m_mat implicit none integer, parameter :: kd = kind(0.0d0) contains subroutine lu2(a, b, x) real(kd), intent(in) :: a(:, :), b(:) real(kd), intent(out) :: x(:) real(kd), allocatable :: c(:, :) c = a ! F20…