Exception for Jmat.Real.lambertw, newton loop step

• could not determine a ground truth

  1. wj = 1.2417946412699573e+183
  2. x = -2.7949007957275274e+253

Specification

Math

\[\begin{array}{l} \\ \begin{array}{l} t_0 := wj \cdot e^{wj}\\ wj - \frac{t\_0 - x}{e^{wj} + t\_0} \end{array} \end{array} \]

FPCore

(FPCore (wj x)
 :precision binary64
 (let* ((t_0 (* wj (exp wj)))) (- wj (/ (- t_0 x) (+ (exp wj) t_0)))))

C

double code(double wj, double x) {
	double t_0 = wj * exp(wj);
	return wj - ((t_0 - x) / (exp(wj) + t_0));
}

Fortran

real(8) function code(wj, x)
    real(8), intent (in) :: wj
    real(8), intent (in) :: x
    real(8) :: t_0
    t_0 = wj * exp(wj)
    code = wj - ((t_0 - x) / (exp(wj) + t_0))
end function

Java

public static double code(double wj, double x) {
	double t_0 = wj * Math.exp(wj);
	return wj - ((t_0 - x) / (Math.exp(wj) + t_0));
}

Python

def code(wj, x):
	t_0 = wj * math.exp(wj)
	return wj - ((t_0 - x) / (math.exp(wj) + t_0))

Julia

function code(wj, x)
	t_0 = Float64(wj * exp(wj))
	return Float64(wj - Float64(Float64(t_0 - x) / Float64(exp(wj) + t_0)))
end

MATLAB

function tmp = code(wj, x)
	t_0 = wj * exp(wj);
	tmp = wj - ((t_0 - x) / (exp(wj) + t_0));
end

Wolfram

code[wj_, x_] := Block[{t$95$0 = N[(wj * N[Exp[wj], $MachinePrecision]), $MachinePrecision]}, N[(wj - N[(N[(t$95$0 - x), $MachinePrecision] / N[(N[Exp[wj], $MachinePrecision] + t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]

Reproduce

herbie shell --seed 2024277 
(FPCore (wj x)
  :name "Jmat.Real.lambertw, newton loop step"
  :precision binary64

  :alt
  (! :herbie-platform default (let ((ew (exp wj))) (- wj (- (/ wj (+ wj 1)) (/ x (+ ew (* wj ew)))))))

  (- wj (/ (- (* wj (exp wj)) x) (+ (exp wj) (* wj (exp wj))))))

Please file a bug report with this information.

Backtrace

match: no matching clause for 'x		L	C
fn->ival-fn	/home/nightlies/.local/share/racket/8.10/pkgs/rival/eval/compile.rkt	17	0
rival-compile	/home/nightlies/.local/share/racket/8.10/pkgs/rival/eval/compile.rkt	253	0
(unnamed)	/usr/share/racket/collects/racket/contract/private/arrow-val-first.rkt	486	18
f115	.../match/compiler.rkt	559	40
compute-true-error	/home/nightlies/herbie/zane-concrete-values/src/core/localize.rkt	131	0
compute-local-errors	/home/nightlies/herbie/zane-concrete-values/src/core/localize.rkt	216	0
actual-errors	/home/nightlies/herbie/zane-concrete-values/src/core/explain.rkt	35	0
generate-timelines	/home/nightlies/herbie/zane-concrete-values/src/core/explain.rkt	498	0
explain!	/home/nightlies/herbie/zane-concrete-values/src/core/mainloop.rkt	323	0
run-improve!	/home/nightlies/herbie/zane-concrete-values/src/core/mainloop.rkt	45	0
get-alternatives/report	/home/nightlies/herbie/zane-concrete-values/src/api/sandbox.rkt	186	0
(unnamed)	/home/nightlies/herbie/zane-concrete-values/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