間違えました。

ようやく有限要素法のソルバーができてデバッグを始めたら、やってはいけないバグに気づいた。係数行列か既知のベクトルがおかしいみたい。で現在書き直し中。