fortran66のブログ

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