[go: up one dir, main page]

Function accurate::two_product_fma [] [src]

pub fn two_product_fma<F>(a: F, b: F) -> (F, F) where F: Float

Product transformation

Transforms a product a * b into the pair (x, y) so that

x = fl(a * b)

is the product of a and b with floating point rounding applied and

y = a * b - x

is the remainder of the product.

References

Based on Ogita, Rump and Oishi 05