yade-dev team mailing list archive

[eudoxos at BerliOS <eudoxos@xxxxxxxxxxxxxxx>]

Author: eudoxos
Date: 2009-08-04 10:46:09 +0200 (Tue, 04 Aug 2009)
New Revision: 1916

Modified:
   trunk/examples/concrete/uniaxial/rb.py
Log:
1. Reset interaction radius after the initial contacts have been created.


Modified: trunk/examples/concrete/uniaxial/rb.py
===================================================================
--- trunk/examples/concrete/uniaxial/rb.py	2009-08-04 06:38:03 UTC (rev 1915)
+++ trunk/examples/concrete/uniaxial/rb.py	2009-08-04 08:46:09 UTC (rev 1916)
@@ -76,7 +76,7 @@
 
 O.engines=[
 	BexResetter(),
-	BoundingVolumeMetaEngine([InteractingSphere2AABB(aabbEnlargeFactor=intRadius),MetaInteractingGeometry2AABB()]),
+	BoundingVolumeMetaEngine([InteractingSphere2AABB(aabbEnlargeFactor=intRadius,label='is2aabb'),MetaInteractingGeometry2AABB()]),
 	InsertionSortCollider(),
 	InteractionDispatchers(
 		[ef2_Sphere_Sphere_Dem3DofGeom(distanceFactor=intRadius,label='ss2d3dg')],
@@ -114,6 +114,10 @@
 		renderer['displacementScale']=(1000,1000,1000) if mode=='tension' else (100,100,100)
 	except ImportError: pass
 	print "init done, will now run."
+	O.step(); O.step(); # to create initial contacts
+	# now reset the interaction radius and go ahead
+	ss2d3dg['distanceFactor']=1.
+	is2aabb['aabbEnlargeFactor']=1.
 	O.run()
 
 def stopIfDamaged():