Re: [sage-release] Re: Sage 8.4.beta7 released

2018-10-03 Thread Kwankyu Lee


On Wednesday, October 3, 2018 at 10:01:40 PM UTC+9, Samuel Lelievre wrote:
>
> Wed 2018-10-03 12:36 UTC+2, Travis Scrimshaw: 
> > 
> > On Wednesday, October 3, 2018 at 7:35:08 PM UTC+10, Kwankyu Lee wrote: 
> >> 
> >> On Wednesday, October 3, 2018 at 3:14:32 PM UTC+9, Travis Scrimshaw 
> wrote: 
> >>> 
> >>> I believe there is a Heisenbug in 
> src/sage/rings/function_field/ideal.py. 
>
> >>> Here is a specific test that can sometimes fail: 
> >>> 
> >>> sage: K. = FunctionField(GF(3^2)); R. = K[] 
> >>> sage: F. = K.extension(t^3 + t^2 - x^4) 
> >>> sage: Oinf = F.maximal_order_infinite() 
> >>> sage: I = Oinf.ideal(1/y) 
> >>> sage: I.factor()   # This will sometimes fail 
> >>> 
> >>> I suspect people without meataxe will not have this failure. 
>
> Is this file added by some optional package, such as meataxe? 
>
> Without meataxe, the directory src/sage/rings/function_field 
> has a file function_field_ideal.py but no file ideal.py --- and 
>
> sage: Oinf = F.maximal_order_infinite() 
>
> gives a NotImplementedError. 
>
> > I created https://trac.sagemath.org/ticket/26389 
> > (Simon, I have cc-ed you). 
>
> Note that the ticket summary has "Hisenbug" for "Heisenbug" 
> and it might make sense to start the ticket description with 
> an indication that this happens with Meataxe installed. 
>

This is on Sage 8.4.beta 7. 

-- 
You received this message because you are subscribed to the Google Groups 
"sage-release" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to sage-release+unsubscr...@googlegroups.com.
To post to this group, send email to sage-release@googlegroups.com.
Visit this group at https://groups.google.com/group/sage-release.
For more options, visit https://groups.google.com/d/optout.


Re: [sage-release] Re: Sage 8.4.beta7 released

2018-10-03 Thread Samuel Lelièvre
Wed 2018-10-03 12:36 UTC+2, Travis Scrimshaw:
>
> On Wednesday, October 3, 2018 at 7:35:08 PM UTC+10, Kwankyu Lee wrote:
>>
>> On Wednesday, October 3, 2018 at 3:14:32 PM UTC+9, Travis Scrimshaw wrote:
>>>
>>> I believe there is a Heisenbug in src/sage/rings/function_field/ideal.py.

>>> Here is a specific test that can sometimes fail:
>>>
>>> sage: K. = FunctionField(GF(3^2)); R. = K[]
>>> sage: F. = K.extension(t^3 + t^2 - x^4)
>>> sage: Oinf = F.maximal_order_infinite()
>>> sage: I = Oinf.ideal(1/y)
>>> sage: I.factor()   # This will sometimes fail
>>>
>>> I suspect people without meataxe will not have this failure.

Is this file added by some optional package, such as meataxe?

Without meataxe, the directory src/sage/rings/function_field
has a file function_field_ideal.py but no file ideal.py --- and

sage: Oinf = F.maximal_order_infinite()

gives a NotImplementedError.

> I created https://trac.sagemath.org/ticket/26389
> (Simon, I have cc-ed you).

Note that the ticket summary has "Hisenbug" for "Heisenbug"
and it might make sense to start the ticket description with
an indication that this happens with Meataxe installed.

-- 
You received this message because you are subscribed to the Google Groups 
"sage-release" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to sage-release+unsubscr...@googlegroups.com.
To post to this group, send email to sage-release@googlegroups.com.
Visit this group at https://groups.google.com/group/sage-release.
For more options, visit https://groups.google.com/d/optout.


[sage-release] Re: Sage 8.4.beta7 released

2018-10-03 Thread tscrim


