Line | Count | Source (jump to first uncovered line) |
1 | #include "cs.h" | |
2 | /* y = A*x+y */ | |
3 | CS_INT cs_gaxpy (const cs *A, const CS_ENTRY *x, CS_ENTRY *y) | |
4 | 0 | { |
5 | 0 | CS_INT p, j, n, *Ap, *Ai ; |
6 | 0 | CS_ENTRY *Ax ; |
7 | 0 | if (!CS_CSC (A) || !x || !y) return (0) ; /* check inputs */ |
8 | 0 | n = A->n ; Ap = A->p ; Ai = A->i ; Ax = A->x ; |
9 | 0 | for (j = 0 ; j < n ; j++) |
10 | 0 | { |
11 | 0 | for (p = Ap [j] ; p < Ap [j+1] ; p++) |
12 | 0 | { |
13 | 0 | y [Ai [p]] += Ax [p] * x [j] ; |
14 | 0 | } |
15 | 0 | } |
16 | 0 | return (1) ; |
17 | 0 | } |