
(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 15 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
(let* ((t_0 (* c (* a -3.0))))
(/
(* c (+ (* -2.0 (* b b)) (* 3.0 (* c a))))
(* (+ t_0 (* (* b b) 2.0)) (+ b (sqrt (+ (* b b) t_0)))))))
double code(double a, double b, double c) {
double t_0 = c * (a * -3.0);
return (c * ((-2.0 * (b * b)) + (3.0 * (c * a)))) / ((t_0 + ((b * b) * 2.0)) * (b + sqrt(((b * b) + t_0))));
}
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
t_0 = c * (a * (-3.0d0))
code = (c * (((-2.0d0) * (b * b)) + (3.0d0 * (c * a)))) / ((t_0 + ((b * b) * 2.0d0)) * (b + sqrt(((b * b) + t_0))))
end function
public static double code(double a, double b, double c) {
double t_0 = c * (a * -3.0);
return (c * ((-2.0 * (b * b)) + (3.0 * (c * a)))) / ((t_0 + ((b * b) * 2.0)) * (b + Math.sqrt(((b * b) + t_0))));
}
def code(a, b, c): t_0 = c * (a * -3.0) return (c * ((-2.0 * (b * b)) + (3.0 * (c * a)))) / ((t_0 + ((b * b) * 2.0)) * (b + math.sqrt(((b * b) + t_0))))
function code(a, b, c) t_0 = Float64(c * Float64(a * -3.0)) return Float64(Float64(c * Float64(Float64(-2.0 * Float64(b * b)) + Float64(3.0 * Float64(c * a)))) / Float64(Float64(t_0 + Float64(Float64(b * b) * 2.0)) * Float64(b + sqrt(Float64(Float64(b * b) + t_0))))) end
function tmp = code(a, b, c) t_0 = c * (a * -3.0); tmp = (c * ((-2.0 * (b * b)) + (3.0 * (c * a)))) / ((t_0 + ((b * b) * 2.0)) * (b + sqrt(((b * b) + t_0)))); end
code[a_, b_, c_] := Block[{t$95$0 = N[(c * N[(a * -3.0), $MachinePrecision]), $MachinePrecision]}, N[(N[(c * N[(N[(-2.0 * N[(b * b), $MachinePrecision]), $MachinePrecision] + N[(3.0 * N[(c * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(t$95$0 + N[(N[(b * b), $MachinePrecision] * 2.0), $MachinePrecision]), $MachinePrecision] * N[(b + N[Sqrt[N[(N[(b * b), $MachinePrecision] + t$95$0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := c \cdot \left(a \cdot -3\right)\\
\frac{c \cdot \left(-2 \cdot \left(b \cdot b\right) + 3 \cdot \left(c \cdot a\right)\right)}{\left(t\_0 + \left(b \cdot b\right) \cdot 2\right) \cdot \left(b + \sqrt{b \cdot b + t\_0}\right)}
\end{array}
\end{array}
Initial program 33.2%
/-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-*.f6433.2%
Simplified33.2%
flip--N/A
associate-/l/N/A
Applied egg-rr34.6%
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
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6499.0%
Simplified99.0%
Applied egg-rr99.0%
Taylor expanded in c around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6499.4%
Simplified99.4%
Final simplification99.4%
(FPCore (a b c)
:precision binary64
(let* ((t_0 (* c (* a -3.0))))
(*
(+ t_0 (- (* b b) (* b b)))
(/ (/ 0.3333333333333333 a) (+ b (sqrt (+ (* b b) t_0)))))))
double code(double a, double b, double c) {
double t_0 = c * (a * -3.0);
return (t_0 + ((b * b) - (b * b))) * ((0.3333333333333333 / a) / (b + sqrt(((b * b) + t_0))));
}
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
t_0 = c * (a * (-3.0d0))
code = (t_0 + ((b * b) - (b * b))) * ((0.3333333333333333d0 / a) / (b + sqrt(((b * b) + t_0))))
end function
public static double code(double a, double b, double c) {
double t_0 = c * (a * -3.0);
return (t_0 + ((b * b) - (b * b))) * ((0.3333333333333333 / a) / (b + Math.sqrt(((b * b) + t_0))));
}
def code(a, b, c): t_0 = c * (a * -3.0) return (t_0 + ((b * b) - (b * b))) * ((0.3333333333333333 / a) / (b + math.sqrt(((b * b) + t_0))))
function code(a, b, c) t_0 = Float64(c * Float64(a * -3.0)) return Float64(Float64(t_0 + Float64(Float64(b * b) - Float64(b * b))) * Float64(Float64(0.3333333333333333 / a) / Float64(b + sqrt(Float64(Float64(b * b) + t_0))))) end
function tmp = code(a, b, c) t_0 = c * (a * -3.0); tmp = (t_0 + ((b * b) - (b * b))) * ((0.3333333333333333 / a) / (b + sqrt(((b * b) + t_0)))); end
code[a_, b_, c_] := Block[{t$95$0 = N[(c * N[(a * -3.0), $MachinePrecision]), $MachinePrecision]}, N[(N[(t$95$0 + N[(N[(b * b), $MachinePrecision] - N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(0.3333333333333333 / a), $MachinePrecision] / N[(b + N[Sqrt[N[(N[(b * b), $MachinePrecision] + t$95$0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := c \cdot \left(a \cdot -3\right)\\
\left(t\_0 + \left(b \cdot b - b \cdot b\right)\right) \cdot \frac{\frac{0.3333333333333333}{a}}{b + \sqrt{b \cdot b + t\_0}}
\end{array}
\end{array}
Initial program 33.2%
/-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-*.f6433.2%
Simplified33.2%
div-invN/A
flip--N/A
associate-*l/N/A
associate-/l*N/A
*-lowering-*.f64N/A
Applied egg-rr34.3%
+-commutativeN/A
associate-+l-N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6499.0%
Applied egg-rr99.0%
Final simplification99.0%
(FPCore (a b c) :precision binary64 (* (/ (/ 0.3333333333333333 a) (+ b (sqrt (+ (* b b) (* c (* a -3.0)))))) (* (* c a) -3.0)))
double code(double a, double b, double c) {
return ((0.3333333333333333 / a) / (b + sqrt(((b * b) + (c * (a * -3.0)))))) * ((c * a) * -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 = ((0.3333333333333333d0 / a) / (b + sqrt(((b * b) + (c * (a * (-3.0d0))))))) * ((c * a) * (-3.0d0))
end function
public static double code(double a, double b, double c) {
return ((0.3333333333333333 / a) / (b + Math.sqrt(((b * b) + (c * (a * -3.0)))))) * ((c * a) * -3.0);
}
def code(a, b, c): return ((0.3333333333333333 / a) / (b + math.sqrt(((b * b) + (c * (a * -3.0)))))) * ((c * a) * -3.0)
function code(a, b, c) return Float64(Float64(Float64(0.3333333333333333 / a) / Float64(b + sqrt(Float64(Float64(b * b) + Float64(c * Float64(a * -3.0)))))) * Float64(Float64(c * a) * -3.0)) end
function tmp = code(a, b, c) tmp = ((0.3333333333333333 / a) / (b + sqrt(((b * b) + (c * (a * -3.0)))))) * ((c * a) * -3.0); end
code[a_, b_, c_] := N[(N[(N[(0.3333333333333333 / a), $MachinePrecision] / N[(b + N[Sqrt[N[(N[(b * b), $MachinePrecision] + N[(c * N[(a * -3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(c * a), $MachinePrecision] * -3.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\frac{0.3333333333333333}{a}}{b + \sqrt{b \cdot b + c \cdot \left(a \cdot -3\right)}} \cdot \left(\left(c \cdot a\right) \cdot -3\right)
\end{array}
Initial program 33.2%
/-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-*.f6433.2%
Simplified33.2%
div-invN/A
flip--N/A
associate-*l/N/A
associate-/l*N/A
*-lowering-*.f64N/A
Applied egg-rr34.3%
Taylor expanded in b around 0
*-lowering-*.f64N/A
*-lowering-*.f6498.9%
Simplified98.9%
Final simplification98.9%
(FPCore (a b c)
:precision binary64
(/
(* (* b b) (+ (* (* c a) -6.0) (/ (* (* 9.0 (* a a)) (* c c)) (* b b))))
(*
a
(+
(* 12.0 (* b (* b b)))
(* a (* 3.0 (+ (* a (* (/ (* c c) b) 2.25)) (* -9.0 (* c b)))))))))
double code(double a, double b, double c) {
return ((b * b) * (((c * a) * -6.0) + (((9.0 * (a * a)) * (c * c)) / (b * b)))) / (a * ((12.0 * (b * (b * b))) + (a * (3.0 * ((a * (((c * c) / b) * 2.25)) + (-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 = ((b * b) * (((c * a) * (-6.0d0)) + (((9.0d0 * (a * a)) * (c * c)) / (b * b)))) / (a * ((12.0d0 * (b * (b * b))) + (a * (3.0d0 * ((a * (((c * c) / b) * 2.25d0)) + ((-9.0d0) * (c * b)))))))
end function
public static double code(double a, double b, double c) {
return ((b * b) * (((c * a) * -6.0) + (((9.0 * (a * a)) * (c * c)) / (b * b)))) / (a * ((12.0 * (b * (b * b))) + (a * (3.0 * ((a * (((c * c) / b) * 2.25)) + (-9.0 * (c * b)))))));
}
def code(a, b, c): return ((b * b) * (((c * a) * -6.0) + (((9.0 * (a * a)) * (c * c)) / (b * b)))) / (a * ((12.0 * (b * (b * b))) + (a * (3.0 * ((a * (((c * c) / b) * 2.25)) + (-9.0 * (c * b)))))))
function code(a, b, c) return Float64(Float64(Float64(b * b) * Float64(Float64(Float64(c * a) * -6.0) + Float64(Float64(Float64(9.0 * Float64(a * a)) * Float64(c * c)) / Float64(b * b)))) / Float64(a * Float64(Float64(12.0 * Float64(b * Float64(b * b))) + Float64(a * Float64(3.0 * Float64(Float64(a * Float64(Float64(Float64(c * c) / b) * 2.25)) + Float64(-9.0 * Float64(c * b)))))))) end
function tmp = code(a, b, c) tmp = ((b * b) * (((c * a) * -6.0) + (((9.0 * (a * a)) * (c * c)) / (b * b)))) / (a * ((12.0 * (b * (b * b))) + (a * (3.0 * ((a * (((c * c) / b) * 2.25)) + (-9.0 * (c * b))))))); end
code[a_, b_, c_] := N[(N[(N[(b * b), $MachinePrecision] * N[(N[(N[(c * a), $MachinePrecision] * -6.0), $MachinePrecision] + N[(N[(N[(9.0 * N[(a * a), $MachinePrecision]), $MachinePrecision] * N[(c * c), $MachinePrecision]), $MachinePrecision] / N[(b * b), $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[(N[(N[(c * c), $MachinePrecision] / b), $MachinePrecision] * 2.25), $MachinePrecision]), $MachinePrecision] + N[(-9.0 * N[(c * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\left(b \cdot b\right) \cdot \left(\left(c \cdot a\right) \cdot -6 + \frac{\left(9 \cdot \left(a \cdot a\right)\right) \cdot \left(c \cdot c\right)}{b \cdot b}\right)}{a \cdot \left(12 \cdot \left(b \cdot \left(b \cdot b\right)\right) + a \cdot \left(3 \cdot \left(a \cdot \left(\frac{c \cdot c}{b} \cdot 2.25\right) + -9 \cdot \left(c \cdot b\right)\right)\right)\right)}
\end{array}
Initial program 33.2%
/-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-*.f6433.2%
Simplified33.2%
flip--N/A
associate-/l/N/A
Applied egg-rr34.6%
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
Simplified31.4%
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-*.f6495.4%
Simplified95.4%
Final simplification95.4%
(FPCore (a b c) :precision binary64 (/ (/ (* c (+ (* 9.0 (* c (* a a))) (* -6.0 (* (* b b) a)))) (* 3.0 a)) (+ (* (* b (* b b)) 4.0) (* c (+ (* -9.0 (* b a)) (* c (* 2.25 (/ (* a a) b))))))))
double code(double a, double b, double c) {
return ((c * ((9.0 * (c * (a * a))) + (-6.0 * ((b * b) * a)))) / (3.0 * a)) / (((b * (b * b)) * 4.0) + (c * ((-9.0 * (b * a)) + (c * (2.25 * ((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 = ((c * ((9.0d0 * (c * (a * a))) + ((-6.0d0) * ((b * b) * a)))) / (3.0d0 * a)) / (((b * (b * b)) * 4.0d0) + (c * (((-9.0d0) * (b * a)) + (c * (2.25d0 * ((a * a) / b))))))
end function
public static double code(double a, double b, double c) {
return ((c * ((9.0 * (c * (a * a))) + (-6.0 * ((b * b) * a)))) / (3.0 * a)) / (((b * (b * b)) * 4.0) + (c * ((-9.0 * (b * a)) + (c * (2.25 * ((a * a) / b))))));
}
def code(a, b, c): return ((c * ((9.0 * (c * (a * a))) + (-6.0 * ((b * b) * a)))) / (3.0 * a)) / (((b * (b * b)) * 4.0) + (c * ((-9.0 * (b * a)) + (c * (2.25 * ((a * a) / b))))))
function code(a, b, c) return Float64(Float64(Float64(c * Float64(Float64(9.0 * Float64(c * Float64(a * a))) + Float64(-6.0 * Float64(Float64(b * b) * a)))) / Float64(3.0 * a)) / Float64(Float64(Float64(b * Float64(b * b)) * 4.0) + Float64(c * Float64(Float64(-9.0 * Float64(b * a)) + Float64(c * Float64(2.25 * Float64(Float64(a * a) / b))))))) end
function tmp = code(a, b, c) tmp = ((c * ((9.0 * (c * (a * a))) + (-6.0 * ((b * b) * a)))) / (3.0 * a)) / (((b * (b * b)) * 4.0) + (c * ((-9.0 * (b * a)) + (c * (2.25 * ((a * a) / b)))))); end
code[a_, b_, c_] := N[(N[(N[(c * N[(N[(9.0 * N[(c * N[(a * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(-6.0 * N[(N[(b * b), $MachinePrecision] * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(3.0 * a), $MachinePrecision]), $MachinePrecision] / N[(N[(N[(b * N[(b * b), $MachinePrecision]), $MachinePrecision] * 4.0), $MachinePrecision] + N[(c * N[(N[(-9.0 * N[(b * a), $MachinePrecision]), $MachinePrecision] + N[(c * N[(2.25 * N[(N[(a * a), $MachinePrecision] / b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\frac{c \cdot \left(9 \cdot \left(c \cdot \left(a \cdot a\right)\right) + -6 \cdot \left(\left(b \cdot b\right) \cdot a\right)\right)}{3 \cdot a}}{\left(b \cdot \left(b \cdot b\right)\right) \cdot 4 + c \cdot \left(-9 \cdot \left(b \cdot a\right) + c \cdot \left(2.25 \cdot \frac{a \cdot a}{b}\right)\right)}
\end{array}
Initial program 33.2%
/-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-*.f6433.2%
Simplified33.2%
flip--N/A
associate-/l/N/A
Applied egg-rr34.6%
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
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6499.0%
Simplified99.0%
Applied egg-rr99.0%
Taylor expanded in c around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-+r+N/A
+-lowering-+.f64N/A
distribute-rgt-outN/A
metadata-evalN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
Simplified95.3%
Final simplification95.3%
(FPCore (a b c)
:precision binary64
(/
(* c (+ (* -6.0 (* (* b b) a)) (* c (* 9.0 (* a a)))))
(*
a
(+
(* 12.0 (* b (* b b)))
(* a (* 3.0 (+ (* a (* (/ (* c c) b) 2.25)) (* -9.0 (* c b)))))))))
double code(double a, double b, double c) {
return (c * ((-6.0 * ((b * b) * a)) + (c * (9.0 * (a * a))))) / (a * ((12.0 * (b * (b * b))) + (a * (3.0 * ((a * (((c * c) / b) * 2.25)) + (-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 * (((-6.0d0) * ((b * b) * a)) + (c * (9.0d0 * (a * a))))) / (a * ((12.0d0 * (b * (b * b))) + (a * (3.0d0 * ((a * (((c * c) / b) * 2.25d0)) + ((-9.0d0) * (c * b)))))))
end function
public static double code(double a, double b, double c) {
return (c * ((-6.0 * ((b * b) * a)) + (c * (9.0 * (a * a))))) / (a * ((12.0 * (b * (b * b))) + (a * (3.0 * ((a * (((c * c) / b) * 2.25)) + (-9.0 * (c * b)))))));
}
def code(a, b, c): return (c * ((-6.0 * ((b * b) * a)) + (c * (9.0 * (a * a))))) / (a * ((12.0 * (b * (b * b))) + (a * (3.0 * ((a * (((c * c) / b) * 2.25)) + (-9.0 * (c * b)))))))
function code(a, b, c) return Float64(Float64(c * Float64(Float64(-6.0 * Float64(Float64(b * b) * a)) + Float64(c * Float64(9.0 * Float64(a * a))))) / Float64(a * Float64(Float64(12.0 * Float64(b * Float64(b * b))) + Float64(a * Float64(3.0 * Float64(Float64(a * Float64(Float64(Float64(c * c) / b) * 2.25)) + Float64(-9.0 * Float64(c * b)))))))) end
function tmp = code(a, b, c) tmp = (c * ((-6.0 * ((b * b) * a)) + (c * (9.0 * (a * a))))) / (a * ((12.0 * (b * (b * b))) + (a * (3.0 * ((a * (((c * c) / b) * 2.25)) + (-9.0 * (c * b))))))); end
code[a_, b_, c_] := N[(N[(c * N[(N[(-6.0 * N[(N[(b * b), $MachinePrecision] * a), $MachinePrecision]), $MachinePrecision] + N[(c * N[(9.0 * N[(a * a), $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[(N[(N[(c * c), $MachinePrecision] / b), $MachinePrecision] * 2.25), $MachinePrecision]), $MachinePrecision] + N[(-9.0 * N[(c * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{c \cdot \left(-6 \cdot \left(\left(b \cdot b\right) \cdot a\right) + c \cdot \left(9 \cdot \left(a \cdot a\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(\frac{c \cdot c}{b} \cdot 2.25\right) + -9 \cdot \left(c \cdot b\right)\right)\right)\right)}
\end{array}
Initial program 33.2%
/-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-*.f6433.2%
Simplified33.2%
flip--N/A
associate-/l/N/A
Applied egg-rr34.6%
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
Simplified31.4%
Taylor expanded in c around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6495.3%
Simplified95.3%
Final simplification95.3%
(FPCore (a b c)
:precision binary64
(/
(/
(* c (+ (* -6.0 (* (* b b) a)) (* 9.0 (* a (* c a)))))
(+
(* 12.0 (* b (* b b)))
(* a (+ (* (/ a (/ b (* c c))) 6.75) (* (* c b) -27.0)))))
a))
double code(double a, double b, double c) {
return ((c * ((-6.0 * ((b * b) * a)) + (9.0 * (a * (c * a))))) / ((12.0 * (b * (b * b))) + (a * (((a / (b / (c * c))) * 6.75) + ((c * b) * -27.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 = ((c * (((-6.0d0) * ((b * b) * a)) + (9.0d0 * (a * (c * a))))) / ((12.0d0 * (b * (b * b))) + (a * (((a / (b / (c * c))) * 6.75d0) + ((c * b) * (-27.0d0)))))) / a
end function
public static double code(double a, double b, double c) {
return ((c * ((-6.0 * ((b * b) * a)) + (9.0 * (a * (c * a))))) / ((12.0 * (b * (b * b))) + (a * (((a / (b / (c * c))) * 6.75) + ((c * b) * -27.0))))) / a;
}
def code(a, b, c): return ((c * ((-6.0 * ((b * b) * a)) + (9.0 * (a * (c * a))))) / ((12.0 * (b * (b * b))) + (a * (((a / (b / (c * c))) * 6.75) + ((c * b) * -27.0))))) / a
function code(a, b, c) return Float64(Float64(Float64(c * Float64(Float64(-6.0 * Float64(Float64(b * b) * a)) + Float64(9.0 * Float64(a * Float64(c * a))))) / Float64(Float64(12.0 * Float64(b * Float64(b * b))) + Float64(a * Float64(Float64(Float64(a / Float64(b / Float64(c * c))) * 6.75) + Float64(Float64(c * b) * -27.0))))) / a) end
function tmp = code(a, b, c) tmp = ((c * ((-6.0 * ((b * b) * a)) + (9.0 * (a * (c * a))))) / ((12.0 * (b * (b * b))) + (a * (((a / (b / (c * c))) * 6.75) + ((c * b) * -27.0))))) / a; end
code[a_, b_, c_] := N[(N[(N[(c * N[(N[(-6.0 * N[(N[(b * b), $MachinePrecision] * a), $MachinePrecision]), $MachinePrecision] + N[(9.0 * N[(a * N[(c * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(12.0 * N[(b * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(a * N[(N[(N[(a / N[(b / N[(c * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * 6.75), $MachinePrecision] + N[(N[(c * b), $MachinePrecision] * -27.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision]
\begin{array}{l}
\\
\frac{\frac{c \cdot \left(-6 \cdot \left(\left(b \cdot b\right) \cdot a\right) + 9 \cdot \left(a \cdot \left(c \cdot a\right)\right)\right)}{12 \cdot \left(b \cdot \left(b \cdot b\right)\right) + a \cdot \left(\frac{a}{\frac{b}{c \cdot c}} \cdot 6.75 + \left(c \cdot b\right) \cdot -27\right)}}{a}
\end{array}
Initial program 33.2%
/-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-*.f6433.2%
Simplified33.2%
flip--N/A
associate-/l/N/A
Applied egg-rr34.6%
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
associate-*r*N/A
*-lowering-*.f64N/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
Simplified95.3%
Applied egg-rr95.3%
Final simplification95.3%
(FPCore (a b c)
:precision binary64
(*
(/ (+ (* -6.0 (* (* b b) a)) (* 9.0 (* a (* c a)))) a)
(/
c
(+
(* 12.0 (* b (* b b)))
(* a (+ (* (/ a (/ b (* c c))) 6.75) (* (* c b) -27.0)))))))
double code(double a, double b, double c) {
return (((-6.0 * ((b * b) * a)) + (9.0 * (a * (c * a)))) / a) * (c / ((12.0 * (b * (b * b))) + (a * (((a / (b / (c * c))) * 6.75) + ((c * b) * -27.0)))));
}
real(8) function code(a, b, c)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
code = ((((-6.0d0) * ((b * b) * a)) + (9.0d0 * (a * (c * a)))) / a) * (c / ((12.0d0 * (b * (b * b))) + (a * (((a / (b / (c * c))) * 6.75d0) + ((c * b) * (-27.0d0))))))
end function
public static double code(double a, double b, double c) {
return (((-6.0 * ((b * b) * a)) + (9.0 * (a * (c * a)))) / a) * (c / ((12.0 * (b * (b * b))) + (a * (((a / (b / (c * c))) * 6.75) + ((c * b) * -27.0)))));
}
def code(a, b, c): return (((-6.0 * ((b * b) * a)) + (9.0 * (a * (c * a)))) / a) * (c / ((12.0 * (b * (b * b))) + (a * (((a / (b / (c * c))) * 6.75) + ((c * b) * -27.0)))))
function code(a, b, c) return Float64(Float64(Float64(Float64(-6.0 * Float64(Float64(b * b) * a)) + Float64(9.0 * Float64(a * Float64(c * a)))) / a) * Float64(c / Float64(Float64(12.0 * Float64(b * Float64(b * b))) + Float64(a * Float64(Float64(Float64(a / Float64(b / Float64(c * c))) * 6.75) + Float64(Float64(c * b) * -27.0)))))) end
function tmp = code(a, b, c) tmp = (((-6.0 * ((b * b) * a)) + (9.0 * (a * (c * a)))) / a) * (c / ((12.0 * (b * (b * b))) + (a * (((a / (b / (c * c))) * 6.75) + ((c * b) * -27.0))))); end
code[a_, b_, c_] := N[(N[(N[(N[(-6.0 * N[(N[(b * b), $MachinePrecision] * a), $MachinePrecision]), $MachinePrecision] + N[(9.0 * N[(a * N[(c * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision] * N[(c / N[(N[(12.0 * N[(b * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(a * N[(N[(N[(a / N[(b / N[(c * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * 6.75), $MachinePrecision] + N[(N[(c * b), $MachinePrecision] * -27.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{-6 \cdot \left(\left(b \cdot b\right) \cdot a\right) + 9 \cdot \left(a \cdot \left(c \cdot a\right)\right)}{a} \cdot \frac{c}{12 \cdot \left(b \cdot \left(b \cdot b\right)\right) + a \cdot \left(\frac{a}{\frac{b}{c \cdot c}} \cdot 6.75 + \left(c \cdot b\right) \cdot -27\right)}
\end{array}
Initial program 33.2%
/-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-*.f6433.2%
Simplified33.2%
flip--N/A
associate-/l/N/A
Applied egg-rr34.6%
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
associate-*r*N/A
*-lowering-*.f64N/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
Simplified95.3%
Applied egg-rr95.3%
Final simplification95.3%
(FPCore (a b c)
:precision binary64
(*
c
(/
(+ (* -6.0 (* (* b b) a)) (* 9.0 (* a (* c a))))
(*
a
(+
(* 12.0 (* b (* b b)))
(* a (+ (* (/ a (/ b (* c c))) 6.75) (* (* c b) -27.0))))))))
double code(double a, double b, double c) {
return c * (((-6.0 * ((b * b) * a)) + (9.0 * (a * (c * a)))) / (a * ((12.0 * (b * (b * b))) + (a * (((a / (b / (c * c))) * 6.75) + ((c * b) * -27.0))))));
}
real(8) function code(a, b, c)
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
code = c * ((((-6.0d0) * ((b * b) * a)) + (9.0d0 * (a * (c * a)))) / (a * ((12.0d0 * (b * (b * b))) + (a * (((a / (b / (c * c))) * 6.75d0) + ((c * b) * (-27.0d0)))))))
end function
public static double code(double a, double b, double c) {
return c * (((-6.0 * ((b * b) * a)) + (9.0 * (a * (c * a)))) / (a * ((12.0 * (b * (b * b))) + (a * (((a / (b / (c * c))) * 6.75) + ((c * b) * -27.0))))));
}
def code(a, b, c): return c * (((-6.0 * ((b * b) * a)) + (9.0 * (a * (c * a)))) / (a * ((12.0 * (b * (b * b))) + (a * (((a / (b / (c * c))) * 6.75) + ((c * b) * -27.0))))))
function code(a, b, c) return Float64(c * Float64(Float64(Float64(-6.0 * Float64(Float64(b * b) * a)) + Float64(9.0 * Float64(a * Float64(c * a)))) / Float64(a * Float64(Float64(12.0 * Float64(b * Float64(b * b))) + Float64(a * Float64(Float64(Float64(a / Float64(b / Float64(c * c))) * 6.75) + Float64(Float64(c * b) * -27.0))))))) end
function tmp = code(a, b, c) tmp = c * (((-6.0 * ((b * b) * a)) + (9.0 * (a * (c * a)))) / (a * ((12.0 * (b * (b * b))) + (a * (((a / (b / (c * c))) * 6.75) + ((c * b) * -27.0)))))); end
code[a_, b_, c_] := N[(c * N[(N[(N[(-6.0 * N[(N[(b * b), $MachinePrecision] * a), $MachinePrecision]), $MachinePrecision] + N[(9.0 * N[(a * N[(c * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(a * N[(N[(12.0 * N[(b * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(a * N[(N[(N[(a / N[(b / N[(c * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * 6.75), $MachinePrecision] + N[(N[(c * b), $MachinePrecision] * -27.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
c \cdot \frac{-6 \cdot \left(\left(b \cdot b\right) \cdot a\right) + 9 \cdot \left(a \cdot \left(c \cdot a\right)\right)}{a \cdot \left(12 \cdot \left(b \cdot \left(b \cdot b\right)\right) + a \cdot \left(\frac{a}{\frac{b}{c \cdot c}} \cdot 6.75 + \left(c \cdot b\right) \cdot -27\right)\right)}
\end{array}
Initial program 33.2%
/-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-*.f6433.2%
Simplified33.2%
flip--N/A
associate-/l/N/A
Applied egg-rr34.6%
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
associate-*r*N/A
*-lowering-*.f64N/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
Simplified95.3%
Applied egg-rr95.3%
Final simplification95.3%
(FPCore (a b c)
:precision binary64
(/
1.0
(/
(+
(* -2.0 b)
(*
c
(+ (* (* c -3.0) (* (/ (* a a) (* b (* b b))) -0.375)) (* 1.5 (/ a b)))))
c)))
double code(double a, double b, double c) {
return 1.0 / (((-2.0 * b) + (c * (((c * -3.0) * (((a * a) / (b * (b * b))) * -0.375)) + (1.5 * (a / b))))) / c);
}
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 / ((((-2.0d0) * b) + (c * (((c * (-3.0d0)) * (((a * a) / (b * (b * b))) * (-0.375d0))) + (1.5d0 * (a / b))))) / c)
end function
public static double code(double a, double b, double c) {
return 1.0 / (((-2.0 * b) + (c * (((c * -3.0) * (((a * a) / (b * (b * b))) * -0.375)) + (1.5 * (a / b))))) / c);
}
def code(a, b, c): return 1.0 / (((-2.0 * b) + (c * (((c * -3.0) * (((a * a) / (b * (b * b))) * -0.375)) + (1.5 * (a / b))))) / c)
function code(a, b, c) return Float64(1.0 / Float64(Float64(Float64(-2.0 * b) + Float64(c * Float64(Float64(Float64(c * -3.0) * Float64(Float64(Float64(a * a) / Float64(b * Float64(b * b))) * -0.375)) + Float64(1.5 * Float64(a / b))))) / c)) end
function tmp = code(a, b, c) tmp = 1.0 / (((-2.0 * b) + (c * (((c * -3.0) * (((a * a) / (b * (b * b))) * -0.375)) + (1.5 * (a / b))))) / c); end
code[a_, b_, c_] := N[(1.0 / N[(N[(N[(-2.0 * b), $MachinePrecision] + N[(c * N[(N[(N[(c * -3.0), $MachinePrecision] * N[(N[(N[(a * a), $MachinePrecision] / N[(b * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * -0.375), $MachinePrecision]), $MachinePrecision] + N[(1.5 * N[(a / b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / c), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{1}{\frac{-2 \cdot b + c \cdot \left(\left(c \cdot -3\right) \cdot \left(\frac{a \cdot a}{b \cdot \left(b \cdot b\right)} \cdot -0.375\right) + 1.5 \cdot \frac{a}{b}\right)}{c}}
\end{array}
Initial program 33.2%
/-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-*.f6433.2%
Simplified33.2%
associate-/l/N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
sqrt-lowering-sqrt.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6433.2%
Applied egg-rr33.2%
Taylor expanded in c around 0
/-lowering-/.f64N/A
Simplified93.5%
Final simplification93.5%
(FPCore (a b c) :precision binary64 (/ 1.0 (+ (/ (* -2.0 b) c) (* a (+ (* (* a -3.0) (* -0.375 (/ c (* b (* b b))))) (/ 1.5 b))))))
double code(double a, double b, double c) {
return 1.0 / (((-2.0 * b) / c) + (a * (((a * -3.0) * (-0.375 * (c / (b * (b * b))))) + (1.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 = 1.0d0 / ((((-2.0d0) * b) / c) + (a * (((a * (-3.0d0)) * ((-0.375d0) * (c / (b * (b * b))))) + (1.5d0 / b))))
end function
public static double code(double a, double b, double c) {
return 1.0 / (((-2.0 * b) / c) + (a * (((a * -3.0) * (-0.375 * (c / (b * (b * b))))) + (1.5 / b))));
}
def code(a, b, c): return 1.0 / (((-2.0 * b) / c) + (a * (((a * -3.0) * (-0.375 * (c / (b * (b * b))))) + (1.5 / b))))
function code(a, b, c) return Float64(1.0 / Float64(Float64(Float64(-2.0 * b) / c) + Float64(a * Float64(Float64(Float64(a * -3.0) * Float64(-0.375 * Float64(c / Float64(b * Float64(b * b))))) + Float64(1.5 / b))))) end
function tmp = code(a, b, c) tmp = 1.0 / (((-2.0 * b) / c) + (a * (((a * -3.0) * (-0.375 * (c / (b * (b * b))))) + (1.5 / b)))); end
code[a_, b_, c_] := N[(1.0 / N[(N[(N[(-2.0 * b), $MachinePrecision] / c), $MachinePrecision] + N[(a * N[(N[(N[(a * -3.0), $MachinePrecision] * N[(-0.375 * N[(c / N[(b * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(1.5 / b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{1}{\frac{-2 \cdot b}{c} + a \cdot \left(\left(a \cdot -3\right) \cdot \left(-0.375 \cdot \frac{c}{b \cdot \left(b \cdot b\right)}\right) + \frac{1.5}{b}\right)}
\end{array}
Initial program 33.2%
/-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-*.f6433.2%
Simplified33.2%
associate-/l/N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
sqrt-lowering-sqrt.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6433.2%
Applied egg-rr33.2%
Taylor expanded in a around 0
+-lowering-+.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified93.5%
Final simplification93.5%
(FPCore (a b c) :precision binary64 (/ 1.0 (/ (+ (* -2.0 b) (/ (* (* c a) 1.5) b)) c)))
double code(double a, double b, double c) {
return 1.0 / (((-2.0 * b) + (((c * a) * 1.5) / b)) / c);
}
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 / ((((-2.0d0) * b) + (((c * a) * 1.5d0) / b)) / c)
end function
public static double code(double a, double b, double c) {
return 1.0 / (((-2.0 * b) + (((c * a) * 1.5) / b)) / c);
}
def code(a, b, c): return 1.0 / (((-2.0 * b) + (((c * a) * 1.5) / b)) / c)
function code(a, b, c) return Float64(1.0 / Float64(Float64(Float64(-2.0 * b) + Float64(Float64(Float64(c * a) * 1.5) / b)) / c)) end
function tmp = code(a, b, c) tmp = 1.0 / (((-2.0 * b) + (((c * a) * 1.5) / b)) / c); end
code[a_, b_, c_] := N[(1.0 / N[(N[(N[(-2.0 * b), $MachinePrecision] + N[(N[(N[(c * a), $MachinePrecision] * 1.5), $MachinePrecision] / b), $MachinePrecision]), $MachinePrecision] / c), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{1}{\frac{-2 \cdot b + \frac{\left(c \cdot a\right) \cdot 1.5}{b}}{c}}
\end{array}
Initial program 33.2%
/-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-*.f6433.2%
Simplified33.2%
associate-/l/N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
sqrt-lowering-sqrt.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6433.2%
Applied egg-rr33.2%
Taylor expanded in c around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6489.8%
Simplified89.8%
Final simplification89.8%
(FPCore (a b c) :precision binary64 (/ 1.0 (+ (* 1.5 (/ a b)) (/ (* -2.0 b) c))))
double code(double a, double b, double c) {
return 1.0 / ((1.5 * (a / b)) + ((-2.0 * b) / c));
}
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 / ((1.5d0 * (a / b)) + (((-2.0d0) * b) / c))
end function
public static double code(double a, double b, double c) {
return 1.0 / ((1.5 * (a / b)) + ((-2.0 * b) / c));
}
def code(a, b, c): return 1.0 / ((1.5 * (a / b)) + ((-2.0 * b) / c))
function code(a, b, c) return Float64(1.0 / Float64(Float64(1.5 * Float64(a / b)) + Float64(Float64(-2.0 * b) / c))) end
function tmp = code(a, b, c) tmp = 1.0 / ((1.5 * (a / b)) + ((-2.0 * b) / c)); end
code[a_, b_, c_] := N[(1.0 / N[(N[(1.5 * N[(a / b), $MachinePrecision]), $MachinePrecision] + N[(N[(-2.0 * b), $MachinePrecision] / c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{1}{1.5 \cdot \frac{a}{b} + \frac{-2 \cdot b}{c}}
\end{array}
Initial program 33.2%
/-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-*.f6433.2%
Simplified33.2%
associate-/l/N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
sqrt-lowering-sqrt.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6433.2%
Applied egg-rr33.2%
Taylor expanded in a around 0
+-lowering-+.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6489.8%
Simplified89.8%
Final simplification89.8%
(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 33.2%
/-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-*.f6433.2%
Simplified33.2%
Taylor expanded in b around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6479.8%
Simplified79.8%
(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 33.2%
/-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-*.f6433.2%
Simplified33.2%
Taylor expanded in b around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6479.8%
Simplified79.8%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f6479.6%
Applied egg-rr79.6%
Final simplification79.6%
herbie shell --seed 2024163
(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)))