eja: fix random_instance() for JordanSpinEJA.

eja: add random_instance() for BilinearFormEJA.

eja: alter BilinearFormEJA to take only a matrix argument.

The BilinearFormEJA constructor originally took both the size "n" and
a bilinear form matrix "B" that acted on the x-bar component of an
algebra element. This was a bit silly: except when n=0, we can
determine the size of the algebra from the bilinear form matrix.

Now, we insist that the entire matrix "B" be passed in, including the
upper-left "1" block. This allows us to pass in a matrix of size zero
to get the trivial algebra, and in all other cases, to infer that the
dimension of the algebra is the size of the matrix.

eja: simplify the element-subalgebra one() method.

eja: add a random comment.

eja: add yet another TODO.

eja: minor improvement to the algebra one() method.

eja: try to speed up super/subalgebra conversion.

eja: drop redundant vector->superalgebra-element conversion.

eja: expand an existing subalgebra test.

eja: add another TODO item.

eja: update a few TODO items.

eja: drop cached superalgebra basis from subalgebras.

eja: drop an unused variable.

eja: echelonize a matrix manually to ensure a decent algorithm is used.

eja: don't compute an unused vector space for the element subalgebra.

eja: add a comment about sub-subalgebra representations.

eja: move the TrivialEJA random_instance() into an override method.

eja: don't orthonormalize the basis when computing minimal polynomials.

eja: drop obsolete _a_regular_element() override method.

eja: make _max_random_instance_size() abstract.

eja: rename _max_test_case_size() -> _max_random_instance_size().

eja: add inner_product() for DirectSumEJA.

eja: add inclusion maps for DirectSumEJA.

eja: add projections for DirectSumEJA.

eja: unindent a large "else" clause.

eja: fix base ring of fast _charpoly_coefficients().

eja: begin generalizing the charpoly-over-QQ optimizations.

As a very simple first step, encode the fact that HadamardEJA and
BilinearFormEJA have rational bases within the type system. This
allows us to compute their charpolys over the rationals instead of
whatever their user-given base ring is.

eja: split "check" args into check_field and check_axioms.

eja: add a few more TODO items.

eja: improve how we get the "real part" of a possibly-complex trace.

eja: enable consistency checks by default in (non-element) subalgebras.

eja: pass check=False for known-good constructions.

eja: pass check=False in the fast charpoly_coefficients() method.

eja: use fuzzy equality test with inexact base rings.

eja: check that the multiplication table is square with check=True.

eja: check EJA properties with check=True.

eja: add DirectSumEJA constructor (and not much else).

eja: add a link to a trac ticket in random_element().

eja: make random_element() more random.

eja: add orthogonality examples/tests for the Peirce decomposition.

eja: add a TODO to fix random_element().

eja: add "of" to the algebra characteristic_polynomial() method name.

eja: add another TODO that I keep forgetting.

eja: remove wayward utf-8 comment.

eja: let random_eja() return zero-dimensional algebras of all types.

eja: fix two more tests to work in dimension zero.

eja: make an inverse() test work in trivial algebras.

eja: make all tests work in trivial algebras.

This lets us drop the "nontrivial" argument to the random_eja()
function, because now we don't need it.

eja: document why an a[0] access is safe.

eja: expound on an existing TODO item.

eja: allow matrix algebras of "size" zero.

eja: remove another completed TODO.

eja: remove a completed TODO.

eja: remove one use of nontrivial=True passed to random_eja().

eja: fix the fast matrix _charpoly_coefficients() method.

eja: speed up _charpoly_coefficients when rank is cached.

eja: fix the determinant in trivial algebras.

eja: remove more vestigial charpoly stuff.

makefile: don't allow tests to time out.

eja: simplify and unify the charpoly/rank stuff.

eja: replace rank computation with length of charpoly coefficients.

If we cache the solution to the charpoly system instead, then the
result is useful in more than one place (for the charpoly itself,
in particular).

eja: get rid of the KnownRankEJA class.

Not needed now that I know how to compute rank.

eja: reimplement random_eja() without KnownRankEJA.

eja: fix tests and pre-cache ranks.

eja: replace the rank accessor with the full computation.

eja: use now somewhat-proven and fast rank computation.

eja: add TODON'T use Cramer's rule.

eja: add two more experimental rank computations.

eja: add voodoo _rank_computation() method for EJAs.

eja: drop the a_jordan_frame() method on EJAs.

I understand this better now, and don't expect this method to work
in general. It works for associative algebras, but we already know
how to find a Jordan frame in associative algebras: the spectral
decomposition, which works on the associative algebra alg({x}).

Revert "eja: start experiment full_spectral_decomposition() method."

This reverts commit 6d32150d3f7c175473a1d590facce1ea7c5ca77a.
I don't think this will ever work reliably.

eja: start experiment full_spectral_decomposition() method.

eja: fix sub-subalgebra element construction.

mjo/eja/eja_element_subalgebra.py: switch two tests back to QQ.

mjo/eja/eja_element.py: test spectral decomposition in subalgebras.

mjo/eja/eja_element.py: remove a dead line.

eja: make AA the default field because everything cool requires it.

mjo/eja/eja_algebra.py: fix unimportant left/right eigenspace terminology.

eja: ignore errors more specifically in a_jordan_frame().

eja: rename orthogonal_idempotents() to a_jordan_frame().

mjo/eja/eja_algebra.py: add WIP orthogonal_idempotents() method.

mjo/ldlt.py: goodbye, use https://trac.sagemath.org/ticket/10332 instead.

mjo/ldlt.py: add fast block-LDLT-based is_positive_semidefinite().

mjo/ldlt.py: move implementation to SageMath proper.

The actual block-LDLT code is now on the u/mjo/ticket/10332 branch.

mjo/ldlt.py: refactor into user-(un)friendly portions.

mjo/ldlt.py: delete naive implementations; fix tests.

mjo/ldlt.py: add examples and tests for block_ldlt().

mjo/ldlt.py: make block_ldlt() work on Hermitian matrices, too.

mjo/ldlt.py: switch some row/column indexing around.

mjo/ldlt.py: rule #1 is never compute the inverse of a matrix.

mjo/ldlt.py: get fast block_ldlt() working.

mjo/ldlt.py: fix a dimension in block_ldlt_naive().

mjo/ldlt.py: fix two bugs in the imperative block_ldlt() function.

mjo/ldlt.py: begin fast block-LDLT implementation.

mjo/ldlt.py: refactor the one-by-one block-LDLT pivot.

mjo/ldlt.py: start a naive Bunch-Kaufman block LDLT.

mjo/ldlt.py: various efficiency improvements to ldlt_fast().

mjo/ldlt.py: avoid some intermediate variables in ldlt_fast().

mjo/ldlt.py: add a vaguely correct non-recursive ldlt_fast().