nonlinear solver works with NLN_GS, derivative computation stuck at 1 with GMRES