
(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 16 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 -2.2e-45)
(atan (- 0.0 (sqrt (/ 1.0 (+ 1.0 (/ (* -19.6 H) (* v v)))))))
(if (<= v 6e-70)
(atan (/ v (pow (/ 1.0 (+ (* -19.6 H) (* v v))) -0.5)))
(atan (/ 1.0 (sqrt (+ 1.0 (* H (/ -19.6 (* v v))))))))))
double code(double v, double H) {
double tmp;
if (v <= -2.2e-45) {
tmp = atan((0.0 - sqrt((1.0 / (1.0 + ((-19.6 * H) / (v * v)))))));
} else if (v <= 6e-70) {
tmp = atan((v / pow((1.0 / ((-19.6 * H) + (v * v))), -0.5)));
} else {
tmp = atan((1.0 / sqrt((1.0 + (H * (-19.6 / (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 <= (-2.2d-45)) then
tmp = atan((0.0d0 - sqrt((1.0d0 / (1.0d0 + (((-19.6d0) * h) / (v * v)))))))
else if (v <= 6d-70) then
tmp = atan((v / ((1.0d0 / (((-19.6d0) * h) + (v * v))) ** (-0.5d0))))
else
tmp = atan((1.0d0 / sqrt((1.0d0 + (h * ((-19.6d0) / (v * v)))))))
end if
code = tmp
end function
public static double code(double v, double H) {
double tmp;
if (v <= -2.2e-45) {
tmp = Math.atan((0.0 - Math.sqrt((1.0 / (1.0 + ((-19.6 * H) / (v * v)))))));
} else if (v <= 6e-70) {
tmp = Math.atan((v / Math.pow((1.0 / ((-19.6 * H) + (v * v))), -0.5)));
} else {
tmp = Math.atan((1.0 / Math.sqrt((1.0 + (H * (-19.6 / (v * v)))))));
}
return tmp;
}
def code(v, H): tmp = 0 if v <= -2.2e-45: tmp = math.atan((0.0 - math.sqrt((1.0 / (1.0 + ((-19.6 * H) / (v * v))))))) elif v <= 6e-70: tmp = math.atan((v / math.pow((1.0 / ((-19.6 * H) + (v * v))), -0.5))) else: tmp = math.atan((1.0 / math.sqrt((1.0 + (H * (-19.6 / (v * v))))))) return tmp
function code(v, H) tmp = 0.0 if (v <= -2.2e-45) tmp = atan(Float64(0.0 - sqrt(Float64(1.0 / Float64(1.0 + Float64(Float64(-19.6 * H) / Float64(v * v))))))); elseif (v <= 6e-70) tmp = atan(Float64(v / (Float64(1.0 / Float64(Float64(-19.6 * H) + Float64(v * v))) ^ -0.5))); else tmp = atan(Float64(1.0 / sqrt(Float64(1.0 + Float64(H * Float64(-19.6 / Float64(v * v))))))); end return tmp end
function tmp_2 = code(v, H) tmp = 0.0; if (v <= -2.2e-45) tmp = atan((0.0 - sqrt((1.0 / (1.0 + ((-19.6 * H) / (v * v))))))); elseif (v <= 6e-70) tmp = atan((v / ((1.0 / ((-19.6 * H) + (v * v))) ^ -0.5))); else tmp = atan((1.0 / sqrt((1.0 + (H * (-19.6 / (v * v))))))); end tmp_2 = tmp; end
code[v_, H_] := If[LessEqual[v, -2.2e-45], N[ArcTan[N[(0.0 - N[Sqrt[N[(1.0 / N[(1.0 + N[(N[(-19.6 * H), $MachinePrecision] / N[(v * v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[v, 6e-70], N[ArcTan[N[(v / N[Power[N[(1.0 / N[(N[(-19.6 * H), $MachinePrecision] + N[(v * v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], -0.5], $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[N[(1.0 / N[Sqrt[N[(1.0 + N[(H * N[(-19.6 / N[(v * v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;v \leq -2.2 \cdot 10^{-45}:\\
\;\;\;\;\tan^{-1} \left(0 - \sqrt{\frac{1}{1 + \frac{-19.6 \cdot H}{v \cdot v}}}\right)\\
\mathbf{elif}\;v \leq 6 \cdot 10^{-70}:\\
\;\;\;\;\tan^{-1} \left(\frac{v}{{\left(\frac{1}{-19.6 \cdot H + v \cdot v}\right)}^{-0.5}}\right)\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1} \left(\frac{1}{\sqrt{1 + H \cdot \frac{-19.6}{v \cdot v}}}\right)\\
\end{array}
\end{array}
if v < -2.19999999999999993e-45Initial program 58.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-eval58.1%
Simplified58.1%
Taylor expanded in v around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6458.1%
Simplified58.1%
Taylor expanded in v around -inf
atan-lowering-atan.f64N/A
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
sqrt-lowering-sqrt.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
if -2.19999999999999993e-45 < v < 6.0000000000000003e-70Initial 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%
pow1/2N/A
flip3-+N/A
clear-numN/A
inv-powN/A
pow-powN/A
metadata-evalN/A
metadata-evalN/A
pow-lowering-pow.f64N/A
Applied egg-rr99.6%
if 6.0000000000000003e-70 < v Initial program 60.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-eval60.5%
Simplified60.5%
Taylor expanded in v around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6460.5%
Simplified60.5%
atan-lowering-atan.f64N/A
pow1/2N/A
unpow-prod-downN/A
associate-/r*N/A
pow1/2N/A
sqrt-prodN/A
rem-square-sqrtN/A
*-inversesN/A
/-lowering-/.f64N/A
pow1/2N/A
sqrt-lowering-sqrt.f64N/A
+-lowering-+.f64N/A
associate-*r/N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64100.0%
Applied egg-rr100.0%
Final simplification99.9%
(FPCore (v H)
:precision binary64
(if (<= v -5e+154)
(atan -1.0)
(if (<= v 6.8e-69)
(atan (/ v (pow (/ 1.0 (+ (* -19.6 H) (* v v))) -0.5)))
(atan (/ 1.0 (sqrt (+ 1.0 (* H (/ -19.6 (* v v))))))))))
double code(double v, double H) {
double tmp;
if (v <= -5e+154) {
tmp = atan(-1.0);
} else if (v <= 6.8e-69) {
tmp = atan((v / pow((1.0 / ((-19.6 * H) + (v * v))), -0.5)));
} else {
tmp = atan((1.0 / sqrt((1.0 + (H * (-19.6 / (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 <= (-5d+154)) then
tmp = atan((-1.0d0))
else if (v <= 6.8d-69) then
tmp = atan((v / ((1.0d0 / (((-19.6d0) * h) + (v * v))) ** (-0.5d0))))
else
tmp = atan((1.0d0 / sqrt((1.0d0 + (h * ((-19.6d0) / (v * v)))))))
end if
code = tmp
end function
public static double code(double v, double H) {
double tmp;
if (v <= -5e+154) {
tmp = Math.atan(-1.0);
} else if (v <= 6.8e-69) {
tmp = Math.atan((v / Math.pow((1.0 / ((-19.6 * H) + (v * v))), -0.5)));
} else {
tmp = Math.atan((1.0 / Math.sqrt((1.0 + (H * (-19.6 / (v * v)))))));
}
return tmp;
}
def code(v, H): tmp = 0 if v <= -5e+154: tmp = math.atan(-1.0) elif v <= 6.8e-69: tmp = math.atan((v / math.pow((1.0 / ((-19.6 * H) + (v * v))), -0.5))) else: tmp = math.atan((1.0 / math.sqrt((1.0 + (H * (-19.6 / (v * v))))))) return tmp
function code(v, H) tmp = 0.0 if (v <= -5e+154) tmp = atan(-1.0); elseif (v <= 6.8e-69) tmp = atan(Float64(v / (Float64(1.0 / Float64(Float64(-19.6 * H) + Float64(v * v))) ^ -0.5))); else tmp = atan(Float64(1.0 / sqrt(Float64(1.0 + Float64(H * Float64(-19.6 / Float64(v * v))))))); end return tmp end
function tmp_2 = code(v, H) tmp = 0.0; if (v <= -5e+154) tmp = atan(-1.0); elseif (v <= 6.8e-69) tmp = atan((v / ((1.0 / ((-19.6 * H) + (v * v))) ^ -0.5))); else tmp = atan((1.0 / sqrt((1.0 + (H * (-19.6 / (v * v))))))); end tmp_2 = tmp; end
code[v_, H_] := If[LessEqual[v, -5e+154], N[ArcTan[-1.0], $MachinePrecision], If[LessEqual[v, 6.8e-69], N[ArcTan[N[(v / N[Power[N[(1.0 / N[(N[(-19.6 * H), $MachinePrecision] + N[(v * v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], -0.5], $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[N[(1.0 / N[Sqrt[N[(1.0 + N[(H * N[(-19.6 / N[(v * v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;v \leq -5 \cdot 10^{+154}:\\
\;\;\;\;\tan^{-1} -1\\
\mathbf{elif}\;v \leq 6.8 \cdot 10^{-69}:\\
\;\;\;\;\tan^{-1} \left(\frac{v}{{\left(\frac{1}{-19.6 \cdot H + v \cdot v}\right)}^{-0.5}}\right)\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1} \left(\frac{1}{\sqrt{1 + H \cdot \frac{-19.6}{v \cdot v}}}\right)\\
\end{array}
\end{array}
if v < -5.00000000000000004e154Initial program 3.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-eval3.1%
Simplified3.1%
Taylor expanded in v around -inf
Simplified100.0%
if -5.00000000000000004e154 < v < 6.80000000000000016e-69Initial program 99.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-eval99.7%
Simplified99.7%
pow1/2N/A
flip3-+N/A
clear-numN/A
inv-powN/A
pow-powN/A
metadata-evalN/A
metadata-evalN/A
pow-lowering-pow.f64N/A
Applied egg-rr99.7%
if 6.80000000000000016e-69 < v Initial program 60.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-eval60.5%
Simplified60.5%
Taylor expanded in v around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6460.5%
Simplified60.5%
atan-lowering-atan.f64N/A
pow1/2N/A
unpow-prod-downN/A
associate-/r*N/A
pow1/2N/A
sqrt-prodN/A
rem-square-sqrtN/A
*-inversesN/A
/-lowering-/.f64N/A
pow1/2N/A
sqrt-lowering-sqrt.f64N/A
+-lowering-+.f64N/A
associate-*r/N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64100.0%
Applied egg-rr100.0%
Final simplification99.8%
(FPCore (v H)
:precision binary64
(if (<= v -5e+154)
(atan -1.0)
(if (<= v 6.2e-69)
(atan (* v (pow (+ (* -19.6 H) (* v v)) -0.5)))
(atan (/ 1.0 (sqrt (+ 1.0 (* H (/ -19.6 (* v v))))))))))
double code(double v, double H) {
double tmp;
if (v <= -5e+154) {
tmp = atan(-1.0);
} else if (v <= 6.2e-69) {
tmp = atan((v * pow(((-19.6 * H) + (v * v)), -0.5)));
} else {
tmp = atan((1.0 / sqrt((1.0 + (H * (-19.6 / (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 <= (-5d+154)) then
tmp = atan((-1.0d0))
else if (v <= 6.2d-69) then
tmp = atan((v * ((((-19.6d0) * h) + (v * v)) ** (-0.5d0))))
else
tmp = atan((1.0d0 / sqrt((1.0d0 + (h * ((-19.6d0) / (v * v)))))))
end if
code = tmp
end function
public static double code(double v, double H) {
double tmp;
if (v <= -5e+154) {
tmp = Math.atan(-1.0);
} else if (v <= 6.2e-69) {
tmp = Math.atan((v * Math.pow(((-19.6 * H) + (v * v)), -0.5)));
} else {
tmp = Math.atan((1.0 / Math.sqrt((1.0 + (H * (-19.6 / (v * v)))))));
}
return tmp;
}
def code(v, H): tmp = 0 if v <= -5e+154: tmp = math.atan(-1.0) elif v <= 6.2e-69: tmp = math.atan((v * math.pow(((-19.6 * H) + (v * v)), -0.5))) else: tmp = math.atan((1.0 / math.sqrt((1.0 + (H * (-19.6 / (v * v))))))) return tmp
function code(v, H) tmp = 0.0 if (v <= -5e+154) tmp = atan(-1.0); elseif (v <= 6.2e-69) tmp = atan(Float64(v * (Float64(Float64(-19.6 * H) + Float64(v * v)) ^ -0.5))); else tmp = atan(Float64(1.0 / sqrt(Float64(1.0 + Float64(H * Float64(-19.6 / Float64(v * v))))))); end return tmp end
function tmp_2 = code(v, H) tmp = 0.0; if (v <= -5e+154) tmp = atan(-1.0); elseif (v <= 6.2e-69) tmp = atan((v * (((-19.6 * H) + (v * v)) ^ -0.5))); else tmp = atan((1.0 / sqrt((1.0 + (H * (-19.6 / (v * v))))))); end tmp_2 = tmp; end
code[v_, H_] := If[LessEqual[v, -5e+154], N[ArcTan[-1.0], $MachinePrecision], If[LessEqual[v, 6.2e-69], N[ArcTan[N[(v * N[Power[N[(N[(-19.6 * H), $MachinePrecision] + N[(v * v), $MachinePrecision]), $MachinePrecision], -0.5], $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[N[(1.0 / N[Sqrt[N[(1.0 + N[(H * N[(-19.6 / N[(v * v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;v \leq -5 \cdot 10^{+154}:\\
\;\;\;\;\tan^{-1} -1\\
\mathbf{elif}\;v \leq 6.2 \cdot 10^{-69}:\\
\;\;\;\;\tan^{-1} \left(v \cdot {\left(-19.6 \cdot H + v \cdot v\right)}^{-0.5}\right)\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1} \left(\frac{1}{\sqrt{1 + H \cdot \frac{-19.6}{v \cdot v}}}\right)\\
\end{array}
\end{array}
if v < -5.00000000000000004e154Initial program 3.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-eval3.1%
Simplified3.1%
Taylor expanded in v around -inf
Simplified100.0%
if -5.00000000000000004e154 < v < 6.1999999999999999e-69Initial program 99.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-eval99.7%
Simplified99.7%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
pow1/2N/A
pow-flipN/A
pow-lowering-pow.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
metadata-eval99.7%
Applied egg-rr99.7%
if 6.1999999999999999e-69 < v Initial program 60.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-eval60.5%
Simplified60.5%
Taylor expanded in v around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6460.5%
Simplified60.5%
atan-lowering-atan.f64N/A
pow1/2N/A
unpow-prod-downN/A
associate-/r*N/A
pow1/2N/A
sqrt-prodN/A
rem-square-sqrtN/A
*-inversesN/A
/-lowering-/.f64N/A
pow1/2N/A
sqrt-lowering-sqrt.f64N/A
+-lowering-+.f64N/A
associate-*r/N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64100.0%
Applied egg-rr100.0%
Final simplification99.8%
(FPCore (v H) :precision binary64 (if (<= v -5e+154) (atan -1.0) (if (<= v 1e+140) (atan (/ v (sqrt (+ (* -19.6 H) (* v v))))) (atan 1.0))))
double code(double v, double H) {
double tmp;
if (v <= -5e+154) {
tmp = atan(-1.0);
} else if (v <= 1e+140) {
tmp = atan((v / sqrt(((-19.6 * H) + (v * v)))));
} 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+154)) then
tmp = atan((-1.0d0))
else if (v <= 1d+140) then
tmp = atan((v / sqrt((((-19.6d0) * h) + (v * v)))))
else
tmp = atan(1.0d0)
end if
code = tmp
end function
public static double code(double v, double H) {
double tmp;
if (v <= -5e+154) {
tmp = Math.atan(-1.0);
} else if (v <= 1e+140) {
tmp = Math.atan((v / Math.sqrt(((-19.6 * H) + (v * v)))));
} else {
tmp = Math.atan(1.0);
}
return tmp;
}
def code(v, H): tmp = 0 if v <= -5e+154: tmp = math.atan(-1.0) elif v <= 1e+140: tmp = math.atan((v / math.sqrt(((-19.6 * H) + (v * v))))) else: tmp = math.atan(1.0) return tmp
function code(v, H) tmp = 0.0 if (v <= -5e+154) tmp = atan(-1.0); elseif (v <= 1e+140) tmp = atan(Float64(v / sqrt(Float64(Float64(-19.6 * H) + Float64(v * v))))); else tmp = atan(1.0); end return tmp end
function tmp_2 = code(v, H) tmp = 0.0; if (v <= -5e+154) tmp = atan(-1.0); elseif (v <= 1e+140) tmp = atan((v / sqrt(((-19.6 * H) + (v * v))))); else tmp = atan(1.0); end tmp_2 = tmp; end
code[v_, H_] := If[LessEqual[v, -5e+154], N[ArcTan[-1.0], $MachinePrecision], If[LessEqual[v, 1e+140], N[ArcTan[N[(v / N[Sqrt[N[(N[(-19.6 * H), $MachinePrecision] + N[(v * v), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[1.0], $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;v \leq -5 \cdot 10^{+154}:\\
\;\;\;\;\tan^{-1} -1\\
\mathbf{elif}\;v \leq 10^{+140}:\\
\;\;\;\;\tan^{-1} \left(\frac{v}{\sqrt{-19.6 \cdot H + v \cdot v}}\right)\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1} 1\\
\end{array}
\end{array}
if v < -5.00000000000000004e154Initial program 3.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-eval3.1%
Simplified3.1%
Taylor expanded in v around -inf
Simplified100.0%
if -5.00000000000000004e154 < v < 1.00000000000000006e140Initial program 99.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-eval99.7%
Simplified99.7%
if 1.00000000000000006e140 < v Initial program 8.4%
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-eval8.4%
Simplified8.4%
Taylor expanded in v around inf
Simplified100.0%
Final simplification99.8%
(FPCore (v H)
:precision binary64
(if (<= v -6.1e-7)
(atan (/ v (* v (- -1.0 (* (/ H (* v v)) -9.8)))))
(if (<= v 2.7e-44)
(atan (/ v (pow (/ -0.05102040816326531 H) -0.5)))
(atan (/ 1.0 (/ (* v (+ 1.0 (/ (/ (* H -9.8) v) v))) v))))))
double code(double v, double H) {
double tmp;
if (v <= -6.1e-7) {
tmp = atan((v / (v * (-1.0 - ((H / (v * v)) * -9.8)))));
} else if (v <= 2.7e-44) {
tmp = atan((v / pow((-0.05102040816326531 / H), -0.5)));
} else {
tmp = atan((1.0 / ((v * (1.0 + (((H * -9.8) / v) / 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 <= (-6.1d-7)) then
tmp = atan((v / (v * ((-1.0d0) - ((h / (v * v)) * (-9.8d0))))))
else if (v <= 2.7d-44) then
tmp = atan((v / (((-0.05102040816326531d0) / h) ** (-0.5d0))))
else
tmp = atan((1.0d0 / ((v * (1.0d0 + (((h * (-9.8d0)) / v) / v))) / v)))
end if
code = tmp
end function
public static double code(double v, double H) {
double tmp;
if (v <= -6.1e-7) {
tmp = Math.atan((v / (v * (-1.0 - ((H / (v * v)) * -9.8)))));
} else if (v <= 2.7e-44) {
tmp = Math.atan((v / Math.pow((-0.05102040816326531 / H), -0.5)));
} else {
tmp = Math.atan((1.0 / ((v * (1.0 + (((H * -9.8) / v) / v))) / v)));
}
return tmp;
}
def code(v, H): tmp = 0 if v <= -6.1e-7: tmp = math.atan((v / (v * (-1.0 - ((H / (v * v)) * -9.8))))) elif v <= 2.7e-44: tmp = math.atan((v / math.pow((-0.05102040816326531 / H), -0.5))) else: tmp = math.atan((1.0 / ((v * (1.0 + (((H * -9.8) / v) / v))) / v))) return tmp
function code(v, H) tmp = 0.0 if (v <= -6.1e-7) tmp = atan(Float64(v / Float64(v * Float64(-1.0 - Float64(Float64(H / Float64(v * v)) * -9.8))))); elseif (v <= 2.7e-44) tmp = atan(Float64(v / (Float64(-0.05102040816326531 / H) ^ -0.5))); else tmp = atan(Float64(1.0 / Float64(Float64(v * Float64(1.0 + Float64(Float64(Float64(H * -9.8) / v) / v))) / v))); end return tmp end
function tmp_2 = code(v, H) tmp = 0.0; if (v <= -6.1e-7) tmp = atan((v / (v * (-1.0 - ((H / (v * v)) * -9.8))))); elseif (v <= 2.7e-44) tmp = atan((v / ((-0.05102040816326531 / H) ^ -0.5))); else tmp = atan((1.0 / ((v * (1.0 + (((H * -9.8) / v) / v))) / v))); end tmp_2 = tmp; end
code[v_, H_] := If[LessEqual[v, -6.1e-7], N[ArcTan[N[(v / N[(v * N[(-1.0 - N[(N[(H / N[(v * v), $MachinePrecision]), $MachinePrecision] * -9.8), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[v, 2.7e-44], N[ArcTan[N[(v / N[Power[N[(-0.05102040816326531 / H), $MachinePrecision], -0.5], $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[N[(1.0 / N[(N[(v * N[(1.0 + N[(N[(N[(H * -9.8), $MachinePrecision] / v), $MachinePrecision] / v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / v), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;v \leq -6.1 \cdot 10^{-7}:\\
\;\;\;\;\tan^{-1} \left(\frac{v}{v \cdot \left(-1 - \frac{H}{v \cdot v} \cdot -9.8\right)}\right)\\
\mathbf{elif}\;v \leq 2.7 \cdot 10^{-44}:\\
\;\;\;\;\tan^{-1} \left(\frac{v}{{\left(\frac{-0.05102040816326531}{H}\right)}^{-0.5}}\right)\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1} \left(\frac{1}{\frac{v \cdot \left(1 + \frac{\frac{H \cdot -9.8}{v}}{v}\right)}{v}}\right)\\
\end{array}
\end{array}
if v < -6.09999999999999983e-7Initial program 56.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-eval56.3%
Simplified56.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--.f6496.0%
Simplified96.0%
if -6.09999999999999983e-7 < v < 2.6999999999999999e-44Initial 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%
pow1/2N/A
flip3-+N/A
clear-numN/A
inv-powN/A
pow-powN/A
metadata-evalN/A
metadata-evalN/A
pow-lowering-pow.f64N/A
Applied egg-rr99.6%
Taylor expanded in v around 0
rem-square-sqrtN/A
unpow2N/A
/-lowering-/.f64N/A
unpow2N/A
rem-square-sqrt90.1%
Simplified90.1%
if 2.6999999999999999e-44 < v Initial program 59.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-eval59.1%
Simplified59.1%
Taylor expanded in v around inf
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6488.3%
Simplified88.3%
associate-/r*N/A
*-inversesN/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
associate-/r*N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f6488.3%
Applied egg-rr88.3%
/-rgt-identityN/A
clear-numN/A
frac-2negN/A
metadata-evalN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
distribute-frac-neg2N/A
metadata-evalN/A
frac-2negN/A
associate-*r/N/A
*-inversesN/A
associate-/r/N/A
clear-numN/A
/-lowering-/.f64N/A
Applied egg-rr88.3%
Final simplification91.6%
(FPCore (v H)
:precision binary64
(if (<= v -7e-7)
(atan (/ v (* v (- -1.0 (* (/ H (* v v)) -9.8)))))
(if (<= v 2.5e-44)
(atan (/ v (sqrt (* -19.6 H))))
(atan (/ 1.0 (/ (* v (+ 1.0 (/ (/ (* H -9.8) v) v))) v))))))
double code(double v, double H) {
double tmp;
if (v <= -7e-7) {
tmp = atan((v / (v * (-1.0 - ((H / (v * v)) * -9.8)))));
} else if (v <= 2.5e-44) {
tmp = atan((v / sqrt((-19.6 * H))));
} else {
tmp = atan((1.0 / ((v * (1.0 + (((H * -9.8) / v) / 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 <= (-7d-7)) then
tmp = atan((v / (v * ((-1.0d0) - ((h / (v * v)) * (-9.8d0))))))
else if (v <= 2.5d-44) then
tmp = atan((v / sqrt(((-19.6d0) * h))))
else
tmp = atan((1.0d0 / ((v * (1.0d0 + (((h * (-9.8d0)) / v) / v))) / v)))
end if
code = tmp
end function
public static double code(double v, double H) {
double tmp;
if (v <= -7e-7) {
tmp = Math.atan((v / (v * (-1.0 - ((H / (v * v)) * -9.8)))));
} else if (v <= 2.5e-44) {
tmp = Math.atan((v / Math.sqrt((-19.6 * H))));
} else {
tmp = Math.atan((1.0 / ((v * (1.0 + (((H * -9.8) / v) / v))) / v)));
}
return tmp;
}
def code(v, H): tmp = 0 if v <= -7e-7: tmp = math.atan((v / (v * (-1.0 - ((H / (v * v)) * -9.8))))) elif v <= 2.5e-44: tmp = math.atan((v / math.sqrt((-19.6 * H)))) else: tmp = math.atan((1.0 / ((v * (1.0 + (((H * -9.8) / v) / v))) / v))) return tmp
function code(v, H) tmp = 0.0 if (v <= -7e-7) tmp = atan(Float64(v / Float64(v * Float64(-1.0 - Float64(Float64(H / Float64(v * v)) * -9.8))))); elseif (v <= 2.5e-44) tmp = atan(Float64(v / sqrt(Float64(-19.6 * H)))); else tmp = atan(Float64(1.0 / Float64(Float64(v * Float64(1.0 + Float64(Float64(Float64(H * -9.8) / v) / v))) / v))); end return tmp end
function tmp_2 = code(v, H) tmp = 0.0; if (v <= -7e-7) tmp = atan((v / (v * (-1.0 - ((H / (v * v)) * -9.8))))); elseif (v <= 2.5e-44) tmp = atan((v / sqrt((-19.6 * H)))); else tmp = atan((1.0 / ((v * (1.0 + (((H * -9.8) / v) / v))) / v))); end tmp_2 = tmp; end
code[v_, H_] := If[LessEqual[v, -7e-7], N[ArcTan[N[(v / N[(v * N[(-1.0 - N[(N[(H / N[(v * v), $MachinePrecision]), $MachinePrecision] * -9.8), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[v, 2.5e-44], N[ArcTan[N[(v / N[Sqrt[N[(-19.6 * H), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[N[(1.0 / N[(N[(v * N[(1.0 + N[(N[(N[(H * -9.8), $MachinePrecision] / v), $MachinePrecision] / v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / v), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;v \leq -7 \cdot 10^{-7}:\\
\;\;\;\;\tan^{-1} \left(\frac{v}{v \cdot \left(-1 - \frac{H}{v \cdot v} \cdot -9.8\right)}\right)\\
\mathbf{elif}\;v \leq 2.5 \cdot 10^{-44}:\\
\;\;\;\;\tan^{-1} \left(\frac{v}{\sqrt{-19.6 \cdot H}}\right)\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1} \left(\frac{1}{\frac{v \cdot \left(1 + \frac{\frac{H \cdot -9.8}{v}}{v}\right)}{v}}\right)\\
\end{array}
\end{array}
if v < -6.99999999999999968e-7Initial program 56.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-eval56.3%
Simplified56.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--.f6496.0%
Simplified96.0%
if -6.99999999999999968e-7 < v < 2.50000000000000019e-44Initial 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-*.f6489.9%
Simplified89.9%
if 2.50000000000000019e-44 < v Initial program 59.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-eval59.1%
Simplified59.1%
Taylor expanded in v around inf
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6488.3%
Simplified88.3%
associate-/r*N/A
*-inversesN/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
associate-/r*N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f6488.3%
Applied egg-rr88.3%
/-rgt-identityN/A
clear-numN/A
frac-2negN/A
metadata-evalN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
distribute-frac-neg2N/A
metadata-evalN/A
frac-2negN/A
associate-*r/N/A
*-inversesN/A
associate-/r/N/A
clear-numN/A
/-lowering-/.f64N/A
Applied egg-rr88.3%
Final simplification91.5%
(FPCore (v H)
:precision binary64
(let* ((t_0 (/ (/ (* H -9.8) v) v)))
(if (<= v -1e-289)
(atan (/ v (* v (- -1.0 t_0))))
(atan (/ 1.0 (/ (* v (+ 1.0 t_0)) v))))))
double code(double v, double H) {
double t_0 = ((H * -9.8) / v) / v;
double tmp;
if (v <= -1e-289) {
tmp = atan((v / (v * (-1.0 - t_0))));
} else {
tmp = atan((1.0 / ((v * (1.0 + t_0)) / v)));
}
return tmp;
}
real(8) function code(v, h)
real(8), intent (in) :: v
real(8), intent (in) :: h
real(8) :: t_0
real(8) :: tmp
t_0 = ((h * (-9.8d0)) / v) / v
if (v <= (-1d-289)) then
tmp = atan((v / (v * ((-1.0d0) - t_0))))
else
tmp = atan((1.0d0 / ((v * (1.0d0 + t_0)) / v)))
end if
code = tmp
end function
public static double code(double v, double H) {
double t_0 = ((H * -9.8) / v) / v;
double tmp;
if (v <= -1e-289) {
tmp = Math.atan((v / (v * (-1.0 - t_0))));
} else {
tmp = Math.atan((1.0 / ((v * (1.0 + t_0)) / v)));
}
return tmp;
}
def code(v, H): t_0 = ((H * -9.8) / v) / v tmp = 0 if v <= -1e-289: tmp = math.atan((v / (v * (-1.0 - t_0)))) else: tmp = math.atan((1.0 / ((v * (1.0 + t_0)) / v))) return tmp
function code(v, H) t_0 = Float64(Float64(Float64(H * -9.8) / v) / v) tmp = 0.0 if (v <= -1e-289) tmp = atan(Float64(v / Float64(v * Float64(-1.0 - t_0)))); else tmp = atan(Float64(1.0 / Float64(Float64(v * Float64(1.0 + t_0)) / v))); end return tmp end
function tmp_2 = code(v, H) t_0 = ((H * -9.8) / v) / v; tmp = 0.0; if (v <= -1e-289) tmp = atan((v / (v * (-1.0 - t_0)))); else tmp = atan((1.0 / ((v * (1.0 + t_0)) / v))); end tmp_2 = tmp; end
code[v_, H_] := Block[{t$95$0 = N[(N[(N[(H * -9.8), $MachinePrecision] / v), $MachinePrecision] / v), $MachinePrecision]}, If[LessEqual[v, -1e-289], N[ArcTan[N[(v / N[(v * N[(-1.0 - t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[N[(1.0 / N[(N[(v * N[(1.0 + t$95$0), $MachinePrecision]), $MachinePrecision] / v), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{\frac{H \cdot -9.8}{v}}{v}\\
\mathbf{if}\;v \leq -1 \cdot 10^{-289}:\\
\;\;\;\;\tan^{-1} \left(\frac{v}{v \cdot \left(-1 - t\_0\right)}\right)\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1} \left(\frac{1}{\frac{v \cdot \left(1 + t\_0\right)}{v}}\right)\\
\end{array}
\end{array}
if v < -1e-289Initial program 70.4%
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-eval70.4%
Simplified70.4%
pow1/2N/A
flip3-+N/A
clear-numN/A
inv-powN/A
pow-powN/A
metadata-evalN/A
metadata-evalN/A
pow-lowering-pow.f64N/A
Applied egg-rr70.4%
Taylor expanded in v around -inf
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r/N/A
unpow2N/A
associate-/r*N/A
associate-*r/N/A
/-lowering-/.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
mul-1-negN/A
neg-sub0N/A
--lowering--.f6473.3%
Simplified73.3%
if -1e-289 < v Initial program 71.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-eval71.8%
Simplified71.8%
Taylor expanded in v around inf
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6472.3%
Simplified72.3%
associate-/r*N/A
*-inversesN/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
associate-/r*N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f6472.3%
Applied egg-rr72.3%
/-rgt-identityN/A
clear-numN/A
frac-2negN/A
metadata-evalN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
distribute-frac-neg2N/A
metadata-evalN/A
frac-2negN/A
associate-*r/N/A
*-inversesN/A
associate-/r/N/A
clear-numN/A
/-lowering-/.f64N/A
Applied egg-rr72.3%
Final simplification72.8%
(FPCore (v H) :precision binary64 (if (<= v 5e-259) (atan (/ v (* v (- -1.0 (* (/ H (* v v)) -9.8))))) (atan (/ 1.0 (/ (* v (+ 1.0 (/ (/ (* H -9.8) v) v))) v)))))
double code(double v, double H) {
double tmp;
if (v <= 5e-259) {
tmp = atan((v / (v * (-1.0 - ((H / (v * v)) * -9.8)))));
} else {
tmp = atan((1.0 / ((v * (1.0 + (((H * -9.8) / v) / 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 <= 5d-259) then
tmp = atan((v / (v * ((-1.0d0) - ((h / (v * v)) * (-9.8d0))))))
else
tmp = atan((1.0d0 / ((v * (1.0d0 + (((h * (-9.8d0)) / v) / v))) / v)))
end if
code = tmp
end function
public static double code(double v, double H) {
double tmp;
if (v <= 5e-259) {
tmp = Math.atan((v / (v * (-1.0 - ((H / (v * v)) * -9.8)))));
} else {
tmp = Math.atan((1.0 / ((v * (1.0 + (((H * -9.8) / v) / v))) / v)));
}
return tmp;
}
def code(v, H): tmp = 0 if v <= 5e-259: tmp = math.atan((v / (v * (-1.0 - ((H / (v * v)) * -9.8))))) else: tmp = math.atan((1.0 / ((v * (1.0 + (((H * -9.8) / v) / v))) / v))) return tmp
function code(v, H) tmp = 0.0 if (v <= 5e-259) tmp = atan(Float64(v / Float64(v * Float64(-1.0 - Float64(Float64(H / Float64(v * v)) * -9.8))))); else tmp = atan(Float64(1.0 / Float64(Float64(v * Float64(1.0 + Float64(Float64(Float64(H * -9.8) / v) / v))) / v))); end return tmp end
function tmp_2 = code(v, H) tmp = 0.0; if (v <= 5e-259) tmp = atan((v / (v * (-1.0 - ((H / (v * v)) * -9.8))))); else tmp = atan((1.0 / ((v * (1.0 + (((H * -9.8) / v) / v))) / v))); end tmp_2 = tmp; end
code[v_, H_] := If[LessEqual[v, 5e-259], N[ArcTan[N[(v / N[(v * N[(-1.0 - N[(N[(H / N[(v * v), $MachinePrecision]), $MachinePrecision] * -9.8), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[N[(1.0 / N[(N[(v * N[(1.0 + N[(N[(N[(H * -9.8), $MachinePrecision] / v), $MachinePrecision] / v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / v), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;v \leq 5 \cdot 10^{-259}:\\
\;\;\;\;\tan^{-1} \left(\frac{v}{v \cdot \left(-1 - \frac{H}{v \cdot v} \cdot -9.8\right)}\right)\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1} \left(\frac{1}{\frac{v \cdot \left(1 + \frac{\frac{H \cdot -9.8}{v}}{v}\right)}{v}}\right)\\
\end{array}
\end{array}
if v < 4.99999999999999977e-259Initial program 72.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-eval72.8%
Simplified72.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.5%
Simplified72.5%
if 4.99999999999999977e-259 < v Initial program 68.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-eval68.7%
Simplified68.7%
Taylor expanded in v around inf
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6473.0%
Simplified73.0%
associate-/r*N/A
*-inversesN/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
associate-/r*N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f6473.1%
Applied egg-rr73.1%
/-rgt-identityN/A
clear-numN/A
frac-2negN/A
metadata-evalN/A
metadata-evalN/A
distribute-neg-inN/A
metadata-evalN/A
distribute-frac-neg2N/A
metadata-evalN/A
frac-2negN/A
associate-*r/N/A
*-inversesN/A
associate-/r/N/A
clear-numN/A
/-lowering-/.f64N/A
Applied egg-rr73.1%
Final simplification72.8%
(FPCore (v H) :precision binary64 (if (<= v -1.4e-100) (atan -1.0) (atan (/ 1.0 (+ 1.0 (/ -1.0 (/ (/ (* v v) H) 9.8)))))))
double code(double v, double H) {
double tmp;
if (v <= -1.4e-100) {
tmp = atan(-1.0);
} else {
tmp = atan((1.0 / (1.0 + (-1.0 / (((v * v) / H) / 9.8)))));
}
return tmp;
}
real(8) function code(v, h)
real(8), intent (in) :: v
real(8), intent (in) :: h
real(8) :: tmp
if (v <= (-1.4d-100)) then
tmp = atan((-1.0d0))
else
tmp = atan((1.0d0 / (1.0d0 + ((-1.0d0) / (((v * v) / h) / 9.8d0)))))
end if
code = tmp
end function
public static double code(double v, double H) {
double tmp;
if (v <= -1.4e-100) {
tmp = Math.atan(-1.0);
} else {
tmp = Math.atan((1.0 / (1.0 + (-1.0 / (((v * v) / H) / 9.8)))));
}
return tmp;
}
def code(v, H): tmp = 0 if v <= -1.4e-100: tmp = math.atan(-1.0) else: tmp = math.atan((1.0 / (1.0 + (-1.0 / (((v * v) / H) / 9.8))))) return tmp
function code(v, H) tmp = 0.0 if (v <= -1.4e-100) tmp = atan(-1.0); else tmp = atan(Float64(1.0 / Float64(1.0 + Float64(-1.0 / Float64(Float64(Float64(v * v) / H) / 9.8))))); end return tmp end
function tmp_2 = code(v, H) tmp = 0.0; if (v <= -1.4e-100) tmp = atan(-1.0); else tmp = atan((1.0 / (1.0 + (-1.0 / (((v * v) / H) / 9.8))))); end tmp_2 = tmp; end
code[v_, H_] := If[LessEqual[v, -1.4e-100], N[ArcTan[-1.0], $MachinePrecision], N[ArcTan[N[(1.0 / N[(1.0 + N[(-1.0 / N[(N[(N[(v * v), $MachinePrecision] / H), $MachinePrecision] / 9.8), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;v \leq -1.4 \cdot 10^{-100}:\\
\;\;\;\;\tan^{-1} -1\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1} \left(\frac{1}{1 + \frac{-1}{\frac{\frac{v \cdot v}{H}}{9.8}}}\right)\\
\end{array}
\end{array}
if v < -1.39999999999999998e-100Initial program 61.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-eval61.7%
Simplified61.7%
Taylor expanded in v around -inf
Simplified87.9%
if -1.39999999999999998e-100 < v Initial program 77.4%
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-eval77.4%
Simplified77.4%
Taylor expanded in v around inf
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6462.1%
Simplified62.1%
associate-/r*N/A
*-inversesN/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
associate-/r*N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f6462.1%
Applied egg-rr62.1%
metadata-evalN/A
associate-*r/N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
associate-/r/N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f6462.1%
Applied egg-rr62.1%
(FPCore (v H)
:precision binary64
(if (<= v -1.4e-100)
(atan -1.0)
(if (<= v 7.2e-140)
(atan (/ -0.10204081632653061 (/ H (* v v))))
(atan 1.0))))
double code(double v, double H) {
double tmp;
if (v <= -1.4e-100) {
tmp = atan(-1.0);
} else if (v <= 7.2e-140) {
tmp = atan((-0.10204081632653061 / (H / (v * v))));
} 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 <= (-1.4d-100)) then
tmp = atan((-1.0d0))
else if (v <= 7.2d-140) then
tmp = atan(((-0.10204081632653061d0) / (h / (v * v))))
else
tmp = atan(1.0d0)
end if
code = tmp
end function
public static double code(double v, double H) {
double tmp;
if (v <= -1.4e-100) {
tmp = Math.atan(-1.0);
} else if (v <= 7.2e-140) {
tmp = Math.atan((-0.10204081632653061 / (H / (v * v))));
} else {
tmp = Math.atan(1.0);
}
return tmp;
}
def code(v, H): tmp = 0 if v <= -1.4e-100: tmp = math.atan(-1.0) elif v <= 7.2e-140: tmp = math.atan((-0.10204081632653061 / (H / (v * v)))) else: tmp = math.atan(1.0) return tmp
function code(v, H) tmp = 0.0 if (v <= -1.4e-100) tmp = atan(-1.0); elseif (v <= 7.2e-140) tmp = atan(Float64(-0.10204081632653061 / Float64(H / Float64(v * v)))); else tmp = atan(1.0); end return tmp end
function tmp_2 = code(v, H) tmp = 0.0; if (v <= -1.4e-100) tmp = atan(-1.0); elseif (v <= 7.2e-140) tmp = atan((-0.10204081632653061 / (H / (v * v)))); else tmp = atan(1.0); end tmp_2 = tmp; end
code[v_, H_] := If[LessEqual[v, -1.4e-100], N[ArcTan[-1.0], $MachinePrecision], If[LessEqual[v, 7.2e-140], N[ArcTan[N[(-0.10204081632653061 / N[(H / N[(v * v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[1.0], $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;v \leq -1.4 \cdot 10^{-100}:\\
\;\;\;\;\tan^{-1} -1\\
\mathbf{elif}\;v \leq 7.2 \cdot 10^{-140}:\\
\;\;\;\;\tan^{-1} \left(\frac{-0.10204081632653061}{\frac{H}{v \cdot v}}\right)\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1} 1\\
\end{array}
\end{array}
if v < -1.39999999999999998e-100Initial program 61.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-eval61.7%
Simplified61.7%
Taylor expanded in v around -inf
Simplified87.9%
if -1.39999999999999998e-100 < v < 7.2000000000000001e-140Initial 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
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6428.9%
Simplified28.9%
associate-/r*N/A
*-inversesN/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
associate-/r*N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f6428.8%
Applied egg-rr28.8%
Taylor expanded in v around 0
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6428.9%
Simplified28.9%
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6428.9%
Applied egg-rr28.9%
if 7.2000000000000001e-140 < v Initial program 65.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-eval65.3%
Simplified65.3%
Taylor expanded in v around inf
Simplified79.6%
(FPCore (v H)
:precision binary64
(if (<= v -1.4e-100)
(atan -1.0)
(if (<= v 2.3e-144)
(atan (* (/ (* v v) H) -0.10204081632653061))
(atan 1.0))))
double code(double v, double H) {
double tmp;
if (v <= -1.4e-100) {
tmp = atan(-1.0);
} else if (v <= 2.3e-144) {
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 <= (-1.4d-100)) then
tmp = atan((-1.0d0))
else if (v <= 2.3d-144) 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 <= -1.4e-100) {
tmp = Math.atan(-1.0);
} else if (v <= 2.3e-144) {
tmp = Math.atan((((v * v) / H) * -0.10204081632653061));
} else {
tmp = Math.atan(1.0);
}
return tmp;
}
def code(v, H): tmp = 0 if v <= -1.4e-100: tmp = math.atan(-1.0) elif v <= 2.3e-144: 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 <= -1.4e-100) tmp = atan(-1.0); elseif (v <= 2.3e-144) tmp = atan(Float64(Float64(Float64(v * v) / H) * -0.10204081632653061)); else tmp = atan(1.0); end return tmp end
function tmp_2 = code(v, H) tmp = 0.0; if (v <= -1.4e-100) tmp = atan(-1.0); elseif (v <= 2.3e-144) tmp = atan((((v * v) / H) * -0.10204081632653061)); else tmp = atan(1.0); end tmp_2 = tmp; end
code[v_, H_] := If[LessEqual[v, -1.4e-100], N[ArcTan[-1.0], $MachinePrecision], If[LessEqual[v, 2.3e-144], N[ArcTan[N[(N[(N[(v * v), $MachinePrecision] / H), $MachinePrecision] * -0.10204081632653061), $MachinePrecision]], $MachinePrecision], N[ArcTan[1.0], $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;v \leq -1.4 \cdot 10^{-100}:\\
\;\;\;\;\tan^{-1} -1\\
\mathbf{elif}\;v \leq 2.3 \cdot 10^{-144}:\\
\;\;\;\;\tan^{-1} \left(\frac{v \cdot v}{H} \cdot -0.10204081632653061\right)\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1} 1\\
\end{array}
\end{array}
if v < -1.39999999999999998e-100Initial program 61.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-eval61.7%
Simplified61.7%
Taylor expanded in v around -inf
Simplified87.9%
if -1.39999999999999998e-100 < v < 2.3e-144Initial 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
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6428.9%
Simplified28.9%
associate-/r*N/A
*-inversesN/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
associate-/r*N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f6428.8%
Applied egg-rr28.8%
Taylor expanded in v around 0
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6428.9%
Simplified28.9%
if 2.3e-144 < v Initial program 65.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-eval65.3%
Simplified65.3%
Taylor expanded in v around inf
Simplified79.6%
Final simplification72.3%
(FPCore (v H)
:precision binary64
(if (<= v -3.6e-177)
(atan -1.0)
(if (<= v 1.12e-141)
(atan (* -0.10204081632653061 (* v (/ v H))))
(atan 1.0))))
double code(double v, double H) {
double tmp;
if (v <= -3.6e-177) {
tmp = atan(-1.0);
} else if (v <= 1.12e-141) {
tmp = atan((-0.10204081632653061 * (v * (v / H))));
} 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 <= (-3.6d-177)) then
tmp = atan((-1.0d0))
else if (v <= 1.12d-141) then
tmp = atan(((-0.10204081632653061d0) * (v * (v / h))))
else
tmp = atan(1.0d0)
end if
code = tmp
end function
public static double code(double v, double H) {
double tmp;
if (v <= -3.6e-177) {
tmp = Math.atan(-1.0);
} else if (v <= 1.12e-141) {
tmp = Math.atan((-0.10204081632653061 * (v * (v / H))));
} else {
tmp = Math.atan(1.0);
}
return tmp;
}
def code(v, H): tmp = 0 if v <= -3.6e-177: tmp = math.atan(-1.0) elif v <= 1.12e-141: tmp = math.atan((-0.10204081632653061 * (v * (v / H)))) else: tmp = math.atan(1.0) return tmp
function code(v, H) tmp = 0.0 if (v <= -3.6e-177) tmp = atan(-1.0); elseif (v <= 1.12e-141) tmp = atan(Float64(-0.10204081632653061 * Float64(v * Float64(v / H)))); else tmp = atan(1.0); end return tmp end
function tmp_2 = code(v, H) tmp = 0.0; if (v <= -3.6e-177) tmp = atan(-1.0); elseif (v <= 1.12e-141) tmp = atan((-0.10204081632653061 * (v * (v / H)))); else tmp = atan(1.0); end tmp_2 = tmp; end
code[v_, H_] := If[LessEqual[v, -3.6e-177], N[ArcTan[-1.0], $MachinePrecision], If[LessEqual[v, 1.12e-141], N[ArcTan[N[(-0.10204081632653061 * N[(v * N[(v / H), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[ArcTan[1.0], $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;v \leq -3.6 \cdot 10^{-177}:\\
\;\;\;\;\tan^{-1} -1\\
\mathbf{elif}\;v \leq 1.12 \cdot 10^{-141}:\\
\;\;\;\;\tan^{-1} \left(-0.10204081632653061 \cdot \left(v \cdot \frac{v}{H}\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1} 1\\
\end{array}
\end{array}
if v < -3.59999999999999983e-177Initial program 65.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-eval65.3%
Simplified65.3%
Taylor expanded in v around -inf
Simplified80.0%
if -3.59999999999999983e-177 < v < 1.12000000000000002e-141Initial program 99.6%
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.6%
Simplified99.6%
Taylor expanded in v around inf
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6435.0%
Simplified35.0%
associate-/r*N/A
*-inversesN/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
associate-/r*N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f6435.0%
Applied egg-rr35.0%
Taylor expanded in v around 0
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6435.0%
Simplified35.0%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f6435.0%
Applied egg-rr35.0%
if 1.12000000000000002e-141 < v Initial program 65.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-eval65.3%
Simplified65.3%
Taylor expanded in v around inf
Simplified79.6%
Final simplification72.3%
(FPCore (v H) :precision binary64 (if (<= v -4.1e-178) (atan -1.0) (atan (/ 1.0 (+ 1.0 (* (/ H v) (/ -9.8 v)))))))
double code(double v, double H) {
double tmp;
if (v <= -4.1e-178) {
tmp = atan(-1.0);
} else {
tmp = atan((1.0 / (1.0 + ((H / v) * (-9.8 / v)))));
}
return tmp;
}
real(8) function code(v, h)
real(8), intent (in) :: v
real(8), intent (in) :: h
real(8) :: tmp
if (v <= (-4.1d-178)) then
tmp = atan((-1.0d0))
else
tmp = atan((1.0d0 / (1.0d0 + ((h / v) * ((-9.8d0) / v)))))
end if
code = tmp
end function
public static double code(double v, double H) {
double tmp;
if (v <= -4.1e-178) {
tmp = Math.atan(-1.0);
} else {
tmp = Math.atan((1.0 / (1.0 + ((H / v) * (-9.8 / v)))));
}
return tmp;
}
def code(v, H): tmp = 0 if v <= -4.1e-178: tmp = math.atan(-1.0) else: tmp = math.atan((1.0 / (1.0 + ((H / v) * (-9.8 / v))))) return tmp
function code(v, H) tmp = 0.0 if (v <= -4.1e-178) tmp = atan(-1.0); else tmp = atan(Float64(1.0 / Float64(1.0 + Float64(Float64(H / v) * Float64(-9.8 / v))))); end return tmp end
function tmp_2 = code(v, H) tmp = 0.0; if (v <= -4.1e-178) tmp = atan(-1.0); else tmp = atan((1.0 / (1.0 + ((H / v) * (-9.8 / v))))); end tmp_2 = tmp; end
code[v_, H_] := If[LessEqual[v, -4.1e-178], N[ArcTan[-1.0], $MachinePrecision], N[ArcTan[N[(1.0 / N[(1.0 + N[(N[(H / v), $MachinePrecision] * N[(-9.8 / v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;v \leq -4.1 \cdot 10^{-178}:\\
\;\;\;\;\tan^{-1} -1\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1} \left(\frac{1}{1 + \frac{H}{v} \cdot \frac{-9.8}{v}}\right)\\
\end{array}
\end{array}
if v < -4.0999999999999999e-178Initial program 65.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-eval65.3%
Simplified65.3%
Taylor expanded in v around -inf
Simplified80.0%
if -4.0999999999999999e-178 < v Initial program 75.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-eval75.7%
Simplified75.7%
Taylor expanded in v around inf
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6466.5%
Simplified66.5%
associate-/r*N/A
*-inversesN/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
clear-numN/A
associate-/r*N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f6466.5%
Applied egg-rr66.5%
associate-/r/N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f6466.5%
Applied egg-rr66.5%
(FPCore (v H) :precision binary64 (if (<= v -4.4e-176) (atan -1.0) (atan (/ v (+ v (/ (* H -9.8) v))))))
double code(double v, double H) {
double tmp;
if (v <= -4.4e-176) {
tmp = atan(-1.0);
} else {
tmp = atan((v / (v + ((H * -9.8) / v))));
}
return tmp;
}
real(8) function code(v, h)
real(8), intent (in) :: v
real(8), intent (in) :: h
real(8) :: tmp
if (v <= (-4.4d-176)) then
tmp = atan((-1.0d0))
else
tmp = atan((v / (v + ((h * (-9.8d0)) / v))))
end if
code = tmp
end function
public static double code(double v, double H) {
double tmp;
if (v <= -4.4e-176) {
tmp = Math.atan(-1.0);
} else {
tmp = Math.atan((v / (v + ((H * -9.8) / v))));
}
return tmp;
}
def code(v, H): tmp = 0 if v <= -4.4e-176: tmp = math.atan(-1.0) else: tmp = math.atan((v / (v + ((H * -9.8) / v)))) return tmp
function code(v, H) tmp = 0.0 if (v <= -4.4e-176) tmp = atan(-1.0); else tmp = atan(Float64(v / Float64(v + Float64(Float64(H * -9.8) / v)))); end return tmp end
function tmp_2 = code(v, H) tmp = 0.0; if (v <= -4.4e-176) tmp = atan(-1.0); else tmp = atan((v / (v + ((H * -9.8) / v)))); end tmp_2 = tmp; end
code[v_, H_] := If[LessEqual[v, -4.4e-176], N[ArcTan[-1.0], $MachinePrecision], N[ArcTan[N[(v / N[(v + N[(N[(H * -9.8), $MachinePrecision] / v), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;v \leq -4.4 \cdot 10^{-176}:\\
\;\;\;\;\tan^{-1} -1\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1} \left(\frac{v}{v + \frac{H \cdot -9.8}{v}}\right)\\
\end{array}
\end{array}
if v < -4.3999999999999997e-176Initial program 65.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-eval65.3%
Simplified65.3%
Taylor expanded in v around -inf
Simplified80.0%
if -4.3999999999999997e-176 < v Initial program 75.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-eval75.7%
Simplified75.7%
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-*.f6466.5%
Simplified66.5%
(FPCore (v H) :precision binary64 (if (<= v -5e-310) (atan -1.0) (atan 1.0)))
double code(double v, double H) {
double tmp;
if (v <= -5e-310) {
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-310)) 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-310) {
tmp = Math.atan(-1.0);
} else {
tmp = Math.atan(1.0);
}
return tmp;
}
def code(v, H): tmp = 0 if v <= -5e-310: tmp = math.atan(-1.0) else: tmp = math.atan(1.0) return tmp
function code(v, H) tmp = 0.0 if (v <= -5e-310) 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-310) tmp = atan(-1.0); else tmp = atan(1.0); end tmp_2 = tmp; end
code[v_, H_] := If[LessEqual[v, -5e-310], N[ArcTan[-1.0], $MachinePrecision], N[ArcTan[1.0], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;v \leq -5 \cdot 10^{-310}:\\
\;\;\;\;\tan^{-1} -1\\
\mathbf{else}:\\
\;\;\;\;\tan^{-1} 1\\
\end{array}
\end{array}
if v < -4.999999999999985e-310Initial program 70.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-eval70.8%
Simplified70.8%
Taylor expanded in v around -inf
Simplified67.8%
if -4.999999999999985e-310 < v Initial program 71.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-eval71.3%
Simplified71.3%
Taylor expanded in v around inf
Simplified66.3%
(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 71.0%
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-eval71.0%
Simplified71.0%
Taylor expanded in v around -inf
Simplified37.1%
herbie shell --seed 2024161
(FPCore (v H)
:name "Optimal throwing angle"
:precision binary64
(atan (/ v (sqrt (- (* v v) (* (* 2.0 9.8) H))))))