eja: add span_of_powers() method.

eja: remove an erroneous comment.

eja: begin implementing an element class.

eja: implement my own algebra subclass.

eja: fix field arguments, add comment.

eja: add eja_minimal_polynomial() function.

eja: add eja_ln() for the rank-two algebra.

eja: add examples for eja_rn.

eja: begin draft implementation of euclidean jordan algebras.

mjo/polynomial.py: remove a redundant line from a test.

mjo/polynomial.py: add another multidiv test based on an exercise.

Revert "mjo/polynomial.py: improve multidiv performance a bit."

This reverts commit 4e4efc9eff5c77a2ca19002b1dfa45598e974c54. Example
4 in the text doesn't work with this implementation, but it does with
my original one. Let's stick with slow and correct, for now.

mjo/polynomial.py: add a (failing!) test for Example 4 in the text.

mjo/polynomial.py: add two more examples from exercises in the text.

mjo/polynomial.py: improve multidiv performance a bit.

We don't need to re-check earlier "denominators" after division has
occurred, even though we *do* need to re-check the current one. This
makes the algorithm closer in practice to the one given in the text.

mjo/polynomial.py: improve tests for new multidiv function.

mjo/orthogonal_polynomials.py: fix doctests.

Add new mjo.polynomial module with the "multidiv" function.

This new module implements the generalized polynomial division
algorithm in Section 2.3 of Cox, Little, and O'Shea.

Revert "... add completely_positive_operators_gens()."

This reverts commit 3d6ec3f4f138a6278be4f393b491aa1ac2bbffa8. The new
function was mathematically wrong, and can't be fixed. Oops.

mjo/cone/completely_positive.py: add completely_positive_operators_gens().

Add a new function to return (as matrices) the generators of the
completely-positive cone of K, where K is some other given cone.

mjo/cone/symmetric_pd: add inverse_symmetric_pd() function.

mjo/cone/symmetric_pd: new module for symmetric positive-definite matrices.

mjo/cone/symmetric_psd: add missing set_random_seed() calls.

mjo/cone/symmetric_psd: don't use recusion; work on trivial matrices.

mjo/cone: improve docstrings the the rearrangement and schur cones.

mjo/cone: rename random_psd() to random_symmetric_psd().

cone/trivial_cone.py: add the trivial_cone() function.

cone/rearrangement.py: simplify the implementation.

Use matrix tricks instead of a local function to construct the
generators.

cone/rearrangement.py: add a "lattice" argument.

cone/rearrangement.py: fix the test for propriety.

We had a test for the rearrangement cone being proper, but it
incorrectly tested the p=0 case. Now that we're sanity-checking p,
that test was throwing a ValueError. Fix it to use p >= 1.

cone/rearrangement.py: add a few more tests.

Add a few more tests to the rearrangement_cone() function. In
particular, delete the has_rearrangement_property() function since it
only served to test the rearrangement cone. The property tests were
moved into the cone doctests.

The cone constructor now also checks that "p" is within bounds, and
throws a ValueError if not. There are tests for this too.

cone/rearrangement.py: add some preliminary references.

cone/schur.py: add "lattice" argument.

cone/nonnegative_orthant.py: add a word to the docs.

cone/nonnegative_orthant.py: add "lattice" argument.

cone/schur.py: add a test for the dual of the Schur cone.

cone/symmetric_psd.py: use two more generator expressions.

cone/symmetric_psd.py: use a generator expression in unit_eigenvectors().

cone/cone.py: use generator expressions instead of lists.

orthogonal_polynomials.py: use generator expressions where applicable.

orthogonal_polynomials.py: use xrange everywhere.

matrix_vector.py: use xrange everywhere.

mjo/interpolation.py: partial conversion to generator expressions.

mjo/interpolation.py: use xrange everywhere.

cone/schur.py: use xrange in one place.

cone/rearrangement.py: use xrange everywhere.

cone/doubly_nonnegative.py: use xrange in two places.

cone/doubly_nonnegative.py: use generator expressions in two places.

cone/faces.py: use xrange in two places.

cone/permutation_invariant.py: use xrange in two places.

symbol_sequence.py: allow _flatten_list_ to work on an iterable.

symbol_sequence.py: use xrange where applicable.

symbol_sequence.py: use generator expressions where applicable.

cone/rearrangement.py: use generator expressions where applicable.

cone/permutation_invariant.py: add test for the nonnegative orthant.

cone/permutation_invariant.py: use generator expressions where applicable.

cone/schur.py: use generator expressions where applicable.

cone/faces.py: use generator expressions where applicable.

cone/schur.py: test that the schur cone induces the majorization ordering.

cone/schur.py: clean up and add some tests.

mjo/cone: add the nonnegative_orthant module.

cone/faces.py: add the missing HilgertHofmannLawson reference that I cited.

cone/faces.py: test that dual-face-of-dual-face is a no-op.

cone/faces.py: test the relationship between dual_face and face_generated_by.

cone/faces.py: add preliminary dual_face() operation.

cone/faces.py: add test for Proposition 3.1 and Corollary 3.9 in Tam.

cone/faces.py: test a lattice join property of face_generated_by().

cone/faces.py: add a relative interior test for face_generated_by().

cone/faces.py: remove a misleading comment.

cone/faces.py: update the wording of two tests.

README: update comments on the test suite.

Kill unnecessary sitedir mangling.

At some point in the past, mangling the sitedir made the test suite
work. Now it doesn't, and we have to do something else. So clean up
the sitedir stuff.

Add the new mjo.cone.faces module with the face_generated_by() function.

makefile: add a "make check" target to run the tests properly.

We need to utter magic words to get the test suite to work, so I've
encoded them in a "make check" target for future reference.

mjo/cone/symmetric_psd.py: fix tests with PYTHONPATH="."

mjo/cone/rearrangement.py: fix tests with PYTHONPATH="."

mjo/cone/permutation_invariant.py: fix tests with PYTHONPATH="."

mjo/cone/doubly_nonnegative.py: fix tests with PYTHONPATH="."

mjo/cone/completely_positive.py: fix tests to work with PYTHONPATH="."

mjo/symbol_sequence.py: fix tests to work with PYTHONPATH="."

mjo/orthogonal_polynomials.py: fix tests to work with PYTHONPATH="."

mjo/matrix_vector.py: make tests work with PYTHONPATH="."

mjo/interpolation.py: allow tests to run with PYTHONPATH=".".

Remove all of the is_foo_on() functions to a sage branch.

Add polynomial ring examples for is_positive_on.

Add another test and an implementation comment for is_lyapunov_like_on.

Add SEEALSO blocks and inexact ring error tests.

Add tests for the "K must be a Cone" TypeError.

Clarify the docstrings for the is_foo_on functions.

Clean up some of the is_foo_on doctests.

Fix the is_positive_on test and give better examples.

Add work-in-progress code for the Schur cone (from David Sossa's work).

Add positive/Z tests and update code for upstream changes.

Add an is_cross_positive() function and implement is_lyapunov_like() using it.

Remove the Pi/Z stuff for inclusion into sage.

Remove references to my unfinished paper.

Implement Z-operators in terms of cross-positive ones.

Add a test for my construction of LL(pi(K,H)).

Begin working on a two-cone pi(K1,K2).

Add LL_cone() convenience function.