On Mo, 2014-07-28 at 14:37 +0200, Eelco Hoogendoorn wrote: > To rephrase my most pressing question: may np.ones((N,2)).mean(0) and > np.ones((2,N)).mean(1) produce different results with the > implementation in the current master? If so, I think that would be > very much regrettable; and if this is a minority opinion, I do hope > that at least this gets documented in a most explicit manner. >
This will always give you different results. Though in master. the difference is more likely to be large, since (often the second one) maybe be less likely to run into bigger numerical issues. > > On Sun, Jul 27, 2014 at 8:26 PM, Sturla Molden > <sturla.mol...@gmail.com> wrote: > Nathaniel Smith <n...@pobox.com> wrote: > > > The problem here is that when summing up the values, the sum > gets > > large enough that after rounding, x + 1 = x and the sum > stops > > increasing. > > > Interesting. That explains why the divide-and-conquer > reduction is much > more robust. > > Thanks :) > > > Sturla > > _______________________________________________ > NumPy-Discussion mailing list > NumPy-Discussion@scipy.org > http://mail.scipy.org/mailman/listinfo/numpy-discussion > > > > _______________________________________________ > NumPy-Discussion mailing list > NumPy-Discussion@scipy.org > http://mail.scipy.org/mailman/listinfo/numpy-discussion _______________________________________________ NumPy-Discussion mailing list NumPy-Discussion@scipy.org http://mail.scipy.org/mailman/listinfo/numpy-discussion
