[SciPy-User] Why does scipy.integrate.odeint still provide a result even after "repeated convergence failures"?
Sturla Molden
sturla.molden at gmail.com
Tue May 12 08:36:54 EDT 2015
Are you trying to fit an ODE or just solve an ODE with known parameters?
Fitting parameters by solving an ODE can be a notoriously difficult
problem with respect to convergence. There are methods such as "multiple
shooting" which can improve the situation. But fitting an ODE is often
an ill-posed problem.
Sturla
On 12/05/15 07:44, Brian Merchant wrote:
> Hi Warren,
>
> Just an update: I was able to "zero" in on the portion of the model that
> causes the issue (I set the ode function to just return 0 for every
> other variable, except the one that causes issues).
>
> I can't find an error in the Python code responsible for calculating the
> derivative for that variable (I stepped through it, it simply implements
> a Hill equation:
> http://en.wikipedia.org/wiki/Hill_equation_%28biochemistry%29), and the
> results provided by odeint in the end don't seem *obviously* out-of-whack.
>
> Brian
>
> On Mon, May 11, 2015 at 9:52 PM, Warren Weckesser
> <warren.weckesser at gmail.com <mailto:warren.weckesser at gmail.com>> wrote:
>
>
>
> On Tue, May 12, 2015 at 12:43 AM, Brian Merchant
> <bhmerchant at gmail.com <mailto:bhmerchant at gmail.com>> wrote:
>
> Hi all,
>
> Why does scipy.integrate.odeint still provide a result even
> after "repeated convergence failures"? Is this result still of use?
>
>
>
> Which version of scipy are you using? Can you give a short example
> that demonstrates the problem?
>
> Warren
>
>
> Kind regards,
> Brian
>
> _______________________________________________
> SciPy-User mailing list
> SciPy-User at scipy.org <mailto:SciPy-User at scipy.org>
> http://mail.scipy.org/mailman/listinfo/scipy-user
>
>
>
> _______________________________________________
> SciPy-User mailing list
> SciPy-User at scipy.org <mailto:SciPy-User at scipy.org>
> http://mail.scipy.org/mailman/listinfo/scipy-user
>
>
>
>
> _______________________________________________
> 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