2011-12-01から1ヶ月間の記事一覧
元々Sunのサイトにあったんですが、Oracleに買収されてリンクが移ったのでメモ。What Every Computer Scientist Should Know About Floating-Point Arithmetic. PDF直リン1 PDF直リン2浮動小数点演算について PDF直リン 付録D
組み込み関数 CMPLX() は REAL(4) がデフォールト値なので、引数が倍精度でも第二引数で KIND を倍精度に指定しないと、単精度で変換されます。 実行結果 ソース・プログラム program Complex implicit none complex (kind = 8) :: c real (kind = 4) :: a r…
FORTRAN66 では DO LOOP 終了後の DO LOOP 変数の値は不定ですが、FORTRAN77 以降は終了時の値です。
昔Householder法とQR法のサブルーチンを外人からもらいましたが、EISPACKの tred と tql で F66 で書かれていました。のちに Numerical Recipes の F77 版に置き換えたけれども、基本的に LAPACK のものだと思います。微妙に変えてあって、tql は旧版で収束…
クリントン(雌の方)がForeign Policy誌10月11日号に発表した論文“America's Pacific Century”がネットで話題だというので暇つぶしにななめ読みしたが、日米安保破棄も何も書いてないw このおっさんの記事『日米安保破棄を真剣に検討し始めた米国』はなん…
子言衛霊公之無道也、康子曰、夫如是、奚而不喪。孔子曰、仲叔圉治賓客、祝鴕治宗廟、王孫賈治軍旅。夫如是、奚其喪。[通釈]孔子が衛の霊公の無道なことを語られたので、季康子が問うて曰うには、「お話のように無道なのに、どうして位を失わないのですか。…
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…
実行結果 ソース・プログラム 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…