yade-dev team mailing list archive

[Bruno Chareyre <bruno.chareyre@xxxxxxxxxxx>]

The analytical way maybe?

for (sphere in clump){
volume+=sphere.volume;
for (othersphere in clump && othersphere.index<sphere.index)
volume-=intersection(sphere,othersphere)
}

This one doesn't work, its a bit more complicated with multiple 
overlaps, but you get the idea.
CGAL has powerfull stuff for intersections, I'll have a look.

Bruno

On 08/07/10 20:29, Václav Šmilauer wrote:
> > This is a "must have", but unfortunately it is not so simple. Most
> > clumps have overlapping spheres, and overlapped volumes should be
> > counted only once. Here you will get spheresVolume>totalVolume, i.e.
> > negative porosity....
> >      
> It can be approximated via sampling the affected volume, though it is
> not trivial to make it efficient for larger number of particles (you'd
> have to use bboxes somehow).
>
>
>
> _______________________________________________
> Mailing list: https://launchpad.net/~yade-dev
> Post to     : yade-dev@xxxxxxxxxxxxxxxxxxx
> Unsubscribe : https://launchpad.net/~yade-dev
> More help   : https://help.launchpad.net/ListHelp
>
>