%
(lattice.rank(), n))
- def d(j):
- v = [1]*n # Create the list of all ones...
- v[j] = 1 - p # Now "fix" the ``j``th entry.
- return v
-
- G = identity_matrix(ZZ,n).rows() + [ d(j) for j in xrange(n) ]
+ I = identity_matrix(ZZ,n)
+ M = matrix.ones(ZZ,n) - p*I
+ G = identity_matrix(ZZ,n).rows() + M.rows()
return Cone(G, lattice=lattice)