好吧!我先自己回复下自己。在稀疏矩阵存储中,通用有两种格式,如下矩阵
A = [ a00 a01 0 0 0 a11 a12 0 a20 0 a22 0 a30 0 0 a33];引用一篇博士论文 CRS scheme stores the non-zero elements the matrix in left-to-right and top-to-bottom order in vector v (row-wise storge). 存储格式如下
v = [a00, a01, a11,a12,a20,a22,a30,a33] cI = [0,1,1,2,0,2,0,3] rI = [0,1,1,2,0,2,0,3]但是在openfoam中,不是这么存储的
openfoam是一种 LDU 格式 对角线 lowerAddr() upperAddr()
具体可参见lduAddressing.H 头文件
因为计算流体力学中,大部分矩阵都是对称(因为一个internal face对应两个体)除了特别的边界条件外,基本上都是对称矩阵。而矩阵对角元素就更好求了,大部分情况直接把该行相加就得到主对角元素。这种sparse matrix设计是符合cfd计算原理的,但是我的问题还是没解决。看到相关资料分享一下。
稀疏矩阵也分好多种,openfoam中的renumberMesh就是尽量保证在是多对角阵。
谢谢!