Jeffrey Sarnoff

You are getting contaminated results. There is a better way to benchmark it. ``` using DoubleFloats, BenchmarkTools dfc = ( df64"3.133239376997663645548490085151484674892E16", df64"-5.479944965767990821079467311839107722107E14", df64"6.290828903904724265980249871997551894090E12", df64"-3.633750176832769659849028554429106299915E10", df64"1.207743757532429576399485415069244807022E8", df64"-2.107485999925074577174305650549367415465E5" ); bfc = (...

You are repeatedly evaluating the coeffs. The string to Double64 takes longer than the string to BigFloat because the first uses the second. establish the coeffs outside of the function...

@heltonmc something like this ``` module BesselJ0 export besselj0 using DoubleFloats const J0_2N = ( df64"3.133239376997663645548490085151484674892E16", df64"-5.479944965767990821079467311839107722107E14", df64"6.290828903904724265980249871997551894090E12", df64"-3.633750176832769659849028554429106299915E10", df64"1.207743757532429576399485415069244807022E8", df64"-2.107485999925074577174305650549367415465E5", df64"1.562826808020631846245296572935547005859E2" ) const J0_2D = ( df64"2.005273201278504733151033654496928968261E18", df64"2.063038558793221244373123294054149790864E16", df64"1.053350447931127971406896594022010524994E14",df64"3.496556557558702583143527876385508882310E11",...

Yes .. the Double32 type was an early artifact of hoping to redouble Double64s to compute Double128 accurately. At the time, Julia did not provide the mechanism for that sort...

There are at least a few approaches to implementing `exp` and `log`. One gives each their own algorithmic specialization and whatever specific computational advantages are to be had. One works...

If you sample deeply through all of the arithmetic and each elementary function, and then do the same with other implementations of double-double, we get at least as many and...

I need to find the version of that number I had decided upon some time ago. I do recall there being some threshold of cooperation throughout the functional realizations. OTOH,...

It looks good, much more useful than was available an hour ago! I need to write a harness of more focused test regions.

super good, smiling bits On Thu, Jan 20, 2022 at 7:27 PM Oscar Smith ***@***.***> wrote: > I just updated the code. I'm seeing about half the error of the...

let me know when you want to merge something