yade-dev team mailing list archive

[Václav Šmilauer <eudoxos@xxxxxxxx>]

> Hi all,
> 
> Maybe I should have report a bug for that, but I do not really get
> where is the issue...
Hi Luc,

thanks for reporting that, it was really a bug. Next time, don't
hesitate to file bug at http://bugs.launchpad.net, such things should
not be forgotten.

The problem was that Ig2_Facet_Sphere_Dem3DofGeom created the
interaction, but the reference (initial) length was computed a little
differently than displacementN() called from in Law2_.._CundallStrack;
it gave displacement of something like 1e-16 (numerical residuum), which
was "positive", hence the interaction was deleted immediately.

It is fixed now (see attachment; r1914) and in the future, the
dispatcher will detect if contact is made real and deleted at the very
same step, reporting an error.

(Perhaps such checks should be done only in debug build?)

Cheers, Vaclav

Attachment: fixed.png
Description: PNG image