yade-dev team mailing list archive

[Bruno Chareyre <bruno.chareyre@xxxxxxxxxxx>]

> 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