- # Now figure out the entries of the right-multiplication
- # matrix for the successive basis elements b0, b1,... of
- # that subspace.
+ # If our superalgebra is a subalgebra of something else, then
+ # these vectors won't have the right coordinates for
+ # V.span_of_basis() unless we use V.from_vector() on them.
+ basis_vectors = map(power_vectors.__getitem__, ind_rows)
+ else:
+ # If we're going to orthonormalize the basis anyway, we
+ # might as well just do Gram-Schmidt on the whole list of
+ # powers. The redundant ones will get zero'd out.
+ superalgebra_basis = gram_schmidt(powers)
+ basis_vectors = [ b.to_vector() for b in superalgebra_basis ]
+
+ W = V.span_of_basis( V.from_vector(v) for v in basis_vectors )