- S = Cone(linspace_gens, K.lattice())
-
- # Since ``S`` is a subspace, its dual is its orthogonal complement
- # (albeit in the wrong lattice).
- S_perp = Cone(S.dual(), K.lattice())
- P = K.intersection(S_perp)
+ # Since ``S`` is a subspace, the rays of its dual generate its
+ # orthogonal complement.
+ P = K.intersection( Cone(S.dual(), K.lattice()) )