Hi, yea, that is how I work now. N = ToricLattice(3) Ns = SomeSubLatticeOfN points = [p for p in P.intersection(plane).integral_points() if p in Ns]
but that seems not to be the nicest way to do it. greatz Am Sa 29 Dez 2012 18:20:57 CET schrieb Volker Braun: > A plane is a polyhedron, too: > > sage: P = polytopes.n_cube(3) > sage: plane = Polyhedron(eqns=[(0,1,1,1)]) > sage: plane.Hrepresentation() > (An equation (1, 1, 1) x + 0 == 0,) > sage: P.intersection(plane).integral_points() > ((-1, 0, 1), (-1, 1, 0), (0, -1, 1), (0, 0, 0), (0, 1, -1), (1, -1, > 0), (1, 0, -1)) > > > > On Saturday, December 29, 2012 3:03:32 PM UTC, Johhannes wrote: > > Hi List, > I need to know all lattice points in a given polytope with respect > to a > sublattice (or considered as submodule) of ZZ^n. > Is there a build-in way to do this? the default 'intersection' if > polyhedrons does not work for modules/lattices, and the > LatticePolytope > does not have a intersection method at all. > > greatz Johannes > > -- > You received this message because you are subscribed to the Google > Groups "sage-support" group. > To post to this group, send email to sage-support@googlegroups.com. > To unsubscribe from this group, send email to > sage-support+unsubscr...@googlegroups.com. > Visit this group at http://groups.google.com/group/sage-support?hl=en. > > -- You received this message because you are subscribed to the Google Groups "sage-support" group. To post to this group, send email to sage-support@googlegroups.com. To unsubscribe from this group, send email to sage-support+unsubscr...@googlegroups.com. Visit this group at http://groups.google.com/group/sage-support?hl=en.
