+       const shared_ptr<Body>&  body1=Body::byId(contact->getId1(),rootBody);
+       const shared_ptr<Body>&  body2=Body::byId(contact->getId2(),rootBody);
+       const shared_ptr<RpmMat>&  
rbp1=YADE_PTR_CAST<RpmMat>(body1->physicalParameters);
+       const shared_ptr<RpmMat>&  
rbp2=YADE_PTR_CAST<RpmMat>(body2->physicalParameters);
+                       rbp1->isDamaged=true;
+                       rbp2->isDamaged=true;
That reminds me (thanks): I know there is a bug in CPM incrementing number of borken contacts per body, which should be protected with #pragma omp critical section {...} so that only one thread at time modifies the rbp1/rbp2. Concurrency could make increments getting lost and other weird things. I will fix it in CPM tonight. Vaclav

_______________________________________________
Mailing list: https://launchpad.net/~yade-dev
Post to     : [email protected]
Unsubscribe : https://launchpad.net/~yade-dev
More help   : https://help.launchpad.net/ListHelp

Reply via email to