PolynomialIdeals
EquidimensionalDecomposition
Decompose an ideal into ideals of distinct dimension
Calling Sequence
Parameters
Description
Examples
References
Compatibility
EquidimensionalDecomposition(J)
J
-
a polynomial ideal
The EquidimensionalDecomposition command computes a sequence of ideals of distinct Hilbert dimension whose intersection is equal to the original ideal. Assuming there are no embedded primes, the prime components of each ideal in the sequence have the same dimension also. In general this decomposition is not unique.
This function is part of the PolynomialIdeals package, and can be used in the form EquidimensionalDecomposition(..) only after executing the command with(PolynomialIdeals). However, it can always be accessed through the long form of the command using PolynomialIdeals[EquidimensionalDecomposition](..).
In the example below, the variety is composed of five points (dimension zero), three curves (dimension one), and one surface (dimension two). The equidimensional decomposition places all of the points, all of the curves, and all of the surfaces into separate ideals.
with⁡PolynomialIdeals:
J≔z4⁢x3+z4⁢y2−z5−z4⁢x−z⁢x4−x⁢y2⁢z+z2⁢x+z⁢x2,x2⁢y2−z⁢x2+z⁢x−y2⁢z−x3⁢z−y⁢z+x3⁢y−x⁢y+x5−x3+y3+z2,z3⁢x4+z3⁢x⁢y2−z3⁢x3−x5−y2⁢z3−z3⁢x2−z4⁢x−x2⁢y2+x4+z3⁢x+z4+x⁢y2+x3+z⁢x2−x2−z⁢x:
E≔EquidimensionalDecomposition⁡J
E≔z2,y2⁢z+y⁢z,y6+2⁢y5+y4,x3⁢z⁢y−x⁢y⁢z+y⁢z,x4+x3⁢y−x3⁢z+x⁢y2+y3−x2−x⁢y,x3⁢y2+x3⁢y−x3⁢z+y4−x⁢y2+y3−x⁢y+x⁢z+y⁢z,y⁢z6+z6+y2−y⁢z+y−z,x3⁢z6−x⁢z6−z7+z6+x3⁢y−x3⁢z−x⁢y+x⁢z−y⁢z+z2+y−z,y⁢z9+z9−x⁢y⁢z6−x⁢z6−x3⁢z3−z4⁢y+x4+x⁢z3+y⁢z3+x⁢y⁢z−x2−x⁢y,x3+y2−x−z
seq⁡HilbertDimension⁡i,i=E
0,1,2
seq⁡map⁡HilbertDimension,PrimaryDecomposition⁡i,x,y,z,i=E
0,0,0,0,0,1,1,1,2
The next example illustrates what happens when embedded primes are present.
K≔x2−y,x3−y⁢z⁢w
K≔x2−y,−y⁢z⁢w+x3
E≔EquidimensionalDecomposition⁡K
E≔w,y2,x⁢y,x2−y,x2−y,y⁢z⁢w−x⁢y
seq⁡HilbertDimension⁡i,w,x,y,z,i=E
1,2
seq⁡map⁡HilbertDimension,PrimaryDecomposition⁡i,w,x,y,z,i=E
1,2,2,1
Becker, T., and Weispfenning, V. Groebner Bases. Springer-Verlag, 1993.
The PolynomialIdeals[EquidimensionalDecomposition] command was updated in Maple 16.
See Also
map
PolynomialIdeals[HilbertDimension]
PolynomialIdeals[Intersect]
PolynomialIdeals[PrimeDecomposition]
PolynomialIdeals[Simplify]
PolynomialIdeals[ZeroDimensionalDecomposition]
Download Help Document