Exception for Linear.Quaternion:$c/ from linear-1.19.1.3, C

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

Specification

Math

\[\begin{array}{l} \\ \left(\left(x \cdot y + y \cdot y\right) - y \cdot z\right) - y \cdot y \end{array} \]

FPCore

(FPCore (x y z)
 :precision binary64
 (- (- (+ (* x y) (* y y)) (* y z)) (* y y)))

C

double code(double x, double y, double z) {
	return (((x * y) + (y * y)) - (y * z)) - (y * y);
}

Fortran

real(8) function code(x, y, z)
    real(8), intent (in) :: x
    real(8), intent (in) :: y
    real(8), intent (in) :: z
    code = (((x * y) + (y * y)) - (y * z)) - (y * y)
end function

Java

public static double code(double x, double y, double z) {
	return (((x * y) + (y * y)) - (y * z)) - (y * y);
}

Python

def code(x, y, z):
	return (((x * y) + (y * y)) - (y * z)) - (y * y)

Julia

function code(x, y, z)
	return Float64(Float64(Float64(Float64(x * y) + Float64(y * y)) - Float64(y * z)) - Float64(y * y))
end

MATLAB

function tmp = code(x, y, z)
	tmp = (((x * y) + (y * y)) - (y * z)) - (y * y);
end

Wolfram

code[x_, y_, z_] := N[(N[(N[(N[(x * y), $MachinePrecision] + N[(y * y), $MachinePrecision]), $MachinePrecision] - N[(y * z), $MachinePrecision]), $MachinePrecision] - N[(y * y), $MachinePrecision]), $MachinePrecision]

Reproduce

herbie shell --seed 2024277 
(FPCore (x y z)
  :name "Linear.Quaternion:$c/ from linear-1.19.1.3, C"
  :precision binary64

  :alt
  (! :herbie-platform default (* (- x z) y))

  (- (- (+ (* x y) (* y y)) (* y z)) (* y y)))

Please file a bug report with this information.

Backtrace

rival-apply: contract violation expected: (or/c ^TYPE? boolean?) given: 0 in: an element of the 2nd argument of (-> rival-machine? (vectorof (or/c ^TYPE? boolean?)) (vectorof any/c)) contract from: <pkgs>/rival/main.rkt blaming: <pkgs>/herbie/core/rival.rkt (assuming the contract is correct) at: <pkgs>/rival/main.rkt:95:24		L	C
raise-blame-error	/usr/share/racket/collects/racket/contract/private/blame.rkt	346	0
loop	.../private/map.rkt	40	19
vector-map/update	/usr/share/racket/collects/racket/vector.rkt	62	0
vector-map	/usr/share/racket/collects/racket/vector.rkt	106	0
rival-apply	/home/nightlies/.local/share/racket/8.10/pkgs/rival/eval/main.rkt	65	0
(unnamed)	/usr/share/racket/collects/racket/contract/private/arrow-val-first.rkt	486	18
(unnamed)	/home/nightlies/herbie/zane-concrete-values/src/core/rival.rkt	76	4
real-apply	/home/nightlies/herbie/zane-concrete-values/src/core/rival.rkt	67	0
(unnamed)	.../private/arrow-val-first.rkt	486	18
<eval-prog-real>	/home/nightlies/herbie/zane-concrete-values/src/core/sampling.rkt	136	2
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