Add tuned CSCMatrix<Scalar>::multiply_with_vector.
Work on nonlinear solvers.