19 size_t A_num_rows,
size_t ,
size_t ,
20 const I* RESTRICT A_pos ,
const I* RESTRICT A_idx,
const V* RESTRICT A_val,
26 #pragma omp for nowait
27 for(
int i = 0; i < (int)A_num_rows; i++)
29 for (
int jj = A_pos[i]; jj < A_pos[i+1]; jj++)
32 y_ptr[i] = DG_FMA(
alpha*A_val[jj], x_ptr[j], y_ptr[i]);
38 #pragma omp for nowait
39 for(
int i = 0; i < (int)A_num_rows; i++)
42 for (
int jj = A_pos[i]; jj < A_pos[i+1]; jj++)
45 temp = DG_FMA(
alpha*A_val[jj], x_ptr[j], temp);
48 y_ptr[i] = DG_FMA(
beta, y_ptr[i], temp);
void spmv_omp_kernel(CSRCache_omp &, size_t A_num_rows, size_t, size_t, const I *RESTRICT A_pos, const I *RESTRICT A_idx, const V *RESTRICT A_val, value_type alpha, value_type beta, const C1 *RESTRICT x_ptr, C2 *RESTRICT y_ptr)
Definition sparsematrix_omp.h:17