yade-dev team mailing list archive

[noreply@xxxxxxxxxxxxx]

------------------------------------------------------------
revno: 2358
committer: bchareyre <bchareyre@dt-rv020>
branch nick: trunk
timestamp: Mon 2010-07-12 21:19:53 +0200
message:
  - Last change in contact plane makes script unstable. Restore stable (but wrong?...) equations.
modified:
  pkg/common/DataClass/Shape/Cylinder.cpp


--
lp:yade
https://code.launchpad.net/~yade-dev/yade/trunk

Your team Yade developers is subscribed to branch lp:yade.
To unsubscribe from this branch go to https://code.launchpad.net/~yade-dev/yade/trunk/+edit-subscription

=== modified file 'pkg/common/DataClass/Shape/Cylinder.cpp'
--- pkg/common/DataClass/Shape/Cylinder.cpp	2010-07-12 18:44:19 +0000
+++ pkg/common/DataClass/Shape/Cylinder.cpp	2010-07-12 19:19:53 +0000
@@ -135,7 +135,10 @@
 	//length only used for display
 	s->length=length;
 	scm->penetrationDepth=s->initLength-length;
-	scm->contactPoint=pChain->pos+segment;
+	
+// 	scm->contactPoint=pChain->pos+segment;//generates instabilites (TODO)
+	scm->contactPoint=pChain->pos+pChain->ori*Vector3r::UnitZ()*length;
+
 #ifdef IGCACHE
 	if (scene->isPeriodic) {
 		Vector3r shiftVel = scene->cell->velGrad*scene->cell->Hsize*c->cellDist.cast<Real>();