[sage-devel] Re: Missing digits in numerical_approx

2011-12-15 Thread Sébastien Labbé
> Somebody needs to open a ticket. http://trac.sagemath.org/sage_trac/ticket/12163 -- To post to this group, send an email to sage-devel@googlegroups.com To unsubscribe from this group, send an email to sage-devel+unsubscr...@googlegroups.com For more options, visit this group at http://groups.

Re: [sage-devel] Re: Missing digits in numerical_approx

2011-12-14 Thread D. S. McNeil
I think problems caused by the log(10,2) ~ 3.32192 assumption occur in several places throughout Sage, and it's the source of at least one bug report (http://trac.sagemath.org/sage_trac/ticket/10164). The following all look suspicious to me: misc/functional.py:prec = int((digits+1) *

[sage-devel] Re: Missing digits in numerical_approx

2011-12-14 Thread Harald Schilly
Via "n??" I looked at the code. There is a line prec = int((digits+1) * 3.32192) + 1 where the 3.32… is an approximation to log_2(10). I highly suspect that's the source of it. Make the trailing 2 a 3 and it should give you more than the requested number of digits. sage: n(log(10,2), prec=100) 3