给定矩阵A，我有一个以以下格式存储的大矩阵；
一个=
1 0 3
5 1 -2
0 0 7
3个向量；
NVPN = [1 3 4 7] - 我在第一列中任意放了一个 1，然后从第二列开始，它是每列非零元素数量的累积总和。
NNVI = [1 2 2 1 2 3] - 每个非零元素的行索引。
CONT = [1 5 1 3 -2 7] - 每个非零元素的值。
我现在需要执行矩阵*矩阵乘法和矩阵*向量乘法。有谁知道是否有任何 FORTRAN 库，我可以修改以适应我的问题，在上面执行此操作？
提前致谢