From: Robert D. <rob...@gm...> - 2014-09-19 16:49:29
|
On 2014-09-19, Volker van Nek <vol...@gm...> wrote: > [ I also implemented exp(x), log(x), sin(x) and friends for bigfloat x. > In contrast to fpe1, .. these implementions are still some kind of raw > material and need more error analysis. And of course it's a lot more > difficult to merge them into Maxima. My idea is to put them into > /share/contrib for now. ] OK by me. > The values of %e, %pi, %gamma and log(2) computed by the described > functions are checked against reference values which were parsed and > rounded to the same precision. For all 0 < $fpprec < 10000 there are no > differences to the reference values. Sounds good. What about differences with respect to the existing implementation? Also, as mentioned by RJF, I would be interested to know if the new algorithms are faster or slower than the existing ones for relatively small values of fpprec. best Robert Dodier |