Passing Additional Arguments To The Jacobian In Scipy.optimize.root

I am solving large systems of equations, but have access to the jacobian matrix. I decided to uses scipy.optimize.root (open to other suggestions), where I can pass the jacobian as

Solution 1:

Use root(..., args=(value,)) to find the root of fun(x,value)==0 at fixed value. Note that args must be a tuple even if a single-element one.