Mathematica 9 is now available
Services & Resources / Wolfram Forums / MathGroup Archive
-----

MathGroup Archive 2008

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Re: ordering of multivariate polynomial terms

  • To: mathgroup at smc.vnet.net
  • Subject: [mg86650] Re: [mg86646] Re: ordering of multivariate polynomial terms
  • From: Andrzej Kozlowski <akoz at mimuw.edu.pl>
  • Date: Sun, 16 Mar 2008 04:39:05 -0500 (EST)
  • References: <frg0d4$evn$1@smc.vnet.net> <200803152246.RAA09977@smc.vnet.net>

Another way to do this is by using the function 
GroebnerBasis`DistributedTermsList, which accepts as an option an 
arbitrary monomial order. So, using the DegreeLexicographic monomial 
order we get:

First[GroebnerBasis`DistributedTermsList[f, {x, y, z},
        MonomialOrder -> DegreeLexicographic]] /.
    {(l_)?VectorQ, m_} :> m*Times @@ ({x, y, z}^l)

{x^4, 4*x^3*y, 6*x^2*y^2, 4*x*y^3, y^4, 5*x^3,
    9*x^2*y, 15*x*y^2, 3*y^3, 7*x^2, 16*x*y, 10*y^2,
    5*x, 4*y, 2}

One can also use this with other monomial orders including user 
defined ones.

Andrzej Kozlowski



On 15 Mar 2008, at 23:46, Szabolcs Horv=E1t wrote:

> Stefan Porubsky wrote:
>> Hello,
>>
>> how it is possible to rearrange the terms of a multivariate 
>> polynomial
>> in Standard or Traditional Form in such a way that the terms are 
>> ordered
>> according to the number of their variables, i.e. first the constant
>> term, then the group of terms depending on one of the variables, then
>> the group of terms depening on two variables, then the group of those
>> with three variables, etc.?
>>
>
> It is not possible to do this because Plus has the Attribute 
> orderless.
>  This means that Mathematica will automatically sort the terms of any
> sum using its built-in ordering rules.
>
> Here's a function that returns a list of the terms sorted by degree:
>
> sortPoly[poly_, vars_] := Module[{x, exp},
>   exp[term_] := Exponent[term /. Alternatives @@ vars -> x, x];
>   Sort[List @@ poly, exp[#1] >= exp[#2] &]
> ]
>
> (The exp[] function returns the degree of a term.)
>
> For example:
>
> f = 2 + 5 x + 7 x^2 + 5 x^3 + x^4 + 4 y + 16 x y + 9 x^2 y + 4 x^3 y =
+
>  10 y^2 + 15 x y^2 + 6 x^2 y^2 + 3 y^3 + 4 x y^3 + y^4
>
> sortPoly[f, {x, y}]
>
> {x^4, 4 x^3 y, 6 x^2 y^2, 4 x y^3, y^4, 5 x^3, 9 x^2 y,
>  15 x y^2, 3 y^3, 7 x^2, 16 x y, 10 y^2, 5 x, 4 y, 2}
>
> Now we can convert the list back to a sum, and use Hold to prevent
> automatic reordering.
>
> Hold@Plus[##] & @@ %
>
> Hold[x^4 + 4 x^3 y + 6 x^2 y^2 + 4 x y^3 + y^4 + 5 x^3 + 9 x^2 y +
>   15 x y^2 + 3 y^3 + 7 x^2 + 16 x y + 10 y^2 + 5 x + 4 y + 2]
>
> ...
>
> I'm not sure if I interpreted your request correctly.  Here's and
> alternative "exp" function that returns the number of different
> variables present in a term:
>
> exp[term_] := Count[FreeQ[term, #] & /@ vars, False]
>
> The full sortPoly:
>
> sortPoly[poly_, vars_] := Module[{x, exp},
>   exp[term_] := Count[FreeQ[term, #] & /@ vars, False];
>   Sort[List @@ poly, exp[#1] <= exp[#2] &]
> ]
>
> Now sortPoly[f, {x,y}] returns
>
> {2, 5 x, 7 x^2, 5 x^3, x^4, 4 y, 10 y^2, 3 y^3, y^4, 16 x y, 9 x^2 y,
>  4 x^3 y, 15 x y^2, 6 x^2 y^2, 4 x y^3}
>



  • Prev by Date: Speeding up FindMinimum setup
  • Next by Date: Re: Operate on Heads
  • Previous by thread: Re: ordering of multivariate polynomial terms
  • Next by thread: RealOnly