> Oh simply cubic. Say that you divide the main cell into 8 cubic parts
> and then you get the average of stress in each of them (I simply want
> to compare them).
Here is the algorithm :
foreach sphere:
cubeStress+=bodyStress(sphere.id)*volumeInsideCube
cubeStress/=cubeVolume
Then if you have N cubes, sum_N(cubeStress)/N will be exactly the stress
in the cell.
Probably similar to what Vincent suggests. I don't think it is
equivalent with weighting by branch length, since volumeInsideCube is a
non linear function of position.
Writing the volumeInsideCube function will be the boring part...
Bruno
_______________________________________________
Mailing list: https://launchpad.net/~yade-dev
Post to : [email protected]
Unsubscribe : https://launchpad.net/~yade-dev
More help : https://help.launchpad.net/ListHelp
