Re: [Numpy-discussion] SFMT (faster mersenne twister)

2014-09-05 Thread Sturla Molden
Robert Kern wrote: > No, that is not what I meant. If the SFMT can be made to output the > same bitstream for the same seed, we can use it (modifying it if > necessary to avoid global state if necessary), but it does not look so > to me. I welcome corrections on that count (in PR form, preferably

Re: [Numpy-discussion] SFMT (faster mersenne twister)

2014-09-05 Thread Robert Kern
On Fri, Sep 5, 2014 at 9:36 PM, Sturla Molden wrote: > On 05/09/14 19:19, Neal Becker wrote: > >> It's a variant of the standard MT rather than just an implementation > >> of it, so we can't just drop it in. You will need to build the > >> infrastructure to support multiple PRNGs first (or rath

Re: [Numpy-discussion] SFMT (faster mersenne twister)

2014-09-05 Thread Sturla Molden
On 05/09/14 19:19, Neal Becker wrote: >> It's a variant of the standard MT rather than just an implementation >> of it, so we can't just drop it in. You will need to build the >> infrastructure to support multiple PRNGs first (or rather, build the >> infrastructure to reuse the non-uniform dist

Re: [Numpy-discussion] SFMT (faster mersenne twister)

2014-09-05 Thread Alan G Isaac
On 9/5/2014 1:19 PM, Neal Becker wrote: > I think it's somewhat debatable whether generating a different sequence of > random numbers counts as breaking backward compatibility. Please: it does. Alan Isaac ___ NumPy-Discussion mailing list NumPy-Discuss

[Numpy-discussion] numpy.vectorize docstrings not shown by help command

2014-09-05 Thread John Haiducek
When I apply numpy.vectorize() to a function, documentation tools behave inconsistently with regard to the new, vectorized function. The function's __doc__ attribute does contain the docstring of the original function as expected, but the built-in help() command displays the documentation of the

Re: [Numpy-discussion] SFMT (faster mersenne twister)

2014-09-05 Thread Phil Hodge
On 09/05/2014 01:19 PM, Neal Becker wrote: > You mean it's not backward compatible because it won't generate exactly the > same > sequence of output for a given seed, and therefore we wouldn't want to make > that > change? > > I think it's somewhat debatable whether generating a different sequenc

Re: [Numpy-discussion] SFMT (faster mersenne twister)

2014-09-05 Thread Robert Kern
On Fri, Sep 5, 2014 at 6:19 PM, Neal Becker wrote: > Robert Kern wrote: > >> On Fri, Sep 5, 2014 at 12:05 PM, Neal Becker wrote: >>> Robert Kern wrote: >>> On Thu, Sep 4, 2014 at 12:32 PM, Neal Becker wrote: > http://www.math.sci.hiroshima-u.ac.jp/~%20m-mat/MT/SFMT/index.html

Re: [Numpy-discussion] SFMT (faster mersenne twister)

2014-09-05 Thread Neal Becker
Robert Kern wrote: > On Fri, Sep 5, 2014 at 12:05 PM, Neal Becker wrote: >> Robert Kern wrote: >> >>> On Thu, Sep 4, 2014 at 12:32 PM, Neal Becker wrote: http://www.math.sci.hiroshima-u.ac.jp/~%20m-mat/MT/SFMT/index.html >>> >>> What would you like to say about it? >>> >> >> If it is faster

Re: [Numpy-discussion] SFMT (faster mersenne twister)

2014-09-05 Thread Robert Kern
On Fri, Sep 5, 2014 at 12:05 PM, Neal Becker wrote: > Robert Kern wrote: > >> On Thu, Sep 4, 2014 at 12:32 PM, Neal Becker wrote: >>> http://www.math.sci.hiroshima-u.ac.jp/~%20m-mat/MT/SFMT/index.html >> >> What would you like to say about it? >> > > If it is faster (and at least as good), maybe

Re: [Numpy-discussion] SFMT (faster mersenne twister)

2014-09-05 Thread Neal Becker
Robert Kern wrote: > On Thu, Sep 4, 2014 at 12:32 PM, Neal Becker wrote: >> http://www.math.sci.hiroshima-u.ac.jp/~%20m-mat/MT/SFMT/index.html > > What would you like to say about it? > If it is faster (and at least as good), maybe we'd like to adopt it to replace that used for mtrand -- --