I'm new to boost. I am trying to use the float128 type for calculating sqrts as a larger part of code, however I am not sure the values obtained are accurate to float128 precision. I have tried running the example for printing log(2) found here, I find the same for sqrt:
However, the value printed log(2) only agrees to 16 decimal places, i.e. double precision. The value I obtain is
0.693147180559945286226763982995180413 rather than 0.693147180559945309417232121458176575
g++ (GCC) 4.8.2 20140120 (Red Hat 4.8.2-15) with boost 1.57.0.
Any advice as to what I am doing incorrectly?
My compile line is:
g++ -I ~/$BOOSTLIB new.cpp -lquadmath -o new
You don’t give your exact code but I suspect that you have fallen into what is an annoying (but unavoidable) pit that is eagerly awaiting the unwary.
(Been there, done that – more than once!)
You need to ensure that *all* values processed really are float128 or other multiprecision types by construction or casting, for example log(float128(2)).
(Otherwise computation will take place in double precision and only then be converted to float128 and printed, as you observe).