[Maxima] Applying vector identities in Maxima
Robert Dodier
robert.dodier at gmail.com
Fri May 19 00:02:11 CDT 2006
Hi Neilen,
> I'm trying to use the vect pacakge to apply some vector identies to an
> expression. Consider
>
> curl(l1*grad(l2)),
>
> where l1 and l2 are scalar functions.
OK, it looks like vect.mac knows only a few identities.
There are probably other problems as well; the demo
script doesn't appear to yield the expected results
(although it's hard to tell what should be expected).
Anyway, I find that the following ...
load (vect);
matchdeclare (xx, scalar_not1p, yy, nonscalarp, [aa, bb], all);
scalar_not1p (e) := e # 1 and scalarp (e);
simp : false;
tellsimpafter (curl (xx*yy), xx*(curl yy) - yy~(grad xx));
tellsimpafter ((xx * yy) ~ aa, xx * (yy ~ aa));
tellsimpafter (aa ~ (xx * yy), xx * (aa ~ yy));
simp : true;
declare (["grad", "laplacian"], nonscalar);
... equips Maxima to handle the specific problem you
mentioned. A couple of identities are represented as
simplification rules. Writing such rules is something of
an art -- let me know if you want to go into details.
Here is an example session.
(%i1) announce_rules_firing : true;
(%o1) true
(%i2) load ("./vect.mac");
By expressrule1 , express(etrue) --> express1(etrue)
(%o2) ./vect.mac
(%i3) load ("./morevect.mac");
yy xx partitions `product'
yy xx partitions `product'
yy xx partitions `product'
(%o3) ./morevect.mac
(%i4) declare ([s, t], scalar);
(%o4) done
(%i5) curl (s * grad t);
By curlrule2 , curl (s*grad t) --> s*curl grad t-(-grad s) ~ grad t
(%o5) s curl grad t - (- grad s) ~ grad t
(%i6) ''%;
By ~rule5 , (-grad s) ~ grad t --> -grad s ~ grad t
By curlrule1 , curl grad t --> 0
(%o6) grad s ~ grad t
Note that not all of the rules were applied to get %o5,
so I reevaluated it (i.e. ''%) to apply some more rules.
The global flag announce_rules_firing controls the
"By foorule, bar --> baz" debug output.
HTH
Robert Dodier
More information about the Maxima
mailing list