矩阵求解
gunslinger
2003年本科,2006年研究生毕业于华中科技大学(原华中理工大学)材料学院。
展开
-
矩阵求解泛谈
矩阵求解AX = b 在工程上使用相当广泛,尤其是在simulation领域。矩阵求解的速度相当程度地决定一款模拟软件的命运。矩阵求解的方法很多,如一般的迭代法,投影法,AMG等。 1)迭代法 不实用,尤其是求解大型矩阵。但可以作为投影法的prediction算子存在。如CGJacobi, CGSSOR,这是把Jacobi, SSOR方法当作CG原创 2011-10-27 21:57:56 · 712 阅读 · 0 评论 -
sparse matrix 存储
由于sparse matrix的特殊性,在存储的时候,一般情况下吃用CSR(Compressed Sparse Row)这种format来实现。通过使用3个数组来描述一个Sparse Matrix. double* mValues; // 保存Sparse Matrix 非0的值。 int* mColumnIds; //保存column ID int* mRowIds; //保存各行第一个非原创 2011-10-27 21:29:56 · 641 阅读 · 0 评论 -
LU decomposition
LU分解,是一种比较有用的求解线性矩阵方程AX = b的方法。当然,它只是在未知数量有限的情况下才能发挥它的优势。怎么求解A = LU呢,方法挺多,这里介绍一种常用的方法。直接上代码: for(int i = 0; i { for(int j = i + 1; j { Aji = Aji / Aii;原创 2011-11-01 21:57:36 · 825 阅读 · 0 评论