CGAL users discussion list

[perrone <>]

I'm an italian student of Salento University (Lecce,Italy) and in my 
mechanical
engineering degree thesis i want to make a model of an aluminium open cell 
foam
using laguerre tessellation.
Some applications show that random Laguerre tessellations generates by systems
of non-overlapping spheres are promising models for foam structures. So,ad
first,i have to make a matlab routine to generete a hard spheres packing.the
algorithm is the seguent:
First,i have to create a set of N spheres with x,y and z center coordinates
uniformly distribuited in a parallelepipedal container and with diameters d
chosen accordingly to a prescribed distribution function(in my case,poisson
distribution). Thus,ovelappings of the spheres in the initial ocnfiguration 
are
typical. Then i have to reduce these overlaps.
After this non-overlapping system of dense spheres is created,i have to find
each y point that satisfy this relationship
   ||x(1)-y||^2-r(1)^2 <= ||x(2)-y||^2-r(2)^2
where
x(1)= x,y,z center coordinates vector of a first sphere
r(1)= radius of the first sphere
x(2)=x,y,z center coordinates vector of a second sphere
r(2)= radius of the second sphere
y= x,y,z coordinates vector of y points satisfying the relationship.

I want to obtain this using CGAL and i have read that you had the same problem
some years ago.
It would be a big pleasure for me if you could answer and help me in such
way,like a code of a piece of it.
I hope in your answer.
With best regards,
 
Francesco Perrone