Average Error: 0.0 → 0.0
Time: 12.4s
Precision: 64
\[-\log \left(\frac{1}{x} - 1\right)\]
↓
\[-\log \left(\mathsf{expm1}\left(-\log x\right)\right)\]
-\log \left(\frac{1}{x} - 1\right)
↓
-\log \left(\mathsf{expm1}\left(-\log x\right)\right)
double f(double x) {
        double r333054 = 1.0;
        double r333055 = x;
        double r333056 = r333054 / r333055;
        double r333057 = r333056 - r333054;
        double r333058 = log(r333057);
        double r333059 = -r333058;
        return r333059;
}
↓
double f(double x) {
        double r333060 = x;
        double r333061 = log(r333060);
        double r333062 = -r333061;
        double r333063 = expm1(r333062);
        double r333064 = log(r333063);
        double r333065 = -r333064;
        return r333065;
}

Error

Bits error versus x

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Derivation

  1. Initial program 0.0

    \[-\log \left(\frac{1}{x} - 1\right)\]
  2. Using strategy rm

  3. Applied add-exp-log0.0

    \[\leadsto -\log \left(\frac{1}{\color{blue}{e^{\log x}}} - 1\right)\]

  4. Applied rec-exp0.0

    \[\leadsto -\log \left(\color{blue}{e^{-\log x}} - 1\right)\]

  5. Applied expm1-def0.0

    \[\leadsto -\log \color{blue}{\left(\mathsf{expm1}\left(-\log x\right)\right)}\]

  6. Final simplification0.0

    \[\leadsto -\log \left(\mathsf{expm1}\left(-\log x\right)\right)\]

Reproduce

herbie shell --seed 2019163 +o rules:numerics
(FPCore (x)
  :name "neg log"
  (- (log (- (/ 1 x) 1))))