But if don't convert `b` to CSC, `b[:, j].toarray()` will fail since it doesn't have the `toarray` member. So we may need extra code to deal with such a situation. In my opinion, `spsolve` is designed for sparse matrices so we may not need to optimize it for non-sparse matrices.