fortran66のブログ

fortran について書きます。

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

密行列と粗行列のCG法を単一のルーチンで記述

密行列と粗行列のCG法(共役勾配法)による線型方程式の解法を、オブジェクト指向機能を用いて単一のルーチンで記述してみました。行列とベクトルの積をそれぞれの行列形式について定義し、それを「*」演算子に総称化させることで可能になります。 abstract…