Targeted Simplify in Mathematica
Posted
by Timo
on Stack Overflow
See other posts from Stack Overflow
or by Timo
Published on 2009-11-17T15:48:38Z
Indexed on
2010/04/07
14:53 UTC
Read the original article
Hit count: 321
mathematica
I generate very long and complex analytic expressions of the general form:
(...something not so complex...)(...ditto...)(...ditto...)...lots...
When I try to use Simplify
, Mathematica grinds to a halt, I am assuming due to the fact that it tries to expand the brackets and or simplify across different brackets. The brackets, while containing long expressions, are easily simplified by Mathematica on their own. Is there some way I can limit the scope of Simplify
to a single bracket at a time?
Edit: Some additional info and progress.
So using the advice from you guys I have now started using something in the vein of
In[1]:= trouble = Log[(x + I y) (x - I y) + Sqrt[(a + I b) (a - I b)]];
In[2]:= Replace[trouble, form_ /; (Head[form] == Times) :> Simplify[form],{3}]
Out[2]= Log[Sqrt[a^2 + b^2] + (x - I y) (x + I y)]
Changing Times
to an appropriate head like Plus
or Power
makes it possible to target the simplification quite accurately. The problem / question that remains, though, is the following: Simplify
will still descend deeper than the level specified to Replace
, e.g.
In[3]:= Replace[trouble, form_ /; (Head[form] == Plus) :> Simplify[form], {1}]
Out[3]= Log[Sqrt[a^2 + b^2] + x^2 + y^2]
simplifies the square root as well.
My plan was to iteratively use Replace
from the bottom up one level at a time, but this clearly will result in vast amount of repeated work by Simplify
and ultimately result in the exact same bogging down of Mathematica I experienced in the outset. Is there a way to restrict Simplify
to a certain level(s)?
I realize that this sort of restriction may not produce optimal results, but the idea here is getting something that is "good enough".
© Stack Overflow or respective owner