0

我想在 Fortran95 中使用 Sparse Blas 来创建矩阵,并且我正在使用点入口构造。使用命令创建矩阵后

调用 duscr_begin(n,n,a,istat)

这里 a 是 n 乘以 n 的矩阵的句柄。在其中插入值后,如何使用其句柄 a 查看最终矩阵?由于我想将矩阵用于其他一些操作,所以我想以三个向量(稀疏)形式(row_index、Col_index、Value)查看矩阵。

关于这个 Sparse Blas 的详细信息在第 3 章中给出,可以在这里看到

http://www.netlib.org/blas/blast-forum/

4

1 回答 1

0

实际上我问的是在 16 天之前,这不仅仅是在屏幕上写一个变量。我正在使用一些称为 Sparse Blas 的库来创建稀疏矩阵。后来通过深入研究图书馆,我找到了解决我的问题的方法,即使用句柄我们如何获得三个向量行、col 和 Val。命令类似于

call accessdata_dsp(mat,a_handle,ierr)
call get_infoa(mat%INFOA,'n',nnz,ierr)
allocate(K0_row(nnz),K0_col(nnz),K0_A(nnz))
K0_row=mat%IA1; K0_col=mat%IA2; K0_A=mat%A

所以这里 nnz 是稀疏矩阵中的非零项,而 K0_row、K0_col 和 K0_A 是我们需要的三个向量,可以在进一步计算中使用。

于 2012-06-15T19:05:59.400 回答