Sparse Example Matrix-Vector Multiply

#include "spblas.h"

/* sparse matrix in coordinate format */

const int N = 4;

const int nnz = 6;

double val[nnz] = { 1.1, 2.2, 2.4, 3.3, 4.1, 4.4 };

double indx[nnz] = { 1, 2, 2, 3, 4, 4};

double jndx[nnz] = { 1, 2, 4, 3, 1, 4};

double x[N] = { 1.0, 1.0, 1.0, 1.0 };

double y[N] = { 0.0, 0.0, 0.0, 0.0 };

int A;

int istat = 0;

double alpha = 0.0;

/* 1: create sparse matrix handle */

c_duscr_coo( &A, N, N, val, indx, jndx, nnz, blas_zero_base, &istat );

/* 2: y = A*x */

c_dusmv( blas_no_trans, alpha, A, x, 1, y, 1 );

/* 3: release matrix handle */

c_usds(A);

Previous slide

Next slide

Back to first slide

View graphic version