[ODE] recent dSetZero() patch wrong?

Nguyen Binh ngbinh at glassegg.com
Thu Oct 23 14:28:53 MST 2003


AP> It did make StepFast working incorrectly. The cause was in stepfast.cpp,
AP> around line 740: "dSetZero (J, 2 *
AP> m * 8);" shouldn't be commented out. I don't know the "meaning" of it, but
AP> it seems from the code that dSetZero is actually needed there...

    Yes, you are right! like *c, J is used as place holder for
    Jacobian sub matrices and it should be set to zero on initialize.


-- 
Best regards,

---------------------------------------------------------------------
   Nguyen Binh
   Software Engineer
   Glass Egg Digital Media
   
   E.Town Building  
   7th Floor, 364 CongHoa Street
   Tan Binh District,
   HoChiMinh City,
   VietNam,

   Phone : +84 8 8109018
   Fax   : +84 8 8109013

     www.glassegg.com
---------------------------------------------------------------------



More information about the ODE mailing list