[Maxima] custom simplification, was: polarform returns a rectangular expression for float argumen

Robert Dodier robert.dodier at gmail.com
Fri Jan 20 00:09:06 CST 2006


On 1/19/06, Barton Willis <willisb at unk.edu> wrote:

> I'd like for the function simplus, for example, to look
> at the property list of non-atom addends to see if there
> was a specific function for adding such objects. If there
> was such a function, simplus would ship the addition
> off to that function. That way simplus could be extended
> without hacking the simplus code. That's all. I don't think
> simplus, simptimes, and ... allow this.

well, the existing tellsimp / tellsimpafter code does something
pretty much like this, except it associates rules with operators
instead of operands.

> (%i3) tellsimp(x * interval(a,b), interval(min(x*a,x*b), max(x*a,x*b)))$

it turns out that tellsimp and tellsimpafter construct different
functions for identical rules -- they don't differ just in putting
new rules before or after the built-in rules. maybe there's
a reason for that.

i like tellsimpafter's approach to + and * -- given a pattern stated
as a two-operand addition, it sweeps up multiple operands which
satisfy the same match predicates as each of the two operands.

so i would get started on interval arithmetic like this:

matchdeclare (ii, intervalp, nn, not_intervalp);
intervalp (e) := not atom (e) and op (e) = 'interval;
not_intervalp (e) := not intervalp (e);

tellsimpafter (ii + nn, FOO (ii, nn));

a + 1;  => a + 1
a + interval (3, 7);  =>  FOO(interval(3, 7), a)
a + interval (3, 7) + b + interval (9, 12) + c;  =>  FOO(interval(9,
12) + interval(3, 7), c + b + a)
interval (3, 7) + interval (9, 12);  =>  FOO(interval(9, 12) +
interval(3, 7), 0)

at this point it's easy to write FOO, since it is known the first
argument is a sum of intervals and the second is a sum of
stuff that isn't intervals. FOO presumably constructs and
returns a new interval.

rules for other operators can be defined similarly.
only rules for + and * make use of the sweep-up matching,
although it makes sense for any n-ary commutative operator.
a slightly burdensome aspect of the rule definition stuff is that
a rule has to be declared for each operator, so a single rule
can't cover all the built-in trig functions, for example.

maybe the notation could be less obscure -- i'd like to see
ii + nn --> FOO (ii, nn) assuming intervalp (ii), not_intervalp (nn);
but that's easy to arrange with macros.

here is a link to the reference manual,
http://maxima.sourceforge.net/docs/manual/en/maxima_38.html
and some additional verbiage.
http://maxima.sourceforge.net/wiki/index.php/User-defined%20simplification
http://maxima.sourceforge.net/wiki/index.php/tellsimpafter%20worked%20example%20%28base%2060%20arithmetic%29

sorry to run on so long -- hth,

robert




More information about the Maxima mailing list