[SciPy-User] lsoda vs. Coulomb friction
Charles R Harris
Wed Feb 3 14:16:38 CST 2010
On Wed, Feb 3, 2010 at 9:28 AM, Ryan Krauss <email@example.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
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the SciPy-User