C++实现共轭梯度法求解对称正定稀疏线性系统

共轭梯度法（Conjugate Gradient Method）是一种迭代算法，用于求解对称正定线性系统。这种方法特别适用于处理稀疏矩阵，因为它只需要存储非零元素，而不需要存储整个矩阵。在本文中，我们将介绍如何使用C++实现共轭梯度法来求解对称正定稀疏线性系统。

首先，我们需要定义问题的输入：对称正定矩阵A、向量b以及初始解x。在这里，我们将使用稀疏矩阵存储格式来表示矩阵A，比如使用压缩列（Compressed Column）格式。这种格式只存储非零元素及其对应的行号，从而减少内存占用。

下面是一个示例的C++代码实现：

#include <iostream>
#include <vector>

// 定义稀疏矩阵的压缩列格式
struct SparseMatrix {