On Wednesday, October 3, 2018 at 7:35:08 PM UTC+10, Kwankyu Lee wrote:
>
>
>
> On Wednesday, October 3, 2018 at 3:14:32 PM UTC+9, tsc...@ucdavis.edu 
> wrote:
>>
>> I believe there is a Heisenbug in src/sage/rings/function_field/ideal.py. 
>> Here is a specific test that can sometimes fail:
>>
>> sage: K. = FunctionField(GF(3^2)); R. = K[]
>> sage: F. = K.extension(t^3 + t^2 - x^4)
>> sage: Oinf = F.maximal_order_infinite()
>> sage: I = Oinf.ideal(1/y)
>> sage: I.factor()   # This will sometimes fail
>>
>> ---
>> TypeError Traceback (most recent call 
>> last)
>>
>> I suspect people without meataxe will not have this failure. 
>>
>
> Perhaps true. I must be the person who ran this test most number of times 
> and I never encountered this failure. 
>

Okay, that is good to know. Do you know offhand (i.e., don't spend much 
time on this) of a simple way I can extract out the relevant matrix 
computation? This is not a minimal example, and it would be useful to get 
it down smaller.

I created https://trac.sagemath.org/ticket/26389 (Simon, I have cc-ed you).

Best,
Travis

-- 
You received this message because you are subscribed to the Google Groups 
"sage-release" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to sage-release+unsubscr...@googlegroups.com.
To post to this group, send email to sage-release@googlegroups.com.
Visit this group at https://groups.google.com/group/sage-release.
For more options, visit https://groups.google.com/d/optout.


[sage-release] Re: Sage 8.4.beta7 released

2018-10-03 Thread Kwankyu Lee


On Wednesday, October 3, 2018 at 3:14:32 PM UTC+9, tsc...@ucdavis.edu wrote:
>
> I believe there is a Heisenbug in src/sage/rings/function_field/ideal.py. 
> Here is a specific test that can sometimes fail:
>
> sage: K. = FunctionField(GF(3^2)); R. = K[]
> sage: F. = K.extension(t^3 + t^2 - x^4)
> sage: Oinf = F.maximal_order_infinite()
> sage: I = Oinf.ideal(1/y)
> sage: I.factor()   # This will sometimes fail
> ---
> TypeError Traceback (most recent call last
> )
>
> I suspect people without meataxe will not have this failure. 
>

Perhaps true. I must be the person who ran this test most number of times 
and I never encountered this failure. 

-- 
You received this message because you are subscribed to the Google Groups 
"sage-release" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to sage-release+unsubscr...@googlegroups.com.
To post to this group, send email to sage-release@googlegroups.com.
Visit this group at https://groups.google.com/group/sage-release.
For more options, visit https://groups.google.com/d/optout.


[sage-release] Re: Sage 8.4.beta7 released

2018-10-03 Thread Simon King
Hi Travis,

On 2018-10-03, tsc...@ucdavis.edu  wrote:
> In the file, I can get anywhere between 0 and 6 tests failing. Note that 
> once this works, it never seems to fail. Also, inserting output into the 
> error, n = 255 and self.order() = 9. I suspect the problem is in 
> Matrix_gfpn_dense as this is not a problem when I go to GF(19^2). I should 
> note that I do have meataxe installed, and I suspect people without meataxe 
> will not have this failure.

When you open a ticket for it, please Cc me.

Best regards,
Simon

-- 
You received this message because you are subscribed to the Google Groups 
"sage-release" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to sage-release+unsubscr...@googlegroups.com.
To post to this group, send email to sage-release@googlegroups.com.
Visit this group at https://groups.google.com/group/sage-release.
For more options, visit https://groups.google.com/d/optout.


[sage-release] Re: Sage 8.4.beta7 released

2018-10-03 Thread tscrim
I believe there is a Heisenbug in src/sage/rings/function_field/ideal.py. 
Here is a specific test that can sometimes fail:

sage: K. = FunctionField(GF(3^2)); R. = K[]
sage: F. = K.extension(t^3 + t^2 - x^4)
sage: Oinf = F.maximal_order_infinite()
sage: I = Oinf.ideal(1/y)
sage: I.factor()   # This will sometimes fail
---
TypeError Traceback (most recent call last)
 in ()
> 1 I.factor()

/home/uqtscrim/sage/local/lib/python2.7/site-packages/sage/rings/
function_field/ideal.pyc in factor(self)
   2671 defined by y^2 + y + (x^2 + 1)/x)^2
   2672 """
-> 2673 return Factorization(self._factor(), cr=True)
   2674 
   2675 def _factor(self):

/home/uqtscrim/sage/local/lib/python2.7/site-packages/sage/rings/function_field/ideal.pyc
 
in _factor(self)
   2693 
   2694 factors = []
-> 2695 for iprime, exp in O._to_iF(self).factor():
   2696 prime = FunctionFieldIdealInfinite_global(O, iprime)
   2697 factors.append((prime, exp))

/home/uqtscrim/sage/local/lib/python2.7/site-packages/sage/rings/function_field/ideal.pyc
 
in factor(self)
   1572 True
   1573 """
