【連立一次方程式】安定化双共役勾配法(Bi-CGSTAB法) C++コード付き

今回は正定値対称行列にしか適用できない共役勾配法(Conjugate Gradient法)をパワーアップした手法である双共役勾配法(Bi-Conjugate Gradient法)を安定化した、安定化双共役勾配法(Bi-conjugate gradient stabilized法)のC++コードを公開します。ややこしいですね。つまりCG法→Bi-CG法→Bi-CGSTAB法のような流れです…