对称矩阵下三角的储存方式

frankgty

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