Currently set to No Index

Just evaluating a polynomial: how hard could it be?

The previous post looked at an example of strange floating point behavior taking from book End of Error. This post looks at another example.
This example, by Siegfried Rump, asks us to evaluate
333.75 y6 + x2 (11 x2y2 – y6 – 121 y4 – 2) + 5.5 y8 + x/(2y)
at x = 77617 and y = 33096.
Here we evaluate Rump’s example in single, double, and quadruple precision.

Read Full Story