2019-06-16 14:18:26 -0500 | commented answer | Looking for stellation of polyhedra I see. As of now, it is now possible to deal with these kinds of stellation of polytopes in Sage. Currently, it can only deal with convex shapes. Eventually, it is a long term goal to be able to deal with collections of polyhedral cells. But currently, it is not really supported as "1 object".... |

2019-02-09 02:53:32 -0500 | received badge | ● Necromancer (source) |

2019-02-08 08:14:26 -0500 | answered a question | No Output Given : [] Since your system is linear, I would use matrices to model this. Perhaps I made a typing mistake while translating from what you wrote, but then I get the following. It seems that your system does not have a solution, the matrix is not full-rank over the polynomials ring over the parameters. |

2019-02-05 15:09:47 -0500 | commented answer | cartesian product polyhedra For a reference on the product of polyhedra, see p.390 in the Handbook of Discrete and Computational Geometry, Chapter 15: https://www.csun.edu/~ctoth/Handbook/... Which is the reference which Sage tries to follow for nomenclature. |

2019-02-05 15:03:52 -0500 | answered a question | Looking for stellation of polyhedra Hi! Could give a precise definition of what you mean by stellation? There is a quick reference page about things you can do with a polyhedron: https://doc.sagemath.org/html/en/them... Perhaps what you mean is related to the following methods. Say `P.barycentric_subdivision()` : constructs a polyhedron out of`P` whose facets correspond to complete flags of faces of`P` .`P.stack(face)` : puts a new vertex a little bit outside of the face`face` .
In general, Sage tries to follow the nomenclature presented in the Handbook of Discrete and Computational Geometry (Chapter 15): https://www.csun.edu/~ctoth/Handbook/... Is the operation described in there? Otherwise, please let us know what is the operation and we will add it to the desired methods. |

2018-12-24 12:20:40 -0500 | received badge | ● Enlightened (source) |

2018-12-24 12:20:38 -0500 | received badge | ● Good Answer (source) |

2018-11-11 11:00:10 -0500 | received badge | ● Necromancer (source) |

2018-11-11 07:06:37 -0500 | answered a question | I.variety() : missing solution values According to the documentation of we read:
Being miscomputed might also include "missing solutions"... This triangular decomposition also highly depends on the ideal and its groebner basis computation. Currently, it seems indeed that getting a 0-dimensional variety out of an ideal is rather hazardous. I am using Bertini to do such computations, eventually is could be interfaced in Sage. Otherwise, maybe the interface to PHC could do the thing: |

2018-07-05 06:47:46 -0500 | received badge | ● Commentator |

2018-07-05 06:47:46 -0500 | commented answer | Can I intersect the boundaries of two polyhedra and display it? You can also do the intersection operation by using "&": |

2018-07-05 06:44:33 -0500 | answered a question | Getting convex hull of a set of points and plotting the polygon Since I see that you have decimal entries, this suggests that you have used "RDF" as a base ring for your polygon. You may check if it is the case by calling P.base_ring() There is a known bug in RDF that pops up when plotting: https://trac.sagemath.org/ticket/24877 Otherwise, if you just want to visualize your polygon your should be able to visualize it with the command ".plot()". More details will be explained in the extended tutorial which will appear in version 8.3: https://trac.sagemath.org/ticket/22572 Under "Geometry", "A longer introduction... " (once version 8.3 is out): |

2018-05-02 17:34:21 -0500 | received badge | ● Nice Answer (source) |

2018-05-02 17:34:03 -0500 | received badge | ● Nice Answer (source) |

2018-05-02 17:33:31 -0500 | received badge | ● Necromancer (source) |

2018-05-02 09:19:17 -0500 | answered a question | find one interior point of a polyhedron If all you want is some point that would be in the interior of the solution set, you can get one by doing: which returns the barycenter of the vertices of the polyhedron and adds the potential generating rays. So it would be in the relative interior of the polyhedron. |

2018-05-02 09:09:50 -0500 | received badge | ● Necromancer (source) |

2018-05-02 09:08:01 -0500 | commented answer | Why is this polar polytope incorrect? Indeed, polarity of polyhedron requires compactness (to get a good duality property, for example that taking twice the dual gives back the original object) and that the origin is in the interior of the (then compact) polyhedron. Hence, in order to allow polarity of polytope that do not contain the origin, it is automatically shifted so that its barycenter is the origin. This allows for a somehow "canonical" polar object. Taking another point in the interior as the origin would also work, but not implemented. |

2018-05-02 08:56:23 -0500 | answered a question | Speed of enumerating integer points in polytope In Sage 8.1, I can reproduce your code as follows (notice that I added "ppl" as a solver since the default behavior changed between your post and this answer). In order to speed things up, I use (You can install it by typing the commands in a terminal: ) and then use the backend which should appear instantaneously. |

2018-05-02 08:35:48 -0500 | commented question | random polytopes Note that generating random d-polytopes is a sensitive thing to do: you end up with probability 1-epsilon with a simplicial polytope, which may or may not be what you want to have. |

2018-05-02 08:32:31 -0500 | received badge | ● Necromancer (source) |

2018-05-02 08:24:44 -0500 | answered a question | Can sage compute the h* vector of a polytope? You can use the You can install them by typing: in a terminal once this ticket has been merged. Then you can type in sage: Note that this requires the latest features of this ticket. |

