yade-dev team mailing list archive

[noreply@xxxxxxxxxxxxx]

------------------------------------------------------------
revno: 2786
committer: Janek Kozicki <cosurgi@xxxxxxxxxx>
branch nick: yade
timestamp: Thu 2011-03-17 10:40:47 +0100
message:
  standalone example almost works. I still have some problem with sample randomness, and exact answers.
modified:
  py/utils.py


--
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 'py/utils.py'
--- py/utils.py	2011-03-16 21:16:04 +0000
+++ py/utils.py	2011-03-17 09:40:47 +0000
@@ -513,17 +513,18 @@
 	:return: the porosity of the sample inside given volume
 
 	Example invocation:
-		>>> O.engines[5].porosity
-		0.3901808
-		>>> utils.voxelPorosityTriaxial(O.engines[5],200,0)
-		0.3908145
-		>>> utils.voxelPorosityTriaxial(O.engines[5],800,0)
-		0.3902078
-		>>> utils.voxelPorosityTriaxial(O.engines[5],200,avg_radius)
-		0.3604238
-		>>> utils.voxelPorosityTriaxial(O.engines[5],800,avg_radius)
-		0.3604175
-	
+		>>> from yade import utils
+		>>> O.reset()
+		>>> rAvg=0.03
+		>>> TriaxialTest(numberOfGrains=200,radiusMean=rAvg).load()
+		>>> O.dt=-1
+		>>> O.run(1000)
+		>>> O.engines[4].porosity
+		0.44007807740143889
+		>>> utils.voxelPorosityTriaxial(O.engines[4],200,0)
+		0.44055412500000002
+		>>> utils.voxelPorosityTriaxial(O.engines[4],200,rAvg)
+		0.36798199999999998
 	"""
 	p_bottom	= O.bodies[triax.wall_bottom_id].state.se3[0]
 	p_top		= O.bodies[triax.wall_top_id   ].state.se3[0]