
(FPCore (x y z t) :precision binary64 (* (/ 1.0 3.0) (acos (* (/ (* 3.0 (/ x (* y 27.0))) (* z 2.0)) (sqrt t)))))
double code(double x, double y, double z, double t) {
return (1.0 / 3.0) * acos((((3.0 * (x / (y * 27.0))) / (z * 2.0)) * sqrt(t)));
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = (1.0d0 / 3.0d0) * acos((((3.0d0 * (x / (y * 27.0d0))) / (z * 2.0d0)) * sqrt(t)))
end function
public static double code(double x, double y, double z, double t) {
return (1.0 / 3.0) * Math.acos((((3.0 * (x / (y * 27.0))) / (z * 2.0)) * Math.sqrt(t)));
}
def code(x, y, z, t): return (1.0 / 3.0) * math.acos((((3.0 * (x / (y * 27.0))) / (z * 2.0)) * math.sqrt(t)))
function code(x, y, z, t) return Float64(Float64(1.0 / 3.0) * acos(Float64(Float64(Float64(3.0 * Float64(x / Float64(y * 27.0))) / Float64(z * 2.0)) * sqrt(t)))) end
function tmp = code(x, y, z, t) tmp = (1.0 / 3.0) * acos((((3.0 * (x / (y * 27.0))) / (z * 2.0)) * sqrt(t))); end
code[x_, y_, z_, t_] := N[(N[(1.0 / 3.0), $MachinePrecision] * N[ArcCos[N[(N[(N[(3.0 * N[(x / N[(y * 27.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(z * 2.0), $MachinePrecision]), $MachinePrecision] * N[Sqrt[t], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{1}{3} \cdot \cos^{-1} \left(\frac{3 \cdot \frac{x}{y \cdot 27}}{z \cdot 2} \cdot \sqrt{t}\right)
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 8 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t) :precision binary64 (* (/ 1.0 3.0) (acos (* (/ (* 3.0 (/ x (* y 27.0))) (* z 2.0)) (sqrt t)))))
double code(double x, double y, double z, double t) {
return (1.0 / 3.0) * acos((((3.0 * (x / (y * 27.0))) / (z * 2.0)) * sqrt(t)));
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = (1.0d0 / 3.0d0) * acos((((3.0d0 * (x / (y * 27.0d0))) / (z * 2.0d0)) * sqrt(t)))
end function
public static double code(double x, double y, double z, double t) {
return (1.0 / 3.0) * Math.acos((((3.0 * (x / (y * 27.0))) / (z * 2.0)) * Math.sqrt(t)));
}
def code(x, y, z, t): return (1.0 / 3.0) * math.acos((((3.0 * (x / (y * 27.0))) / (z * 2.0)) * math.sqrt(t)))
function code(x, y, z, t) return Float64(Float64(1.0 / 3.0) * acos(Float64(Float64(Float64(3.0 * Float64(x / Float64(y * 27.0))) / Float64(z * 2.0)) * sqrt(t)))) end
function tmp = code(x, y, z, t) tmp = (1.0 / 3.0) * acos((((3.0 * (x / (y * 27.0))) / (z * 2.0)) * sqrt(t))); end
code[x_, y_, z_, t_] := N[(N[(1.0 / 3.0), $MachinePrecision] * N[ArcCos[N[(N[(N[(3.0 * N[(x / N[(y * 27.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(z * 2.0), $MachinePrecision]), $MachinePrecision] * N[Sqrt[t], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{1}{3} \cdot \cos^{-1} \left(\frac{3 \cdot \frac{x}{y \cdot 27}}{z \cdot 2} \cdot \sqrt{t}\right)
\end{array}
(FPCore (x y z t) :precision binary64 (+ -1.0 (fma 0.3333333333333333 (acos (* 0.05555555555555555 (* (sqrt t) (/ (/ x y) z)))) 1.0)))
double code(double x, double y, double z, double t) {
return -1.0 + fma(0.3333333333333333, acos((0.05555555555555555 * (sqrt(t) * ((x / y) / z)))), 1.0);
}
function code(x, y, z, t) return Float64(-1.0 + fma(0.3333333333333333, acos(Float64(0.05555555555555555 * Float64(sqrt(t) * Float64(Float64(x / y) / z)))), 1.0)) end
code[x_, y_, z_, t_] := N[(-1.0 + N[(0.3333333333333333 * N[ArcCos[N[(0.05555555555555555 * N[(N[Sqrt[t], $MachinePrecision] * N[(N[(x / y), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
-1 + \mathsf{fma}\left(0.3333333333333333, \cos^{-1} \left(0.05555555555555555 \cdot \left(\sqrt{t} \cdot \frac{\frac{x}{y}}{z}\right)\right), 1\right)
\end{array}
Initial program 98.5%
Simplified98.5%
*-commutative98.5%
associate-/l/98.1%
associate-*l/98.1%
*-commutative98.1%
Applied egg-rr98.1%
associate-*l/98.1%
associate-*r*98.1%
expm1-log1p-u98.1%
expm1-undefine99.6%
Applied egg-rr99.6%
sub-neg99.6%
metadata-eval99.6%
+-commutative99.6%
log1p-undefine97.2%
rem-exp-log97.2%
+-commutative97.2%
fma-define99.6%
Simplified99.6%
Taylor expanded in t around -inf 0.0%
mul-1-neg0.0%
distribute-rgt-neg-in0.0%
*-commutative0.0%
unpow20.0%
rem-square-sqrt99.6%
neg-mul-199.6%
distribute-frac-neg299.6%
distribute-lft-neg-out99.6%
associate-/r*100.0%
neg-mul-1100.0%
neg-mul-1100.0%
times-frac100.0%
metadata-eval100.0%
Simplified100.0%
Final simplification100.0%
(FPCore (x y z t) :precision binary64 (+ -1.0 (fma 0.3333333333333333 (acos (* 0.05555555555555555 (* x (/ (sqrt t) (* y z))))) 1.0)))
double code(double x, double y, double z, double t) {
return -1.0 + fma(0.3333333333333333, acos((0.05555555555555555 * (x * (sqrt(t) / (y * z))))), 1.0);
}
function code(x, y, z, t) return Float64(-1.0 + fma(0.3333333333333333, acos(Float64(0.05555555555555555 * Float64(x * Float64(sqrt(t) / Float64(y * z))))), 1.0)) end
code[x_, y_, z_, t_] := N[(-1.0 + N[(0.3333333333333333 * N[ArcCos[N[(0.05555555555555555 * N[(x * N[(N[Sqrt[t], $MachinePrecision] / N[(y * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
-1 + \mathsf{fma}\left(0.3333333333333333, \cos^{-1} \left(0.05555555555555555 \cdot \left(x \cdot \frac{\sqrt{t}}{y \cdot z}\right)\right), 1\right)
\end{array}
Initial program 98.5%
Simplified98.5%
*-commutative98.5%
associate-/l/98.1%
associate-*l/98.1%
*-commutative98.1%
Applied egg-rr98.1%
associate-*l/98.1%
associate-*r*98.1%
expm1-log1p-u98.1%
expm1-undefine99.6%
Applied egg-rr99.6%
sub-neg99.6%
metadata-eval99.6%
+-commutative99.6%
log1p-undefine97.2%
rem-exp-log97.2%
+-commutative97.2%
fma-define99.6%
Simplified99.6%
Final simplification99.6%
(FPCore (x y z t) :precision binary64 (* 0.3333333333333333 (cbrt (pow (acos 0.0) 3.0))))
double code(double x, double y, double z, double t) {
return 0.3333333333333333 * cbrt(pow(acos(0.0), 3.0));
}
public static double code(double x, double y, double z, double t) {
return 0.3333333333333333 * Math.cbrt(Math.pow(Math.acos(0.0), 3.0));
}
function code(x, y, z, t) return Float64(0.3333333333333333 * cbrt((acos(0.0) ^ 3.0))) end
code[x_, y_, z_, t_] := N[(0.3333333333333333 * N[Power[N[Power[N[ArcCos[0.0], $MachinePrecision], 3.0], $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
0.3333333333333333 \cdot \sqrt[3]{{\cos^{-1} 0}^{3}}
\end{array}
Initial program 98.5%
Simplified98.5%
*-commutative98.5%
associate-/l/98.1%
associate-*l/98.1%
*-commutative98.1%
Applied egg-rr98.1%
associate-*l/98.1%
expm1-log1p-u98.1%
associate-*r*98.1%
expm1-undefine98.1%
*-commutative98.1%
*-commutative98.1%
associate-*r*98.1%
associate-*l/98.1%
associate-/l*98.1%
associate-*l*98.1%
*-commutative98.1%
associate-/r*98.1%
Applied egg-rr98.1%
sub-neg98.1%
metadata-eval98.1%
+-commutative98.1%
log1p-undefine98.1%
rem-exp-log98.1%
+-commutative98.1%
associate-*l/98.1%
associate-*r/98.1%
associate-*l/98.1%
associate-/l*97.3%
*-commutative97.3%
associate-/l*98.1%
associate-*r/98.1%
*-commutative98.1%
associate-*r*98.1%
Simplified98.1%
Taylor expanded in x around 0 97.0%
add-cbrt-cube98.5%
pow398.5%
metadata-eval98.5%
Applied egg-rr98.5%
(FPCore (x y z t)
:precision binary64
(*
0.3333333333333333
(acos
(+
-1.0
(*
t
(+
(/ 1.0 t)
(* 0.05555555555555555 (* (sqrt (/ 1.0 t)) (/ x (* y z))))))))))
double code(double x, double y, double z, double t) {
return 0.3333333333333333 * acos((-1.0 + (t * ((1.0 / t) + (0.05555555555555555 * (sqrt((1.0 / t)) * (x / (y * z))))))));
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = 0.3333333333333333d0 * acos(((-1.0d0) + (t * ((1.0d0 / t) + (0.05555555555555555d0 * (sqrt((1.0d0 / t)) * (x / (y * z))))))))
end function
public static double code(double x, double y, double z, double t) {
return 0.3333333333333333 * Math.acos((-1.0 + (t * ((1.0 / t) + (0.05555555555555555 * (Math.sqrt((1.0 / t)) * (x / (y * z))))))));
}
def code(x, y, z, t): return 0.3333333333333333 * math.acos((-1.0 + (t * ((1.0 / t) + (0.05555555555555555 * (math.sqrt((1.0 / t)) * (x / (y * z))))))))
function code(x, y, z, t) return Float64(0.3333333333333333 * acos(Float64(-1.0 + Float64(t * Float64(Float64(1.0 / t) + Float64(0.05555555555555555 * Float64(sqrt(Float64(1.0 / t)) * Float64(x / Float64(y * z))))))))) end
function tmp = code(x, y, z, t) tmp = 0.3333333333333333 * acos((-1.0 + (t * ((1.0 / t) + (0.05555555555555555 * (sqrt((1.0 / t)) * (x / (y * z)))))))); end
code[x_, y_, z_, t_] := N[(0.3333333333333333 * N[ArcCos[N[(-1.0 + N[(t * N[(N[(1.0 / t), $MachinePrecision] + N[(0.05555555555555555 * N[(N[Sqrt[N[(1.0 / t), $MachinePrecision]], $MachinePrecision] * N[(x / N[(y * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
0.3333333333333333 \cdot \cos^{-1} \left(-1 + t \cdot \left(\frac{1}{t} + 0.05555555555555555 \cdot \left(\sqrt{\frac{1}{t}} \cdot \frac{x}{y \cdot z}\right)\right)\right)
\end{array}
Initial program 98.5%
Simplified98.5%
*-commutative98.5%
associate-/l/98.1%
associate-*l/98.1%
*-commutative98.1%
Applied egg-rr98.1%
associate-*l/98.1%
expm1-log1p-u98.1%
associate-*r*98.1%
expm1-undefine98.1%
*-commutative98.1%
*-commutative98.1%
associate-*r*98.1%
associate-*l/98.1%
associate-/l*98.1%
associate-*l*98.1%
*-commutative98.1%
associate-/r*98.1%
Applied egg-rr98.1%
sub-neg98.1%
metadata-eval98.1%
+-commutative98.1%
log1p-undefine98.1%
rem-exp-log98.1%
+-commutative98.1%
associate-*l/98.1%
associate-*r/98.1%
associate-*l/98.1%
associate-/l*97.3%
*-commutative97.3%
associate-/l*98.1%
associate-*r/98.1%
*-commutative98.1%
associate-*r*98.1%
Simplified98.1%
Taylor expanded in t around inf 98.3%
Final simplification98.3%
(FPCore (x y z t) :precision binary64 (* 0.3333333333333333 (acos (+ -1.0 (* x (+ (* 0.05555555555555555 (/ (sqrt t) (* y z))) (/ 1.0 x)))))))
double code(double x, double y, double z, double t) {
return 0.3333333333333333 * acos((-1.0 + (x * ((0.05555555555555555 * (sqrt(t) / (y * z))) + (1.0 / x)))));
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = 0.3333333333333333d0 * acos(((-1.0d0) + (x * ((0.05555555555555555d0 * (sqrt(t) / (y * z))) + (1.0d0 / x)))))
end function
public static double code(double x, double y, double z, double t) {
return 0.3333333333333333 * Math.acos((-1.0 + (x * ((0.05555555555555555 * (Math.sqrt(t) / (y * z))) + (1.0 / x)))));
}
def code(x, y, z, t): return 0.3333333333333333 * math.acos((-1.0 + (x * ((0.05555555555555555 * (math.sqrt(t) / (y * z))) + (1.0 / x)))))
function code(x, y, z, t) return Float64(0.3333333333333333 * acos(Float64(-1.0 + Float64(x * Float64(Float64(0.05555555555555555 * Float64(sqrt(t) / Float64(y * z))) + Float64(1.0 / x)))))) end
function tmp = code(x, y, z, t) tmp = 0.3333333333333333 * acos((-1.0 + (x * ((0.05555555555555555 * (sqrt(t) / (y * z))) + (1.0 / x))))); end
code[x_, y_, z_, t_] := N[(0.3333333333333333 * N[ArcCos[N[(-1.0 + N[(x * N[(N[(0.05555555555555555 * N[(N[Sqrt[t], $MachinePrecision] / N[(y * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(1.0 / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
0.3333333333333333 \cdot \cos^{-1} \left(-1 + x \cdot \left(0.05555555555555555 \cdot \frac{\sqrt{t}}{y \cdot z} + \frac{1}{x}\right)\right)
\end{array}
Initial program 98.5%
Simplified98.5%
*-commutative98.5%
associate-/l/98.1%
associate-*l/98.1%
*-commutative98.1%
Applied egg-rr98.1%
associate-*l/98.1%
expm1-log1p-u98.1%
associate-*r*98.1%
expm1-undefine98.1%
*-commutative98.1%
*-commutative98.1%
associate-*r*98.1%
associate-*l/98.1%
associate-/l*98.1%
associate-*l*98.1%
*-commutative98.1%
associate-/r*98.1%
Applied egg-rr98.1%
sub-neg98.1%
metadata-eval98.1%
+-commutative98.1%
log1p-undefine98.1%
rem-exp-log98.1%
+-commutative98.1%
associate-*l/98.1%
associate-*r/98.1%
associate-*l/98.1%
associate-/l*97.3%
*-commutative97.3%
associate-/l*98.1%
associate-*r/98.1%
*-commutative98.1%
associate-*r*98.1%
Simplified98.1%
Taylor expanded in x around inf 98.3%
Taylor expanded in t around -inf 0.0%
associate-*r/0.0%
associate-*r/0.0%
associate-*r*0.0%
unpow20.0%
*-commutative0.0%
rem-square-sqrt98.3%
associate-*l*98.3%
metadata-eval98.3%
*-rgt-identity98.3%
Simplified98.3%
(FPCore (x y z t) :precision binary64 (* 0.3333333333333333 (acos (+ -1.0 (/ (+ z (* (/ x y) (* 0.05555555555555555 (sqrt t)))) z)))))
double code(double x, double y, double z, double t) {
return 0.3333333333333333 * acos((-1.0 + ((z + ((x / y) * (0.05555555555555555 * sqrt(t)))) / z)));
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = 0.3333333333333333d0 * acos(((-1.0d0) + ((z + ((x / y) * (0.05555555555555555d0 * sqrt(t)))) / z)))
end function
public static double code(double x, double y, double z, double t) {
return 0.3333333333333333 * Math.acos((-1.0 + ((z + ((x / y) * (0.05555555555555555 * Math.sqrt(t)))) / z)));
}
def code(x, y, z, t): return 0.3333333333333333 * math.acos((-1.0 + ((z + ((x / y) * (0.05555555555555555 * math.sqrt(t)))) / z)))
function code(x, y, z, t) return Float64(0.3333333333333333 * acos(Float64(-1.0 + Float64(Float64(z + Float64(Float64(x / y) * Float64(0.05555555555555555 * sqrt(t)))) / z)))) end
function tmp = code(x, y, z, t) tmp = 0.3333333333333333 * acos((-1.0 + ((z + ((x / y) * (0.05555555555555555 * sqrt(t)))) / z))); end
code[x_, y_, z_, t_] := N[(0.3333333333333333 * N[ArcCos[N[(-1.0 + N[(N[(z + N[(N[(x / y), $MachinePrecision] * N[(0.05555555555555555 * N[Sqrt[t], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
0.3333333333333333 \cdot \cos^{-1} \left(-1 + \frac{z + \frac{x}{y} \cdot \left(0.05555555555555555 \cdot \sqrt{t}\right)}{z}\right)
\end{array}
Initial program 98.5%
Simplified98.5%
*-commutative98.5%
associate-/l/98.1%
associate-*l/98.1%
*-commutative98.1%
Applied egg-rr98.1%
associate-*l/98.1%
expm1-log1p-u98.1%
associate-*r*98.1%
expm1-undefine98.1%
*-commutative98.1%
*-commutative98.1%
associate-*r*98.1%
associate-*l/98.1%
associate-/l*98.1%
associate-*l*98.1%
*-commutative98.1%
associate-/r*98.1%
Applied egg-rr98.1%
sub-neg98.1%
metadata-eval98.1%
+-commutative98.1%
log1p-undefine98.1%
rem-exp-log98.1%
+-commutative98.1%
associate-*l/98.1%
associate-*r/98.1%
associate-*l/98.1%
associate-/l*97.3%
*-commutative97.3%
associate-/l*98.1%
associate-*r/98.1%
*-commutative98.1%
associate-*r*98.1%
Simplified98.1%
Taylor expanded in z around 0 98.5%
associate-*r*98.5%
Simplified98.5%
Final simplification98.5%
(FPCore (x y z t) :precision binary64 (* 0.3333333333333333 (acos (* (sqrt t) (* 0.05555555555555555 (/ (/ x y) z))))))
double code(double x, double y, double z, double t) {
return 0.3333333333333333 * acos((sqrt(t) * (0.05555555555555555 * ((x / y) / z))));
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = 0.3333333333333333d0 * acos((sqrt(t) * (0.05555555555555555d0 * ((x / y) / z))))
end function
public static double code(double x, double y, double z, double t) {
return 0.3333333333333333 * Math.acos((Math.sqrt(t) * (0.05555555555555555 * ((x / y) / z))));
}
def code(x, y, z, t): return 0.3333333333333333 * math.acos((math.sqrt(t) * (0.05555555555555555 * ((x / y) / z))))
function code(x, y, z, t) return Float64(0.3333333333333333 * acos(Float64(sqrt(t) * Float64(0.05555555555555555 * Float64(Float64(x / y) / z))))) end
function tmp = code(x, y, z, t) tmp = 0.3333333333333333 * acos((sqrt(t) * (0.05555555555555555 * ((x / y) / z)))); end
code[x_, y_, z_, t_] := N[(0.3333333333333333 * N[ArcCos[N[(N[Sqrt[t], $MachinePrecision] * N[(0.05555555555555555 * N[(N[(x / y), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
0.3333333333333333 \cdot \cos^{-1} \left(\sqrt{t} \cdot \left(0.05555555555555555 \cdot \frac{\frac{x}{y}}{z}\right)\right)
\end{array}
Initial program 98.5%
Simplified98.5%
Final simplification98.5%
(FPCore (x y z t) :precision binary64 (* 0.3333333333333333 (acos 0.0)))
double code(double x, double y, double z, double t) {
return 0.3333333333333333 * acos(0.0);
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = 0.3333333333333333d0 * acos(0.0d0)
end function
public static double code(double x, double y, double z, double t) {
return 0.3333333333333333 * Math.acos(0.0);
}
def code(x, y, z, t): return 0.3333333333333333 * math.acos(0.0)
function code(x, y, z, t) return Float64(0.3333333333333333 * acos(0.0)) end
function tmp = code(x, y, z, t) tmp = 0.3333333333333333 * acos(0.0); end
code[x_, y_, z_, t_] := N[(0.3333333333333333 * N[ArcCos[0.0], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
0.3333333333333333 \cdot \cos^{-1} 0
\end{array}
Initial program 98.5%
Simplified98.5%
*-commutative98.5%
associate-/l/98.1%
associate-*l/98.1%
*-commutative98.1%
Applied egg-rr98.1%
associate-*l/98.1%
expm1-log1p-u98.1%
associate-*r*98.1%
expm1-undefine98.1%
*-commutative98.1%
*-commutative98.1%
associate-*r*98.1%
associate-*l/98.1%
associate-/l*98.1%
associate-*l*98.1%
*-commutative98.1%
associate-/r*98.1%
Applied egg-rr98.1%
sub-neg98.1%
metadata-eval98.1%
+-commutative98.1%
log1p-undefine98.1%
rem-exp-log98.1%
+-commutative98.1%
associate-*l/98.1%
associate-*r/98.1%
associate-*l/98.1%
associate-/l*97.3%
*-commutative97.3%
associate-/l*98.1%
associate-*r/98.1%
*-commutative98.1%
associate-*r*98.1%
Simplified98.1%
Taylor expanded in x around 0 97.0%
pow197.0%
metadata-eval97.0%
Applied egg-rr97.0%
unpow197.0%
Simplified97.0%
(FPCore (x y z t) :precision binary64 (/ (acos (* (/ (/ x 27.0) (* y z)) (/ (sqrt t) (/ 2.0 3.0)))) 3.0))
double code(double x, double y, double z, double t) {
return acos((((x / 27.0) / (y * z)) * (sqrt(t) / (2.0 / 3.0)))) / 3.0;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = acos((((x / 27.0d0) / (y * z)) * (sqrt(t) / (2.0d0 / 3.0d0)))) / 3.0d0
end function
public static double code(double x, double y, double z, double t) {
return Math.acos((((x / 27.0) / (y * z)) * (Math.sqrt(t) / (2.0 / 3.0)))) / 3.0;
}
def code(x, y, z, t): return math.acos((((x / 27.0) / (y * z)) * (math.sqrt(t) / (2.0 / 3.0)))) / 3.0
function code(x, y, z, t) return Float64(acos(Float64(Float64(Float64(x / 27.0) / Float64(y * z)) * Float64(sqrt(t) / Float64(2.0 / 3.0)))) / 3.0) end
function tmp = code(x, y, z, t) tmp = acos((((x / 27.0) / (y * z)) * (sqrt(t) / (2.0 / 3.0)))) / 3.0; end
code[x_, y_, z_, t_] := N[(N[ArcCos[N[(N[(N[(x / 27.0), $MachinePrecision] / N[(y * z), $MachinePrecision]), $MachinePrecision] * N[(N[Sqrt[t], $MachinePrecision] / N[(2.0 / 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / 3.0), $MachinePrecision]
\begin{array}{l}
\\
\frac{\cos^{-1} \left(\frac{\frac{x}{27}}{y \cdot z} \cdot \frac{\sqrt{t}}{\frac{2}{3}}\right)}{3}
\end{array}
herbie shell --seed 2024123
(FPCore (x y z t)
:name "Diagrams.Solve.Polynomial:cubForm from diagrams-solve-0.1, D"
:precision binary64
:alt
(! :herbie-platform default (/ (acos (* (/ (/ x 27) (* y z)) (/ (sqrt t) (/ 2 3)))) 3))
(* (/ 1.0 3.0) (acos (* (/ (* 3.0 (/ x (* y 27.0))) (* z 2.0)) (sqrt t)))))