@wds15 yep, that’s better.
@charlesm93 Since KINSetMaxSetupCalls is set to 1, it’s exact Newton where Jacobian gets updated every step. How’s that going to affect the performance of a large system? Will modified newton achieve the same, say, for Laplace approximation, with better efficiency?