[SciPy-User] optimize.leastsq does not converge with full Jacobian?
Sebastian Walter
sebastian.walter at gmail.com
Fri Oct 2 04:10:32 EDT 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 at 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 at scipy.org
> http://mail.scipy.org/mailman/listinfo/scipy-user
>
More information about the SciPy-User
mailing list