Re: [sage-devel] potential bug for calculating Heegner points

2022-07-06 Thread John Cremona
On Wed, 6 Jul 2022 at 10:28, John Cremona wrote: > > The problem is in lines 3373-3374 and 3382 of > src/sage/schemes/elliptic_curves/heegner.py. The floating point > approximation of the point has x-coordinate 0.500 - > 1.32287565553229*I and the code tries to find an algebraic numbe

Re: [sage-devel] potential bug for calculating Heegner points

2022-07-06 Thread John Cremona
The problem is in lines 3373-3374 and 3382 of src/sage/schemes/elliptic_curves/heegner.py. The floating point approximation of the point has x-coordinate 0.500 - 1.32287565553229*I and the code tries to find an algebraic number to match (i.e. a polynomial over QQ with it as approximate

Re: [sage-devel] potential bug for calculating Heegner points

2022-07-06 Thread John Cremona
On Tue, 5 Jul 2022, 23:21 'Justin C. Walker' via sage-devel, < sage-devel@googlegroups.com> wrote: > > > > On Jul 5, 2022, at 07:00 , Debanjana wrote: > > > > sage: E = EllipticCurve('11a1') > > sage: P = E.heegner_point(-7) > > sage: t = P.point_exact() > > sage: t.domain() > > Spectrum of Numbe

Re: [sage-devel] potential bug for calculating Heegner points

2022-07-05 Thread 'Justin C. Walker' via sage-devel
> On Jul 5, 2022, at 07:00 , Debanjana wrote: > > sage: E = EllipticCurve('11a1') > sage: P = E.heegner_point(-7) > sage: t = P.point_exact() > sage: t.domain() > Spectrum of Number Field in a with defining polynomial x^2 + x + 20 > sage: t.domain().base_ring().discriminant() > -79 > > The an

[sage-devel] potential bug for calculating Heegner points

2022-07-05 Thread Debanjana
sage: E = EllipticCurve('11a1') sage: P = E.heegner_point(-7) sage: t = P.point_exact() sage: t.domain() Spectrum of Number Field in a with defining polynomial x^2 + x + 20 sage: t.domain().base_ring().discriminant() -79 The answer should be -7 but Sage gives -79 -- A smile is a curve that can s