On Tue, 26 Apr 2016, Kiran Kedlaya wrote:
Yes, but the docstring itself still needs to be corrected: it says 'for
example "Matrix.random(GF(3), 2)" never generates "[[2,0],[0,2]]".' without
the flag.
OK, I added that to http://trac.sagemath.org/ticket/20465 .
--
Jori Mäntysalo
On Sunday, April 24, 2016 at 4:13:44 AM UTC-7, Jori Mäntysalo wrote:
>
> On Sun, 24 Apr 2016, Volker Braun wrote:
>
> >Warning: Matrices generated are not uniformly distributed. For
> > unimodular matrices over finite field this function does not even
> >
> > At least its not THAT
On Sun, 24 Apr 2016, Volker Braun wrote:
Warning: Matrices generated are not uniformly distributed. For
unimodular matrices over finite field this function does not even
At least its not THAT predictable ;-)
... after some attempts ...
sage: MatrixSpace(GF(3), 2, 2).random_element()
On Sunday, April 24, 2016 at 1:01:39 AM UTC+2, Kiran Kedlaya wrote:
>
> This seems to be filling out each row using successive random 32-bit words
> from gmp, lopping off the extra bits at the end. Is that a bad idea for
> some reason?
>
Its only lopping off the excess bits at the last column
On Sunday, April 24, 2016 at 12:36:28 AM UTC+2, Dima Pasechnik wrote:
>
> sage: random_matrix?
> ...
>Warning: Matrices generated are not uniformly distributed. For
> unimodular matrices over finite field this function does not even
> generate all of them: for example
On Sunday, April 24, 2016 at 12:04:27 AM UTC+1, Kiran Kedlaya wrote:
>
> On Saturday, April 23, 2016 at 3:36:28 PM UTC-7, Dima Pasechnik wrote:
>>
>>
>>
>> On Saturday, April 23, 2016 at 10:16:33 PM UTC+1, Kiran Kedlaya wrote:
>>>
>>> I just tried the following experiment on several different
On Saturday, April 23, 2016 at 3:36:28 PM UTC-7, Dima Pasechnik wrote:
>
>
>
> On Saturday, April 23, 2016 at 10:16:33 PM UTC+1, Kiran Kedlaya wrote:
>>
>> I just tried the following experiment on several different machines;
>> several different recent versions of Sage; algorithm = ple, m4ri; and
On Saturday, April 23, 2016 at 3:20:32 PM UTC-7, Volker Braun wrote:
>
> I'm guessing that its not a coincidence that Mersenne twister has a period
> of 2^19937 ;-)
>
> Probably not!
> The relevant code is
> in sage.matrix.matrix_mod2_dense.Matrix_mod2_dense.randomize
>
> Not sure exactly
On Saturday, April 23, 2016 at 10:16:33 PM UTC+1, Kiran Kedlaya wrote:
>
> I just tried the following experiment on several different machines;
> several different recent versions of Sage; algorithm = ple, m4ri; and n =
> 2, 25000, 3. In all cases, the result is identical:
>
> sage: M
I'm guessing that its not a coincidence that Mersenne twister has a period
of 2^19937 ;-)
The relevant code is
in sage.matrix.matrix_mod2_dense.Matrix_mod2_dense.randomize
Not sure exactly whats wrong but taking rstate.c_random() % (number of
columns) at the very least leads to modulo bias
10 matches
Mail list logo