MathGroup Archive 1998

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

Search the Archive

Bug in PolynomialReduce

  • To: mathgroup at smc.vnet.net
  • Subject: [mg15005] Bug in PolynomialReduce
  • From: "Greg Arnold" <Gregory.Arnold at sensors.wpafb.af.mil>
  • Date: Sat, 5 Dec 1998 01:30:10 -0500
  • Sender: owner-wri-mathgroup at wolfram.com

Hello all,

I'm trying to use GroebnerBasis and PolynomialReduce to simplify
equations which I'm trying to integrate.  My general problem is that I
want to be able to substitute for complex functions.  However, I'm not
getting the expected result even for this fairly simple case:

In[46]:=
gb=GroebnerBasis[{r1==z/y,r2==x/y},{x,y,z,r1,r2}]; gb // InputForm

Out[46]//InputForm=
{r1*y - z, -r1 + z/y, x - r2*y}

In[48]:=
Length[gb]

Out[48]=
3

In[53]:=
pr=PolynomialReduce[ z/y+x/y,gb,{x,y,z,r1,r2}]; pr // InputForm

Out[53]//InputForm=
{{0, 0, y^(-1), 0}, (r2*y + z)/y}


(1) I expect to have only 3 coefficients (not 4) in the first part, and
(2) the residual (the last part) should reduce to r1+r2

Am I doing something wrong?  I'm running 3.0.0.0 WinNT.  Is there a bug
that has been fixed in 3.0.1?

P.S.  I did notice that if I choose the order {x,z,y,r1,r2}, then I will
get the expected result for both (1) and (2)!.  However, I cannot think
of a way to decide which way(s) will work and which will not.

Thanks in advance!

Greg



  • Prev by Date: Re: Contour Integration
  • Next by Date: Re: Re: Way to evaluate D[(1-x^2)y''[x],{x,n
  • Previous by thread: Re: error from MathLink, plotting results from a c++ simulation
  • Next by thread: Re: Bug in PolynomialReduce