PARI/GP Bug report logs - #1151
rnfequation(,,1) is broken

Package: pari; Maintainer for pari is Aurel Page <aurel.page@normalesup.org>; Source for pari is src:pari.

Reported by: Bill Allombert <Bill.Allombert@math.u-bordeaux1.fr>

Date: Sat, 18 Dec 2010 16:03:01 UTC

Severity: normal

Done: Karim Belabas <Karim.Belabas@math.u-bordeaux1.fr>

Bug is archived. No further changes may be made.

Full log

Message #20 received at 1151-close@pari.math.u-bordeaux.fr (full text, mbox, reply):

Received: (at 1151-close) by pari.math.u-bordeaux.fr; 8 Feb 2011 20:45:46 +0000
From kb@math.u-bordeaux1.fr Tue Feb 08 21:45:46 2011
Received: from smail.math.u-bordeaux1.fr ([147.210.16.22])
	by pari.math.u-bordeaux1.fr with esmtp (Exim 4.69)
	(envelope-from <kb@math.u-bordeaux1.fr>)
	id 1PmuRW-0003tN-Dh
	for 1151-close@pari.math.u-bordeaux.fr; Tue, 08 Feb 2011 21:45:46 +0100
Received: from localhost (localhost [127.0.0.1])
	by smail.math.u-bordeaux1.fr (Postfix) with ESMTP id 805F1512B1;
	Tue,  8 Feb 2011 21:45:41 +0100 (CET)
X-Virus-Scanned: amavisd-new at math.u-bordeaux1.fr
X-Spam-Score: 2.603
X-Spam-Level: **
X-Spam-Status: No, score=2.603 tagged_above=2 required=6 tests=[BAYES_00=-1,
	DNS_FROM_OPENWHOIS=1.13, NO_DNS_FOR_FROM=1.496,
	RCVD_IN_SORBS_DUL=0.877, RDNS_DYNAMIC=0.1]
Received: from smail.math.u-bordeaux1.fr ([127.0.0.1])
	by localhost (smail.math.u-bordeaux1.fr [127.0.0.1]) (amavisd-new, port 10025)
	with LMTP id Rt3o6AkW8qUl; Tue,  8 Feb 2011 21:45:37 +0100 (CET)
Received: from lfant (tal33-1-82-226-197-154.fbx.proxad.net [82.226.197.154])
	(using TLSv1 with cipher AES256-SHA (256/256 bits))
	(No client certificate requested)
	by smail.math.u-bordeaux1.fr (Postfix) with ESMTP id AE167512AF;
	Tue,  8 Feb 2011 21:45:37 +0100 (CET)
Received: from kb by lfant with local (Exim 4.72)
	(envelope-from <kb@koala.org>)
	id 1PmuTI-0004xW-K9; Tue, 08 Feb 2011 21:47:36 +0100
Date: Tue, 8 Feb 2011 21:47:36 +0100
From: Karim Belabas <Karim.Belabas@math.u-bordeaux1.fr>
To: Bill Allombert <Bill.Allombert@math.u-bordeaux1.fr>,
	1151-close@pari.math.u-bordeaux.fr
Subject: Re: Bug#1151: rnfequation(,,1) is broken
Message-ID: <20110208204736.GA18741@math.u-bordeaux1.fr>
References: <20101218155503.GC9570@yellowpig>
 <20110131235925.GE3266@yellowpig>
 <20110201162149.GF3266@yellowpig>
MIME-Version: 1.0
Content-Type: text/plain; charset=us-ascii
Content-Disposition: inline
In-Reply-To: <20110201162149.GF3266@yellowpig>
User-Agent: Mutt/1.5.20 (2009-06-14)

* Bill Allombert [2011-02-01 17:57]:
> On Tue, Feb 01, 2011 at 12:59:25AM +0100, Bill Allombert wrote:
> > On Sat, Dec 18, 2010 at 04:55:03PM +0100, Bill Allombert wrote:
> > > Package: pari
> > > Version: 2.4.3, 2.3.5
> > > 
> > > Hello Karim,
> > > This issue was reported on the pari-users list
> > > <http://pari.math.u-bordeaux.fr/archives/pari-users-1012/msg00004.html>.
> > > I report it for reference.
> > > 
> > > ? f = (y^17-1)/(y-1);
> > > ? g = Mod(1,f)*x+Mod(y^6,f);
> > > ? rnfequation(f,g,1);
> > >   ***   at top-level: rnfequation(f,g,1)
> > >   ***                 ^------------------
> > >   *** rnfequation: the PARI stack overflows !
> > 
> > Well, according to my calculation the correct LERS is:
> > (-x^4 + x^3 - x^2)*y + (-x^7 + x^6 - x^5)
> > 
> > Furthermore, the values returned by Flx_resultant_all for x=2 are 144 times too large.
> > For x=3 the values are random but the quotient H0/H1 seems correct
> 
> The attached patch is probably incorrect, but it allows both 'make test-all' and 
> this test
> for(i=1,16,print(i);rnfequation(polcyclo(17,y),-x+y^i,1))
> 
> to work. Strange...

Fixed in svn, using quite different formulas (cf Algorithm 3.3.7 in GTM138).

Thanks !

    K.B.
--
Karim Belabas, IMB (UMR 5251)  Tel: (+33) (0)5 40 00 26 17
Universite Bordeaux 1          Fax: (+33) (0)5 40 00 69 50
351, cours de la Liberation    http://www.math.u-bordeaux1.fr/~belabas/
F-33405 Talence (France)       http://pari.math.u-bordeaux1.fr/  [PARI/GP]
`

Send a report that this bug log contains spam.

Bill Allombert <allomber@math.u-bordeaux.fr>. Last modified: Sat Aug 26 17:27:12 2023; Machine Name: pari

PARI/GP Bug tracking system

Debbugs is free software and licensed under the terms of the GNU Public License version 2. The current version can be obtained from https://bugs.debian.org/debbugs-source/.

Copyright © 1999 Darren O. Benham, 1997,2003 nCipher Corporation Ltd, 1994-97 Ian Jackson, 2005-2017 Don Armstrong, and many other contributors.