C++ 实现复数 GMRES 算法详解

复数的 GMRES 算法与实数 GMRES 算法基本相同，只是需要对复数进行适当的处理。以下是复数 GMRES 算法的基本步骤：

1. 定义复数向量和矩阵类型。可以使用 C++ 标准库中的 'std::complex' 类型表示复数。

2. 实现复数向量的加法、减法、数乘、内积等基本运算。

3. 实现复数向量的范数计算函数。

4. 实现复数矩阵向复数向量的乘法函数。

5. 实现复数 GMRES 算法的主函数。在每次迭代中，需要进行以下步骤：

   a. 计算残差向量 'r'。

   b. 对残差向量进行正交化，得到一个新的 Krylov 子空间。

   c. 在新的 Krylov 子空间中求解最小二乘问题，得到新的解向量 'x'。

   d. 计算新的残差向量 'r'，并判断是否达到收敛条件。

6. 在实现过程中，需要注意复数的运算规则和精度问题。例如，对于复数的判断相等，应该比较其实部和虚部的差的绝对值是否小于某个阈值。