[ODE] ODE crashes if fdir1 vector collinear with vector of normal

Dmitry Medvedev pegasi at inbox.ru
Thu May 3 06:52:09 MST 2007

 Hello, everybody.

 Lets see the following situation for example:

  fdir = {1.0, 0.0, 0.0}
  contacts normal vector = { a, 0.0, 0.0}, where a - some value

  division by zero occurs in the quickstep.cpp file in function "static void SOR_LCP(....." in line 
    Ad[i] = sor_w / (sum + cfm[i]);

 Perhaps some special handling should be implemented for such cases.

 Thank you.

More information about the ODE mailing list