Accelerated Sparse Matrix Operations in Nonlinear Least Squares Solvers