- # First, compute the basis B...
- x0 = self.zero()
- c = 1
- for g in self.gens():
- x0 += c*g
- c +=1
- if not x0.is_regular():
- raise ValueError("don't know a regular element")
-
- V = x0.vector().parent().ambient_vector_space()
- V1 = V.span_of_basis( (x0**k).vector() for k in range(self.rank()) )
- B = (V1.basis() + V1.complement().basis())
-