[SciPy-User] lsoda vs. Coulomb friction
Charles R Harris
charlesr.harris at gmail.com
Wed Feb 3 15:16:38 EST 2010
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>
More information about the SciPy-User
mailing list