使用CUSP打印矩阵的csr格式
Posted
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了使用CUSP打印矩阵的csr格式相关的知识,希望对你有一定的参考价值。
我正在尝试使用CUSP和CUDA将稀疏矩阵转换为csr格式。我稍微修改了CUSP文档中显示的代码:
#include <cusp/array2d.h>
#include <cusp/coo_matrix.h>
#include <cusp/csr_matrix.h>
#include <cusp/print.h>
int main(void)
{
// create a simple example
cusp::array2d<float, cusp::host_memory> A(3,4);
A(0,0) = 10; A(0,1) = 0; A(0,2) = 20; A(0,3) = 0;
A(1,0) = 0; A(1,1) = 30; A(1,2) = 0; A(1,3) = 40;
A(2,0) = 50; A(2,1) = 60; A(2,2) = 70; A(2,3) = 80;
print(A);
// save A to disk in MatrixMarket format
cusp::io::write_matrix_market_file(A, "A.mtx");
// load A from disk into a coo_matrix
cusp::csr_matrix<int, float, cusp::device_memory> B;
cusp::io::read_matrix_market_file(B, "A.mtx");
cusp::io::write_matrix_market_file(B,"B.mtx");
// print B
cusp::print(B);
return 0;
}
但我得到的结果是使用CU格式的矩阵,索引移动了1:
%%MatrixMarket matrix coordinate real general
3 4 8
1 1 10
1 3 20
2 2 30
2 4 40
3 1 50
3 2 60
3 3 70
3 4 80
有帮助吗?
谢谢。
答案
谢谢你的回复!事实上,我必须亲自编写csr格式,尽管它并不像我最初想的那么难。我的代码结果如下:
// Imprimimos la matriz en formato csr:
cout<<"
=======================
";
cout<<"Matriz de adyacencia en formato CSR:
";
cout<<"Row-offset :";
for (int j=0 ; j<10 ; j++) cout<<B.row_offsets[j]<<" ";
cout<<"
"; cout<<"Column index :";
for (int j=0 ; j<17 ; j++) cout<<B.column_indices[j]<<" ";
cout<<"
"; cout<<"Values:";
for (int j=0 ; j<17 ; j++) cout<<B.values[j]<<" ";
cout<<"
==================
";
以上是关于使用CUSP打印矩阵的csr格式的主要内容,如果未能解决你的问题,请参考以下文章
以压缩稀疏行格式 (csr_matrix) 对矩阵中的值取对数