2018-05-02 08:09:19 -0500 | commented question | How to solve a large list of simultaneous equations Which question is it a duplicate of? |

2018-01-10 02:58:18 -0500 | commented question | Obtaining a certain list of posets via SAGE Since what you are looking for are exactly distributive lattices, perhaps it would be good to change the title of the question to indicate that this is what is looked for. Something like "Obtaining distributive lattices from a list of posets via SAGE" |

2017-10-19 01:10:05 -0500 | received badge | ● Necromancer (source) |

2017-10-11 10:23:50 -0500 | received badge | ● Nice Answer (source) |

2017-10-11 02:06:00 -0500 | answered a question | Hi, there: I am computing the eigenvalues and eigenvectors of the following complex matrix. It looks that Sage does not give the right answer? The ring of complex floats So while doing computation in floating point entries (i.e. in An alternative is to use an algebraic extension of |

2017-10-11 01:47:59 -0500 | commented question | Hi, there: I am computing the eigenvalues and eigenvectors of the following complex matrix. It looks that Sage does not give the right answer? Could you add the values that you obtain and perhaps the version of Sage that you are using? What are the expected eigenvalues? In Sage version 8.1beta7 I get sage: B.eigenvalues()
[3.917412364954441 + 0.054291888902286756 Is that correct? |

2017-10-10 15:12:28 -0500 | commented answer | Obtaining lattices quickly in SAGE This is true. I appended the answer. |

2017-10-10 13:51:09 -0500 | answered a question | Projection along affine hull Here is an example. It requires to load the Let's say we have the following polytope: and we are interested in the affine space spanned by the points we first get a linear subspace and create the projection matrix: Then, we create a function that will emulate the affine projection (notice the addition of the first element of the affine basis). Then, we project the polytope and use the transformed coordinates of the vertices to create a new polytope. Now, we can see the difference between and the resulting polytope: |

2017-10-10 12:55:27 -0500 | commented question | Projection along affine hull Could we have a precise definition of what is meant by "projection of |

2017-10-09 10:32:59 -0500 | answered a question | Obtaining lattices quickly in SAGE It is not usual to bundle a database of combinatorial objects coming with Sage by default. It would be nice to have a (fast) enumerated set of all lattices on a given number of elements. Although my impression is that the set of lattices on a given number of vertices is not easy to generate recursively. Indeed, the lattice property of a poset does not behave well with respect to inclusion, so a recursive enumeration would not work as smoothly as it does for posets. That said, depending on what you want to do, you could use an iterator: to generate all connected lattices one by one, but it does not increase the speed of generation... |

2017-10-09 09:14:57 -0500 | answered a question | force 0/1 in V-representation of a set of inequalities From what I understand of your question, you could scale the polyhedron by the lcm of the denominators of the coordinates of the vertices. This will give you a polyhedron with integer vertices. Nevertheless, in general, it is not true that every (combinatorial type of) polyhedron object has a 0/1 V-representation. So your example should have a special structure that allows this to be possible. Unfortunately, we can not help further since we don't have an explicit example... |

2017-10-09 08:59:14 -0500 | answered a question | combinatorial equivalence for Polyhedra / isomorphism for lattices Hi, Since version 7.6 of Sagemath, there is exactly the required method for Polyhedron objects: It does not requires the computation of the full face lattice (checks vertex/facet adjacency graph isomorphism), hence should be faster than asking for isomorphic face lattices. The |

2017-10-09 08:59:09 -0500 | received badge | ● Necromancer (source) |

2017-10-09 08:59:09 -0500 | received badge | ● Teacher (source) |

2017-10-09 08:48:47 -0500 | answered a question | Visualizing Polyhedron with specified range First, to remove the arrows, one possible way is to use Second, to get a bigger picture, there is a work around by multiplying the rays by a certain scalar. This is not as clean as specifying the ranges on each axis, but is a start. Currently, to plot the polyhedron object, it uses the data of the polyhedron and it is not (yet) possible to set a certain range for the axis. |

2017-10-09 07:34:58 -0500 | received badge | ● Editor (source) |

2017-10-09 07:34:08 -0500 | answered a question | Can the base ring of a polyhedron be restricted? The current difference between the base ring A polytope (a compact polyhedron) may have all vertices with integer coordinates, so that it could be called a "Lattice polytope", but its The tickets Add .change_ring() method for polyhedra and Add .change_backend() method for polyhedra will then make it possible to apply the desired changes on the base ring and backend, if possible. |

2017-10-09 07:19:15 -0500 | answered a question | Is there a function to map faces of a cone its dual? As far as I know, there is no such function in Sage yet. But the following ticket contains a discussion about this issue: Ticket 17215. It may eventually be implemented... |

2017-10-09 06:45:57 -0500 | commented answer | Generating Tikz codes from Sage for drawing scattering points in 3D You can change the vertex, edge and facet color directly in the tikz method using the parameters "vertex_color", "edge_color" and "facet_color". They can take any string argument which tikz can interpret as a color. One slight issue with the above loop on integral points is that it may not iterate the interior lattice points according to the projection, i.e. from the furthest to closest to the screen, and maybe some point nodes intersect in a weird way. You could sort them as follows: |

Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.