
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (/ (* (exp (- (/ (* sinTheta_i sinTheta_O) v))) (/ (* cosTheta_i cosTheta_O) v)) (* (* (sinh (/ 1.0 v)) 2.0) v)))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return (expf(-((sinTheta_i * sinTheta_O) / v)) * ((cosTheta_i * cosTheta_O) / v)) / ((sinhf((1.0f / v)) * 2.0f) * v);
}
real(4) function code(costheta_i, costheta_o, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = (exp(-((sintheta_i * sintheta_o) / v)) * ((costheta_i * costheta_o) / v)) / ((sinh((1.0e0 / v)) * 2.0e0) * v)
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(exp(Float32(-Float32(Float32(sinTheta_i * sinTheta_O) / v))) * Float32(Float32(cosTheta_i * cosTheta_O) / v)) / Float32(Float32(sinh(Float32(Float32(1.0) / v)) * Float32(2.0)) * v)) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = (exp(-((sinTheta_i * sinTheta_O) / v)) * ((cosTheta_i * cosTheta_O) / v)) / ((sinh((single(1.0) / v)) * single(2.0)) * v); end
\begin{array}{l}
\\
\frac{e^{-\frac{sinTheta_i \cdot sinTheta_O}{v}} \cdot \frac{cosTheta_i \cdot cosTheta_O}{v}}{\left(\sinh \left(\frac{1}{v}\right) \cdot 2\right) \cdot v}
\end{array}
Sampling outcomes in binary32 precision:
Herbie found 18 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (/ (* (exp (- (/ (* sinTheta_i sinTheta_O) v))) (/ (* cosTheta_i cosTheta_O) v)) (* (* (sinh (/ 1.0 v)) 2.0) v)))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return (expf(-((sinTheta_i * sinTheta_O) / v)) * ((cosTheta_i * cosTheta_O) / v)) / ((sinhf((1.0f / v)) * 2.0f) * v);
}
real(4) function code(costheta_i, costheta_o, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = (exp(-((sintheta_i * sintheta_o) / v)) * ((costheta_i * costheta_o) / v)) / ((sinh((1.0e0 / v)) * 2.0e0) * v)
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(exp(Float32(-Float32(Float32(sinTheta_i * sinTheta_O) / v))) * Float32(Float32(cosTheta_i * cosTheta_O) / v)) / Float32(Float32(sinh(Float32(Float32(1.0) / v)) * Float32(2.0)) * v)) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = (exp(-((sinTheta_i * sinTheta_O) / v)) * ((cosTheta_i * cosTheta_O) / v)) / ((sinh((single(1.0) / v)) * single(2.0)) * v); end
\begin{array}{l}
\\
\frac{e^{-\frac{sinTheta_i \cdot sinTheta_O}{v}} \cdot \frac{cosTheta_i \cdot cosTheta_O}{v}}{\left(\sinh \left(\frac{1}{v}\right) \cdot 2\right) \cdot v}
\end{array}
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* (* cosTheta_O (* cosTheta_i (/ 1.0 v))) (* (exp (/ sinTheta_O (/ v (- sinTheta_i)))) (* (/ 1.0 v) (/ 0.5 (sinh (/ 1.0 v)))))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return (cosTheta_O * (cosTheta_i * (1.0f / v))) * (expf((sinTheta_O / (v / -sinTheta_i))) * ((1.0f / v) * (0.5f / sinhf((1.0f / v)))));
}
real(4) function code(costheta_i, costheta_o, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = (costheta_o * (costheta_i * (1.0e0 / v))) * (exp((sintheta_o / (v / -sintheta_i))) * ((1.0e0 / v) * (0.5e0 / sinh((1.0e0 / v)))))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(cosTheta_O * Float32(cosTheta_i * Float32(Float32(1.0) / v))) * Float32(exp(Float32(sinTheta_O / Float32(v / Float32(-sinTheta_i)))) * Float32(Float32(Float32(1.0) / v) * Float32(Float32(0.5) / sinh(Float32(Float32(1.0) / v)))))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = (cosTheta_O * (cosTheta_i * (single(1.0) / v))) * (exp((sinTheta_O / (v / -sinTheta_i))) * ((single(1.0) / v) * (single(0.5) / sinh((single(1.0) / v))))); end
\begin{array}{l}
\\
\left(cosTheta_O \cdot \left(cosTheta_i \cdot \frac{1}{v}\right)\right) \cdot \left(e^{\frac{sinTheta_O}{\frac{v}{-sinTheta_i}}} \cdot \left(\frac{1}{v} \cdot \frac{0.5}{\sinh \left(\frac{1}{v}\right)}\right)\right)
\end{array}
Initial program 98.4%
distribute-neg-frac98.4%
*-commutative98.4%
distribute-rgt-neg-in98.4%
associate-*l/98.4%
*-commutative98.4%
Simplified98.4%
div-inv98.4%
associate-*l/98.5%
associate-*r/98.5%
*-commutative98.5%
Applied egg-rr98.5%
associate-*l*98.5%
*-commutative98.5%
associate-*l/98.5%
associate-*r/98.4%
associate-/l*98.4%
associate-/r*98.7%
*-commutative98.7%
Simplified98.7%
div-inv98.8%
Applied egg-rr98.8%
div-inv98.8%
Applied egg-rr98.8%
expm1-log1p-u98.7%
expm1-udef94.0%
associate-/r*94.0%
metadata-eval94.0%
Applied egg-rr94.0%
expm1-def98.7%
expm1-log1p98.8%
Simplified98.8%
Final simplification98.8%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* (* cosTheta_O (* cosTheta_i (/ 1.0 v))) (* (exp (/ sinTheta_O (/ v (- sinTheta_i)))) (/ (/ 0.5 v) (sinh (/ 1.0 v))))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return (cosTheta_O * (cosTheta_i * (1.0f / v))) * (expf((sinTheta_O / (v / -sinTheta_i))) * ((0.5f / v) / sinhf((1.0f / v))));
}
real(4) function code(costheta_i, costheta_o, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = (costheta_o * (costheta_i * (1.0e0 / v))) * (exp((sintheta_o / (v / -sintheta_i))) * ((0.5e0 / v) / sinh((1.0e0 / v))))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(cosTheta_O * Float32(cosTheta_i * Float32(Float32(1.0) / v))) * Float32(exp(Float32(sinTheta_O / Float32(v / Float32(-sinTheta_i)))) * Float32(Float32(Float32(0.5) / v) / sinh(Float32(Float32(1.0) / v))))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = (cosTheta_O * (cosTheta_i * (single(1.0) / v))) * (exp((sinTheta_O / (v / -sinTheta_i))) * ((single(0.5) / v) / sinh((single(1.0) / v)))); end
\begin{array}{l}
\\
\left(cosTheta_O \cdot \left(cosTheta_i \cdot \frac{1}{v}\right)\right) \cdot \left(e^{\frac{sinTheta_O}{\frac{v}{-sinTheta_i}}} \cdot \frac{\frac{0.5}{v}}{\sinh \left(\frac{1}{v}\right)}\right)
\end{array}
Initial program 98.4%
distribute-neg-frac98.4%
*-commutative98.4%
distribute-rgt-neg-in98.4%
associate-*l/98.4%
*-commutative98.4%
Simplified98.4%
div-inv98.4%
associate-*l/98.5%
associate-*r/98.5%
*-commutative98.5%
Applied egg-rr98.5%
associate-*l*98.5%
*-commutative98.5%
associate-*l/98.5%
associate-*r/98.4%
associate-/l*98.4%
associate-/r*98.7%
*-commutative98.7%
Simplified98.7%
div-inv98.8%
Applied egg-rr98.8%
*-un-lft-identity98.7%
Applied egg-rr98.8%
associate-*r/98.7%
times-frac98.7%
metadata-eval98.7%
associate-/r*98.4%
associate-*r/98.4%
metadata-eval98.4%
associate-/r*98.7%
Simplified98.8%
Final simplification98.8%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* (* (exp (/ sinTheta_O (/ v (- sinTheta_i)))) (/ (/ 0.5 v) (sinh (/ 1.0 v)))) (* cosTheta_O (/ cosTheta_i v))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return (expf((sinTheta_O / (v / -sinTheta_i))) * ((0.5f / v) / sinhf((1.0f / v)))) * (cosTheta_O * (cosTheta_i / v));
}
real(4) function code(costheta_i, costheta_o, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = (exp((sintheta_o / (v / -sintheta_i))) * ((0.5e0 / v) / sinh((1.0e0 / v)))) * (costheta_o * (costheta_i / v))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(exp(Float32(sinTheta_O / Float32(v / Float32(-sinTheta_i)))) * Float32(Float32(Float32(0.5) / v) / sinh(Float32(Float32(1.0) / v)))) * Float32(cosTheta_O * Float32(cosTheta_i / v))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = (exp((sinTheta_O / (v / -sinTheta_i))) * ((single(0.5) / v) / sinh((single(1.0) / v)))) * (cosTheta_O * (cosTheta_i / v)); end
\begin{array}{l}
\\
\left(e^{\frac{sinTheta_O}{\frac{v}{-sinTheta_i}}} \cdot \frac{\frac{0.5}{v}}{\sinh \left(\frac{1}{v}\right)}\right) \cdot \left(cosTheta_O \cdot \frac{cosTheta_i}{v}\right)
\end{array}
Initial program 98.4%
distribute-neg-frac98.4%
*-commutative98.4%
distribute-rgt-neg-in98.4%
associate-*l/98.4%
*-commutative98.4%
Simplified98.4%
div-inv98.4%
associate-*l/98.5%
associate-*r/98.5%
*-commutative98.5%
Applied egg-rr98.5%
associate-*l*98.5%
*-commutative98.5%
associate-*l/98.5%
associate-*r/98.4%
associate-/l*98.4%
associate-/r*98.7%
*-commutative98.7%
Simplified98.7%
*-un-lft-identity98.7%
Applied egg-rr98.7%
associate-*r/98.7%
times-frac98.7%
metadata-eval98.7%
associate-/r*98.4%
associate-*r/98.4%
metadata-eval98.4%
associate-/r*98.7%
Simplified98.7%
Final simplification98.7%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (/ (* cosTheta_O cosTheta_i) (* (* (sinh (/ 1.0 v)) (* v (* v 2.0))) (exp (* sinTheta_i (/ sinTheta_O v))))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return (cosTheta_O * cosTheta_i) / ((sinhf((1.0f / v)) * (v * (v * 2.0f))) * expf((sinTheta_i * (sinTheta_O / v))));
}
real(4) function code(costheta_i, costheta_o, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = (costheta_o * costheta_i) / ((sinh((1.0e0 / v)) * (v * (v * 2.0e0))) * exp((sintheta_i * (sintheta_o / v))))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(cosTheta_O * cosTheta_i) / Float32(Float32(sinh(Float32(Float32(1.0) / v)) * Float32(v * Float32(v * Float32(2.0)))) * exp(Float32(sinTheta_i * Float32(sinTheta_O / v))))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = (cosTheta_O * cosTheta_i) / ((sinh((single(1.0) / v)) * (v * (v * single(2.0)))) * exp((sinTheta_i * (sinTheta_O / v)))); end
\begin{array}{l}
\\
\frac{cosTheta_O \cdot cosTheta_i}{\left(\sinh \left(\frac{1}{v}\right) \cdot \left(v \cdot \left(v \cdot 2\right)\right)\right) \cdot e^{sinTheta_i \cdot \frac{sinTheta_O}{v}}}
\end{array}
Initial program 98.4%
associate-*l/98.5%
times-frac98.4%
exp-neg98.4%
associate-*l/98.4%
*-lft-identity98.4%
associate-/l/98.4%
associate-*l*98.4%
associate-*l*98.5%
*-commutative98.5%
*-commutative98.5%
associate-*l/98.5%
Simplified98.5%
Final simplification98.5%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (/ (* cosTheta_O (/ cosTheta_i (* (* (sinh (/ 1.0 v)) (* v v)) 2.0))) (exp (/ sinTheta_i (/ v sinTheta_O)))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return (cosTheta_O * (cosTheta_i / ((sinhf((1.0f / v)) * (v * v)) * 2.0f))) / expf((sinTheta_i / (v / sinTheta_O)));
}
real(4) function code(costheta_i, costheta_o, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = (costheta_o * (costheta_i / ((sinh((1.0e0 / v)) * (v * v)) * 2.0e0))) / exp((sintheta_i / (v / sintheta_o)))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(cosTheta_O * Float32(cosTheta_i / Float32(Float32(sinh(Float32(Float32(1.0) / v)) * Float32(v * v)) * Float32(2.0)))) / exp(Float32(sinTheta_i / Float32(v / sinTheta_O)))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = (cosTheta_O * (cosTheta_i / ((sinh((single(1.0) / v)) * (v * v)) * single(2.0)))) / exp((sinTheta_i / (v / sinTheta_O))); end
\begin{array}{l}
\\
\frac{cosTheta_O \cdot \frac{cosTheta_i}{\left(\sinh \left(\frac{1}{v}\right) \cdot \left(v \cdot v\right)\right) \cdot 2}}{e^{\frac{sinTheta_i}{\frac{v}{sinTheta_O}}}}
\end{array}
Initial program 98.4%
associate-*l/98.5%
times-frac98.4%
exp-neg98.4%
associate-*l/98.4%
*-lft-identity98.4%
associate-/l/98.4%
associate-*l*98.4%
associate-*l*98.5%
*-commutative98.5%
*-commutative98.5%
associate-*l/98.5%
Simplified98.5%
expm1-log1p-u98.5%
expm1-udef48.5%
times-frac48.5%
*-commutative48.5%
*-commutative48.5%
Applied egg-rr48.5%
expm1-def98.5%
expm1-log1p98.5%
associate-*r/98.5%
associate-*r*98.5%
unpow298.5%
associate-*r*98.5%
unpow298.5%
associate-*r/98.5%
associate-/l*98.5%
Simplified98.5%
Final simplification98.5%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* (* cosTheta_i (/ cosTheta_O v)) (/ (/ 1.0 v) (- (exp (/ 1.0 v)) (exp (/ -1.0 v))))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return (cosTheta_i * (cosTheta_O / v)) * ((1.0f / v) / (expf((1.0f / v)) - expf((-1.0f / v))));
}
real(4) function code(costheta_i, costheta_o, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = (costheta_i * (costheta_o / v)) * ((1.0e0 / v) / (exp((1.0e0 / v)) - exp(((-1.0e0) / v))))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(cosTheta_i * Float32(cosTheta_O / v)) * Float32(Float32(Float32(1.0) / v) / Float32(exp(Float32(Float32(1.0) / v)) - exp(Float32(Float32(-1.0) / v))))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = (cosTheta_i * (cosTheta_O / v)) * ((single(1.0) / v) / (exp((single(1.0) / v)) - exp((single(-1.0) / v)))); end
\begin{array}{l}
\\
\left(cosTheta_i \cdot \frac{cosTheta_O}{v}\right) \cdot \frac{\frac{1}{v}}{e^{\frac{1}{v}} - e^{\frac{-1}{v}}}
\end{array}
Initial program 98.4%
*-commutative98.4%
associate-*r/98.5%
*-commutative98.5%
associate-*l/98.5%
*-commutative98.5%
*-commutative98.5%
associate-*r*98.5%
associate-/l/98.5%
exp-neg98.5%
associate-/l/98.5%
associate-/r*98.5%
metadata-eval98.5%
associate-*l/98.5%
*-commutative98.5%
exp-prod98.5%
Simplified98.5%
Taylor expanded in sinTheta_O around 0 98.2%
rec-exp98.2%
distribute-neg-frac98.2%
metadata-eval98.2%
Simplified98.2%
expm1-log1p-u97.8%
Applied egg-rr97.8%
Taylor expanded in v around 0 98.2%
associate-/r*98.4%
Simplified98.4%
Final simplification98.4%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* (/ cosTheta_i (* v v)) (/ cosTheta_O (- (exp (/ 1.0 v)) (exp (/ -1.0 v))))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return (cosTheta_i / (v * v)) * (cosTheta_O / (expf((1.0f / v)) - expf((-1.0f / v))));
}
real(4) function code(costheta_i, costheta_o, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = (costheta_i / (v * v)) * (costheta_o / (exp((1.0e0 / v)) - exp(((-1.0e0) / v))))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(cosTheta_i / Float32(v * v)) * Float32(cosTheta_O / Float32(exp(Float32(Float32(1.0) / v)) - exp(Float32(Float32(-1.0) / v))))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = (cosTheta_i / (v * v)) * (cosTheta_O / (exp((single(1.0) / v)) - exp((single(-1.0) / v)))); end
\begin{array}{l}
\\
\frac{cosTheta_i}{v \cdot v} \cdot \frac{cosTheta_O}{e^{\frac{1}{v}} - e^{\frac{-1}{v}}}
\end{array}
Initial program 98.4%
*-commutative98.4%
associate-*r/98.5%
*-commutative98.5%
associate-*l/98.5%
*-commutative98.5%
*-commutative98.5%
associate-*r*98.5%
associate-/l/98.5%
exp-neg98.5%
associate-/l/98.5%
associate-/r*98.5%
metadata-eval98.5%
associate-*l/98.5%
*-commutative98.5%
exp-prod98.5%
Simplified98.5%
Taylor expanded in sinTheta_O around 0 98.2%
times-frac98.4%
unpow298.4%
rec-exp98.4%
distribute-neg-frac98.4%
metadata-eval98.4%
Simplified98.4%
Final simplification98.4%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v)
:precision binary32
(let* ((t_0 (* cosTheta_i (/ cosTheta_O v))))
(if (<= v 0.5099999904632568)
(* t_0 (/ 1.0 (* v (+ (exp (/ 1.0 v)) -1.0))))
(+
(* t_0 (/ 0.009722222222222222 (pow v 4.0)))
(* t_0 (+ 0.5 (/ -0.08333333333333333 (* v v))))))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
float t_0 = cosTheta_i * (cosTheta_O / v);
float tmp;
if (v <= 0.5099999904632568f) {
tmp = t_0 * (1.0f / (v * (expf((1.0f / v)) + -1.0f)));
} else {
tmp = (t_0 * (0.009722222222222222f / powf(v, 4.0f))) + (t_0 * (0.5f + (-0.08333333333333333f / (v * v))));
}
return tmp;
}
real(4) function code(costheta_i, costheta_o, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
real(4) :: t_0
real(4) :: tmp
t_0 = costheta_i * (costheta_o / v)
if (v <= 0.5099999904632568e0) then
tmp = t_0 * (1.0e0 / (v * (exp((1.0e0 / v)) + (-1.0e0))))
else
tmp = (t_0 * (0.009722222222222222e0 / (v ** 4.0e0))) + (t_0 * (0.5e0 + ((-0.08333333333333333e0) / (v * v))))
end if
code = tmp
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) t_0 = Float32(cosTheta_i * Float32(cosTheta_O / v)) tmp = Float32(0.0) if (v <= Float32(0.5099999904632568)) tmp = Float32(t_0 * Float32(Float32(1.0) / Float32(v * Float32(exp(Float32(Float32(1.0) / v)) + Float32(-1.0))))); else tmp = Float32(Float32(t_0 * Float32(Float32(0.009722222222222222) / (v ^ Float32(4.0)))) + Float32(t_0 * Float32(Float32(0.5) + Float32(Float32(-0.08333333333333333) / Float32(v * v))))); end return tmp end
function tmp_2 = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) t_0 = cosTheta_i * (cosTheta_O / v); tmp = single(0.0); if (v <= single(0.5099999904632568)) tmp = t_0 * (single(1.0) / (v * (exp((single(1.0) / v)) + single(-1.0)))); else tmp = (t_0 * (single(0.009722222222222222) / (v ^ single(4.0)))) + (t_0 * (single(0.5) + (single(-0.08333333333333333) / (v * v)))); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := cosTheta_i \cdot \frac{cosTheta_O}{v}\\
\mathbf{if}\;v \leq 0.5099999904632568:\\
\;\;\;\;t_0 \cdot \frac{1}{v \cdot \left(e^{\frac{1}{v}} + -1\right)}\\
\mathbf{else}:\\
\;\;\;\;t_0 \cdot \frac{0.009722222222222222}{{v}^{4}} + t_0 \cdot \left(0.5 + \frac{-0.08333333333333333}{v \cdot v}\right)\\
\end{array}
\end{array}
if v < 0.50999999Initial program 98.1%
*-commutative98.1%
associate-*r/98.1%
*-commutative98.1%
associate-*l/98.1%
*-commutative98.1%
*-commutative98.1%
associate-*r*98.1%
associate-/l/98.1%
exp-neg98.1%
associate-/l/98.1%
associate-/r*98.1%
metadata-eval98.1%
associate-*l/98.1%
*-commutative98.1%
exp-prod98.1%
Simplified98.1%
Taylor expanded in sinTheta_O around 0 97.9%
rec-exp97.9%
distribute-neg-frac97.9%
metadata-eval97.9%
Simplified97.9%
Taylor expanded in v around inf 68.1%
if 0.50999999 < v Initial program 98.8%
*-commutative98.8%
associate-*r/99.0%
*-commutative99.0%
associate-*l/99.1%
*-commutative99.1%
*-commutative99.1%
associate-*r*99.1%
associate-/l/99.1%
exp-neg99.1%
associate-/l/99.1%
associate-/r*99.1%
metadata-eval99.1%
associate-*l/99.1%
*-commutative99.1%
exp-prod99.1%
Simplified99.1%
Taylor expanded in sinTheta_O around 0 98.7%
rec-exp98.7%
distribute-neg-frac98.7%
metadata-eval98.7%
Simplified98.7%
expm1-log1p-u98.2%
Applied egg-rr98.2%
Taylor expanded in v around inf 71.0%
+-commutative71.0%
associate--l+71.0%
associate-*r/71.0%
metadata-eval71.0%
sub-neg71.0%
associate-*r/71.0%
metadata-eval71.0%
distribute-neg-frac71.0%
metadata-eval71.0%
unpow271.0%
Simplified71.0%
distribute-rgt-in71.1%
Applied egg-rr71.1%
Final simplification69.3%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v)
:precision binary32
(let* ((t_0 (* cosTheta_i (/ cosTheta_O v))))
(if (<= v 0.5099999904632568)
(* t_0 (/ 1.0 (* v (+ (exp (/ 1.0 v)) -1.0))))
(*
t_0
(+
(/ 0.009722222222222222 (pow v 4.0))
(+ 0.5 (/ -0.08333333333333333 (* v v))))))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
float t_0 = cosTheta_i * (cosTheta_O / v);
float tmp;
if (v <= 0.5099999904632568f) {
tmp = t_0 * (1.0f / (v * (expf((1.0f / v)) + -1.0f)));
} else {
tmp = t_0 * ((0.009722222222222222f / powf(v, 4.0f)) + (0.5f + (-0.08333333333333333f / (v * v))));
}
return tmp;
}
real(4) function code(costheta_i, costheta_o, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
real(4) :: t_0
real(4) :: tmp
t_0 = costheta_i * (costheta_o / v)
if (v <= 0.5099999904632568e0) then
tmp = t_0 * (1.0e0 / (v * (exp((1.0e0 / v)) + (-1.0e0))))
else
tmp = t_0 * ((0.009722222222222222e0 / (v ** 4.0e0)) + (0.5e0 + ((-0.08333333333333333e0) / (v * v))))
end if
code = tmp
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) t_0 = Float32(cosTheta_i * Float32(cosTheta_O / v)) tmp = Float32(0.0) if (v <= Float32(0.5099999904632568)) tmp = Float32(t_0 * Float32(Float32(1.0) / Float32(v * Float32(exp(Float32(Float32(1.0) / v)) + Float32(-1.0))))); else tmp = Float32(t_0 * Float32(Float32(Float32(0.009722222222222222) / (v ^ Float32(4.0))) + Float32(Float32(0.5) + Float32(Float32(-0.08333333333333333) / Float32(v * v))))); end return tmp end
function tmp_2 = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) t_0 = cosTheta_i * (cosTheta_O / v); tmp = single(0.0); if (v <= single(0.5099999904632568)) tmp = t_0 * (single(1.0) / (v * (exp((single(1.0) / v)) + single(-1.0)))); else tmp = t_0 * ((single(0.009722222222222222) / (v ^ single(4.0))) + (single(0.5) + (single(-0.08333333333333333) / (v * v)))); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := cosTheta_i \cdot \frac{cosTheta_O}{v}\\
\mathbf{if}\;v \leq 0.5099999904632568:\\
\;\;\;\;t_0 \cdot \frac{1}{v \cdot \left(e^{\frac{1}{v}} + -1\right)}\\
\mathbf{else}:\\
\;\;\;\;t_0 \cdot \left(\frac{0.009722222222222222}{{v}^{4}} + \left(0.5 + \frac{-0.08333333333333333}{v \cdot v}\right)\right)\\
\end{array}
\end{array}
if v < 0.50999999Initial program 98.1%
*-commutative98.1%
associate-*r/98.1%
*-commutative98.1%
associate-*l/98.1%
*-commutative98.1%
*-commutative98.1%
associate-*r*98.1%
associate-/l/98.1%
exp-neg98.1%
associate-/l/98.1%
associate-/r*98.1%
metadata-eval98.1%
associate-*l/98.1%
*-commutative98.1%
exp-prod98.1%
Simplified98.1%
Taylor expanded in sinTheta_O around 0 97.9%
rec-exp97.9%
distribute-neg-frac97.9%
metadata-eval97.9%
Simplified97.9%
Taylor expanded in v around inf 68.1%
if 0.50999999 < v Initial program 98.8%
*-commutative98.8%
associate-*r/99.0%
*-commutative99.0%
associate-*l/99.1%
*-commutative99.1%
*-commutative99.1%
associate-*r*99.1%
associate-/l/99.1%
exp-neg99.1%
associate-/l/99.1%
associate-/r*99.1%
metadata-eval99.1%
associate-*l/99.1%
*-commutative99.1%
exp-prod99.1%
Simplified99.1%
Taylor expanded in sinTheta_O around 0 98.7%
rec-exp98.7%
distribute-neg-frac98.7%
metadata-eval98.7%
Simplified98.7%
expm1-log1p-u98.2%
Applied egg-rr98.2%
Taylor expanded in v around inf 71.0%
+-commutative71.0%
associate--l+71.0%
associate-*r/71.0%
metadata-eval71.0%
sub-neg71.0%
associate-*r/71.0%
metadata-eval71.0%
distribute-neg-frac71.0%
metadata-eval71.0%
unpow271.0%
Simplified71.0%
Final simplification69.2%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v)
:precision binary32
(let* ((t_0 (* cosTheta_i (/ cosTheta_O v))))
(if (<= v 0.5099999904632568)
(* t_0 (/ 1.0 (* v (+ (exp (/ 1.0 v)) -1.0))))
(*
t_0
(-
(+ 0.5 (/ 0.009722222222222222 (pow v 4.0)))
(/ 0.08333333333333333 (* v v)))))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
float t_0 = cosTheta_i * (cosTheta_O / v);
float tmp;
if (v <= 0.5099999904632568f) {
tmp = t_0 * (1.0f / (v * (expf((1.0f / v)) + -1.0f)));
} else {
tmp = t_0 * ((0.5f + (0.009722222222222222f / powf(v, 4.0f))) - (0.08333333333333333f / (v * v)));
}
return tmp;
}
real(4) function code(costheta_i, costheta_o, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
real(4) :: t_0
real(4) :: tmp
t_0 = costheta_i * (costheta_o / v)
if (v <= 0.5099999904632568e0) then
tmp = t_0 * (1.0e0 / (v * (exp((1.0e0 / v)) + (-1.0e0))))
else
tmp = t_0 * ((0.5e0 + (0.009722222222222222e0 / (v ** 4.0e0))) - (0.08333333333333333e0 / (v * v)))
end if
code = tmp
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) t_0 = Float32(cosTheta_i * Float32(cosTheta_O / v)) tmp = Float32(0.0) if (v <= Float32(0.5099999904632568)) tmp = Float32(t_0 * Float32(Float32(1.0) / Float32(v * Float32(exp(Float32(Float32(1.0) / v)) + Float32(-1.0))))); else tmp = Float32(t_0 * Float32(Float32(Float32(0.5) + Float32(Float32(0.009722222222222222) / (v ^ Float32(4.0)))) - Float32(Float32(0.08333333333333333) / Float32(v * v)))); end return tmp end
function tmp_2 = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) t_0 = cosTheta_i * (cosTheta_O / v); tmp = single(0.0); if (v <= single(0.5099999904632568)) tmp = t_0 * (single(1.0) / (v * (exp((single(1.0) / v)) + single(-1.0)))); else tmp = t_0 * ((single(0.5) + (single(0.009722222222222222) / (v ^ single(4.0)))) - (single(0.08333333333333333) / (v * v))); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := cosTheta_i \cdot \frac{cosTheta_O}{v}\\
\mathbf{if}\;v \leq 0.5099999904632568:\\
\;\;\;\;t_0 \cdot \frac{1}{v \cdot \left(e^{\frac{1}{v}} + -1\right)}\\
\mathbf{else}:\\
\;\;\;\;t_0 \cdot \left(\left(0.5 + \frac{0.009722222222222222}{{v}^{4}}\right) - \frac{0.08333333333333333}{v \cdot v}\right)\\
\end{array}
\end{array}
if v < 0.50999999Initial program 98.1%
*-commutative98.1%
associate-*r/98.1%
*-commutative98.1%
associate-*l/98.1%
*-commutative98.1%
*-commutative98.1%
associate-*r*98.1%
associate-/l/98.1%
exp-neg98.1%
associate-/l/98.1%
associate-/r*98.1%
metadata-eval98.1%
associate-*l/98.1%
*-commutative98.1%
exp-prod98.1%
Simplified98.1%
Taylor expanded in sinTheta_O around 0 97.9%
rec-exp97.9%
distribute-neg-frac97.9%
metadata-eval97.9%
Simplified97.9%
Taylor expanded in v around inf 68.1%
if 0.50999999 < v Initial program 98.8%
*-commutative98.8%
associate-*r/99.0%
*-commutative99.0%
associate-*l/99.1%
*-commutative99.1%
*-commutative99.1%
associate-*r*99.1%
associate-/l/99.1%
exp-neg99.1%
associate-/l/99.1%
associate-/r*99.1%
metadata-eval99.1%
associate-*l/99.1%
*-commutative99.1%
exp-prod99.1%
Simplified99.1%
Taylor expanded in sinTheta_O around 0 98.7%
rec-exp98.7%
distribute-neg-frac98.7%
metadata-eval98.7%
Simplified98.7%
expm1-log1p-u98.2%
Applied egg-rr98.2%
Taylor expanded in v around inf 71.0%
associate-*r/71.0%
metadata-eval71.0%
associate-*r/71.0%
metadata-eval71.0%
unpow271.0%
Simplified71.0%
Final simplification69.2%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v)
:precision binary32
(if (<= v 0.4650000035762787)
(* (* cosTheta_i (/ cosTheta_O v)) (/ 1.0 (* v (+ (exp (/ 1.0 v)) -1.0))))
(/
(* cosTheta_O (/ cosTheta_i (* 2.0 (+ v (/ 0.16666666666666666 v)))))
(exp (/ sinTheta_i (/ v sinTheta_O))))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
float tmp;
if (v <= 0.4650000035762787f) {
tmp = (cosTheta_i * (cosTheta_O / v)) * (1.0f / (v * (expf((1.0f / v)) + -1.0f)));
} else {
tmp = (cosTheta_O * (cosTheta_i / (2.0f * (v + (0.16666666666666666f / v))))) / expf((sinTheta_i / (v / sinTheta_O)));
}
return tmp;
}
real(4) function code(costheta_i, costheta_o, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
real(4) :: tmp
if (v <= 0.4650000035762787e0) then
tmp = (costheta_i * (costheta_o / v)) * (1.0e0 / (v * (exp((1.0e0 / v)) + (-1.0e0))))
else
tmp = (costheta_o * (costheta_i / (2.0e0 * (v + (0.16666666666666666e0 / v))))) / exp((sintheta_i / (v / sintheta_o)))
end if
code = tmp
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = Float32(0.0) if (v <= Float32(0.4650000035762787)) tmp = Float32(Float32(cosTheta_i * Float32(cosTheta_O / v)) * Float32(Float32(1.0) / Float32(v * Float32(exp(Float32(Float32(1.0) / v)) + Float32(-1.0))))); else tmp = Float32(Float32(cosTheta_O * Float32(cosTheta_i / Float32(Float32(2.0) * Float32(v + Float32(Float32(0.16666666666666666) / v))))) / exp(Float32(sinTheta_i / Float32(v / sinTheta_O)))); end return tmp end
function tmp_2 = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = single(0.0); if (v <= single(0.4650000035762787)) tmp = (cosTheta_i * (cosTheta_O / v)) * (single(1.0) / (v * (exp((single(1.0) / v)) + single(-1.0)))); else tmp = (cosTheta_O * (cosTheta_i / (single(2.0) * (v + (single(0.16666666666666666) / v))))) / exp((sinTheta_i / (v / sinTheta_O))); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;v \leq 0.4650000035762787:\\
\;\;\;\;\left(cosTheta_i \cdot \frac{cosTheta_O}{v}\right) \cdot \frac{1}{v \cdot \left(e^{\frac{1}{v}} + -1\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{cosTheta_O \cdot \frac{cosTheta_i}{2 \cdot \left(v + \frac{0.16666666666666666}{v}\right)}}{e^{\frac{sinTheta_i}{\frac{v}{sinTheta_O}}}}\\
\end{array}
\end{array}
if v < 0.465000004Initial program 98.1%
*-commutative98.1%
associate-*r/98.1%
*-commutative98.1%
associate-*l/98.1%
*-commutative98.1%
*-commutative98.1%
associate-*r*98.1%
associate-/l/98.1%
exp-neg98.1%
associate-/l/98.1%
associate-/r*98.1%
metadata-eval98.1%
associate-*l/98.1%
*-commutative98.1%
exp-prod98.1%
Simplified98.1%
Taylor expanded in sinTheta_O around 0 98.0%
rec-exp98.0%
distribute-neg-frac98.0%
metadata-eval98.0%
Simplified98.0%
Taylor expanded in v around inf 69.4%
if 0.465000004 < v Initial program 98.7%
associate-*l/98.9%
times-frac98.8%
exp-neg98.8%
associate-*l/98.8%
*-lft-identity98.8%
associate-/l/98.8%
associate-*l*98.8%
associate-*l*98.8%
*-commutative98.8%
*-commutative98.8%
associate-*l/98.8%
Simplified98.8%
expm1-log1p-u98.8%
expm1-udef48.3%
times-frac48.3%
*-commutative48.3%
*-commutative48.3%
Applied egg-rr48.3%
expm1-def98.9%
expm1-log1p98.9%
associate-*r/98.9%
associate-*r*98.9%
unpow298.9%
associate-*r*98.9%
unpow298.9%
associate-*r/98.9%
associate-/l*98.9%
Simplified98.9%
Taylor expanded in v around inf 67.0%
associate-*r/67.0%
metadata-eval67.0%
Simplified67.0%
Final simplification68.4%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v)
:precision binary32
(let* ((t_0 (+ (* v 2.0) (/ 0.3333333333333333 v))))
(if (<= v 0.4650000035762787)
(* (* cosTheta_i (/ cosTheta_O v)) (/ 1.0 (* v (+ (exp (/ 1.0 v)) -1.0))))
(-
(/ cosTheta_i (/ t_0 cosTheta_O))
(/ (* sinTheta_i (* cosTheta_i (* cosTheta_O sinTheta_O))) (* v t_0))))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
float t_0 = (v * 2.0f) + (0.3333333333333333f / v);
float tmp;
if (v <= 0.4650000035762787f) {
tmp = (cosTheta_i * (cosTheta_O / v)) * (1.0f / (v * (expf((1.0f / v)) + -1.0f)));
} else {
tmp = (cosTheta_i / (t_0 / cosTheta_O)) - ((sinTheta_i * (cosTheta_i * (cosTheta_O * sinTheta_O))) / (v * t_0));
}
return tmp;
}
real(4) function code(costheta_i, costheta_o, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
real(4) :: t_0
real(4) :: tmp
t_0 = (v * 2.0e0) + (0.3333333333333333e0 / v)
if (v <= 0.4650000035762787e0) then
tmp = (costheta_i * (costheta_o / v)) * (1.0e0 / (v * (exp((1.0e0 / v)) + (-1.0e0))))
else
tmp = (costheta_i / (t_0 / costheta_o)) - ((sintheta_i * (costheta_i * (costheta_o * sintheta_o))) / (v * t_0))
end if
code = tmp
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) t_0 = Float32(Float32(v * Float32(2.0)) + Float32(Float32(0.3333333333333333) / v)) tmp = Float32(0.0) if (v <= Float32(0.4650000035762787)) tmp = Float32(Float32(cosTheta_i * Float32(cosTheta_O / v)) * Float32(Float32(1.0) / Float32(v * Float32(exp(Float32(Float32(1.0) / v)) + Float32(-1.0))))); else tmp = Float32(Float32(cosTheta_i / Float32(t_0 / cosTheta_O)) - Float32(Float32(sinTheta_i * Float32(cosTheta_i * Float32(cosTheta_O * sinTheta_O))) / Float32(v * t_0))); end return tmp end
function tmp_2 = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) t_0 = (v * single(2.0)) + (single(0.3333333333333333) / v); tmp = single(0.0); if (v <= single(0.4650000035762787)) tmp = (cosTheta_i * (cosTheta_O / v)) * (single(1.0) / (v * (exp((single(1.0) / v)) + single(-1.0)))); else tmp = (cosTheta_i / (t_0 / cosTheta_O)) - ((sinTheta_i * (cosTheta_i * (cosTheta_O * sinTheta_O))) / (v * t_0)); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := v \cdot 2 + \frac{0.3333333333333333}{v}\\
\mathbf{if}\;v \leq 0.4650000035762787:\\
\;\;\;\;\left(cosTheta_i \cdot \frac{cosTheta_O}{v}\right) \cdot \frac{1}{v \cdot \left(e^{\frac{1}{v}} + -1\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{cosTheta_i}{\frac{t_0}{cosTheta_O}} - \frac{sinTheta_i \cdot \left(cosTheta_i \cdot \left(cosTheta_O \cdot sinTheta_O\right)\right)}{v \cdot t_0}\\
\end{array}
\end{array}
if v < 0.465000004Initial program 98.1%
*-commutative98.1%
associate-*r/98.1%
*-commutative98.1%
associate-*l/98.1%
*-commutative98.1%
*-commutative98.1%
associate-*r*98.1%
associate-/l/98.1%
exp-neg98.1%
associate-/l/98.1%
associate-/r*98.1%
metadata-eval98.1%
associate-*l/98.1%
*-commutative98.1%
exp-prod98.1%
Simplified98.1%
Taylor expanded in sinTheta_O around 0 98.0%
rec-exp98.0%
distribute-neg-frac98.0%
metadata-eval98.0%
Simplified98.0%
Taylor expanded in v around inf 69.4%
if 0.465000004 < v Initial program 98.7%
associate-*l/98.9%
times-frac98.8%
exp-neg98.8%
associate-*l/98.8%
*-lft-identity98.8%
associate-/l/98.8%
associate-*l*98.8%
associate-*l*98.8%
*-commutative98.8%
*-commutative98.8%
associate-*l/98.8%
Simplified98.8%
Taylor expanded in v around inf 67.0%
Taylor expanded in sinTheta_O around 0 67.0%
associate-/l*67.0%
associate-*r/67.0%
metadata-eval67.0%
*-commutative67.0%
mul-1-neg67.0%
associate-*r/67.0%
metadata-eval67.0%
*-commutative67.0%
Simplified67.0%
Final simplification68.4%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v)
:precision binary32
(let* ((t_0 (+ (* v 2.0) (/ 0.3333333333333333 v))))
(-
(/ cosTheta_i (/ t_0 cosTheta_O))
(/ (* sinTheta_i (* cosTheta_i (* cosTheta_O sinTheta_O))) (* v t_0)))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
float t_0 = (v * 2.0f) + (0.3333333333333333f / v);
return (cosTheta_i / (t_0 / cosTheta_O)) - ((sinTheta_i * (cosTheta_i * (cosTheta_O * sinTheta_O))) / (v * t_0));
}
real(4) function code(costheta_i, costheta_o, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
real(4) :: t_0
t_0 = (v * 2.0e0) + (0.3333333333333333e0 / v)
code = (costheta_i / (t_0 / costheta_o)) - ((sintheta_i * (costheta_i * (costheta_o * sintheta_o))) / (v * t_0))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) t_0 = Float32(Float32(v * Float32(2.0)) + Float32(Float32(0.3333333333333333) / v)) return Float32(Float32(cosTheta_i / Float32(t_0 / cosTheta_O)) - Float32(Float32(sinTheta_i * Float32(cosTheta_i * Float32(cosTheta_O * sinTheta_O))) / Float32(v * t_0))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) t_0 = (v * single(2.0)) + (single(0.3333333333333333) / v); tmp = (cosTheta_i / (t_0 / cosTheta_O)) - ((sinTheta_i * (cosTheta_i * (cosTheta_O * sinTheta_O))) / (v * t_0)); end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := v \cdot 2 + \frac{0.3333333333333333}{v}\\
\frac{cosTheta_i}{\frac{t_0}{cosTheta_O}} - \frac{sinTheta_i \cdot \left(cosTheta_i \cdot \left(cosTheta_O \cdot sinTheta_O\right)\right)}{v \cdot t_0}
\end{array}
\end{array}
Initial program 98.4%
associate-*l/98.5%
times-frac98.4%
exp-neg98.4%
associate-*l/98.4%
*-lft-identity98.4%
associate-/l/98.4%
associate-*l*98.4%
associate-*l*98.5%
*-commutative98.5%
*-commutative98.5%
associate-*l/98.5%
Simplified98.5%
Taylor expanded in v around inf 60.4%
Taylor expanded in sinTheta_O around 0 60.4%
associate-/l*60.4%
associate-*r/60.4%
metadata-eval60.4%
*-commutative60.4%
mul-1-neg60.4%
associate-*r/60.4%
metadata-eval60.4%
*-commutative60.4%
Simplified60.4%
Final simplification60.4%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (/ cosTheta_i (/ (+ (* v 2.0) (/ 0.3333333333333333 v)) cosTheta_O)))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return cosTheta_i / (((v * 2.0f) + (0.3333333333333333f / v)) / cosTheta_O);
}
real(4) function code(costheta_i, costheta_o, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = costheta_i / (((v * 2.0e0) + (0.3333333333333333e0 / v)) / costheta_o)
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(cosTheta_i / Float32(Float32(Float32(v * Float32(2.0)) + Float32(Float32(0.3333333333333333) / v)) / cosTheta_O)) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = cosTheta_i / (((v * single(2.0)) + (single(0.3333333333333333) / v)) / cosTheta_O); end
\begin{array}{l}
\\
\frac{cosTheta_i}{\frac{v \cdot 2 + \frac{0.3333333333333333}{v}}{cosTheta_O}}
\end{array}
Initial program 98.4%
associate-*l/98.5%
times-frac98.4%
exp-neg98.4%
associate-*l/98.4%
*-lft-identity98.4%
associate-/l/98.4%
associate-*l*98.4%
associate-*l*98.5%
*-commutative98.5%
*-commutative98.5%
associate-*l/98.5%
Simplified98.5%
Taylor expanded in v around inf 60.4%
Taylor expanded in sinTheta_O around 0 60.4%
associate-/l*60.4%
associate-*r/60.4%
metadata-eval60.4%
*-commutative60.4%
Simplified60.4%
Final simplification60.4%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* 0.5 (* cosTheta_O (/ cosTheta_i v))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return 0.5f * (cosTheta_O * (cosTheta_i / v));
}
real(4) function code(costheta_i, costheta_o, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = 0.5e0 * (costheta_o * (costheta_i / v))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(0.5) * Float32(cosTheta_O * Float32(cosTheta_i / v))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = single(0.5) * (cosTheta_O * (cosTheta_i / v)); end
\begin{array}{l}
\\
0.5 \cdot \left(cosTheta_O \cdot \frac{cosTheta_i}{v}\right)
\end{array}
Initial program 98.4%
*-commutative98.4%
associate-*r/98.5%
*-commutative98.5%
associate-*l/98.5%
*-commutative98.5%
*-commutative98.5%
associate-*r*98.5%
associate-/l/98.5%
exp-neg98.5%
associate-/l/98.5%
associate-/r*98.5%
metadata-eval98.5%
associate-*l/98.5%
*-commutative98.5%
exp-prod98.5%
Simplified98.5%
Taylor expanded in v around inf 54.4%
associate-*l/54.4%
*-commutative54.4%
Simplified54.4%
Final simplification54.4%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* 0.5 (/ cosTheta_i (/ v cosTheta_O))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return 0.5f * (cosTheta_i / (v / cosTheta_O));
}
real(4) function code(costheta_i, costheta_o, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = 0.5e0 * (costheta_i / (v / costheta_o))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(0.5) * Float32(cosTheta_i / Float32(v / cosTheta_O))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = single(0.5) * (cosTheta_i / (v / cosTheta_O)); end
\begin{array}{l}
\\
0.5 \cdot \frac{cosTheta_i}{\frac{v}{cosTheta_O}}
\end{array}
Initial program 98.4%
distribute-neg-frac98.4%
*-commutative98.4%
distribute-rgt-neg-in98.4%
associate-*l/98.4%
*-commutative98.4%
Simplified98.4%
div-inv98.4%
associate-*l/98.5%
associate-*r/98.5%
*-commutative98.5%
Applied egg-rr98.5%
associate-*l*98.5%
*-commutative98.5%
associate-*l/98.5%
associate-*r/98.4%
associate-/l*98.4%
associate-/r*98.7%
*-commutative98.7%
Simplified98.7%
Taylor expanded in v around inf 54.4%
associate-/l*54.4%
Simplified54.4%
Final simplification54.4%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* 0.5 (* cosTheta_i (/ cosTheta_O v))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return 0.5f * (cosTheta_i * (cosTheta_O / v));
}
real(4) function code(costheta_i, costheta_o, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = 0.5e0 * (costheta_i * (costheta_o / v))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(0.5) * Float32(cosTheta_i * Float32(cosTheta_O / v))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = single(0.5) * (cosTheta_i * (cosTheta_O / v)); end
\begin{array}{l}
\\
0.5 \cdot \left(cosTheta_i \cdot \frac{cosTheta_O}{v}\right)
\end{array}
Initial program 98.4%
*-commutative98.4%
associate-*r/98.5%
*-commutative98.5%
associate-*l/98.5%
*-commutative98.5%
*-commutative98.5%
associate-*r*98.5%
associate-/l/98.5%
exp-neg98.5%
associate-/l/98.5%
associate-/r*98.5%
metadata-eval98.5%
associate-*l/98.5%
*-commutative98.5%
exp-prod98.5%
Simplified98.5%
Taylor expanded in sinTheta_O around 0 98.2%
rec-exp98.2%
distribute-neg-frac98.2%
metadata-eval98.2%
Simplified98.2%
expm1-log1p-u97.8%
Applied egg-rr97.8%
Taylor expanded in v around inf 54.4%
Final simplification54.4%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (/ 0.5 (/ v (* cosTheta_O cosTheta_i))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return 0.5f / (v / (cosTheta_O * cosTheta_i));
}
real(4) function code(costheta_i, costheta_o, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = 0.5e0 / (v / (costheta_o * costheta_i))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(0.5) / Float32(v / Float32(cosTheta_O * cosTheta_i))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = single(0.5) / (v / (cosTheta_O * cosTheta_i)); end
\begin{array}{l}
\\
\frac{0.5}{\frac{v}{cosTheta_O \cdot cosTheta_i}}
\end{array}
Initial program 98.4%
distribute-neg-frac98.4%
*-commutative98.4%
distribute-rgt-neg-in98.4%
associate-*l/98.4%
*-commutative98.4%
Simplified98.4%
div-inv98.4%
associate-*l/98.5%
associate-*r/98.5%
*-commutative98.5%
Applied egg-rr98.5%
associate-*l*98.5%
*-commutative98.5%
associate-*l/98.5%
associate-*r/98.4%
associate-/l*98.4%
associate-/r*98.7%
*-commutative98.7%
Simplified98.7%
Taylor expanded in v around inf 54.4%
associate-*r/54.4%
associate-/l*54.6%
Simplified54.6%
Final simplification54.6%
herbie shell --seed 2023230
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v)
:name "HairBSDF, Mp, upper"
:precision binary32
:pre (and (and (and (and (and (and (<= -1.0 cosTheta_i) (<= cosTheta_i 1.0)) (and (<= -1.0 cosTheta_O) (<= cosTheta_O 1.0))) (and (<= -1.0 sinTheta_i) (<= sinTheta_i 1.0))) (and (<= -1.0 sinTheta_O) (<= sinTheta_O 1.0))) (< 0.1 v)) (<= v 1.5707964))
(/ (* (exp (- (/ (* sinTheta_i sinTheta_O) v))) (/ (* cosTheta_i cosTheta_O) v)) (* (* (sinh (/ 1.0 v)) 2.0) v)))