Accurate Evaluation of Polynomials in Legendre Basis

Algorithm 4

Product of three floating point numbers.

function [𝑥,𝑦]=𝑇ℎ𝑟𝑒𝑒𝑃𝑟𝑜𝑑(𝑎,𝑏,𝑐)

[𝑡,ℎ]=𝑇𝑤𝑜𝑃𝑟𝑜𝑑(𝑎,𝑏)

[𝑥,𝑒]=𝑇𝑤𝑜𝑃𝑟𝑜𝑑(𝑡,𝑐)

𝑦=𝑐⊗ℎ⊕𝑒