我正在做无网格粒子法MPS的程序

logxing

先不求效率，第一个目标能计算就行。
影响效率的一个是临近粒子的搜索，一个是线性方程组的求解（大型稀疏矩阵）。
欢迎一起讨论下。

logxing

今天集中精力考虑解线性方程组的方法，都说是用ICCG方法，还没怎么看懂。
不过要处理的是100万个粒子，也就是100万阶的稀疏矩阵，每行非0元素得有10到150左右，先别说计算，连内存都不够啊。