The code is a copy the one in Operation_on_polyhedra that was only working with Polyhedron_3. For the moment only default parameters are used. NamedParameters will be introduced later.