Use CSC matrices