pan>;
175 for (IT j = 0; j< A.n; ++j)
177 for(IT k= A.jc[j]; k < A.jc[j+1]; ++k)
179 for(
int i=0; i<D; ++i)
181 y[j][i] += A.num[k] * x[A.ir[k]][i];
void csc_gaxpy_mm(const Csc< NT, IT > &A, array< NT, D > *x, array< NT, D > *y)
void csc_gaxpy_mm_trans(const Csc< NT, IT > &A, array< NT, D > *x, array< NT, D > *y)
ITYPE getlogicalnnz() const