[SciPy-User] optimize.leastsq does not converge with full Jacobian?

Sebastian Walter sebastian.walter@gmail....
Fri Oct 2 03:10:32 CDT 2009


I've added your example to the unit test of pyadolc:
http://github.com/b45ch1/pyadolc/blob/master/tests/complicated_tests.py
starting at  line 248.
Let me know if you have objections.

On Thu, Oct 1, 2009 at 3:01 PM, Sturla Molden <sturla@molden.no> wrote:
> Pauli Virtanen skrev:
>> def jacobian(t, x, y):
>>>     j = np.zeros((2,x.shape[0]))
>>>     Vmax, Km = t
>>>     j[0,:] = x/(Km + x)
>>>     j[1,:] = -Vmax*x/((Km + x)**2)
>>>     return j
>>>
>>
>> Sign error.
>>
>>
>
> Thanks Pauli.
>
> I differentiated the Michaelis-Menten function instead of the residuals.
> That goes for "incredibly stupid" then.
>
> Sturla
>
> _______________________________________________
> SciPy-User mailing list
> SciPy-User@scipy.org
> http://mail.scipy.org/mailman/listinfo/scipy-user
>


More information about the SciPy-User mailing list