Exception for Kahan's exp quotient

• 24.4% of points produce a very large (infinite) output. You may want to add a precondition.

Specification

Math

\[\begin{array}{l} \\ \frac{e^{x} - 1}{x} \end{array} \]

FPCore

(FPCore (x) :precision binary64 (/ (- (exp x) 1.0) x))

C

double code(double x) {
	return (exp(x) - 1.0) / x;
}

Fortran

real(8) function code(x)
    real(8), intent (in) :: x
    code = (exp(x) - 1.0d0) / x
end function

Java

public static double code(double x) {
	return (Math.exp(x) - 1.0) / x;
}

Python

def code(x):
	return (math.exp(x) - 1.0) / x

Julia

function code(x)
	return Float64(Float64(exp(x) - 1.0) / x)
end

MATLAB

function tmp = code(x)
	tmp = (exp(x) - 1.0) / x;
end

Wolfram

code[x_] := N[(N[(N[Exp[x], $MachinePrecision] - 1.0), $MachinePrecision] / x), $MachinePrecision]

Reproduce

herbie shell --seed 2024276 
(FPCore (x)
  :name "Kahan's exp quotient"
  :precision binary64

  :alt
  (! :herbie-platform default (if (and (< x 1) (> x -1)) (/ (- (exp x) 1) (log (exp x))) (/ (- (exp x) 1) x)))

  (/ (- (exp x) 1.0) x))

Please file a bug report with this information.

Backtrace

match-define: no matching clause for (ddlf +nan.0 +nan.0 #t 3.5689737640740505e+33 187701457679585500.0)		L	C
logabs	/home/nightlies/herbie/bhargav-log-space/src/core/logspace.rkt	189	0
representable?	/home/nightlies/herbie/bhargav-log-space/src/core/logspace.rkt	65	0
f149	.../match/compiler.rkt	559	40
predict-errors	/home/nightlies/herbie/bhargav-log-space/src/core/explain.rkt	102	0
explain	/home/nightlies/herbie/bhargav-log-space/src/core/explain.rkt	712	0
explain!	/home/nightlies/herbie/bhargav-log-space/src/core/mainloop.rkt	323	0
run-improve!	/home/nightlies/herbie/bhargav-log-space/src/core/mainloop.rkt	45	0
get-alternatives/report	/home/nightlies/herbie/bhargav-log-space/src/api/sandbox.rkt	186	0
(unnamed)	/home/nightlies/herbie/bhargav-log-space/src/api/sandbox.rkt	282	6
(unnamed)	/usr/share/racket/pkgs/profile-lib/main.rkt	40	10
profile-thunk	/usr/share/racket/pkgs/profile-lib/main.rkt	9	0
(unnamed)	/usr/share/racket/collects/racket/engine.rkt	42	24