#10767: Fractional ideals intersection gives wrong awnsers in some cases.
-----------------------------+----------------------------------------------
Reporter: mderickx | Owner: was
Type: defect | Status: new
Priority: major | Milestone: sage-4.7
Component: number theory | Keywords:
Work_issues: | Upstream: N/A
Reviewer: | Author:
Merged: | Dependencies:
-----------------------------+----------------------------------------------
Comment(by johanbosman):
This seems to be a bug in Pari:
{{{
GP/PARI CALCULATOR Version 2.4.3 (development svn-12623)
i386 running darwin (x86-64/GMP-4.2.1 kernel) 64-bit version
compiled: Mar 4 2011, gcc-4.2.1 (Apple Inc. build 5664)
(readline v6.1 enabled, extended help enabled)
Copyright (C) 2000-2008 The PARI Group
PARI/GP is free software, covered by the GNU General Public License, and
comes
WITHOUT ANY WARRANTY WHATSOEVER.
Type ? for help, \q to quit.
Type ?12 for how to get moral (and possibly technical) support.
parisize = 8000000, primelimit = 500509
? K = nfinit(x^2-2)
%1 = [x^2 - 2, [2, 0], 8, 1, [[1,
-1.4142135623730950488016887242096980786; 1,
1.4142135623730950488016887242096980786], [1,
-1.4142135623730950488016887242096980786; 1,
1.4142135623730950488016887242096980786], [1, -1; 1, 1], [2, 0; 0, 4], [4,
0; 0, 2], [2, 0; 0, 1], [2, [0, 2; 1, 0]]],
[-1.4142135623730950488016887242096980786,
1.4142135623730950488016887242096980786], [1, x], [1, 0; 0, 1], [1, 0, 0,
2; 0, 1, 1, 0]]
? a = idealhnf(K, 1/2)
%2 =
[1/2 0]
[0 1/2]
? idealintersect(K, a, a)
%3 =
[1/4 0]
[0 1/4]
}}}
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/10767#comment:1>
Sage <http://www.sagemath.org>
Sage: Creating a Viable Open Source Alternative to Magma, Maple, Mathematica,
and MATLAB
--
You received this message because you are subscribed to the Google Groups
"sage-trac" group.
To post to this group, send email to [email protected].
To unsubscribe from this group, send email to
[email protected].
For more options, visit this group at
http://groups.google.com/group/sage-trac?hl=en.
