Exception for Toniolo and Linder, Equation (2)

Math

\[\sin^{-1} \left(\sqrt{\frac{1 - {\left(\frac{Om}{Omc}\right)}^{2}}{1 + 2 \cdot {\left(\frac{t}{\ell}\right)}^{2}}}\right) \]

FPCore

(FPCore (t l Om Omc)
 :precision binary64
 (asin
  (sqrt (/ (- 1.0 (pow (/ Om Omc) 2.0)) (+ 1.0 (* 2.0 (pow (/ t l) 2.0)))))))

C

double code(double t, double l, double Om, double Omc) {
	return asin(sqrt(((1.0 - pow((Om / Omc), 2.0)) / (1.0 + (2.0 * pow((t / l), 2.0))))));
}

Fortran

real(8) function code(t, l, om, omc)
    real(8), intent (in) :: t
    real(8), intent (in) :: l
    real(8), intent (in) :: om
    real(8), intent (in) :: omc
    code = asin(sqrt(((1.0d0 - ((om / omc) ** 2.0d0)) / (1.0d0 + (2.0d0 * ((t / l) ** 2.0d0))))))
end function

Java

public static double code(double t, double l, double Om, double Omc) {
	return Math.asin(Math.sqrt(((1.0 - Math.pow((Om / Omc), 2.0)) / (1.0 + (2.0 * Math.pow((t / l), 2.0))))));
}

Python

def code(t, l, Om, Omc):
	return math.asin(math.sqrt(((1.0 - math.pow((Om / Omc), 2.0)) / (1.0 + (2.0 * math.pow((t / l), 2.0))))))

Julia

function code(t, l, Om, Omc)
	return asin(sqrt(Float64(Float64(1.0 - (Float64(Om / Omc) ^ 2.0)) / Float64(1.0 + Float64(2.0 * (Float64(t / l) ^ 2.0))))))
end

MATLAB

function tmp = code(t, l, Om, Omc)
	tmp = asin(sqrt(((1.0 - ((Om / Omc) ^ 2.0)) / (1.0 + (2.0 * ((t / l) ^ 2.0))))));
end

Wolfram

code[t_, l_, Om_, Omc_] := N[ArcSin[N[Sqrt[N[(N[(1.0 - N[Power[N[(Om / Omc), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] / N[(1.0 + N[(2.0 * N[Power[N[(t / l), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]

Reproduce

Please include this information when filing a bug report:

herbie shell --seed 2022203 
(FPCore (t l Om Omc)
  :name "Toniolo and Linder, Equation (2)"
  :precision binary64
  (asin (sqrt (/ (- 1.0 (pow (/ Om Omc) 2.0)) (+ 1.0 (* 2.0 (pow (/ t l) 2.0)))))))

Backtrace

hash-ref: no value found for key key: '~		L	C
get-parametric-operator	/data/pavpan/nightlies/herbie/enumo-rat/src/syntax/syntax.rkt	218	0
f336	.../match/compiler.rkt	548	40
expand-parametric	/data/pavpan/nightlies/herbie/enumo-rat/src/syntax/sugar.rkt	62	0
(unnamed)	/data/pavpan/nightlies/herbie/enumo-rat/src/syntax/rules.rkt	176	8
generate-rules-for	/data/pavpan/nightlies/herbie/enumo-rat/src/syntax/rules.rkt	161	0
generate-missing-rules	/data/pavpan/nightlies/herbie/enumo-rat/src/syntax/rules.rkt	187	0
(unnamed)	/data/pavpan/nightlies/herbie/enumo-rat/src/syntax/rules.rkt	29	26
get-swaps	/data/pavpan/nightlies/herbie/enumo-rat/src/symmetry.rkt	6	0
connected-components	/data/pavpan/nightlies/herbie/enumo-rat/src/symmetry.rkt	17	0
run-improve!	/data/pavpan/nightlies/herbie/enumo-rat/src/mainloop.rkt	296	0
(unnamed)	/data/pavpan/nightlies/herbie/enumo-rat/src/sandbox.rkt	55	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