yade-dev team mailing list archive

[noreply@xxxxxxxxxxxxx]

------------------------------------------------------------
revno: 2239
committer: Václav Šmilauer <eudoxos@xxxxxxxx>
branch nick: trunk
timestamp: Tue 2010-05-18 23:22:14 +0200
message:
  Fix compilation error (Bruno, you must have been terribly tired today)
modified:
  pkg/dem/Engine/GlobalEngine/PeriIsoCompressor.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/dem/Engine/GlobalEngine/PeriIsoCompressor.cpp'
--- pkg/dem/Engine/GlobalEngine/PeriIsoCompressor.cpp	2010-05-18 16:29:56 +0000
+++ pkg/dem/Engine/GlobalEngine/PeriIsoCompressor.cpp	2010-05-18 21:22:14 +0000
@@ -127,9 +127,9 @@
 		//   (reversedForces?-1.:1.)*(Body::byId(I->getId1())->state->pos-Body::byId(I->getId2())->state->pos);
 		Vector3r branch= ( reversedForces?-1.:1. ) *gsc->normal* ( gsc->refR1+gsc->refR2 );
 		// tensorial product f*branch (hand-write the tensor product to prevent matrix instanciation inside the loop by makeTensorProduct)
-		stressTensor(0,0)+=f(0)*branch(0); stressTensor(1,0)+=f(1)*branch(0); stressTensor(2,0)+=f(2)*branch(0);
-		stressTensor(0,1)+=f(0)*branch(1); stressTensor(1,1)+=f(1)*branch(1); stressTensor(2,1)+=f(2)*branch(1);
-		stressTensor(0,2)+=f(0)*branch(2); stressTensor(1,2)+=f(1)*branch(2); stressTensor(2,2)+=f(2)*branch(2);
+		stressTensor(0,0)+=f[0]*branch[0]; stressTensor(1,0)+=f[1]*branch[0]; stressTensor(2,0)+=f[2]*branch[0];
+		stressTensor(0,1)+=f[0]*branch[1]; stressTensor(1,1)+=f[1]*branch[1]; stressTensor(2,1)+=f[2]*branch[1];
+		stressTensor(0,2)+=f[0]*branch[2]; stressTensor(1,2)+=f[1]*branch[2]; stressTensor(2,2)+=f[2]*branch[2];
 		//stressTensor+=makeTensorProduct( f, branch );
 		if( !dynCell )
 		{