Following your instructions exactly, with my current development
branch being 6.4.beta1 (commit 79d9c2f) all I get is:
sage: %runfile test_polynomials.py
sage: test1()
---------------------------------------------------------------------------
ValueError Traceback (most recent call last)
<ipython-input-2-e3097d5bf3e6> in <module>()
----> 1 test1()
/home/jec/sage/test_polynomials.py in test1()
16 # create an extension with at least one root
17 # and convert the polynomial
---> 18 K2 = GF(3**d, 'z')
19 p2 = p1.change_ring(K2)
20
/home/jec/sage/local/lib/python2.7/site-packages/sage/structure/factory.so
in sage.structure.factory.UniqueFactory.__call__
(build/cythonized/sage/structure/factory.c:1244)()
/home/jec/sage/local/lib/python2.7/site-packages/sage/rings/finite_rings/constructor.pyc
in create_key_and_extra_args(self, order, name, modulus, names, impl,
proof, **kwds)
368 order = int(order)
369 if order <= 1:
--> 370 raise ValueError("the order of a finite field
must be > 1.")
371
372 if arith.is_prime(order):
ValueError: the order of a finite field must be > 1.
and (in a separate run)
age: %runfile test_polynomials.py
sage: test2()
---------------------------------------------------------------------------
NameError Traceback (most recent call last)
<ipython-input-2-35ebc1c532fc> in <module>()
----> 1 test2()
/home/jec/sage/test_polynomials.py in test2()
47 """
48 R = PolynomialRing(GF(3), 'T')
---> 49 K2 = GF(3**degree, 'z')
50
51 N = 0
NameError: global name 'degree' is not defined
-- is your test file correct?
John
On 26 August 2014 23:25, Vincent Delecroix <20100.delecr...@gmail.com> wrote:
> Hello,
>
> Working on #16681, it appears that there is a problem with
> factorization of polynomials over finite fields... (despite the ticket
> name, the bug has nothing to do with the recently implemented
> algebraic closure of finite fields).
>
> * if you want to help (i.e. run the tests on your computer) *
>
> In the attached file there are two functions test1 and test2. Just
> start a fresh Sage, and do
> {{{
> sage: %runfile test_polynomials.py
> sage: test1()
> }}}
> you should get a lot of complains (let say one each 10secs/20secs).
>
> Then start an other fresh Sage and run
> {{{
> sage: %runfile test_polynomials.py
> sage: test2()
> }}}
> and you should just see message saying that everything is fine (let
> say each 2min).
>
> Please confirm me that you got the same output (and if not, tell me
> the Sage version, how did you obtain it, and what is your operating
> system).
>
> * if you want to debug *
>
> I do not exactly understand where the bug comes from since it fails
> once but not twice (calling test1() you should see a lot of ""calling
> once more p.factor().prod() gives back the polynomial!!"). The
> factorization is done through pari at lines 3370-3373 of
> sage/rings/polynomial/polynomial_element.pyx
> {{{
> elif is_FiniteField(R):
> v = [x._pari_("a") for x in self.list()]
> f = pari(v).Polrev()
> G = list(f.factor())
> }}}
> So I suspect something is wrong in between Sage and pari about
> conversions between finite fields.
>
> Please help!
> Vincent
>
> --
> You received this message because you are subscribed to the Google Groups
> "sage-devel" group.
> To unsubscribe from this group and stop receiving emails from it, send an
> email to sage-devel+unsubscr...@googlegroups.com.
> To post to this group, send email to sage-devel@googlegroups.com.
> Visit this group at http://groups.google.com/group/sage-devel.
> For more options, visit https://groups.google.com/d/optout.
--
You received this message because you are subscribed to the Google Groups
"sage-devel" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to sage-devel+unsubscr...@googlegroups.com.
To post to this group, send email to sage-devel@googlegroups.com.
Visit this group at http://groups.google.com/group/sage-devel.
For more options, visit https://groups.google.com/d/optout.