-> 1574 return Factorization(self._factor(), cr=True)
   1575 
   1576 def _factor(self):

/home/uqtscrim/sage/local/lib/python2.7/site-packages/sage/rings/
function_field/ideal.pyc in _factor(self)
   1598 primes = set([o.ideal(p) for p,_ in d.factor()] + [p for p,_ 
in i.ideal_below().factor()])
   1599 for prime in primes:
-> 1600 qs = [q[0] for q in O.decomposition(prime)]
   1601 for q in qs:
   1602 exp = q.valuation(self)

/home/uqtscrim/sage/local/lib/python2.7/site-packages/sage/misc/cachefunc.pyx 
in sage.misc.cachefunc.CachedMethodCaller.__call__ (build/cythonized/sage/
misc/cachefunc.c:10824)()
   1951 return cache[k]
   1952 except KeyError:
-> 1953 w = self._instance_call(*args, **kwds)
   1954 cache[k] = w
   1955 return w

/home/uqtscrim/sage/local/lib/python2.7/site-packages/sage/misc/cachefunc.pyx 
in sage.misc.cachefunc.CachedMethodCaller._instance_call (build/cythonized/
sage/misc/cachefunc.c:10280)()
   1827 True
   1828 """
-> 1829 return self.f(self._instance, *args, **kwds)
   1830 
   1831 cdef fix_args_kwds(self, tuple args, dict kwds):

/home/uqtscrim/sage/local/lib/python2.7/site-packages/sage/rings/function_field/order.pyc
 
in decomposition(self, ideal)
   1838 
   1839 # J_1, J_2, ...
-> 1840 Jb =[Kb[0]] + [div(Kb[j],Kb[j-1]) for j in 
range(1,len(Kb))]
   1841 
   1842 # H_1, H_2, ...

/home/uqtscrim/sage/local/lib/python2.7/site-packages/sage/rings/function_field/order.pyc
 
in div(Ib, Jb)
   1805 sJb = Jb.row_space()
   1806 sIb = Ib.row_space()
-> 1807 sJbsIb,proj_sJbsIb,lift_sJbsIb = 
sJb.quotient_abstract(sIb)
   1808 supplement_basis = [lift_sJbsIb(v) for v in 
sJbsIb.basis()]
   1809 

/home/uqtscrim/sage/local/lib/python2.7/site-packages/sage/modules/free_module.pyc
 
in quotient_abstract(self, sub, check)
   4546 # Calling is_subspace may be way too slow and repeat work 
done below.
   4547 # It will be very desirable to somehow do this step better.
-> 4548 if check and (not is_FreeModule(sub) or not 
sub.is_subspace(self)):
   4549 try:
   4550 sub = self.subspace(sub)

/home/uqtscrim/sage/local/lib/python2.7/site-packages/sage/modules/free_module.pyc
 
in is_subspace(self, other)
   3790 True
   3791 """
-> 3792 return self.is_submodule(other)
   3793 
   3794 def span(self, gens, base_ring=None, check=True, 
already_echelonized=False):

/home/uqtscrim/sage/local/lib/python2.7/site-packages/sage/modules/
free_module.pyc in is_submodule(self, other)
   1458 from sage.misc.flatten import flatten
   1459 return all(x in S for x in flatten(M))
-> 1460 return all(x in S for x in M.list())
   1461 
   1462 def __iter__(self):

/home/uqtscrim/sage/local/lib/python2.7/site-packages/sage/matrix/matrix0.pyx 
in sage.matrix.matrix0.Matrix.list (build/cythonized/sage/matrix/matrix0.c:
4165)()
154 [  y   x 2*x + y]
155 """
--> 156 return list(self._list())
157 
158 def _list(self):

/home/uqtscrim/sage/local/lib/python2.7/site-packages/sage/matrix/matrix_gfpn_dense.pyx
 
in sage.matrix.matrix_gfpn_dense.Matrix_gfpn_dense._list 
(build/cythonized/sage/matrix/matrix_gfpn_dense.c:9038)()
908 FfStepPtr(&(p))
909 sig_check()
--> 910 
x.extend([self._converter.int_to_field(FfToInt(FfExtract(p,j))) for j in 
range(self.Data.Noc)])
911 self.cache('list', x)