]> gitweb.michael.orlitzky.com - sage.d.git/commit
eja: drop custom gram_schmidt() routine that isn't noticeably faster.
authorMichael Orlitzky <michael@orlitzky.com>
Tue, 29 Oct 2019 01:57:11 +0000 (21:57 -0400)
committerMichael Orlitzky <michael@orlitzky.com>
Tue, 29 Oct 2019 01:57:11 +0000 (21:57 -0400)
commit1e9700cdd04434465ffcad148d078f7fa361e426
tree92166e1cc90aff18a42685b2bc22b3495141c959
parentc16c11d5b4cd015b6d7516a175a92f46e03bd2d7
eja: drop custom gram_schmidt() routine that isn't noticeably faster.

Using Gram-Schmidt in SageMath is a bit of a pain, since you have to
put your vector space elements into a matrix, which means you have to
convert them to vectors first. And then you wind up transposing the
matrix twice, and then you have to convert everything back to vectors
and then back to algebra elements, dropping the zero vectors/elements
in the process. Sounds slow, right?

Well, testing shows that jumping through all those hoops and using the
default implementation isn't much slower than the custom routine I
wrote. So, for now, let's delete it.

Beware that this breaks the construction of subalgebras when the field
is not exact because SageMath won't orthonormalize over an inexact
field.
mjo/eja/eja_subalgebra.py
mjo/eja/eja_utils.py