
(FPCore (a b c) :precision binary64 (/ (+ (- b) (sqrt (- (* b b) (* (* 4.0 a) c)))) (* 2.0 a)))
double code(double a, double b, double c) {
return (-b + sqrt(((b * b) - ((4.0 * a) * c)))) / (2.0 * a);
}
real(8) function code(a, b, c)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
code = (-b + sqrt(((b * b) - ((4.0d0 * a) * c)))) / (2.0d0 * a)
end function
public static double code(double a, double b, double c) {
return (-b + Math.sqrt(((b * b) - ((4.0 * a) * c)))) / (2.0 * a);
}
def code(a, b, c): return (-b + math.sqrt(((b * b) - ((4.0 * a) * c)))) / (2.0 * a)
function code(a, b, c) return Float64(Float64(Float64(-b) + sqrt(Float64(Float64(b * b) - Float64(Float64(4.0 * a) * c)))) / Float64(2.0 * a)) end
function tmp = code(a, b, c) tmp = (-b + sqrt(((b * b) - ((4.0 * a) * c)))) / (2.0 * a); end
code[a_, b_, c_] := N[(N[((-b) + N[Sqrt[N[(N[(b * b), $MachinePrecision] - N[(N[(4.0 * a), $MachinePrecision] * c), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(2.0 * a), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\left(-b\right) + \sqrt{b \cdot b - \left(4 \cdot a\right) \cdot c}}{2 \cdot a}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 11 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (a b c) :precision binary64 (/ (+ (- b) (sqrt (- (* b b) (* (* 4.0 a) c)))) (* 2.0 a)))
double code(double a, double b, double c) {
return (-b + sqrt(((b * b) - ((4.0 * a) * c)))) / (2.0 * a);
}
real(8) function code(a, b, c)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
code = (-b + sqrt(((b * b) - ((4.0d0 * a) * c)))) / (2.0d0 * a)
end function
public static double code(double a, double b, double c) {
return (-b + Math.sqrt(((b * b) - ((4.0 * a) * c)))) / (2.0 * a);
}
def code(a, b, c): return (-b + math.sqrt(((b * b) - ((4.0 * a) * c)))) / (2.0 * a)
function code(a, b, c) return Float64(Float64(Float64(-b) + sqrt(Float64(Float64(b * b) - Float64(Float64(4.0 * a) * c)))) / Float64(2.0 * a)) end
function tmp = code(a, b, c) tmp = (-b + sqrt(((b * b) - ((4.0 * a) * c)))) / (2.0 * a); end
code[a_, b_, c_] := N[(N[((-b) + N[Sqrt[N[(N[(b * b), $MachinePrecision] - N[(N[(4.0 * a), $MachinePrecision] * c), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(2.0 * a), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\left(-b\right) + \sqrt{b \cdot b - \left(4 \cdot a\right) \cdot c}}{2 \cdot a}
\end{array}
(FPCore (a b c) :precision binary64 (/ (* -2.0 c) (+ b (sqrt (+ (* a (* c -4.0)) (* b b))))))
double code(double a, double b, double c) {
return (-2.0 * c) / (b + sqrt(((a * (c * -4.0)) + (b * b))));
}
real(8) function code(a, b, c)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
code = ((-2.0d0) * c) / (b + sqrt(((a * (c * (-4.0d0))) + (b * b))))
end function
public static double code(double a, double b, double c) {
return (-2.0 * c) / (b + Math.sqrt(((a * (c * -4.0)) + (b * b))));
}
def code(a, b, c): return (-2.0 * c) / (b + math.sqrt(((a * (c * -4.0)) + (b * b))))
function code(a, b, c) return Float64(Float64(-2.0 * c) / Float64(b + sqrt(Float64(Float64(a * Float64(c * -4.0)) + Float64(b * b))))) end
function tmp = code(a, b, c) tmp = (-2.0 * c) / (b + sqrt(((a * (c * -4.0)) + (b * b)))); end
code[a_, b_, c_] := N[(N[(-2.0 * c), $MachinePrecision] / N[(b + N[Sqrt[N[(N[(a * N[(c * -4.0), $MachinePrecision]), $MachinePrecision] + N[(b * b), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{-2 \cdot c}{b + \sqrt{a \cdot \left(c \cdot -4\right) + b \cdot b}}
\end{array}
Initial program 29.8%
+-commutative29.8%
flip-+29.8%
add-sqr-sqrt30.8%
associate-*l*30.8%
cancel-sign-sub-inv30.8%
metadata-eval30.8%
sqr-neg30.8%
add-sqr-sqrt0.0%
sqrt-unprod0.8%
sqr-neg0.8%
sqrt-prod1.5%
add-sqr-sqrt0.8%
unsub-neg0.8%
+-commutative0.8%
Applied egg-rr30.8%
+-commutative30.8%
associate--l+99.4%
*-commutative99.4%
+-inverses99.4%
mul0-rgt99.4%
distribute-lft-out99.4%
metadata-eval99.4%
+-commutative99.4%
*-commutative99.4%
Simplified99.4%
div-inv99.2%
associate-/l*99.2%
*-commutative99.2%
Applied egg-rr99.2%
*-commutative99.2%
associate-*r/99.2%
times-frac99.4%
*-lft-identity99.4%
associate-/r*99.7%
*-commutative99.7%
*-commutative99.7%
times-frac99.7%
metadata-eval99.7%
associate-*l*99.7%
Simplified99.7%
Taylor expanded in a around 0 99.8%
unpow1/299.8%
Applied egg-rr99.8%
(FPCore (a b c)
:precision binary64
(let* ((t_0 (* a (* c c))) (t_1 (* b (* b b))))
(/
(* -2.0 c)
(*
b
(+
2.0
(+
(* -2.0 (+ (* a (/ c (* b b))) (/ (* a t_0) (* (* b b) (* b b)))))
(* -4.0 (/ (* a (* t_0 (* c a))) (* t_1 t_1)))))))))
double code(double a, double b, double c) {
double t_0 = a * (c * c);
double t_1 = b * (b * b);
return (-2.0 * c) / (b * (2.0 + ((-2.0 * ((a * (c / (b * b))) + ((a * t_0) / ((b * b) * (b * b))))) + (-4.0 * ((a * (t_0 * (c * a))) / (t_1 * t_1))))));
}
real(8) function code(a, b, c)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8) :: t_0
real(8) :: t_1
t_0 = a * (c * c)
t_1 = b * (b * b)
code = ((-2.0d0) * c) / (b * (2.0d0 + (((-2.0d0) * ((a * (c / (b * b))) + ((a * t_0) / ((b * b) * (b * b))))) + ((-4.0d0) * ((a * (t_0 * (c * a))) / (t_1 * t_1))))))
end function
public static double code(double a, double b, double c) {
double t_0 = a * (c * c);
double t_1 = b * (b * b);
return (-2.0 * c) / (b * (2.0 + ((-2.0 * ((a * (c / (b * b))) + ((a * t_0) / ((b * b) * (b * b))))) + (-4.0 * ((a * (t_0 * (c * a))) / (t_1 * t_1))))));
}
def code(a, b, c): t_0 = a * (c * c) t_1 = b * (b * b) return (-2.0 * c) / (b * (2.0 + ((-2.0 * ((a * (c / (b * b))) + ((a * t_0) / ((b * b) * (b * b))))) + (-4.0 * ((a * (t_0 * (c * a))) / (t_1 * t_1))))))
function code(a, b, c) t_0 = Float64(a * Float64(c * c)) t_1 = Float64(b * Float64(b * b)) return Float64(Float64(-2.0 * c) / Float64(b * Float64(2.0 + Float64(Float64(-2.0 * Float64(Float64(a * Float64(c / Float64(b * b))) + Float64(Float64(a * t_0) / Float64(Float64(b * b) * Float64(b * b))))) + Float64(-4.0 * Float64(Float64(a * Float64(t_0 * Float64(c * a))) / Float64(t_1 * t_1))))))) end
function tmp = code(a, b, c) t_0 = a * (c * c); t_1 = b * (b * b); tmp = (-2.0 * c) / (b * (2.0 + ((-2.0 * ((a * (c / (b * b))) + ((a * t_0) / ((b * b) * (b * b))))) + (-4.0 * ((a * (t_0 * (c * a))) / (t_1 * t_1)))))); end
code[a_, b_, c_] := Block[{t$95$0 = N[(a * N[(c * c), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(b * N[(b * b), $MachinePrecision]), $MachinePrecision]}, N[(N[(-2.0 * c), $MachinePrecision] / N[(b * N[(2.0 + N[(N[(-2.0 * N[(N[(a * N[(c / N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(a * t$95$0), $MachinePrecision] / N[(N[(b * b), $MachinePrecision] * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(-4.0 * N[(N[(a * N[(t$95$0 * N[(c * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t$95$1 * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := a \cdot \left(c \cdot c\right)\\
t_1 := b \cdot \left(b \cdot b\right)\\
\frac{-2 \cdot c}{b \cdot \left(2 + \left(-2 \cdot \left(a \cdot \frac{c}{b \cdot b} + \frac{a \cdot t\_0}{\left(b \cdot b\right) \cdot \left(b \cdot b\right)}\right) + -4 \cdot \frac{a \cdot \left(t\_0 \cdot \left(c \cdot a\right)\right)}{t\_1 \cdot t\_1}\right)\right)}
\end{array}
\end{array}
Initial program 29.8%
+-commutative29.8%
flip-+29.8%
add-sqr-sqrt30.8%
associate-*l*30.8%
cancel-sign-sub-inv30.8%
metadata-eval30.8%
sqr-neg30.8%
add-sqr-sqrt0.0%
sqrt-unprod0.8%
sqr-neg0.8%
sqrt-prod1.5%
add-sqr-sqrt0.8%
unsub-neg0.8%
+-commutative0.8%
Applied egg-rr30.8%
+-commutative30.8%
associate--l+99.4%
*-commutative99.4%
+-inverses99.4%
mul0-rgt99.4%
distribute-lft-out99.4%
metadata-eval99.4%
+-commutative99.4%
*-commutative99.4%
Simplified99.4%
div-inv99.2%
associate-/l*99.2%
*-commutative99.2%
Applied egg-rr99.2%
*-commutative99.2%
associate-*r/99.2%
times-frac99.4%
*-lft-identity99.4%
associate-/r*99.7%
*-commutative99.7%
*-commutative99.7%
times-frac99.7%
metadata-eval99.7%
associate-*l*99.7%
Simplified99.7%
Taylor expanded in a around 0 99.8%
Taylor expanded in b around inf 95.3%
Simplified95.3%
Final simplification95.3%
(FPCore (a b c) :precision binary64 (/ (* -2.0 c) (+ (* b 2.0) (* c (* -2.0 (+ (/ a b) (* (* a a) (/ c (* b (* b b))))))))))
double code(double a, double b, double c) {
return (-2.0 * c) / ((b * 2.0) + (c * (-2.0 * ((a / b) + ((a * a) * (c / (b * (b * b))))))));
}
real(8) function code(a, b, c)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
code = ((-2.0d0) * c) / ((b * 2.0d0) + (c * ((-2.0d0) * ((a / b) + ((a * a) * (c / (b * (b * b))))))))
end function
public static double code(double a, double b, double c) {
return (-2.0 * c) / ((b * 2.0) + (c * (-2.0 * ((a / b) + ((a * a) * (c / (b * (b * b))))))));
}
def code(a, b, c): return (-2.0 * c) / ((b * 2.0) + (c * (-2.0 * ((a / b) + ((a * a) * (c / (b * (b * b))))))))
function code(a, b, c) return Float64(Float64(-2.0 * c) / Float64(Float64(b * 2.0) + Float64(c * Float64(-2.0 * Float64(Float64(a / b) + Float64(Float64(a * a) * Float64(c / Float64(b * Float64(b * b))))))))) end
function tmp = code(a, b, c) tmp = (-2.0 * c) / ((b * 2.0) + (c * (-2.0 * ((a / b) + ((a * a) * (c / (b * (b * b)))))))); end
code[a_, b_, c_] := N[(N[(-2.0 * c), $MachinePrecision] / N[(N[(b * 2.0), $MachinePrecision] + N[(c * N[(-2.0 * N[(N[(a / b), $MachinePrecision] + N[(N[(a * a), $MachinePrecision] * N[(c / N[(b * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{-2 \cdot c}{b \cdot 2 + c \cdot \left(-2 \cdot \left(\frac{a}{b} + \left(a \cdot a\right) \cdot \frac{c}{b \cdot \left(b \cdot b\right)}\right)\right)}
\end{array}
Initial program 29.8%
+-commutative29.8%
flip-+29.8%
add-sqr-sqrt30.8%
associate-*l*30.8%
cancel-sign-sub-inv30.8%
metadata-eval30.8%
sqr-neg30.8%
add-sqr-sqrt0.0%
sqrt-unprod0.8%
sqr-neg0.8%
sqrt-prod1.5%
add-sqr-sqrt0.8%
unsub-neg0.8%
+-commutative0.8%
Applied egg-rr30.8%
+-commutative30.8%
associate--l+99.4%
*-commutative99.4%
+-inverses99.4%
mul0-rgt99.4%
distribute-lft-out99.4%
metadata-eval99.4%
+-commutative99.4%
*-commutative99.4%
Simplified99.4%
div-inv99.2%
associate-/l*99.2%
*-commutative99.2%
Applied egg-rr99.2%
*-commutative99.2%
associate-*r/99.2%
times-frac99.4%
*-lft-identity99.4%
associate-/r*99.7%
*-commutative99.7%
*-commutative99.7%
times-frac99.7%
metadata-eval99.7%
associate-*l*99.7%
Simplified99.7%
Taylor expanded in a around 0 99.8%
Taylor expanded in c around 0 94.1%
*-commutative94.1%
distribute-lft-out94.1%
associate-/l*94.1%
unpow294.1%
cube-unmult94.1%
Simplified94.1%
(FPCore (a b c) :precision binary64 (/ (* -2.0 c) (+ b (+ b (* c (* -2.0 (+ (/ a b) (* (/ c (* b b)) (/ (* a a) b)))))))))
double code(double a, double b, double c) {
return (-2.0 * c) / (b + (b + (c * (-2.0 * ((a / b) + ((c / (b * b)) * ((a * a) / b)))))));
}
real(8) function code(a, b, c)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
code = ((-2.0d0) * c) / (b + (b + (c * ((-2.0d0) * ((a / b) + ((c / (b * b)) * ((a * a) / b)))))))
end function
public static double code(double a, double b, double c) {
return (-2.0 * c) / (b + (b + (c * (-2.0 * ((a / b) + ((c / (b * b)) * ((a * a) / b)))))));
}
def code(a, b, c): return (-2.0 * c) / (b + (b + (c * (-2.0 * ((a / b) + ((c / (b * b)) * ((a * a) / b)))))))
function code(a, b, c) return Float64(Float64(-2.0 * c) / Float64(b + Float64(b + Float64(c * Float64(-2.0 * Float64(Float64(a / b) + Float64(Float64(c / Float64(b * b)) * Float64(Float64(a * a) / b)))))))) end
function tmp = code(a, b, c) tmp = (-2.0 * c) / (b + (b + (c * (-2.0 * ((a / b) + ((c / (b * b)) * ((a * a) / b))))))); end
code[a_, b_, c_] := N[(N[(-2.0 * c), $MachinePrecision] / N[(b + N[(b + N[(c * N[(-2.0 * N[(N[(a / b), $MachinePrecision] + N[(N[(c / N[(b * b), $MachinePrecision]), $MachinePrecision] * N[(N[(a * a), $MachinePrecision] / b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{-2 \cdot c}{b + \left(b + c \cdot \left(-2 \cdot \left(\frac{a}{b} + \frac{c}{b \cdot b} \cdot \frac{a \cdot a}{b}\right)\right)\right)}
\end{array}
Initial program 29.8%
+-commutative29.8%
flip-+29.8%
add-sqr-sqrt30.8%
associate-*l*30.8%
cancel-sign-sub-inv30.8%
metadata-eval30.8%
sqr-neg30.8%
add-sqr-sqrt0.0%
sqrt-unprod0.8%
sqr-neg0.8%
sqrt-prod1.5%
add-sqr-sqrt0.8%
unsub-neg0.8%
+-commutative0.8%
Applied egg-rr30.8%
+-commutative30.8%
associate--l+99.4%
*-commutative99.4%
+-inverses99.4%
mul0-rgt99.4%
distribute-lft-out99.4%
metadata-eval99.4%
+-commutative99.4%
*-commutative99.4%
Simplified99.4%
div-inv99.2%
associate-/l*99.2%
*-commutative99.2%
Applied egg-rr99.2%
*-commutative99.2%
associate-*r/99.2%
times-frac99.4%
*-lft-identity99.4%
associate-/r*99.7%
*-commutative99.7%
*-commutative99.7%
times-frac99.7%
metadata-eval99.7%
associate-*l*99.7%
Simplified99.7%
Taylor expanded in a around 0 99.8%
Taylor expanded in c around 0 94.0%
distribute-lft-out94.0%
*-commutative94.0%
unpow394.0%
unpow294.0%
times-frac94.0%
unpow294.0%
unpow294.0%
Simplified94.0%
(FPCore (a b c) :precision binary64 (/ (* -2.0 c) (+ (* b 2.0) (* c (* -2.0 (/ a b))))))
double code(double a, double b, double c) {
return (-2.0 * c) / ((b * 2.0) + (c * (-2.0 * (a / b))));
}
real(8) function code(a, b, c)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
code = ((-2.0d0) * c) / ((b * 2.0d0) + (c * ((-2.0d0) * (a / b))))
end function
public static double code(double a, double b, double c) {
return (-2.0 * c) / ((b * 2.0) + (c * (-2.0 * (a / b))));
}
def code(a, b, c): return (-2.0 * c) / ((b * 2.0) + (c * (-2.0 * (a / b))))
function code(a, b, c) return Float64(Float64(-2.0 * c) / Float64(Float64(b * 2.0) + Float64(c * Float64(-2.0 * Float64(a / b))))) end
function tmp = code(a, b, c) tmp = (-2.0 * c) / ((b * 2.0) + (c * (-2.0 * (a / b)))); end
code[a_, b_, c_] := N[(N[(-2.0 * c), $MachinePrecision] / N[(N[(b * 2.0), $MachinePrecision] + N[(c * N[(-2.0 * N[(a / b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{-2 \cdot c}{b \cdot 2 + c \cdot \left(-2 \cdot \frac{a}{b}\right)}
\end{array}
Initial program 29.8%
+-commutative29.8%
flip-+29.8%
add-sqr-sqrt30.8%
associate-*l*30.8%
cancel-sign-sub-inv30.8%
metadata-eval30.8%
sqr-neg30.8%
add-sqr-sqrt0.0%
sqrt-unprod0.8%
sqr-neg0.8%
sqrt-prod1.5%
add-sqr-sqrt0.8%
unsub-neg0.8%
+-commutative0.8%
Applied egg-rr30.8%
+-commutative30.8%
associate--l+99.4%
*-commutative99.4%
+-inverses99.4%
mul0-rgt99.4%
distribute-lft-out99.4%
metadata-eval99.4%
+-commutative99.4%
*-commutative99.4%
Simplified99.4%
div-inv99.2%
associate-/l*99.2%
*-commutative99.2%
Applied egg-rr99.2%
*-commutative99.2%
associate-*r/99.2%
times-frac99.4%
*-lft-identity99.4%
associate-/r*99.7%
*-commutative99.7%
*-commutative99.7%
times-frac99.7%
metadata-eval99.7%
associate-*l*99.7%
Simplified99.7%
Taylor expanded in a around 0 99.8%
Taylor expanded in a around 0 91.3%
+-commutative91.3%
*-commutative91.3%
associate-*r/91.3%
*-commutative91.3%
*-commutative91.3%
associate-*r*91.3%
associate-/l*91.3%
*-commutative91.3%
associate-*r/91.3%
Simplified91.3%
(FPCore (a b c) :precision binary64 (/ (* -2.0 c) (+ b (+ b (* a (* -2.0 (/ c b)))))))
double code(double a, double b, double c) {
return (-2.0 * c) / (b + (b + (a * (-2.0 * (c / b)))));
}
real(8) function code(a, b, c)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
code = ((-2.0d0) * c) / (b + (b + (a * ((-2.0d0) * (c / b)))))
end function
public static double code(double a, double b, double c) {
return (-2.0 * c) / (b + (b + (a * (-2.0 * (c / b)))));
}
def code(a, b, c): return (-2.0 * c) / (b + (b + (a * (-2.0 * (c / b)))))
function code(a, b, c) return Float64(Float64(-2.0 * c) / Float64(b + Float64(b + Float64(a * Float64(-2.0 * Float64(c / b)))))) end
function tmp = code(a, b, c) tmp = (-2.0 * c) / (b + (b + (a * (-2.0 * (c / b))))); end
code[a_, b_, c_] := N[(N[(-2.0 * c), $MachinePrecision] / N[(b + N[(b + N[(a * N[(-2.0 * N[(c / b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{-2 \cdot c}{b + \left(b + a \cdot \left(-2 \cdot \frac{c}{b}\right)\right)}
\end{array}
Initial program 29.8%
+-commutative29.8%
flip-+29.8%
add-sqr-sqrt30.8%
associate-*l*30.8%
cancel-sign-sub-inv30.8%
metadata-eval30.8%
sqr-neg30.8%
add-sqr-sqrt0.0%
sqrt-unprod0.8%
sqr-neg0.8%
sqrt-prod1.5%
add-sqr-sqrt0.8%
unsub-neg0.8%
+-commutative0.8%
Applied egg-rr30.8%
+-commutative30.8%
associate--l+99.4%
*-commutative99.4%
+-inverses99.4%
mul0-rgt99.4%
distribute-lft-out99.4%
metadata-eval99.4%
+-commutative99.4%
*-commutative99.4%
Simplified99.4%
div-inv99.2%
associate-/l*99.2%
*-commutative99.2%
Applied egg-rr99.2%
*-commutative99.2%
associate-*r/99.2%
times-frac99.4%
*-lft-identity99.4%
associate-/r*99.7%
*-commutative99.7%
*-commutative99.7%
times-frac99.7%
metadata-eval99.7%
associate-*l*99.7%
Simplified99.7%
Taylor expanded in a around 0 99.8%
Taylor expanded in a around 0 91.3%
associate-/l*91.3%
associate-*l*91.3%
*-commutative91.3%
associate-*r*91.3%
*-commutative91.3%
Simplified91.3%
Final simplification91.3%
(FPCore (a b c) :precision binary64 (* (/ -1.0 b) (+ c (* a (/ (* c c) (* b b))))))
double code(double a, double b, double c) {
return (-1.0 / b) * (c + (a * ((c * c) / (b * b))));
}
real(8) function code(a, b, c)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
code = ((-1.0d0) / b) * (c + (a * ((c * c) / (b * b))))
end function
public static double code(double a, double b, double c) {
return (-1.0 / b) * (c + (a * ((c * c) / (b * b))));
}
def code(a, b, c): return (-1.0 / b) * (c + (a * ((c * c) / (b * b))))
function code(a, b, c) return Float64(Float64(-1.0 / b) * Float64(c + Float64(a * Float64(Float64(c * c) / Float64(b * b))))) end
function tmp = code(a, b, c) tmp = (-1.0 / b) * (c + (a * ((c * c) / (b * b)))); end
code[a_, b_, c_] := N[(N[(-1.0 / b), $MachinePrecision] * N[(c + N[(a * N[(N[(c * c), $MachinePrecision] / N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{-1}{b} \cdot \left(c + a \cdot \frac{c \cdot c}{b \cdot b}\right)
\end{array}
Initial program 29.8%
Taylor expanded in a around 0 91.0%
+-commutative91.0%
mul-1-neg91.0%
unsub-neg91.0%
mul-1-neg91.0%
associate-/l*91.0%
distribute-rgt-neg-in91.0%
neg-sub091.0%
unpow291.0%
cube-unmult91.0%
times-frac91.0%
Simplified91.0%
Taylor expanded in b around 0 90.5%
distribute-lft-out90.5%
cube-unmult90.7%
times-frac90.6%
unpow290.6%
unpow290.6%
*-commutative90.6%
Simplified90.6%
Taylor expanded in a around 0 90.7%
unpow290.7%
associate-/l*90.7%
unpow290.7%
Simplified90.7%
(FPCore (a b c) :precision binary64 (* (/ -1.0 b) (* c (+ (* a (/ c (* b b))) 1.0))))
double code(double a, double b, double c) {
return (-1.0 / b) * (c * ((a * (c / (b * b))) + 1.0));
}
real(8) function code(a, b, c)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
code = ((-1.0d0) / b) * (c * ((a * (c / (b * b))) + 1.0d0))
end function
public static double code(double a, double b, double c) {
return (-1.0 / b) * (c * ((a * (c / (b * b))) + 1.0));
}
def code(a, b, c): return (-1.0 / b) * (c * ((a * (c / (b * b))) + 1.0))
function code(a, b, c) return Float64(Float64(-1.0 / b) * Float64(c * Float64(Float64(a * Float64(c / Float64(b * b))) + 1.0))) end
function tmp = code(a, b, c) tmp = (-1.0 / b) * (c * ((a * (c / (b * b))) + 1.0)); end
code[a_, b_, c_] := N[(N[(-1.0 / b), $MachinePrecision] * N[(c * N[(N[(a * N[(c / N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{-1}{b} \cdot \left(c \cdot \left(a \cdot \frac{c}{b \cdot b} + 1\right)\right)
\end{array}
Initial program 29.8%
Taylor expanded in a around 0 91.0%
+-commutative91.0%
mul-1-neg91.0%
unsub-neg91.0%
mul-1-neg91.0%
associate-/l*91.0%
distribute-rgt-neg-in91.0%
neg-sub091.0%
unpow291.0%
cube-unmult91.0%
times-frac91.0%
Simplified91.0%
Taylor expanded in b around 0 90.5%
distribute-lft-out90.5%
cube-unmult90.7%
times-frac90.6%
unpow290.6%
unpow290.6%
*-commutative90.6%
Simplified90.6%
Taylor expanded in c around 0 90.7%
unpow290.7%
associate-/l*90.7%
Simplified90.7%
Final simplification90.7%
(FPCore (a b c) :precision binary64 (* c (/ -1.0 b)))
double code(double a, double b, double c) {
return c * (-1.0 / b);
}
real(8) function code(a, b, c)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
code = c * ((-1.0d0) / b)
end function
public static double code(double a, double b, double c) {
return c * (-1.0 / b);
}
def code(a, b, c): return c * (-1.0 / b)
function code(a, b, c) return Float64(c * Float64(-1.0 / b)) end
function tmp = code(a, b, c) tmp = c * (-1.0 / b); end
code[a_, b_, c_] := N[(c * N[(-1.0 / b), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
c \cdot \frac{-1}{b}
\end{array}
Initial program 29.8%
Taylor expanded in a around 0 91.0%
+-commutative91.0%
mul-1-neg91.0%
unsub-neg91.0%
mul-1-neg91.0%
associate-/l*91.0%
distribute-rgt-neg-in91.0%
neg-sub091.0%
unpow291.0%
cube-unmult91.0%
times-frac91.0%
Simplified91.0%
Taylor expanded in b around 0 90.5%
distribute-lft-out90.5%
cube-unmult90.7%
times-frac90.6%
unpow290.6%
unpow290.6%
*-commutative90.6%
Simplified90.6%
Taylor expanded in a around 0 82.1%
Final simplification82.1%
(FPCore (a b c) :precision binary64 (/ c (- b)))
double code(double a, double b, double c) {
return c / -b;
}
real(8) function code(a, b, c)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
code = c / -b
end function
public static double code(double a, double b, double c) {
return c / -b;
}
def code(a, b, c): return c / -b
function code(a, b, c) return Float64(c / Float64(-b)) end
function tmp = code(a, b, c) tmp = c / -b; end
code[a_, b_, c_] := N[(c / (-b)), $MachinePrecision]
\begin{array}{l}
\\
\frac{c}{-b}
\end{array}
Initial program 29.8%
Taylor expanded in b around inf 82.3%
mul-1-neg82.3%
neg-sub082.3%
Simplified82.3%
Final simplification82.3%
(FPCore (a b c) :precision binary64 (/ b a))
double code(double a, double b, double c) {
return b / a;
}
real(8) function code(a, b, c)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
code = b / a
end function
public static double code(double a, double b, double c) {
return b / a;
}
def code(a, b, c): return b / a
function code(a, b, c) return Float64(b / a) end
function tmp = code(a, b, c) tmp = b / a; end
code[a_, b_, c_] := N[(b / a), $MachinePrecision]
\begin{array}{l}
\\
\frac{b}{a}
\end{array}
Initial program 29.8%
+-commutative29.8%
flip-+29.8%
add-sqr-sqrt30.8%
associate-*l*30.8%
cancel-sign-sub-inv30.8%
metadata-eval30.8%
sqr-neg30.8%
add-sqr-sqrt0.0%
sqrt-unprod0.8%
sqr-neg0.8%
sqrt-prod1.5%
add-sqr-sqrt0.8%
unsub-neg0.8%
+-commutative0.8%
Applied egg-rr30.8%
+-commutative30.8%
associate--l+99.4%
*-commutative99.4%
+-inverses99.4%
mul0-rgt99.4%
distribute-lft-out99.4%
metadata-eval99.4%
+-commutative99.4%
*-commutative99.4%
Simplified99.4%
div-inv99.2%
associate-/l*99.2%
*-commutative99.2%
Applied egg-rr99.2%
*-commutative99.2%
associate-*r/99.2%
times-frac99.4%
*-lft-identity99.4%
associate-/r*99.7%
*-commutative99.7%
*-commutative99.7%
times-frac99.7%
metadata-eval99.7%
associate-*l*99.7%
Simplified99.7%
Taylor expanded in a around 0 91.2%
+-commutative91.3%
*-commutative91.3%
associate-*r/91.3%
*-commutative91.3%
*-commutative91.3%
associate-*r*91.3%
associate-/l*91.3%
*-commutative91.3%
associate-*r/91.3%
Simplified91.2%
Taylor expanded in a around inf 1.6%
herbie shell --seed 2024098
(FPCore (a b c)
:name "Quadratic roots, medium range"
:precision binary64
:pre (and (and (and (< 1.1102230246251565e-16 a) (< a 9007199254740992.0)) (and (< 1.1102230246251565e-16 b) (< b 9007199254740992.0))) (and (< 1.1102230246251565e-16 c) (< c 9007199254740992.0)))
(/ (+ (- b) (sqrt (- (* b b) (* (* 4.0 a) c)))) (* 2.0 a)))