Dear Sage developers, the current version 2.5.2 of the flint library doesn't use enough floating point precision in _arith_bernoulli_number_zeta(). This causes bernoulli(250408) to be incorrect.
Since Sage uses flint as the default algorithm for n<=300000 the error affects Sage and can be observed, e.g., by comparing bernoulli(250408,'flint') and bernoulli(250408,'bernmm'). The bug is fixed in a development branch of flint, see this issue: https://github.com/wbhart/flint2/issues/288 Best, Alex -- 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 [email protected]. To post to this group, send email to [email protected]. Visit this group at https://groups.google.com/group/sage-devel. For more options, visit https://groups.google.com/d/optout.
