Fortran Coder

标题: 对称矩阵下三角的储存方式 [打印本页]

作者: frankgty 时间: 2019-11-16 23:46
标题: 对称矩阵下三角的储存方式
有一个对称的稀疏矩阵，在mkl中进行很多操作时（如矩阵向量相乘，稀疏储存方式转换等) 会遇到选择输入矩阵是uplo的情况，此时若想仅保存矩阵下三角，应该以何种方式对它进行储存? 是将下三角元素以一维数组方式储存呢，还是以二维数组（第一行1个元素，第二行2个元素......第n行n个元素）进行储存? 如果是用二维数组储存的话，输出的循环该怎么写?

作者: vvt 时间: 2019-11-17 09:12
MKL提供构建稀疏矩阵的函数，用户不需要自行构建，也不必关心内部构建的结构细节问题。

稀疏矩阵存储有很多格式，诸如 DOK，LIL，CSC，CSR 什么的，有兴趣可以了解~

欢迎光临 Fortran Coder (http://bbs.fcode.cn/)