稀疏BLAS支持
稀疏基本线性代数子程序(BLAS)定义了向量和稀疏矩阵的一组基本操作，可用于创建优化的高级线性代数功能。GSL为稀疏矩阵支持有限数量的BLAS操作。
头文件gsl_spblas.h包含稀疏BLAS函数的原型和相关声明。
44.1 稀疏BLAS操作
int gsl_spblas_dgemv(const CBLAS_TRANSPOSE_t TransA, const double alpha,
const gsl_spmatrix * A, const gsl_vector * x, const double beta,
gsl_vector * y)
本函数计算矩阵向量乘积和y←αop(A)x+βy，其中op(A)=A，AT表示TransA = CblasNoTrans, cblastans。不支持计算，因此x和y必须是不同的向量。矩阵A可以是三元组或压缩格式。
int gsl_spblas_dgemm(const double alpha, const gsl_spmatrix * A,
const gsl_spmatrix * B, gsl_spmatrix * C)
本函数计算稀疏矩阵-矩阵乘积C=αAB。矩阵必须是压缩格式的。
44.2 参考文献和深入阅读
这些函数使用的算法，可以在以下来源找到描述：
• Davis, T. A., Direct Methods for Sparse Linear Systems, SIAM, 2006.
• CSparse software library, https://www.cise.ufl.edu/research/sparse/CSparse