Bouland and Aaronson, Equation (24)

Specification

?
\[\begin{array}{l} \\ \left({\left(a \cdot a + b \cdot b\right)}^{2} + 4 \cdot \left(\left(a \cdot a\right) \cdot \left(1 - a\right) + \left(b \cdot b\right) \cdot \left(3 + a\right)\right)\right) - 1 \end{array} \]
(FPCore (a b)
 :precision binary64
 (-
  (+
   (pow (+ (* a a) (* b b)) 2.0)
   (* 4.0 (+ (* (* a a) (- 1.0 a)) (* (* b b) (+ 3.0 a)))))
  1.0))
double code(double a, double b) {
	return (pow(((a * a) + (b * b)), 2.0) + (4.0 * (((a * a) * (1.0 - a)) + ((b * b) * (3.0 + a))))) - 1.0;
}

real(8) function code(a, b)
    real(8), intent (in) :: a
    real(8), intent (in) :: b
    code = ((((a * a) + (b * b)) ** 2.0d0) + (4.0d0 * (((a * a) * (1.0d0 - a)) + ((b * b) * (3.0d0 + a))))) - 1.0d0
end function

public static double code(double a, double b) {
	return (Math.pow(((a * a) + (b * b)), 2.0) + (4.0 * (((a * a) * (1.0 - a)) + ((b * b) * (3.0 + a))))) - 1.0;
}

def code(a, b):
	return (math.pow(((a * a) + (b * b)), 2.0) + (4.0 * (((a * a) * (1.0 - a)) + ((b * b) * (3.0 + a))))) - 1.0

function code(a, b)
	return Float64(Float64((Float64(Float64(a * a) + Float64(b * b)) ^ 2.0) + Float64(4.0 * Float64(Float64(Float64(a * a) * Float64(1.0 - a)) + Float64(Float64(b * b) * Float64(3.0 + a))))) - 1.0)
end

function tmp = code(a, b)
	tmp = ((((a * a) + (b * b)) ^ 2.0) + (4.0 * (((a * a) * (1.0 - a)) + ((b * b) * (3.0 + a))))) - 1.0;
end

code[a_, b_] := N[(N[(N[Power[N[(N[(a * a), $MachinePrecision] + N[(b * b), $MachinePrecision]), $MachinePrecision], 2.0], $MachinePrecision] + N[(4.0 * N[(N[(N[(a * a), $MachinePrecision] * N[(1.0 - a), $MachinePrecision]), $MachinePrecision] + N[(N[(b * b), $MachinePrecision] * N[(3.0 + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - 1.0), $MachinePrecision]

\begin{array}{l}

\\
\left({\left(a \cdot a + b \cdot b\right)}^{2} + 4 \cdot \left(\left(a \cdot a\right) \cdot \left(1 - a\right) + \left(b \cdot b\right) \cdot \left(3 + a\right)\right)\right) - 1
\end{array}

Reproduce

?
herbie shell --seed 2024196 
(FPCore (a b)
  :name "Bouland and Aaronson, Equation (24)"
  :precision binary64
  (- (+ (pow (+ (* a a) (* b b)) 2.0) (* 4.0 (+ (* (* a a) (- 1.0 a)) (* (* b b) (+ 3.0 a))))) 1.0))

Please file a bug report with this information.

Backtrace

vector-set!: contract violation expected: exact-nonnegative-integer? given: -1 argument position: 2nd first argument...: '#(#f #t #t #f #f #f #f #f #f #f #f #f #f #f #f #f #f #f #f #f #t #f #t #f #f #f #f #f #f #f #f #f #f #f #f #f #f #f #f #f #f #f #f #f)LC
backward-pass/home/nightlies/rival/main/eval/adjust.rkt80
rival-machine-adjust/home/nightlies/rival/main/eval/run.rkt980
rival-machine-full/home/nightlies/rival/main/eval/main.rkt160
rival-apply/home/nightlies/rival/main/eval/main.rkt560
(unnamed)/usr/local/racket/collects/racket/contract/private/arrow-val-first.rkt48618
(unnamed)/home/nightlies/herbie/zane-odyssey-timeline/src/sampling.rkt1594
ival-eval/home/nightlies/herbie/zane-odyssey-timeline/src/sampling.rkt1530
<eval-prog-real>/home/nightlies/herbie/zane-odyssey-timeline/src/sampling.rkt1882
compute-local-errors/home/nightlies/herbie/zane-odyssey-timeline/src/core/localize.rkt1050
batch-localize-errors/home/nightlies/herbie/zane-odyssey-timeline/src/core/localize.rkt910
localize!/home/nightlies/herbie/zane-odyssey-timeline/src/mainloop.rkt1870
run-iter!/home/nightlies/herbie/zane-odyssey-timeline/src/mainloop.rkt700
mutate!/home/nightlies/herbie/zane-odyssey-timeline/src/mainloop.rkt620
run-improve!/home/nightlies/herbie/zane-odyssey-timeline/src/mainloop.rkt480
get-alternatives/report/home/nightlies/herbie/zane-odyssey-timeline/src/sandbox.rkt1670
(unnamed)/home/nightlies/herbie/zane-odyssey-timeline/src/sandbox.rkt2616
(unnamed)/usr/local/racket/share/pkgs/profile-lib/main.rkt4010
profile-thunk/usr/local/racket/share/pkgs/profile-lib/main.rkt90
(unnamed)/usr/local/racket/collects/racket/engine.rkt4224