- A = matrix(field, vector_basis)
- # uh oh, this is only the "usual" inner product
- Q,R = A.gram_schmidt(orthonormal=True)
- self._deorthonormalization_matrix = R.inverse().transpose()
- vector_basis = Q.rows()
+ from mjo.eja.eja_utils import gram_schmidt
+ vector_basis = gram_schmidt(vector_basis, inner_product)