
(FPCore (a b c) :precision binary64 (/ (+ (- b) (sqrt (- (* b b) (* (* 3.0 a) c)))) (* 3.0 a)))
double code(double a, double b, double c) {
return (-b + sqrt(((b * b) - ((3.0 * a) * c)))) / (3.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) - ((3.0d0 * a) * c)))) / (3.0d0 * a)
end function
public static double code(double a, double b, double c) {
return (-b + Math.sqrt(((b * b) - ((3.0 * a) * c)))) / (3.0 * a);
}
def code(a, b, c): return (-b + math.sqrt(((b * b) - ((3.0 * a) * c)))) / (3.0 * a)
function code(a, b, c) return Float64(Float64(Float64(-b) + sqrt(Float64(Float64(b * b) - Float64(Float64(3.0 * a) * c)))) / Float64(3.0 * a)) end
function tmp = code(a, b, c) tmp = (-b + sqrt(((b * b) - ((3.0 * a) * c)))) / (3.0 * a); end
code[a_, b_, c_] := N[(N[((-b) + N[Sqrt[N[(N[(b * b), $MachinePrecision] - N[(N[(3.0 * a), $MachinePrecision] * c), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(3.0 * a), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\left(-b\right) + \sqrt{b \cdot b - \left(3 \cdot a\right) \cdot c}}{3 \cdot a}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 10 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (a b c) :precision binary64 (/ (+ (- b) (sqrt (- (* b b) (* (* 3.0 a) c)))) (* 3.0 a)))
double code(double a, double b, double c) {
return (-b + sqrt(((b * b) - ((3.0 * a) * c)))) / (3.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) - ((3.0d0 * a) * c)))) / (3.0d0 * a)
end function
public static double code(double a, double b, double c) {
return (-b + Math.sqrt(((b * b) - ((3.0 * a) * c)))) / (3.0 * a);
}
def code(a, b, c): return (-b + math.sqrt(((b * b) - ((3.0 * a) * c)))) / (3.0 * a)
function code(a, b, c) return Float64(Float64(Float64(-b) + sqrt(Float64(Float64(b * b) - Float64(Float64(3.0 * a) * c)))) / Float64(3.0 * a)) end
function tmp = code(a, b, c) tmp = (-b + sqrt(((b * b) - ((3.0 * a) * c)))) / (3.0 * a); end
code[a_, b_, c_] := N[(N[((-b) + N[Sqrt[N[(N[(b * b), $MachinePrecision] - N[(N[(3.0 * a), $MachinePrecision] * c), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(3.0 * a), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\left(-b\right) + \sqrt{b \cdot b - \left(3 \cdot a\right) \cdot c}}{3 \cdot a}
\end{array}
(FPCore (a b c) :precision binary64 (/ (* (* a c) -3.0) (* a (* 3.0 (+ b (sqrt (+ (* b b) (* a (* c -3.0)))))))))
double code(double a, double b, double c) {
return ((a * c) * -3.0) / (a * (3.0 * (b + sqrt(((b * b) + (a * (c * -3.0)))))));
}
real(8) function code(a, b, c)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
code = ((a * c) * (-3.0d0)) / (a * (3.0d0 * (b + sqrt(((b * b) + (a * (c * (-3.0d0))))))))
end function
public static double code(double a, double b, double c) {
return ((a * c) * -3.0) / (a * (3.0 * (b + Math.sqrt(((b * b) + (a * (c * -3.0)))))));
}
def code(a, b, c): return ((a * c) * -3.0) / (a * (3.0 * (b + math.sqrt(((b * b) + (a * (c * -3.0)))))))
function code(a, b, c) return Float64(Float64(Float64(a * c) * -3.0) / Float64(a * Float64(3.0 * Float64(b + sqrt(Float64(Float64(b * b) + Float64(a * Float64(c * -3.0)))))))) end
function tmp = code(a, b, c) tmp = ((a * c) * -3.0) / (a * (3.0 * (b + sqrt(((b * b) + (a * (c * -3.0))))))); end
code[a_, b_, c_] := N[(N[(N[(a * c), $MachinePrecision] * -3.0), $MachinePrecision] / N[(a * N[(3.0 * N[(b + N[Sqrt[N[(N[(b * b), $MachinePrecision] + N[(a * N[(c * -3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\left(a \cdot c\right) \cdot -3}{a \cdot \left(3 \cdot \left(b + \sqrt{b \cdot b + a \cdot \left(c \cdot -3\right)}\right)\right)}
\end{array}
Initial program 30.4%
/-lowering-/.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
sqrt-lowering-sqrt.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
*-commutativeN/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
metadata-evalN/A
*-lowering-*.f6430.4%
Simplified30.4%
flip--N/A
associate-/l/N/A
Applied egg-rr31.7%
Taylor expanded in b around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6499.0%
Simplified99.0%
Applied egg-rr99.1%
Taylor expanded in b around 0
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f6499.2%
Simplified99.2%
(FPCore (a b c)
:precision binary64
(/
(* a (+ (* 9.0 (* a (* c c))) (* -6.0 (* c (* b b)))))
(+
(* (* a 12.0) (* b (* b b)))
(*
c
(* 3.0 (+ (* (* a c) (* (/ (* a a) b) 2.25)) (* a (* (* a b) -9.0))))))))
double code(double a, double b, double c) {
return (a * ((9.0 * (a * (c * c))) + (-6.0 * (c * (b * b))))) / (((a * 12.0) * (b * (b * b))) + (c * (3.0 * (((a * c) * (((a * a) / b) * 2.25)) + (a * ((a * b) * -9.0))))));
}
real(8) function code(a, b, c)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
code = (a * ((9.0d0 * (a * (c * c))) + ((-6.0d0) * (c * (b * b))))) / (((a * 12.0d0) * (b * (b * b))) + (c * (3.0d0 * (((a * c) * (((a * a) / b) * 2.25d0)) + (a * ((a * b) * (-9.0d0)))))))
end function
public static double code(double a, double b, double c) {
return (a * ((9.0 * (a * (c * c))) + (-6.0 * (c * (b * b))))) / (((a * 12.0) * (b * (b * b))) + (c * (3.0 * (((a * c) * (((a * a) / b) * 2.25)) + (a * ((a * b) * -9.0))))));
}
def code(a, b, c): return (a * ((9.0 * (a * (c * c))) + (-6.0 * (c * (b * b))))) / (((a * 12.0) * (b * (b * b))) + (c * (3.0 * (((a * c) * (((a * a) / b) * 2.25)) + (a * ((a * b) * -9.0))))))
function code(a, b, c) return Float64(Float64(a * Float64(Float64(9.0 * Float64(a * Float64(c * c))) + Float64(-6.0 * Float64(c * Float64(b * b))))) / Float64(Float64(Float64(a * 12.0) * Float64(b * Float64(b * b))) + Float64(c * Float64(3.0 * Float64(Float64(Float64(a * c) * Float64(Float64(Float64(a * a) / b) * 2.25)) + Float64(a * Float64(Float64(a * b) * -9.0))))))) end
function tmp = code(a, b, c) tmp = (a * ((9.0 * (a * (c * c))) + (-6.0 * (c * (b * b))))) / (((a * 12.0) * (b * (b * b))) + (c * (3.0 * (((a * c) * (((a * a) / b) * 2.25)) + (a * ((a * b) * -9.0)))))); end
code[a_, b_, c_] := N[(N[(a * N[(N[(9.0 * N[(a * N[(c * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(-6.0 * N[(c * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(N[(a * 12.0), $MachinePrecision] * N[(b * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(c * N[(3.0 * N[(N[(N[(a * c), $MachinePrecision] * N[(N[(N[(a * a), $MachinePrecision] / b), $MachinePrecision] * 2.25), $MachinePrecision]), $MachinePrecision] + N[(a * N[(N[(a * b), $MachinePrecision] * -9.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{a \cdot \left(9 \cdot \left(a \cdot \left(c \cdot c\right)\right) + -6 \cdot \left(c \cdot \left(b \cdot b\right)\right)\right)}{\left(a \cdot 12\right) \cdot \left(b \cdot \left(b \cdot b\right)\right) + c \cdot \left(3 \cdot \left(\left(a \cdot c\right) \cdot \left(\frac{a \cdot a}{b} \cdot 2.25\right) + a \cdot \left(\left(a \cdot b\right) \cdot -9\right)\right)\right)}
\end{array}
Initial program 30.4%
/-lowering-/.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
sqrt-lowering-sqrt.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
*-commutativeN/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
metadata-evalN/A
*-lowering-*.f6430.4%
Simplified30.4%
flip--N/A
associate-/l/N/A
Applied egg-rr31.7%
Taylor expanded in c around 0
+-lowering-+.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
distribute-lft-outN/A
*-lowering-*.f64N/A
Simplified27.4%
Taylor expanded in a around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6494.4%
Simplified94.4%
Final simplification94.4%
(FPCore (a b c)
:precision binary64
(/
(+ (* (* c c) (* 9.0 (* a a))) (* -6.0 (* a (* c (* b b)))))
(*
a
(+
(* 12.0 (* b (* b b)))
(* a (* 3.0 (+ (* a (* 2.25 (/ (* c c) b))) (* -9.0 (* c b)))))))))
double code(double a, double b, double c) {
return (((c * c) * (9.0 * (a * a))) + (-6.0 * (a * (c * (b * b))))) / (a * ((12.0 * (b * (b * b))) + (a * (3.0 * ((a * (2.25 * ((c * c) / b))) + (-9.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 = (((c * c) * (9.0d0 * (a * a))) + ((-6.0d0) * (a * (c * (b * b))))) / (a * ((12.0d0 * (b * (b * b))) + (a * (3.0d0 * ((a * (2.25d0 * ((c * c) / b))) + ((-9.0d0) * (c * b)))))))
end function
public static double code(double a, double b, double c) {
return (((c * c) * (9.0 * (a * a))) + (-6.0 * (a * (c * (b * b))))) / (a * ((12.0 * (b * (b * b))) + (a * (3.0 * ((a * (2.25 * ((c * c) / b))) + (-9.0 * (c * b)))))));
}
def code(a, b, c): return (((c * c) * (9.0 * (a * a))) + (-6.0 * (a * (c * (b * b))))) / (a * ((12.0 * (b * (b * b))) + (a * (3.0 * ((a * (2.25 * ((c * c) / b))) + (-9.0 * (c * b)))))))
function code(a, b, c) return Float64(Float64(Float64(Float64(c * c) * Float64(9.0 * Float64(a * a))) + Float64(-6.0 * Float64(a * Float64(c * Float64(b * b))))) / Float64(a * Float64(Float64(12.0 * Float64(b * Float64(b * b))) + Float64(a * Float64(3.0 * Float64(Float64(a * Float64(2.25 * Float64(Float64(c * c) / b))) + Float64(-9.0 * Float64(c * b)))))))) end
function tmp = code(a, b, c) tmp = (((c * c) * (9.0 * (a * a))) + (-6.0 * (a * (c * (b * b))))) / (a * ((12.0 * (b * (b * b))) + (a * (3.0 * ((a * (2.25 * ((c * c) / b))) + (-9.0 * (c * b))))))); end
code[a_, b_, c_] := N[(N[(N[(N[(c * c), $MachinePrecision] * N[(9.0 * N[(a * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(-6.0 * N[(a * N[(c * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(a * N[(N[(12.0 * N[(b * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(a * N[(3.0 * N[(N[(a * N[(2.25 * N[(N[(c * c), $MachinePrecision] / b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(-9.0 * N[(c * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\left(c \cdot c\right) \cdot \left(9 \cdot \left(a \cdot a\right)\right) + -6 \cdot \left(a \cdot \left(c \cdot \left(b \cdot b\right)\right)\right)}{a \cdot \left(12 \cdot \left(b \cdot \left(b \cdot b\right)\right) + a \cdot \left(3 \cdot \left(a \cdot \left(2.25 \cdot \frac{c \cdot c}{b}\right) + -9 \cdot \left(c \cdot b\right)\right)\right)\right)}
\end{array}
Initial program 30.4%
/-lowering-/.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
sqrt-lowering-sqrt.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
*-commutativeN/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
metadata-evalN/A
*-lowering-*.f6430.4%
Simplified30.4%
flip--N/A
associate-/l/N/A
Applied egg-rr31.7%
Taylor expanded in b around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6499.0%
Simplified99.0%
Taylor expanded in a around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
distribute-lft-outN/A
Simplified94.4%
Taylor expanded in b around 0
+-commutativeN/A
+-lowering-+.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6494.4%
Simplified94.4%
Final simplification94.4%
(FPCore (a b c)
:precision binary64
(/
(* a (+ (* 9.0 (* a (* c c))) (* -6.0 (* c (* b b)))))
(*
a
(+
(* 12.0 (* b (* b b)))
(* a (* 3.0 (+ (* a (* 2.25 (/ (* c c) b))) (* -9.0 (* c b)))))))))
double code(double a, double b, double c) {
return (a * ((9.0 * (a * (c * c))) + (-6.0 * (c * (b * b))))) / (a * ((12.0 * (b * (b * b))) + (a * (3.0 * ((a * (2.25 * ((c * c) / b))) + (-9.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 = (a * ((9.0d0 * (a * (c * c))) + ((-6.0d0) * (c * (b * b))))) / (a * ((12.0d0 * (b * (b * b))) + (a * (3.0d0 * ((a * (2.25d0 * ((c * c) / b))) + ((-9.0d0) * (c * b)))))))
end function
public static double code(double a, double b, double c) {
return (a * ((9.0 * (a * (c * c))) + (-6.0 * (c * (b * b))))) / (a * ((12.0 * (b * (b * b))) + (a * (3.0 * ((a * (2.25 * ((c * c) / b))) + (-9.0 * (c * b)))))));
}
def code(a, b, c): return (a * ((9.0 * (a * (c * c))) + (-6.0 * (c * (b * b))))) / (a * ((12.0 * (b * (b * b))) + (a * (3.0 * ((a * (2.25 * ((c * c) / b))) + (-9.0 * (c * b)))))))
function code(a, b, c) return Float64(Float64(a * Float64(Float64(9.0 * Float64(a * Float64(c * c))) + Float64(-6.0 * Float64(c * Float64(b * b))))) / Float64(a * Float64(Float64(12.0 * Float64(b * Float64(b * b))) + Float64(a * Float64(3.0 * Float64(Float64(a * Float64(2.25 * Float64(Float64(c * c) / b))) + Float64(-9.0 * Float64(c * b)))))))) end
function tmp = code(a, b, c) tmp = (a * ((9.0 * (a * (c * c))) + (-6.0 * (c * (b * b))))) / (a * ((12.0 * (b * (b * b))) + (a * (3.0 * ((a * (2.25 * ((c * c) / b))) + (-9.0 * (c * b))))))); end
code[a_, b_, c_] := N[(N[(a * N[(N[(9.0 * N[(a * N[(c * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(-6.0 * N[(c * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(a * N[(N[(12.0 * N[(b * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(a * N[(3.0 * N[(N[(a * N[(2.25 * N[(N[(c * c), $MachinePrecision] / b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(-9.0 * N[(c * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{a \cdot \left(9 \cdot \left(a \cdot \left(c \cdot c\right)\right) + -6 \cdot \left(c \cdot \left(b \cdot b\right)\right)\right)}{a \cdot \left(12 \cdot \left(b \cdot \left(b \cdot b\right)\right) + a \cdot \left(3 \cdot \left(a \cdot \left(2.25 \cdot \frac{c \cdot c}{b}\right) + -9 \cdot \left(c \cdot b\right)\right)\right)\right)}
\end{array}
Initial program 30.4%
/-lowering-/.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
sqrt-lowering-sqrt.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
*-commutativeN/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
metadata-evalN/A
*-lowering-*.f6430.4%
Simplified30.4%
flip--N/A
associate-/l/N/A
Applied egg-rr31.7%
Taylor expanded in b around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6499.0%
Simplified99.0%
Taylor expanded in a around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
distribute-lft-outN/A
Simplified94.4%
Taylor expanded in a around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6494.3%
Simplified94.3%
Final simplification94.3%
(FPCore (a b c) :precision binary64 (/ (+ (+ (/ (* (* a (* c c)) -0.375) (* b b)) (* c -0.5)) (/ (* (* (* a a) -0.5625) (* c (* c c))) (* (* b b) (* b b)))) b))
double code(double a, double b, double c) {
return (((((a * (c * c)) * -0.375) / (b * b)) + (c * -0.5)) + ((((a * a) * -0.5625) * (c * (c * c))) / ((b * b) * (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 = (((((a * (c * c)) * (-0.375d0)) / (b * b)) + (c * (-0.5d0))) + ((((a * a) * (-0.5625d0)) * (c * (c * c))) / ((b * b) * (b * b)))) / b
end function
public static double code(double a, double b, double c) {
return (((((a * (c * c)) * -0.375) / (b * b)) + (c * -0.5)) + ((((a * a) * -0.5625) * (c * (c * c))) / ((b * b) * (b * b)))) / b;
}
def code(a, b, c): return (((((a * (c * c)) * -0.375) / (b * b)) + (c * -0.5)) + ((((a * a) * -0.5625) * (c * (c * c))) / ((b * b) * (b * b)))) / b
function code(a, b, c) return Float64(Float64(Float64(Float64(Float64(Float64(a * Float64(c * c)) * -0.375) / Float64(b * b)) + Float64(c * -0.5)) + Float64(Float64(Float64(Float64(a * a) * -0.5625) * Float64(c * Float64(c * c))) / Float64(Float64(b * b) * Float64(b * b)))) / b) end
function tmp = code(a, b, c) tmp = (((((a * (c * c)) * -0.375) / (b * b)) + (c * -0.5)) + ((((a * a) * -0.5625) * (c * (c * c))) / ((b * b) * (b * b)))) / b; end
code[a_, b_, c_] := N[(N[(N[(N[(N[(N[(a * N[(c * c), $MachinePrecision]), $MachinePrecision] * -0.375), $MachinePrecision] / N[(b * b), $MachinePrecision]), $MachinePrecision] + N[(c * -0.5), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(N[(a * a), $MachinePrecision] * -0.5625), $MachinePrecision] * N[(c * N[(c * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(b * b), $MachinePrecision] * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / b), $MachinePrecision]
\begin{array}{l}
\\
\frac{\left(\frac{\left(a \cdot \left(c \cdot c\right)\right) \cdot -0.375}{b \cdot b} + c \cdot -0.5\right) + \frac{\left(\left(a \cdot a\right) \cdot -0.5625\right) \cdot \left(c \cdot \left(c \cdot c\right)\right)}{\left(b \cdot b\right) \cdot \left(b \cdot b\right)}}{b}
\end{array}
Initial program 30.4%
/-lowering-/.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
sqrt-lowering-sqrt.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
*-commutativeN/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
metadata-evalN/A
*-lowering-*.f6430.4%
Simplified30.4%
Taylor expanded in a around 0
Simplified94.6%
Taylor expanded in b around inf
/-lowering-/.f64N/A
Simplified93.2%
Final simplification93.2%
(FPCore (a b c)
:precision binary64
(+
(/ (* c -0.5) b)
(*
a
(/
(+ (/ (* -0.5625 (* a (* c (* c c)))) (* b b)) (* (* c c) -0.375))
(* b (* b b))))))
double code(double a, double b, double c) {
return ((c * -0.5) / b) + (a * ((((-0.5625 * (a * (c * (c * c)))) / (b * b)) + ((c * c) * -0.375)) / (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 = ((c * (-0.5d0)) / b) + (a * (((((-0.5625d0) * (a * (c * (c * c)))) / (b * b)) + ((c * c) * (-0.375d0))) / (b * (b * b))))
end function
public static double code(double a, double b, double c) {
return ((c * -0.5) / b) + (a * ((((-0.5625 * (a * (c * (c * c)))) / (b * b)) + ((c * c) * -0.375)) / (b * (b * b))));
}
def code(a, b, c): return ((c * -0.5) / b) + (a * ((((-0.5625 * (a * (c * (c * c)))) / (b * b)) + ((c * c) * -0.375)) / (b * (b * b))))
function code(a, b, c) return Float64(Float64(Float64(c * -0.5) / b) + Float64(a * Float64(Float64(Float64(Float64(-0.5625 * Float64(a * Float64(c * Float64(c * c)))) / Float64(b * b)) + Float64(Float64(c * c) * -0.375)) / Float64(b * Float64(b * b))))) end
function tmp = code(a, b, c) tmp = ((c * -0.5) / b) + (a * ((((-0.5625 * (a * (c * (c * c)))) / (b * b)) + ((c * c) * -0.375)) / (b * (b * b)))); end
code[a_, b_, c_] := N[(N[(N[(c * -0.5), $MachinePrecision] / b), $MachinePrecision] + N[(a * N[(N[(N[(N[(-0.5625 * N[(a * N[(c * N[(c * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(b * b), $MachinePrecision]), $MachinePrecision] + N[(N[(c * c), $MachinePrecision] * -0.375), $MachinePrecision]), $MachinePrecision] / N[(b * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{c \cdot -0.5}{b} + a \cdot \frac{\frac{-0.5625 \cdot \left(a \cdot \left(c \cdot \left(c \cdot c\right)\right)\right)}{b \cdot b} + \left(c \cdot c\right) \cdot -0.375}{b \cdot \left(b \cdot b\right)}
\end{array}
Initial program 30.4%
/-lowering-/.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
sqrt-lowering-sqrt.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
*-commutativeN/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
metadata-evalN/A
*-lowering-*.f6430.4%
Simplified30.4%
Taylor expanded in a around 0
Simplified94.6%
Taylor expanded in b around inf
/-lowering-/.f64N/A
+-lowering-+.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6493.2%
Simplified93.2%
Final simplification93.2%
(FPCore (a b c) :precision binary64 (/ (+ (* c -0.5) (/ (* -0.375 (* c (* a c))) (* b b))) b))
double code(double a, double b, double c) {
return ((c * -0.5) + ((-0.375 * (c * (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 = ((c * (-0.5d0)) + (((-0.375d0) * (c * (a * c))) / (b * b))) / b
end function
public static double code(double a, double b, double c) {
return ((c * -0.5) + ((-0.375 * (c * (a * c))) / (b * b))) / b;
}
def code(a, b, c): return ((c * -0.5) + ((-0.375 * (c * (a * c))) / (b * b))) / b
function code(a, b, c) return Float64(Float64(Float64(c * -0.5) + Float64(Float64(-0.375 * Float64(c * Float64(a * c))) / Float64(b * b))) / b) end
function tmp = code(a, b, c) tmp = ((c * -0.5) + ((-0.375 * (c * (a * c))) / (b * b))) / b; end
code[a_, b_, c_] := N[(N[(N[(c * -0.5), $MachinePrecision] + N[(N[(-0.375 * N[(c * N[(a * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / b), $MachinePrecision]
\begin{array}{l}
\\
\frac{c \cdot -0.5 + \frac{-0.375 \cdot \left(c \cdot \left(a \cdot c\right)\right)}{b \cdot b}}{b}
\end{array}
Initial program 30.4%
/-lowering-/.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
sqrt-lowering-sqrt.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
*-commutativeN/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
metadata-evalN/A
*-lowering-*.f6430.4%
Simplified30.4%
Taylor expanded in b around inf
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6490.7%
Simplified90.7%
Final simplification90.7%
(FPCore (a b c) :precision binary64 (* c (+ (/ -0.5 b) (/ (* (* a c) -0.375) (* b (* b b))))))
double code(double a, double b, double c) {
return c * ((-0.5 / b) + (((a * c) * -0.375) / (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 = c * (((-0.5d0) / b) + (((a * c) * (-0.375d0)) / (b * (b * b))))
end function
public static double code(double a, double b, double c) {
return c * ((-0.5 / b) + (((a * c) * -0.375) / (b * (b * b))));
}
def code(a, b, c): return c * ((-0.5 / b) + (((a * c) * -0.375) / (b * (b * b))))
function code(a, b, c) return Float64(c * Float64(Float64(-0.5 / b) + Float64(Float64(Float64(a * c) * -0.375) / Float64(b * Float64(b * b))))) end
function tmp = code(a, b, c) tmp = c * ((-0.5 / b) + (((a * c) * -0.375) / (b * (b * b)))); end
code[a_, b_, c_] := N[(c * N[(N[(-0.5 / b), $MachinePrecision] + N[(N[(N[(a * c), $MachinePrecision] * -0.375), $MachinePrecision] / N[(b * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
c \cdot \left(\frac{-0.5}{b} + \frac{\left(a \cdot c\right) \cdot -0.375}{b \cdot \left(b \cdot b\right)}\right)
\end{array}
Initial program 30.4%
/-lowering-/.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
sqrt-lowering-sqrt.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
*-commutativeN/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
metadata-evalN/A
*-lowering-*.f6430.4%
Simplified30.4%
Taylor expanded in c around 0
sub-negN/A
associate-*r/N/A
associate-*r*N/A
associate-*l/N/A
associate-*r/N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f64N/A
associate-*r/N/A
associate-*l/N/A
associate-*r*N/A
/-lowering-/.f64N/A
Simplified90.4%
Final simplification90.4%
(FPCore (a b c) :precision binary64 (/ (* c -0.5) b))
double code(double a, double b, double c) {
return (c * -0.5) / 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 * (-0.5d0)) / b
end function
public static double code(double a, double b, double c) {
return (c * -0.5) / b;
}
def code(a, b, c): return (c * -0.5) / b
function code(a, b, c) return Float64(Float64(c * -0.5) / b) end
function tmp = code(a, b, c) tmp = (c * -0.5) / b; end
code[a_, b_, c_] := N[(N[(c * -0.5), $MachinePrecision] / b), $MachinePrecision]
\begin{array}{l}
\\
\frac{c \cdot -0.5}{b}
\end{array}
Initial program 30.4%
/-lowering-/.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
sqrt-lowering-sqrt.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
*-commutativeN/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
metadata-evalN/A
*-lowering-*.f6430.4%
Simplified30.4%
Taylor expanded in b around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6482.1%
Simplified82.1%
(FPCore (a b c) :precision binary64 (* c (/ -0.5 b)))
double code(double a, double b, double c) {
return c * (-0.5 / 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 * ((-0.5d0) / b)
end function
public static double code(double a, double b, double c) {
return c * (-0.5 / b);
}
def code(a, b, c): return c * (-0.5 / b)
function code(a, b, c) return Float64(c * Float64(-0.5 / b)) end
function tmp = code(a, b, c) tmp = c * (-0.5 / b); end
code[a_, b_, c_] := N[(c * N[(-0.5 / b), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
c \cdot \frac{-0.5}{b}
\end{array}
Initial program 30.4%
/-lowering-/.f64N/A
+-commutativeN/A
unsub-negN/A
--lowering--.f64N/A
sqrt-lowering-sqrt.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
*-commutativeN/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
metadata-evalN/A
*-lowering-*.f6430.4%
Simplified30.4%
Taylor expanded in b around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6482.1%
Simplified82.1%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f6481.9%
Applied egg-rr81.9%
Final simplification81.9%
herbie shell --seed 2024161
(FPCore (a b c)
:name "Cubic critical, 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) (* (* 3.0 a) c)))) (* 3.0 a)))