迭代方法的美妙之处在于它们需要您做的只是矩阵向量乘法。在您的情况下，矩阵的乘积$A$带有矢量$y$可以写成$z=Ay = (B^TCB)y= B^T(C(By))$这表明您只需要三个矩阵向量乘积，但没有矩阵矩阵乘积。

如果$C$恰好是对称的和半正定的，那么也是$A$您可以使用共轭梯度方法。否则，您可能想考虑另一种迭代方法 GMRES。