Fortran 95中的稀疏Blas

时间:2012-05-29 23:21:24

标签: fortran sparse-matrix fortran90 blas finite-element-analysis

我想在Fortran95中使用Sparse Blas只是为了创建矩阵而我正在使用点入口构造。使用命令

创建矩阵后

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

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

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

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

1 个答案:

答案 0 :(得分:0)

实际上我所问的是在16天之前,它不只是在屏幕上写一个变量。我正在使用一些称为Sparse Blas的库来创建稀疏矩阵。后来通过挖掘到库我找到了我的问题的解决方案,使用句柄我们如何得到三个向量row,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是我们所需的三个向量,可以用于进一步的计算。