43 void scscmv (
int isym,
int indexbase,
int m,
int n,
44 double *val,
int *indx,
int *pntr,
63 for (i = 0; i < m; i++)
68 for (j = 0; j <n ; j++)
70 jbgn = pntr[j]-indexbase;
71 jend = pntr[j + 1]-indexbase;
73 for (i = jbgn; i < jend; i++)
75 y[indx[i]-indexbase] += val[i] * x[j];
76 if (indx[i]-indexbase != j && isym) y[j] += val[i]*x[indx[i]-indexbase];
void scscmv(int isym, int indexbase, int m, int n, double *val, int *indx, int *pntr, double *x, double *b)