Error intersecting polyhedra

asked 2013-06-06 10:13:30 -0600

Bill gravatar image

updated 2015-01-13 13:35:46 -0600

FrédéricC gravatar image

I have a dictionary of polyhedra which I want to consider various intersections and unions of. For now, my dictionary has four polyhedra in it. When running the following for loop:

v={} # contains coordinates of the tetrahedron, centered at (0,0,0)
v[1] = vector([1,0,-1/sqrt(2)])
v[2] = vector([-1,0,-1/sqrt(2)])
v[3] = vector([0,1,1/sqrt(2)])
v[4] = vector([0,-1,1/sqrt(2)])

V = {1,2,3,4}

b={} # will contain coordinates of vertices of the subdivision of the tetrahedron

for X in powerset(V):
  q = len(X)
  if q > 0:
    b[tuple(X)] = (1/q)*sum([v[i] for i in X])

poly_int = Polyhedron([b[key] for key in b.keys()]) # this is a tetrahedron containing
# all of my polyhedra
for i in [1,2,3,4]:
  poly_int = poly[i] & poly_int

I get a very nasty error:

AttributeError: 'Polyhedra_RDF_cdd_with_category.element_class'
object has no attribute '_Vrepresentation'

(The traceback for the error is quite long. I can include it if needed.)

Manually computing

poly[1] & poly[2] & poly[3] & poly[4]

gives the same error. HOWEVER, if instead I compute

(poly[1] & poly[2]) & (poly[3] & poly[4])

I get no errors. While this is a fine workaround for this particular case, I'm scripting these intersections and need some way to actually compute these intersections automatically. Using something like

poly[1].intersection(poly[2]) # etc.

doesn't help either. Any thoughts?

edit retag flag offensive close merge delete

Comments

If you could edit your question to include the definition of `b` that would probably help; I'm trying some basic polyhedra and just getting the empty polyhedron even with multiple `&`s, and the diagnosis you give indicates it's probably a very specific edge case that is going wrong. (Also note that `intersection` *should* (?) be the same as `&`, since this particular binary operator turns into `__and__`, I believe.)

kcrisman gravatar imagekcrisman ( 2013-06-06 14:33:42 -0600 )edit

@kcrisman OK. I included a definition of b. It's a dictionary containing the coordinates of the barycentric subdivision of a 3-simplex.

Bill gravatar imageBill ( 2013-06-06 14:40:15 -0600 )edit

Sorry, one more question - what is `poly`? I thought maybe `v` is `poly`, but that doesn't make sense.

kcrisman gravatar imagekcrisman ( 2013-06-06 14:49:00 -0600 )edit

Right. Sorry. Each poly[i] is a polyhedron defined as a convex hull on some of the points from b.

Bill gravatar imageBill ( 2013-06-06 14:50:29 -0600 )edit

With code, please... :-) Again, my first guesses at *which* points didn't work. Sorry for asking for *all* of it, but a "minimal working example" is key for debugging this. We don't want to try to guess what you mean.

kcrisman gravatar imagekcrisman ( 2013-06-06 14:57:51 -0600 )edit