
(FPCore (v H) :precision binary64 (atan (/ v (sqrt (- (* v v) (* (* 2.0 9.8) H))))))
double code(double v, double H) {
return atan((v / sqrt(((v * v) - ((2.0 * 9.8) * H)))));
}
real(8) function code(v, h)
real(8), intent (in) :: v
real(8), intent (in) :: h
code = atan((v / sqrt(((v * v) - ((2.0d0 * 9.8d0) * h)))))
end function
public static double code(double v, double H) {
return Math.atan((v / Math.sqrt(((v * v) - ((2.0 * 9.8) * H)))));
}
def code(v, H): return math.atan((v / math.sqrt(((v * v) - ((2.0 * 9.8) * H)))))
function code(v, H) return atan(Float64(v / sqrt(Float64(Float64(v * v) - Float64(Float64(2.0 * 9.8) * H))))) end
function tmp = code(v, H) tmp = atan((v / sqrt(((v * v) - ((2.0 * 9.8) * H))))); end
code[v_, H_] := N[ArcTan[N[(v / N[Sqrt[N[(N[(v * v), $MachinePrecision] - N[(N[(2.0 * 9.8), $MachinePrecision] * H), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1} \left(\frac{v}{\sqrt{v \cdot v - \left(2 \cdot 9.8\right) \cdot H}}\right)
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 11 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (v H) :precision binary64 (atan (/ v (sqrt (- (* v v) (* (* 2.0 9.8) H))))))
double code(double v, double H) {
return atan((v / sqrt(((v * v) - ((2.0 * 9.8) * H)))));
}
real(8) function code(v, h)
real(8), intent (in) :: v
real(8), intent (in) :: h
code = atan((v / sqrt(((v * v) - ((2.0d0 * 9.8d0) * h)))))
end function
public static double code(double v, double H) {
return Math.atan((v / Math.sqrt(((v * v) - ((2.0 * 9.8) * H)))));
}
def code(v, H): return math.atan((v / math.sqrt(((v * v) - ((2.0 * 9.8) * H)))))
function code(v, H) return atan(Float64(v / sqrt(Float64(Float64(v * v) - Float64(Float64(2.0 * 9.8) * H))))) end
function tmp = code(v, H) tmp = atan((v / sqrt(((v * v) - ((2.0 * 9.8) * H))))); end
code[v_, H_] := N[ArcTan[N[(v / N[Sqrt[N[(N[(v * v), $MachinePrecision] - N[(N[(2.0 * 9.8), $MachinePrecision] * H), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1} \left(\frac{v}{\sqrt{v \cdot v - \left(2 \cdot 9.8\right) \cdot H}}\right)
\end{array}
(FPCore (v H) :precision binary64 (if (<= v -4e+114) (- 0.0 (atan (/ (/ v (+ 1.0 (/ -9.8 (/ (* v v) H)))) v))) (if (<= v 2e+112) (atan (/ v (sqrt (+ (* v v) (* H -19.6))))) (atan 1.0))))
double code(double v, double H) {
double tmp;
if (v <= -4e+114) {
tmp = 0.0 - atan(((v / (1.0 + (-9.8 / ((v * v) / H)))) / v));
} else if (v <= 2e+112) {
tmp = atan((v / sqrt(((v * v) + (H * -19.6)))));
} else {
tmp = atan(1.0);
}
return tmp;
}
real(8) function code(v, h)
real(8), intent (in) :: v
real(8), intent (in) :: h
real(8) :: tmp
if (v <= (-4d+114)) then
tmp = 0.0d0 - atan(((v / (1.0d0 + ((-9.8d0) / ((v * v) / h)))) / v))
else if (v <= 2d+112) then
tmp = atan((v / sqrt(((v * v) + (h * (-19.6d0))))))
else
tmp = atan(1.0d0)
end if
code = tmp
end function
public static double code(double v, double H) {
double tmp;
if (v <= -4e+114) {
tmp = 0.0 - Math.atan(((v / (1.0 + (-9.8 / ((v * v) / H)))) / v));
} else if (v <= 2e+112) {
tmp = Math.atan((v / Math.sqrt(((v * v) + (H * -19.6)))));
} else {
tmp = Math.atan(1.0);
}
return tmp;
}
def code(v, H): tmp = 0 if v <= -4e+114: tmp = 0.0 - math.atan(((v / (1.0 + (-9.8 / ((v * v) / H)))) / v)) elif v <= 2e+112: tmp = math.atan((v / math.sqrt(((v * v) + (H * -19.6))))) else: tmp = math.atan(1.0) return tmp
function code(v, H) tmp = 0.0 if (v <= -4e+114) tmp = Float64(0.0 - atan(Float64(Float64(v / Float64(1.0 + Float64(-9.8 / Float64(Float64(v * v) / H)))) / v))); elseif (v <= 2e+112) tmp = atan(Float64(v / sqrt(Float64(Float64(v * v) + Float64(H * -19.6))))); else tmp = atan(1.0); end return tmp end
function tmp_2 = code(v, H) tmp = 0.0; if (v <= -4e+114) tmp = 0.0 - atan(((v / (1.0 + (-9.8 / ((v * v) / H)))) / v)); elseif (v <= 2e+112) tmp = atan((v / sqrt(((v * v) + (H * -19.6))))); else tmp = atan(1.0); end tmp_2 = tmp; end
code[v_, H_] := If[LessEqual[v, -4e+114], N[(0.0 - N[ArcTan[N[(N[(v / N[(1.0 + N[(-9.8 / N[(N[(v * v), $MachinePrecision] / H), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / v), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[v, 2e+112], N[ArcTan[N[(v / N[Sqrt[N[(N[(v * v), $MachinePrecision] + N[(H * -19.6), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[1.0], $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;v \leq -4 \cdot 10^{+114}:\\
\;\;\;\;0 - \tan^{-1} \left(\frac{\frac{v}{1 + \frac{-9.8}{\frac{v \cdot v}{H}}}}{v}\right)\\
\mathbf{elif}\;v \leq 2 \cdot 10^{+112}:\\
\;\;\;\;\tan^{-1} \left(\frac{v}{\sqrt{v \cdot v + H \cdot -19.6}}\right)\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1} 1\\
\end{array}
\end{array}
if v < -4e114Initial program 21.8%
atan-lowering-atan.f64N/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
metadata-evalN/A
metadata-eval21.8%
Simplified21.8%
Taylor expanded in v around -inf
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
mul-1-negN/A
neg-sub0N/A
--lowering--.f6498.4%
Simplified98.4%
sub0-negN/A
associate-/r*N/A
distribute-frac-neg2N/A
atan-negN/A
neg-lowering-neg.f64N/A
atan-lowering-atan.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6498.4%
Applied egg-rr98.4%
if -4e114 < v < 1.9999999999999999e112Initial program 99.1%
atan-lowering-atan.f64N/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
metadata-evalN/A
metadata-eval99.1%
Simplified99.1%
if 1.9999999999999999e112 < v Initial program 17.3%
atan-lowering-atan.f64N/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
metadata-evalN/A
metadata-eval17.3%
Simplified17.3%
Taylor expanded in v around inf
Simplified100.0%
Final simplification99.1%
(FPCore (v H)
:precision binary64
(if (<= v -3.8e-16)
(- 0.0 (atan (/ (/ v (+ 1.0 (/ -9.8 (/ (* v v) H)))) v)))
(if (<= v 6.8e-29)
(atan (/ v (sqrt (* H -19.6))))
(atan (+ 1.0 (* H (/ (+ 9.8 (* (/ H (* v v)) 144.06)) (* v v))))))))
double code(double v, double H) {
double tmp;
if (v <= -3.8e-16) {
tmp = 0.0 - atan(((v / (1.0 + (-9.8 / ((v * v) / H)))) / v));
} else if (v <= 6.8e-29) {
tmp = atan((v / sqrt((H * -19.6))));
} else {
tmp = atan((1.0 + (H * ((9.8 + ((H / (v * v)) * 144.06)) / (v * v)))));
}
return tmp;
}
real(8) function code(v, h)
real(8), intent (in) :: v
real(8), intent (in) :: h
real(8) :: tmp
if (v <= (-3.8d-16)) then
tmp = 0.0d0 - atan(((v / (1.0d0 + ((-9.8d0) / ((v * v) / h)))) / v))
else if (v <= 6.8d-29) then
tmp = atan((v / sqrt((h * (-19.6d0)))))
else
tmp = atan((1.0d0 + (h * ((9.8d0 + ((h / (v * v)) * 144.06d0)) / (v * v)))))
end if
code = tmp
end function
public static double code(double v, double H) {
double tmp;
if (v <= -3.8e-16) {
tmp = 0.0 - Math.atan(((v / (1.0 + (-9.8 / ((v * v) / H)))) / v));
} else if (v <= 6.8e-29) {
tmp = Math.atan((v / Math.sqrt((H * -19.6))));
} else {
tmp = Math.atan((1.0 + (H * ((9.8 + ((H / (v * v)) * 144.06)) / (v * v)))));
}
return tmp;
}
def code(v, H): tmp = 0 if v <= -3.8e-16: tmp = 0.0 - math.atan(((v / (1.0 + (-9.8 / ((v * v) / H)))) / v)) elif v <= 6.8e-29: tmp = math.atan((v / math.sqrt((H * -19.6)))) else: tmp = math.atan((1.0 + (H * ((9.8 + ((H / (v * v)) * 144.06)) / (v * v))))) return tmp
function code(v, H) tmp = 0.0 if (v <= -3.8e-16) tmp = Float64(0.0 - atan(Float64(Float64(v / Float64(1.0 + Float64(-9.8 / Float64(Float64(v * v) / H)))) / v))); elseif (v <= 6.8e-29) tmp = atan(Float64(v / sqrt(Float64(H * -19.6)))); else tmp = atan(Float64(1.0 + Float64(H * Float64(Float64(9.8 + Float64(Float64(H / Float64(v * v)) * 144.06)) / Float64(v * v))))); end return tmp end
function tmp_2 = code(v, H) tmp = 0.0; if (v <= -3.8e-16) tmp = 0.0 - atan(((v / (1.0 + (-9.8 / ((v * v) / H)))) / v)); elseif (v <= 6.8e-29) tmp = atan((v / sqrt((H * -19.6)))); else tmp = atan((1.0 + (H * ((9.8 + ((H / (v * v)) * 144.06)) / (v * v))))); end tmp_2 = tmp; end
code[v_, H_] := If[LessEqual[v, -3.8e-16], N[(0.0 - N[ArcTan[N[(N[(v / N[(1.0 + N[(-9.8 / N[(N[(v * v), $MachinePrecision] / H), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / v), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[v, 6.8e-29], N[ArcTan[N[(v / N[Sqrt[N[(H * -19.6), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[N[(1.0 + N[(H * N[(N[(9.8 + N[(N[(H / N[(v * v), $MachinePrecision]), $MachinePrecision] * 144.06), $MachinePrecision]), $MachinePrecision] / N[(v * v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;v \leq -3.8 \cdot 10^{-16}:\\
\;\;\;\;0 - \tan^{-1} \left(\frac{\frac{v}{1 + \frac{-9.8}{\frac{v \cdot v}{H}}}}{v}\right)\\
\mathbf{elif}\;v \leq 6.8 \cdot 10^{-29}:\\
\;\;\;\;\tan^{-1} \left(\frac{v}{\sqrt{H \cdot -19.6}}\right)\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1} \left(1 + H \cdot \frac{9.8 + \frac{H}{v \cdot v} \cdot 144.06}{v \cdot v}\right)\\
\end{array}
\end{array}
if v < -3.80000000000000012e-16Initial program 49.9%
atan-lowering-atan.f64N/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
metadata-evalN/A
metadata-eval49.9%
Simplified49.9%
Taylor expanded in v around -inf
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
mul-1-negN/A
neg-sub0N/A
--lowering--.f6493.9%
Simplified93.9%
sub0-negN/A
associate-/r*N/A
distribute-frac-neg2N/A
atan-negN/A
neg-lowering-neg.f64N/A
atan-lowering-atan.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6493.9%
Applied egg-rr93.9%
if -3.80000000000000012e-16 < v < 6.79999999999999945e-29Initial program 99.5%
atan-lowering-atan.f64N/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
metadata-evalN/A
metadata-eval99.5%
Simplified99.5%
Taylor expanded in v around 0
*-lowering-*.f6488.7%
Simplified88.7%
if 6.79999999999999945e-29 < v Initial program 52.7%
atan-lowering-atan.f64N/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
metadata-evalN/A
metadata-eval52.7%
Simplified52.7%
Taylor expanded in H around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
associate-*r/N/A
associate-*l/N/A
metadata-evalN/A
associate-*r/N/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate-*r/N/A
metadata-evalN/A
associate-*l/N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
Simplified91.0%
Taylor expanded in v around inf
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6491.2%
Simplified91.2%
Final simplification91.3%
(FPCore (v H)
:precision binary64
(if (<= v -1.6e-16)
(- 0.0 (atan (/ (/ v (+ 1.0 (/ -9.8 (/ (* v v) H)))) v)))
(if (<= v 3.15e-28)
(atan (* v (sqrt (/ -0.05102040816326531 H))))
(atan (+ 1.0 (* H (/ (+ 9.8 (* (/ H (* v v)) 144.06)) (* v v))))))))
double code(double v, double H) {
double tmp;
if (v <= -1.6e-16) {
tmp = 0.0 - atan(((v / (1.0 + (-9.8 / ((v * v) / H)))) / v));
} else if (v <= 3.15e-28) {
tmp = atan((v * sqrt((-0.05102040816326531 / H))));
} else {
tmp = atan((1.0 + (H * ((9.8 + ((H / (v * v)) * 144.06)) / (v * v)))));
}
return tmp;
}
real(8) function code(v, h)
real(8), intent (in) :: v
real(8), intent (in) :: h
real(8) :: tmp
if (v <= (-1.6d-16)) then
tmp = 0.0d0 - atan(((v / (1.0d0 + ((-9.8d0) / ((v * v) / h)))) / v))
else if (v <= 3.15d-28) then
tmp = atan((v * sqrt(((-0.05102040816326531d0) / h))))
else
tmp = atan((1.0d0 + (h * ((9.8d0 + ((h / (v * v)) * 144.06d0)) / (v * v)))))
end if
code = tmp
end function
public static double code(double v, double H) {
double tmp;
if (v <= -1.6e-16) {
tmp = 0.0 - Math.atan(((v / (1.0 + (-9.8 / ((v * v) / H)))) / v));
} else if (v <= 3.15e-28) {
tmp = Math.atan((v * Math.sqrt((-0.05102040816326531 / H))));
} else {
tmp = Math.atan((1.0 + (H * ((9.8 + ((H / (v * v)) * 144.06)) / (v * v)))));
}
return tmp;
}
def code(v, H): tmp = 0 if v <= -1.6e-16: tmp = 0.0 - math.atan(((v / (1.0 + (-9.8 / ((v * v) / H)))) / v)) elif v <= 3.15e-28: tmp = math.atan((v * math.sqrt((-0.05102040816326531 / H)))) else: tmp = math.atan((1.0 + (H * ((9.8 + ((H / (v * v)) * 144.06)) / (v * v))))) return tmp
function code(v, H) tmp = 0.0 if (v <= -1.6e-16) tmp = Float64(0.0 - atan(Float64(Float64(v / Float64(1.0 + Float64(-9.8 / Float64(Float64(v * v) / H)))) / v))); elseif (v <= 3.15e-28) tmp = atan(Float64(v * sqrt(Float64(-0.05102040816326531 / H)))); else tmp = atan(Float64(1.0 + Float64(H * Float64(Float64(9.8 + Float64(Float64(H / Float64(v * v)) * 144.06)) / Float64(v * v))))); end return tmp end
function tmp_2 = code(v, H) tmp = 0.0; if (v <= -1.6e-16) tmp = 0.0 - atan(((v / (1.0 + (-9.8 / ((v * v) / H)))) / v)); elseif (v <= 3.15e-28) tmp = atan((v * sqrt((-0.05102040816326531 / H)))); else tmp = atan((1.0 + (H * ((9.8 + ((H / (v * v)) * 144.06)) / (v * v))))); end tmp_2 = tmp; end
code[v_, H_] := If[LessEqual[v, -1.6e-16], N[(0.0 - N[ArcTan[N[(N[(v / N[(1.0 + N[(-9.8 / N[(N[(v * v), $MachinePrecision] / H), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / v), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[v, 3.15e-28], N[ArcTan[N[(v * N[Sqrt[N[(-0.05102040816326531 / H), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[N[(1.0 + N[(H * N[(N[(9.8 + N[(N[(H / N[(v * v), $MachinePrecision]), $MachinePrecision] * 144.06), $MachinePrecision]), $MachinePrecision] / N[(v * v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;v \leq -1.6 \cdot 10^{-16}:\\
\;\;\;\;0 - \tan^{-1} \left(\frac{\frac{v}{1 + \frac{-9.8}{\frac{v \cdot v}{H}}}}{v}\right)\\
\mathbf{elif}\;v \leq 3.15 \cdot 10^{-28}:\\
\;\;\;\;\tan^{-1} \left(v \cdot \sqrt{\frac{-0.05102040816326531}{H}}\right)\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1} \left(1 + H \cdot \frac{9.8 + \frac{H}{v \cdot v} \cdot 144.06}{v \cdot v}\right)\\
\end{array}
\end{array}
if v < -1.60000000000000011e-16Initial program 49.9%
atan-lowering-atan.f64N/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
metadata-evalN/A
metadata-eval49.9%
Simplified49.9%
Taylor expanded in v around -inf
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
mul-1-negN/A
neg-sub0N/A
--lowering--.f6493.9%
Simplified93.9%
sub0-negN/A
associate-/r*N/A
distribute-frac-neg2N/A
atan-negN/A
neg-lowering-neg.f64N/A
atan-lowering-atan.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6493.9%
Applied egg-rr93.9%
if -1.60000000000000011e-16 < v < 3.1499999999999999e-28Initial program 99.5%
atan-lowering-atan.f64N/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
metadata-evalN/A
metadata-eval99.5%
Simplified99.5%
Taylor expanded in v around 0
atan-lowering-atan.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6499.5%
Simplified99.5%
Taylor expanded in H around inf
/-lowering-/.f6488.6%
Simplified88.6%
if 3.1499999999999999e-28 < v Initial program 52.7%
atan-lowering-atan.f64N/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
metadata-evalN/A
metadata-eval52.7%
Simplified52.7%
Taylor expanded in H around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
associate-*r/N/A
associate-*l/N/A
metadata-evalN/A
associate-*r/N/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate-*r/N/A
metadata-evalN/A
associate-*l/N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
Simplified91.0%
Taylor expanded in v around inf
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6491.2%
Simplified91.2%
Final simplification91.3%
(FPCore (v H) :precision binary64 (if (<= v 7.5e-187) (atan (/ v (* v (- -1.0 (/ (* -9.8 (/ H v)) v))))) (atan (+ 1.0 (* H (/ (+ 9.8 (* (/ H (* v v)) 144.06)) (* v v)))))))
double code(double v, double H) {
double tmp;
if (v <= 7.5e-187) {
tmp = atan((v / (v * (-1.0 - ((-9.8 * (H / v)) / v)))));
} else {
tmp = atan((1.0 + (H * ((9.8 + ((H / (v * v)) * 144.06)) / (v * v)))));
}
return tmp;
}
real(8) function code(v, h)
real(8), intent (in) :: v
real(8), intent (in) :: h
real(8) :: tmp
if (v <= 7.5d-187) then
tmp = atan((v / (v * ((-1.0d0) - (((-9.8d0) * (h / v)) / v)))))
else
tmp = atan((1.0d0 + (h * ((9.8d0 + ((h / (v * v)) * 144.06d0)) / (v * v)))))
end if
code = tmp
end function
public static double code(double v, double H) {
double tmp;
if (v <= 7.5e-187) {
tmp = Math.atan((v / (v * (-1.0 - ((-9.8 * (H / v)) / v)))));
} else {
tmp = Math.atan((1.0 + (H * ((9.8 + ((H / (v * v)) * 144.06)) / (v * v)))));
}
return tmp;
}
def code(v, H): tmp = 0 if v <= 7.5e-187: tmp = math.atan((v / (v * (-1.0 - ((-9.8 * (H / v)) / v))))) else: tmp = math.atan((1.0 + (H * ((9.8 + ((H / (v * v)) * 144.06)) / (v * v))))) return tmp
function code(v, H) tmp = 0.0 if (v <= 7.5e-187) tmp = atan(Float64(v / Float64(v * Float64(-1.0 - Float64(Float64(-9.8 * Float64(H / v)) / v))))); else tmp = atan(Float64(1.0 + Float64(H * Float64(Float64(9.8 + Float64(Float64(H / Float64(v * v)) * 144.06)) / Float64(v * v))))); end return tmp end
function tmp_2 = code(v, H) tmp = 0.0; if (v <= 7.5e-187) tmp = atan((v / (v * (-1.0 - ((-9.8 * (H / v)) / v))))); else tmp = atan((1.0 + (H * ((9.8 + ((H / (v * v)) * 144.06)) / (v * v))))); end tmp_2 = tmp; end
code[v_, H_] := If[LessEqual[v, 7.5e-187], N[ArcTan[N[(v / N[(v * N[(-1.0 - N[(N[(-9.8 * N[(H / v), $MachinePrecision]), $MachinePrecision] / v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[N[(1.0 + N[(H * N[(N[(9.8 + N[(N[(H / N[(v * v), $MachinePrecision]), $MachinePrecision] * 144.06), $MachinePrecision]), $MachinePrecision] / N[(v * v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;v \leq 7.5 \cdot 10^{-187}:\\
\;\;\;\;\tan^{-1} \left(\frac{v}{v \cdot \left(-1 - \frac{-9.8 \cdot \frac{H}{v}}{v}\right)}\right)\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1} \left(1 + H \cdot \frac{9.8 + \frac{H}{v \cdot v} \cdot 144.06}{v \cdot v}\right)\\
\end{array}
\end{array}
if v < 7.5000000000000004e-187Initial program 69.9%
atan-lowering-atan.f64N/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
metadata-evalN/A
metadata-eval69.9%
Simplified69.9%
Taylor expanded in v around -inf
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
mul-1-negN/A
neg-sub0N/A
--lowering--.f6466.0%
Simplified66.0%
associate-/r*N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6466.0%
Applied egg-rr66.0%
if 7.5000000000000004e-187 < v Initial program 61.8%
atan-lowering-atan.f64N/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
metadata-evalN/A
metadata-eval61.8%
Simplified61.8%
Taylor expanded in H around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
associate-*r/N/A
associate-*l/N/A
metadata-evalN/A
associate-*r/N/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate-*r/N/A
metadata-evalN/A
associate-*l/N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
Simplified75.9%
Taylor expanded in v around inf
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6477.9%
Simplified77.9%
Final simplification71.0%
(FPCore (v H) :precision binary64 (if (<= v -9.5e-277) (- 0.0 (atan (/ (/ v (+ 1.0 (/ -9.8 (/ (* v v) H)))) v))) (atan (/ v (+ v (/ (* -9.8 H) v))))))
double code(double v, double H) {
double tmp;
if (v <= -9.5e-277) {
tmp = 0.0 - atan(((v / (1.0 + (-9.8 / ((v * v) / H)))) / v));
} else {
tmp = atan((v / (v + ((-9.8 * H) / v))));
}
return tmp;
}
real(8) function code(v, h)
real(8), intent (in) :: v
real(8), intent (in) :: h
real(8) :: tmp
if (v <= (-9.5d-277)) then
tmp = 0.0d0 - atan(((v / (1.0d0 + ((-9.8d0) / ((v * v) / h)))) / v))
else
tmp = atan((v / (v + (((-9.8d0) * h) / v))))
end if
code = tmp
end function
public static double code(double v, double H) {
double tmp;
if (v <= -9.5e-277) {
tmp = 0.0 - Math.atan(((v / (1.0 + (-9.8 / ((v * v) / H)))) / v));
} else {
tmp = Math.atan((v / (v + ((-9.8 * H) / v))));
}
return tmp;
}
def code(v, H): tmp = 0 if v <= -9.5e-277: tmp = 0.0 - math.atan(((v / (1.0 + (-9.8 / ((v * v) / H)))) / v)) else: tmp = math.atan((v / (v + ((-9.8 * H) / v)))) return tmp
function code(v, H) tmp = 0.0 if (v <= -9.5e-277) tmp = Float64(0.0 - atan(Float64(Float64(v / Float64(1.0 + Float64(-9.8 / Float64(Float64(v * v) / H)))) / v))); else tmp = atan(Float64(v / Float64(v + Float64(Float64(-9.8 * H) / v)))); end return tmp end
function tmp_2 = code(v, H) tmp = 0.0; if (v <= -9.5e-277) tmp = 0.0 - atan(((v / (1.0 + (-9.8 / ((v * v) / H)))) / v)); else tmp = atan((v / (v + ((-9.8 * H) / v)))); end tmp_2 = tmp; end
code[v_, H_] := If[LessEqual[v, -9.5e-277], N[(0.0 - N[ArcTan[N[(N[(v / N[(1.0 + N[(-9.8 / N[(N[(v * v), $MachinePrecision] / H), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / v), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[ArcTan[N[(v / N[(v + N[(N[(-9.8 * H), $MachinePrecision] / v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;v \leq -9.5 \cdot 10^{-277}:\\
\;\;\;\;0 - \tan^{-1} \left(\frac{\frac{v}{1 + \frac{-9.8}{\frac{v \cdot v}{H}}}}{v}\right)\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1} \left(\frac{v}{v + \frac{-9.8 \cdot H}{v}}\right)\\
\end{array}
\end{array}
if v < -9.5e-277Initial program 65.8%
atan-lowering-atan.f64N/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
metadata-evalN/A
metadata-eval65.8%
Simplified65.8%
Taylor expanded in v around -inf
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
mul-1-negN/A
neg-sub0N/A
--lowering--.f6472.7%
Simplified72.7%
sub0-negN/A
associate-/r*N/A
distribute-frac-neg2N/A
atan-negN/A
neg-lowering-neg.f64N/A
atan-lowering-atan.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6472.7%
Applied egg-rr72.7%
if -9.5e-277 < v Initial program 67.3%
atan-lowering-atan.f64N/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
metadata-evalN/A
metadata-eval67.3%
Simplified67.3%
Taylor expanded in H around 0
*-commutativeN/A
associate-*l/N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
associate-*r/N/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
associate-*r/N/A
*-commutativeN/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6468.8%
Simplified68.8%
Final simplification70.8%
(FPCore (v H) :precision binary64 (if (<= v -2e-280) (atan (/ v (* v (- -1.0 (/ (* -9.8 (/ H v)) v))))) (atan (/ v (+ v (/ (* -9.8 H) v))))))
double code(double v, double H) {
double tmp;
if (v <= -2e-280) {
tmp = atan((v / (v * (-1.0 - ((-9.8 * (H / v)) / v)))));
} else {
tmp = atan((v / (v + ((-9.8 * H) / v))));
}
return tmp;
}
real(8) function code(v, h)
real(8), intent (in) :: v
real(8), intent (in) :: h
real(8) :: tmp
if (v <= (-2d-280)) then
tmp = atan((v / (v * ((-1.0d0) - (((-9.8d0) * (h / v)) / v)))))
else
tmp = atan((v / (v + (((-9.8d0) * h) / v))))
end if
code = tmp
end function
public static double code(double v, double H) {
double tmp;
if (v <= -2e-280) {
tmp = Math.atan((v / (v * (-1.0 - ((-9.8 * (H / v)) / v)))));
} else {
tmp = Math.atan((v / (v + ((-9.8 * H) / v))));
}
return tmp;
}
def code(v, H): tmp = 0 if v <= -2e-280: tmp = math.atan((v / (v * (-1.0 - ((-9.8 * (H / v)) / v))))) else: tmp = math.atan((v / (v + ((-9.8 * H) / v)))) return tmp
function code(v, H) tmp = 0.0 if (v <= -2e-280) tmp = atan(Float64(v / Float64(v * Float64(-1.0 - Float64(Float64(-9.8 * Float64(H / v)) / v))))); else tmp = atan(Float64(v / Float64(v + Float64(Float64(-9.8 * H) / v)))); end return tmp end
function tmp_2 = code(v, H) tmp = 0.0; if (v <= -2e-280) tmp = atan((v / (v * (-1.0 - ((-9.8 * (H / v)) / v))))); else tmp = atan((v / (v + ((-9.8 * H) / v)))); end tmp_2 = tmp; end
code[v_, H_] := If[LessEqual[v, -2e-280], N[ArcTan[N[(v / N[(v * N[(-1.0 - N[(N[(-9.8 * N[(H / v), $MachinePrecision]), $MachinePrecision] / v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[N[(v / N[(v + N[(N[(-9.8 * H), $MachinePrecision] / v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;v \leq -2 \cdot 10^{-280}:\\
\;\;\;\;\tan^{-1} \left(\frac{v}{v \cdot \left(-1 - \frac{-9.8 \cdot \frac{H}{v}}{v}\right)}\right)\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1} \left(\frac{v}{v + \frac{-9.8 \cdot H}{v}}\right)\\
\end{array}
\end{array}
if v < -1.9999999999999999e-280Initial program 66.3%
atan-lowering-atan.f64N/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
metadata-evalN/A
metadata-eval66.3%
Simplified66.3%
Taylor expanded in v around -inf
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
mul-1-negN/A
neg-sub0N/A
--lowering--.f6472.4%
Simplified72.4%
associate-/r*N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6472.5%
Applied egg-rr72.5%
if -1.9999999999999999e-280 < v Initial program 66.8%
atan-lowering-atan.f64N/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
metadata-evalN/A
metadata-eval66.8%
Simplified66.8%
Taylor expanded in H around 0
*-commutativeN/A
associate-*l/N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
associate-*r/N/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
associate-*r/N/A
*-commutativeN/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6469.1%
Simplified69.1%
Final simplification70.8%
(FPCore (v H) :precision binary64 (if (<= v -2e-280) (atan (/ v (* v (- -1.0 (* -9.8 (/ H (* v v))))))) (atan (/ v (+ v (/ (* -9.8 H) v))))))
double code(double v, double H) {
double tmp;
if (v <= -2e-280) {
tmp = atan((v / (v * (-1.0 - (-9.8 * (H / (v * v)))))));
} else {
tmp = atan((v / (v + ((-9.8 * H) / v))));
}
return tmp;
}
real(8) function code(v, h)
real(8), intent (in) :: v
real(8), intent (in) :: h
real(8) :: tmp
if (v <= (-2d-280)) then
tmp = atan((v / (v * ((-1.0d0) - ((-9.8d0) * (h / (v * v)))))))
else
tmp = atan((v / (v + (((-9.8d0) * h) / v))))
end if
code = tmp
end function
public static double code(double v, double H) {
double tmp;
if (v <= -2e-280) {
tmp = Math.atan((v / (v * (-1.0 - (-9.8 * (H / (v * v)))))));
} else {
tmp = Math.atan((v / (v + ((-9.8 * H) / v))));
}
return tmp;
}
def code(v, H): tmp = 0 if v <= -2e-280: tmp = math.atan((v / (v * (-1.0 - (-9.8 * (H / (v * v))))))) else: tmp = math.atan((v / (v + ((-9.8 * H) / v)))) return tmp
function code(v, H) tmp = 0.0 if (v <= -2e-280) tmp = atan(Float64(v / Float64(v * Float64(-1.0 - Float64(-9.8 * Float64(H / Float64(v * v))))))); else tmp = atan(Float64(v / Float64(v + Float64(Float64(-9.8 * H) / v)))); end return tmp end
function tmp_2 = code(v, H) tmp = 0.0; if (v <= -2e-280) tmp = atan((v / (v * (-1.0 - (-9.8 * (H / (v * v))))))); else tmp = atan((v / (v + ((-9.8 * H) / v)))); end tmp_2 = tmp; end
code[v_, H_] := If[LessEqual[v, -2e-280], N[ArcTan[N[(v / N[(v * N[(-1.0 - N[(-9.8 * N[(H / N[(v * v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[N[(v / N[(v + N[(N[(-9.8 * H), $MachinePrecision] / v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;v \leq -2 \cdot 10^{-280}:\\
\;\;\;\;\tan^{-1} \left(\frac{v}{v \cdot \left(-1 - -9.8 \cdot \frac{H}{v \cdot v}\right)}\right)\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1} \left(\frac{v}{v + \frac{-9.8 \cdot H}{v}}\right)\\
\end{array}
\end{array}
if v < -1.9999999999999999e-280Initial program 66.3%
atan-lowering-atan.f64N/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
metadata-evalN/A
metadata-eval66.3%
Simplified66.3%
Taylor expanded in v around -inf
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
mul-1-negN/A
neg-sub0N/A
--lowering--.f6472.4%
Simplified72.4%
if -1.9999999999999999e-280 < v Initial program 66.8%
atan-lowering-atan.f64N/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
metadata-evalN/A
metadata-eval66.8%
Simplified66.8%
Taylor expanded in H around 0
*-commutativeN/A
associate-*l/N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
associate-*r/N/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
associate-*r/N/A
*-commutativeN/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6469.1%
Simplified69.1%
Final simplification70.8%
(FPCore (v H) :precision binary64 (if (<= v -5.2e-94) (atan -1.0) (if (<= v 4e-187) (atan (* v (* (/ v H) 0.10204081632653061))) (atan 1.0))))
double code(double v, double H) {
double tmp;
if (v <= -5.2e-94) {
tmp = atan(-1.0);
} else if (v <= 4e-187) {
tmp = atan((v * ((v / H) * 0.10204081632653061)));
} else {
tmp = atan(1.0);
}
return tmp;
}
real(8) function code(v, h)
real(8), intent (in) :: v
real(8), intent (in) :: h
real(8) :: tmp
if (v <= (-5.2d-94)) then
tmp = atan((-1.0d0))
else if (v <= 4d-187) then
tmp = atan((v * ((v / h) * 0.10204081632653061d0)))
else
tmp = atan(1.0d0)
end if
code = tmp
end function
public static double code(double v, double H) {
double tmp;
if (v <= -5.2e-94) {
tmp = Math.atan(-1.0);
} else if (v <= 4e-187) {
tmp = Math.atan((v * ((v / H) * 0.10204081632653061)));
} else {
tmp = Math.atan(1.0);
}
return tmp;
}
def code(v, H): tmp = 0 if v <= -5.2e-94: tmp = math.atan(-1.0) elif v <= 4e-187: tmp = math.atan((v * ((v / H) * 0.10204081632653061))) else: tmp = math.atan(1.0) return tmp
function code(v, H) tmp = 0.0 if (v <= -5.2e-94) tmp = atan(-1.0); elseif (v <= 4e-187) tmp = atan(Float64(v * Float64(Float64(v / H) * 0.10204081632653061))); else tmp = atan(1.0); end return tmp end
function tmp_2 = code(v, H) tmp = 0.0; if (v <= -5.2e-94) tmp = atan(-1.0); elseif (v <= 4e-187) tmp = atan((v * ((v / H) * 0.10204081632653061))); else tmp = atan(1.0); end tmp_2 = tmp; end
code[v_, H_] := If[LessEqual[v, -5.2e-94], N[ArcTan[-1.0], $MachinePrecision], If[LessEqual[v, 4e-187], N[ArcTan[N[(v * N[(N[(v / H), $MachinePrecision] * 0.10204081632653061), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[1.0], $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;v \leq -5.2 \cdot 10^{-94}:\\
\;\;\;\;\tan^{-1} -1\\
\mathbf{elif}\;v \leq 4 \cdot 10^{-187}:\\
\;\;\;\;\tan^{-1} \left(v \cdot \left(\frac{v}{H} \cdot 0.10204081632653061\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1} 1\\
\end{array}
\end{array}
if v < -5.19999999999999988e-94Initial program 57.9%
atan-lowering-atan.f64N/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
metadata-evalN/A
metadata-eval57.9%
Simplified57.9%
Taylor expanded in v around -inf
Simplified84.3%
if -5.19999999999999988e-94 < v < 4.0000000000000001e-187Initial program 99.5%
atan-lowering-atan.f64N/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
metadata-evalN/A
metadata-eval99.5%
Simplified99.5%
Taylor expanded in v around -inf
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
mul-1-negN/A
neg-sub0N/A
--lowering--.f6420.8%
Simplified20.8%
Taylor expanded in H around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6420.8%
Simplified20.8%
Taylor expanded in v around 0
associate-*r/N/A
*-commutativeN/A
associate-*r/N/A
metadata-evalN/A
associate-*r/N/A
atan-lowering-atan.f64N/A
associate-*r/N/A
metadata-evalN/A
associate-*r/N/A
associate-*l/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6420.8%
Simplified20.8%
associate-/l*N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6420.9%
Applied egg-rr20.9%
if 4.0000000000000001e-187 < v Initial program 61.8%
atan-lowering-atan.f64N/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
metadata-evalN/A
metadata-eval61.8%
Simplified61.8%
Taylor expanded in v around inf
Simplified76.3%
(FPCore (v H) :precision binary64 (if (<= v -1.02e-103) (atan -1.0) (atan (/ v (+ v (/ (* -9.8 H) v))))))
double code(double v, double H) {
double tmp;
if (v <= -1.02e-103) {
tmp = atan(-1.0);
} else {
tmp = atan((v / (v + ((-9.8 * H) / v))));
}
return tmp;
}
real(8) function code(v, h)
real(8), intent (in) :: v
real(8), intent (in) :: h
real(8) :: tmp
if (v <= (-1.02d-103)) then
tmp = atan((-1.0d0))
else
tmp = atan((v / (v + (((-9.8d0) * h) / v))))
end if
code = tmp
end function
public static double code(double v, double H) {
double tmp;
if (v <= -1.02e-103) {
tmp = Math.atan(-1.0);
} else {
tmp = Math.atan((v / (v + ((-9.8 * H) / v))));
}
return tmp;
}
def code(v, H): tmp = 0 if v <= -1.02e-103: tmp = math.atan(-1.0) else: tmp = math.atan((v / (v + ((-9.8 * H) / v)))) return tmp
function code(v, H) tmp = 0.0 if (v <= -1.02e-103) tmp = atan(-1.0); else tmp = atan(Float64(v / Float64(v + Float64(Float64(-9.8 * H) / v)))); end return tmp end
function tmp_2 = code(v, H) tmp = 0.0; if (v <= -1.02e-103) tmp = atan(-1.0); else tmp = atan((v / (v + ((-9.8 * H) / v)))); end tmp_2 = tmp; end
code[v_, H_] := If[LessEqual[v, -1.02e-103], N[ArcTan[-1.0], $MachinePrecision], N[ArcTan[N[(v / N[(v + N[(N[(-9.8 * H), $MachinePrecision] / v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;v \leq -1.02 \cdot 10^{-103}:\\
\;\;\;\;\tan^{-1} -1\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1} \left(\frac{v}{v + \frac{-9.8 \cdot H}{v}}\right)\\
\end{array}
\end{array}
if v < -1.01999999999999998e-103Initial program 58.2%
atan-lowering-atan.f64N/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
metadata-evalN/A
metadata-eval58.2%
Simplified58.2%
Taylor expanded in v around -inf
Simplified83.5%
if -1.01999999999999998e-103 < v Initial program 72.5%
atan-lowering-atan.f64N/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
metadata-evalN/A
metadata-eval72.5%
Simplified72.5%
Taylor expanded in H around 0
*-commutativeN/A
associate-*l/N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
associate-*r/N/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
associate-*r/N/A
*-commutativeN/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6461.6%
Simplified61.6%
Final simplification70.8%
(FPCore (v H) :precision binary64 (if (<= v -5e-311) (atan -1.0) (atan 1.0)))
double code(double v, double H) {
double tmp;
if (v <= -5e-311) {
tmp = atan(-1.0);
} else {
tmp = atan(1.0);
}
return tmp;
}
real(8) function code(v, h)
real(8), intent (in) :: v
real(8), intent (in) :: h
real(8) :: tmp
if (v <= (-5d-311)) then
tmp = atan((-1.0d0))
else
tmp = atan(1.0d0)
end if
code = tmp
end function
public static double code(double v, double H) {
double tmp;
if (v <= -5e-311) {
tmp = Math.atan(-1.0);
} else {
tmp = Math.atan(1.0);
}
return tmp;
}
def code(v, H): tmp = 0 if v <= -5e-311: tmp = math.atan(-1.0) else: tmp = math.atan(1.0) return tmp
function code(v, H) tmp = 0.0 if (v <= -5e-311) tmp = atan(-1.0); else tmp = atan(1.0); end return tmp end
function tmp_2 = code(v, H) tmp = 0.0; if (v <= -5e-311) tmp = atan(-1.0); else tmp = atan(1.0); end tmp_2 = tmp; end
code[v_, H_] := If[LessEqual[v, -5e-311], N[ArcTan[-1.0], $MachinePrecision], N[ArcTan[1.0], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;v \leq -5 \cdot 10^{-311}:\\
\;\;\;\;\tan^{-1} -1\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1} 1\\
\end{array}
\end{array}
if v < -5.00000000000023e-311Initial program 66.5%
atan-lowering-atan.f64N/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
metadata-evalN/A
metadata-eval66.5%
Simplified66.5%
Taylor expanded in v around -inf
Simplified67.5%
if -5.00000000000023e-311 < v Initial program 66.5%
atan-lowering-atan.f64N/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
metadata-evalN/A
metadata-eval66.5%
Simplified66.5%
Taylor expanded in v around inf
Simplified67.4%
(FPCore (v H) :precision binary64 (atan -1.0))
double code(double v, double H) {
return atan(-1.0);
}
real(8) function code(v, h)
real(8), intent (in) :: v
real(8), intent (in) :: h
code = atan((-1.0d0))
end function
public static double code(double v, double H) {
return Math.atan(-1.0);
}
def code(v, H): return math.atan(-1.0)
function code(v, H) return atan(-1.0) end
function tmp = code(v, H) tmp = atan(-1.0); end
code[v_, H_] := N[ArcTan[-1.0], $MachinePrecision]
\begin{array}{l}
\\
\tan^{-1} -1
\end{array}
Initial program 66.5%
atan-lowering-atan.f64N/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
metadata-evalN/A
metadata-eval66.5%
Simplified66.5%
Taylor expanded in v around -inf
Simplified36.2%
herbie shell --seed 2024160
(FPCore (v H)
:name "Optimal throwing angle"
:precision binary64
(atan (/ v (sqrt (- (* v v) (* (* 2.0 9.8) H))))))