pan>;
173  }
174 
175  for (IT j = 0; j< A.n; ++j)
176  {
177  for(IT k= A.jc[j]; k < A.jc[j+1]; ++k)
178  {
179  for(int i=0; i<D; ++i)
180  {
181  y[j][i] += A.num[k] * x[A.ir[k]][i];
182  }
183  }
184  }
185 }
186 
187 
188 #include "csc.cpp" // Template member function definitions need to be known to the compiler
189 #endif
190 
Csc::getir
ITYPE * getir() const
Definition: csc.h:48
csc_gaxpy_mm
void csc_gaxpy_mm(const Csc< NT, IT > &A, array< NT, D > *x, array< NT, D > *y)
Definition: csc.h:142
BmCsb
Definition: bmcsb.h:21
csc_gaxpy_mm_trans
void csc_gaxpy_mm_trans(const Csc< NT, IT > &A, array< NT, D > *x, array< NT, D > *y)
Definition: csc.h:166
Csc::rowsize
ITYPE rowsize() const
Definition: csc.h:46
Csc::Csc
Csc()
Definition: csc.h:18
triple.h
csc.cpp
Csc::getnum
T * getnum() const
Definition: csc.h:49
BmSym
Definition: bmsym.h:50
CsbSym
Definition: csbsym.h:43
Csc::getlogicalnnz
ITYPE getlogicalnnz() const
Definition: csc.h:50
Triple
Definition: csc.h:12