Just wanted to check before I modified my code. If my jacobian matrices are in fact sparse, will any of the nonlinear solvers be able to take advantage of that structure if I were to build up the jacobians as a sparse data structure? -gideon