← Back to team overview

yade-dev team mailing list archive

[svn] r1716 - trunk/pkg/dem/Engine/StandAloneEngine

 

Author: eudoxos
Date: 2009-03-09 22:47:54 +0100 (Mon, 09 Mar 2009)
New Revision: 1716

Modified:
   trunk/pkg/dem/Engine/StandAloneEngine/ElasticContactLaw.cpp
Log:
hotfix related to last commit


Modified: trunk/pkg/dem/Engine/StandAloneEngine/ElasticContactLaw.cpp
===================================================================
--- trunk/pkg/dem/Engine/StandAloneEngine/ElasticContactLaw.cpp	2009-03-09 20:57:46 UTC (rev 1715)
+++ trunk/pkg/dem/Engine/StandAloneEngine/ElasticContactLaw.cpp	2009-03-09 21:47:54 UTC (rev 1716)
@@ -180,18 +180,18 @@
 	////////// PFC3d SlipModel
 	
 			Vector3r f=currentContactPhysics->normalForce + shearForce;
-			Vector3r c1x(currentContactGeometry->contactPoint-de1->se3.position),
-				c2x(currentContactGeometry->contactPoint-de2->se3.position);
+			Vector3r _c1x(currentContactGeometry->contactPoint-de1->se3.position),
+				_c2x(currentContactGeometry->contactPoint-de2->se3.position);
 			#ifdef BEX_CONTAINER
 				ncb->bex.addForce (id1,-f);
 				ncb->bex.addForce (id2,+f);
-				ncb->bex.addTorque(id1,-c1x.Cross(f));
-				ncb->bex.addTorque(id2, c2x.Cross(f));
+				ncb->bex.addTorque(id1,-_c1x.Cross(f));
+				ncb->bex.addTorque(id2, _c2x.Cross(f));
 			#else
 				static_cast<Force*>   ( ncb->physicalActions->find( id1 , actionForceIndex).get() )->force    -= f;
 				static_cast<Force*>   ( ncb->physicalActions->find( id2 , actionForceIndex ).get() )->force    += f;
-				static_cast<Momentum*>( ncb->physicalActions->find( id1 , actionMomentumIndex ).get() )->momentum -= c1x.Cross(f);
-				static_cast<Momentum*>( ncb->physicalActions->find( id2 , actionMomentumIndex ).get() )->momentum += c2x.Cross(f);
+				static_cast<Momentum*>( ncb->physicalActions->find( id1 , actionMomentumIndex ).get() )->momentum -= _c1x.Cross(f);
+				static_cast<Momentum*>( ncb->physicalActions->find( id2 , actionMomentumIndex ).get() )->momentum += _c2x.Cross(f);
 			#endif
 			
 			currentContactPhysics->prevNormal = currentContactGeometry->normal;