[SciPy-User] lsoda vs. Coulomb friction

[Charles R Harris]

On Wed, Feb 3, 2010 at 9:28 AM, Ryan Krauss <ryanlists at gmail.com> wrote:

> I am trying to use odeint (i.e. lsoda) on a mechanical system that
> involves a mass and friction modeled as a viscous term plus Coulomb
> friction.  The discontinuity near 0 velocity makes lsoda mad.  It
> complains that it has to do excess work (lazy algorithm :).  In spite
> of its complaining, the result leads to fairly good agreement between
> model and experiment.  But is there a better way to handle this?
>
>
You might try some of the other stiff solvers. The algorithms for stiff
equations tend to be implicit, i.e., they have to solve non-linear equations
and are inherently slower for that reason. Other than that I haven't any
suggestions.

Chuck
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.scipy.org/pipermail/scipy-user/attachments/20100203/aa47a9f6/attachment.html>