稀疏矩阵的转置_基于动力自由度缩聚的模型修正（二）稀疏矩阵的Harwell-Boeing压缩储存格式及重构...

由于个人对ANSYS采用的Harwell-Boeing稀疏矩阵储存格式不了解，特地查阅了相关资料进行学习。紧接着上一篇文章，整理了关于该储存方法的资料，以便后续查用。

---

一、稀疏矩阵的存储需求

由于大型稀疏矩阵仅含有少量的非零元素，有效信息的密度低。相比之下，数量庞大的零元素在计算机储存时会占用大量不必要的存储空间。为避免这一问题，仅对非零元素进行存储和操作，可有效降低内存开销，同时提高后续计算的效率。

二、稀疏矩阵的存储方法

1. COO格式

基于非零元素的储存思想，首先考虑采用三组一维数组对稀疏矩阵非零元素的行指标、列指标以及数值分别进行存放，通过行列坐标来描述元素位置，即常规的Coordinate Format (COO)。该方法对矩阵元素按顺序存放，便于依照顺序处理的矩阵运算。

COO格式储存示意

采用matlab生成图示的稀疏矩阵mat的COO格式:

spamat = [3,0,4,0,0;0,2,0,0,0;0,0,0,0,0;0,0,8,3,5;0,0,5,0,0;0,0,0,6,3]; % 生成示例矩阵
coo.numelem = length(find(spamat~=0)); % 获取非零元素数目
[coo.numrow,coo.numcol] = s