simplify/siderels
simplify with respect to side relations
Calling Sequence
Parameters
Description
Examples
References
simplify(expr, eqns)
simplify(expr, eqns, vars)
simplify(expr, eqns, tord)
simplify(expr, eqns, 'mindeg')
expr
-
expression
eqns
set or list of equations (an expression e is understood as the equation e=0)
vars
(optional) set or list of variables appearing in eqns
tord
(optional) a monomial order in the variables appearing in eqns
Simplification of expr with respect to the side relations eqns is performed. The result is an expression which is mathematically equivalent to expr but which is in normal form with respect to the specified side relations.
In all cases simplify computes a Groebner basis for the set of side relations. This basis contains all the algebraic consequences of eqns, so that if an expression can be rewritten in terms of the relations simplify will reduce it to zero. Otherwise, simplify returns a remainder in which none of the terms are reducible.
To ensure a finite division process simplify imposes a total order on terms which can appear. This order, called a monomial order, is used for both the Groebner basis computation and the division. The monomial order can be specified by an optional third argument, and different orders can drastically affect a non-zero result:
If a list of variables vars is specified then pure lexicographic (plex) order is used and the leftmost variables are eliminated first.
If a set of variables vars is given then a total degree (tdeg) order is used and terms with the highest degree in all the variables are eliminated first. Terms with the same total degree are ordered heuristically.
If no third argument was specified then lexicographic order is used and variables appearing on the left hand side of any equation in eqns are eliminated first.
If a monomial order tord is specified then that order is used for all computations. See Groebner[MonomialOrders] for a list of the available monomial orders.
If expr is a rational expression then one of three algorithms can be used. If a list or set of variables vars is specified or if the third argument is omitted then simplification is applied to the numerator and denominator separately. This method is very fast, however no attempt is made to simplify the ratio of the numerator and denominator so simpler representations of the fraction may be missed.
If a monomial order tord is specified then the largest monomial appearing in the numerator and denominator is minimized with respect to tord. The result is a simplified canonical form with respect to the monomial order. Note that in some cases a common factor may be introduced to satisfy this requirement, as in example 2 below. If the denominator is invertible or if it divides the numerator modulo the side relations then a quotient is computed and it is returned instead.
The optional argument mindeg minimizes the sum of the degrees of the resulting fraction. It uses a global search to find an equivalent fraction ab with dega+degb minimal. The result is not canonical, but in some cases the output may be preferable to a canonical form, as in example 3 below. In general the search is exponential time in dega+degb.
If the expression is neither a polynomial nor a rational function in the variables, simplify will recurse on the input expr and apply the normalization to any subexpressions which are rational functions in the variables.
Example 1: trigonometric polynomials
eqns≔sinx2+cosx2=1:
e≔sinx3−11sinx2cosx+3cosx3−sinxcosx+2
simplifye,eqns
14cosx3−cosx2sinx+−sinx−11cosx+sinx+2
simplifye,eqns,sinx,cosx
simplifye,eqns,cosx,sinx
−14sinx2−sinx+3cosx+sinx3+2
f≔8sinx4cosx+15sinx2cosx3−15sinx2cosx+7cosx5−14cosx3+7cosx
simplifyf,eqns
0
r≔1−cosx2+sinxcosxsinxcosx+cosx2
simplifyr,eqns
1−cosx2+sinxcosxcosxcosx+sinx
simplifyr,eqns,mindeg
sinxcosx
Example 2: reducing fractions to a canonical form
eqns2≔xy5−x−y
r≔y5+x+yx−y
simplifyr,eqns2,tdegx,y
x+1x+yxx−y
simplifyr,eqns2,plexx,y
−y4y5+y−1y5−2
simplifyr,eqns2,mindeg
r2≔−x3−yx2−x4−yx3+2xy4+2y4+y5+x+yyx3−x4+2xy4−4+x−y
r2≔2xy4+y5−x4−yx3+2y4−x3−yx2+x+y2xy4−x4+yx3+x−y−4
simplifyr2,eqns2
−2xy4−y5+x4+yx3−2y4+x3+yx2−x−y−2xy4+x4−yx3−x+y+4
simplifyr2,eqns2,tdegx,y
Example 3: the canonical form algorithm may increase total degree
eqns3≔x5+xy−1
r≔x3y3−x4+x−1x2−y2+1
simplifyr,eqns3,tdegx,y
x3+y−1x2+−y4−1x+y3x2x2−y2+1
simplifyr,eqns3,mindeg
x3y3−x4+x−1x2−y2+1
Example 4: the side relations define a field
eqns4≔x+y=z,z3−z2−zy+2y2=1,z3+y2=1,z2+zy−y2=0
eqns4≔x+y=z,z3+y2=1,−y2+zy+z2=0,z3+2y2−zy−z2=1
h1≔36z4y2+36zy4−36zy2−12z2+zy−12y2+12xz−12xy+23z4+43z3y−23z2y2−43zy3+23y4
simplifyh1,eqns4
h2≔zy2+z3y3−3z2y+233+56z2+73zy−116y2+12xz−12xy+23z4+43z3y−23z2y2−43zy3+23y4
h2≔zy2+y3z3−3z2y+233+56z2+73zy−116y2+12xz−12xy+23z4+43z3y−23z2y2−43zy3+23y4
simplifyh2,eqns4,x,y,z
−343+28−17yz2+24+7yz+13y
simplifyh2,eqns4,lexdegx,z,y
293−y5−2y2
h2inv≔simplify1h2,eqns4,tdegx,z,y
h2inv≔447375252342772565971+252548433y342772565971+815384160y242772565971+935711723+4097948zy42772565971+236183994z42772565971
simplifyh2h2inv,eqns4
1
Example 5: a difficult example that can not be handled by other methods
eqns5≔y2−x3+xy−1
eqns5≔−x3+xy+y2−1
algcurvesparametrizationopeqns5,x,y,t
Error, (in `algcurves/ratpar`) genus is not zero
r≔−y3+1−2yx2−x2y5−2xy2−y6x−y7+y6y3−2y2+1
r≔−x2y5−xy6−y7+y6−2yx2−2y2x−y3+1y3−2y2+1
simplifyr,eqns5,tdegx,y
y5+x+yx−y
Cox, D.; Little, J.; and O'Shea, D. Ideals, Varieties, and Algorithms. 2nd ed. Springer-Verlag, 1997.
Monagan, M., and Pearce, R. "Rational Simplification Modulo a Polynomial Ideal." ISSAC 2006 proceedings.
See Also
algcurves[parametrization]
Groebner[Basis]
Groebner[MonomialOrders]
Groebner[NormalForm]
Download Help Document