I do not understand the Output of gen_coeffs.gen_two_diode, which is:
fjac: array([[nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,nan, nan,