Извините, что разочаровал, но, насколько я могу судить, я боюсь, что ответ - нет.Я много работал с разреженными данными и не нашел глубоких библиотек (например, уровня BLAS) для манипуляций с разреженными матрицами.На более высоком уровне, чем BLAS, существует обширная поддержка в Matlab, R, Python и других языках.
Из того, что я видел, это возникает из-за различий в типах разреженных матриц (разбросанных,симметричный, полосатый или трехдиагональный, блочная диагональ), их содержимое (двоичное, целое, вещественное) и их математические свойства (например, положительно определенное - не гарантировано, полный ранг - не распространено) и т. д. имеет тенденцию усложнять оптимизацию.
Вместо этого я, как правило, пишу свой собственный код для обработки разреженных матриц, повторно используя все, что я могу для хранения и простых вычислений, таких как умножение.
Со временем появится хорошая библиотека низкого уровня,но я еще не видел это.