稀疏矩阵的保存和右乘向量,简单的实现。同样暂时不关心速度/效率。
为了能够方便地用 A*b 的操作,矩阵使用Compressed Row Storage(CRS)格式或者称为 Compressed sparse row(CSR)格式储存。
使用Eigen的结果作为比对测试。
程序需要读取文件格式说明如下,如果需要diy测试可以自行准备:
矩阵数据格式
spA: 对称稀疏方阵
rows cols
r0 c0 v0
r1 c1 v1
...
rn cn vn
b: 列向量
v1
v2
...
vn