UNPKG

double-double

Version:

Pure double-double precision functions *with strict error bounds*.

15 lines (14 loc) 584 B
/** * Returns the product of two double-double-precision floating point numbers. * * * relative error bound: 7u^2, i.e. fl(a+b) = (a+b)(1+ϵ), * where ϵ <= 7u^2, u = 0.5 * Number.EPSILON * the error bound is not sharp - the worst case that could be found by the * authors were 5u^2 * * * ALGORITHM 10 of https://hal.archives-ouvertes.fr/hal-01351529v3/document * @param x a double-double precision floating point number * @param y another double-double precision floating point number */ declare function ddMultDd(x: number[], y: number[]): number[]; export { ddMultDd };