fortran66のブログ

fortran について書きます。

2011-12-02から1日間の記事一覧

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…