Maxima function ordmexpt

This project aims to develop a new version of the Maxima function ordmexpt, which is a subroutine for Maxima's ordering predicate great. This rewrite attempts to clarify the logic of the function and to fix bug #4383: great not transitive (so simplifya not idempotent). Specifically, bug #4383 causes the following expression to not simplify correctly

(%i1)	(declare(z,complex), domain : complex)$
(%i2)	exp(z) + sqrt(exp(z)) + exp(z);
(%o2)	%e^z+sqrt(%e^z)+%e^z

Maxima's ordering predicate great is closely linked to simplification, so tweaking great can cause syntactic testsuite failures. To help distinguish syntactic from semantic failures, this project includes a modified function approx-alike that applies several transformations to both the expected and actual outputs. This modified function is not intended as a replacement to approx-alike.

The last I tried, my alternative ordmexpt function along with the modified function approx-alike runs the testsuite, including the share testsuite, with thirty-one unexpected failures. A typical failure is (test done with 'domain : 'complex')

********************** Problem 537 (line 2960) ***************
Input:
                         - 3
                         ───
               2 c z      2    c z
integrate((b %e      + a)    %e   , z)


Result:
            2 c z
     sqrt(%e     )
───────────────────────
             2 c z
a c sqrt(b %e      + a)

This differed from the expected result:
           c z
         %e
───────────────────────
             2 c z
a c sqrt(b %e      + a)

The file rtest_great tests the function great. There is another test file (rtest_shame) of serious bugs that I have collected.

Running the testsuite calls ordmexpt about 2.9 million times, so we need to be concerned with its efficiency. In addition to fixing bug #4383, the modified ordmexpt function fixes two testsuite failures-they are rtest1, #183 and rtest_limit_extra #259; the limit bug is

Running tests in rtest_limit_extra:
********************** Problem 259 (line 909) ***************
Input:
                 - a
                 ───
              a   2
limit((%i + 1)  2   , a, inf)


Result:
ind

... Which was correct, but was expected to be wrong due to a known bug in Maxima or SBCL.

This bug fix is not related to the extra simplifications in approx-alike. Another bug that this code fixes is (see #4455 'great' errors)

(%i2)	t^c + sqrt(t^c) + 42;
(%o2)	t^c+sqrt(t^c)+42

(%i3)	expand(%,0,0);
(%o3)	sqrt(t^c)+t^c+42

(%i4)	expand(%,0,0);
(%o4)	t^c+sqrt(t^c)+42