Fascinating but way over my head

For implementing it, I would need a file with a few thousand x/y pairs in REAL10 format. If somebody can create one with the boost library, I can supply the Assembly proxy for comparisons and benchmarks, using
FastMath. Below an example using a Log10 function.
FastMath FastLog10 ; define a math function
For_ fct=0.0 To 10.0 Step 0.5 ; max 10,000 iterations
fld fct ; X
fstp REAL10 ptr [edi]
void Log10(fct) ; Y (built-in MasmBasic function)
fstp REAL10 ptr [edi+REAL10]
add edi, 2*REAL10
Next
FastMath
Print Str$("Log(5.0)=%Jf", FastLog10(5)v)