本发明公开了一种带多级Cache的SIMD众核处理器上的稀疏矩阵存储方法，包括：(1)获取矩阵A中行非零元个数最大值a，处理器SIMD单元可同时计算的非零元个数b，计算大于a且为b的倍数的最小值作为临时行宽度；(2)对矩阵A，数组Value和Colidx分别顺序存储每一行的非零元值和列坐标，每行非零元个数不足临时行宽度的，在该行尾分别补0和‑1；(3)对Colidx和Value按b列分块；(4)每个列块按行压缩，使列块中有非零元的行集中于列块上部；(6)列块按b行分块，得到子块；(7)去除全零子块，将子块按行存储。本发明的存储方法将稀疏矩阵分成稠密子块，在保留非零元局部性的基础上，提高了对处理器SIMD处理单元以及寄存器的利用率，提高了稀疏矩阵向量乘的性能。
