| 60 | // h(x) = x^3 * a^2 - 1, |
| 61 | #ifdef LIBC_TARGET_CPU_HAS_FMA_DOUBLE |
| 62 | double get_error(const DoubleDouble &x_3, const DoubleDouble &a_sq) { |
| 63 | return fputil::multiply_add(x_3.hi, a_sq.hi, -1.0) + |
| 64 | fputil::multiply_add(x_3.lo, a_sq.hi, x_3.hi * a_sq.lo); |
| 65 | } |
| 66 | #else |
| 67 | double get_error(const DoubleDouble &x_3, const DoubleDouble &a_sq) { |
| 68 | DoubleDouble x_3_a_sq = fputil::quick_mult(a_sq, x_3); |
no test coverage detected