% compute yCSC = A*x with A in CSC format function yCSC = SpMV_CSC(valACSC,rowindexACSC,colbeginACSC,n,x); % yCSC = zeros(n,1); % allocate space for answer for j = 1:length(colbeginACSC)-1, for i=colbeginACSC(j):colbeginACSC(j+1)-1, yCSC(rowindexACSC(i)) = yCSC(rowindexACSC(i)) + valACSC(i)*x(j); end end