
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v)
:precision binary32
(exp
(+
(+
(-
(- (/ (* cosTheta_i cosTheta_O) v) (/ (* sinTheta_i sinTheta_O) v))
(/ 1.0 v))
0.6931)
(log (/ 1.0 (* 2.0 v))))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return expf(((((((cosTheta_i * cosTheta_O) / v) - ((sinTheta_i * sinTheta_O) / v)) - (1.0f / v)) + 0.6931f) + logf((1.0f / (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(((((((costheta_i * costheta_o) / v) - ((sintheta_i * sintheta_o) / v)) - (1.0e0 / v)) + 0.6931e0) + log((1.0e0 / (2.0e0 * v)))))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return exp(Float32(Float32(Float32(Float32(Float32(Float32(cosTheta_i * cosTheta_O) / v) - Float32(Float32(sinTheta_i * sinTheta_O) / v)) - Float32(Float32(1.0) / v)) + Float32(0.6931)) + log(Float32(Float32(1.0) / Float32(Float32(2.0) * v))))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = exp(((((((cosTheta_i * cosTheta_O) / v) - ((sinTheta_i * sinTheta_O) / v)) - (single(1.0) / v)) + single(0.6931)) + log((single(1.0) / (single(2.0) * v))))); end
\begin{array}{l}
\\
e^{\left(\left(\left(\frac{cosTheta\_i \cdot cosTheta\_O}{v} - \frac{sinTheta\_i \cdot sinTheta\_O}{v}\right) - \frac{1}{v}\right) + 0.6931\right) + \log \left(\frac{1}{2 \cdot v}\right)}
\end{array}
Sampling outcomes in binary32 precision:
Herbie found 11 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v)
:precision binary32
(exp
(+
(+
(-
(- (/ (* cosTheta_i cosTheta_O) v) (/ (* sinTheta_i sinTheta_O) v))
(/ 1.0 v))
0.6931)
(log (/ 1.0 (* 2.0 v))))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return expf(((((((cosTheta_i * cosTheta_O) / v) - ((sinTheta_i * sinTheta_O) / v)) - (1.0f / v)) + 0.6931f) + logf((1.0f / (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(((((((costheta_i * costheta_o) / v) - ((sintheta_i * sintheta_o) / v)) - (1.0e0 / v)) + 0.6931e0) + log((1.0e0 / (2.0e0 * v)))))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return exp(Float32(Float32(Float32(Float32(Float32(Float32(cosTheta_i * cosTheta_O) / v) - Float32(Float32(sinTheta_i * sinTheta_O) / v)) - Float32(Float32(1.0) / v)) + Float32(0.6931)) + log(Float32(Float32(1.0) / Float32(Float32(2.0) * v))))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = exp(((((((cosTheta_i * cosTheta_O) / v) - ((sinTheta_i * sinTheta_O) / v)) - (single(1.0) / v)) + single(0.6931)) + log((single(1.0) / (single(2.0) * v))))); end
\begin{array}{l}
\\
e^{\left(\left(\left(\frac{cosTheta\_i \cdot cosTheta\_O}{v} - \frac{sinTheta\_i \cdot sinTheta\_O}{v}\right) - \frac{1}{v}\right) + 0.6931\right) + \log \left(\frac{1}{2 \cdot v}\right)}
\end{array}
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v)
:precision binary32
(let* ((t_0 (/ (fma cosTheta_i cosTheta_O -1.0) v)))
(*
(* (/ 0.5 v) (pow (exp (cbrt (pow t_0 2.0))) (cbrt t_0)))
(pow (cbrt (exp 0.6931)) 3.0))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
float t_0 = fmaf(cosTheta_i, cosTheta_O, -1.0f) / v;
return ((0.5f / v) * powf(expf(cbrtf(powf(t_0, 2.0f))), cbrtf(t_0))) * powf(cbrtf(expf(0.6931f)), 3.0f);
}
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) t_0 = Float32(fma(cosTheta_i, cosTheta_O, Float32(-1.0)) / v) return Float32(Float32(Float32(Float32(0.5) / v) * (exp(cbrt((t_0 ^ Float32(2.0)))) ^ cbrt(t_0))) * (cbrt(exp(Float32(0.6931))) ^ Float32(3.0))) end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{\mathsf{fma}\left(cosTheta\_i, cosTheta\_O, -1\right)}{v}\\
\left(\frac{0.5}{v} \cdot {\left(e^{\sqrt[3]{{t\_0}^{2}}}\right)}^{\left(\sqrt[3]{t\_0}\right)}\right) \cdot {\left(\sqrt[3]{e^{0.6931}}\right)}^{3}
\end{array}
\end{array}
Initial program 99.8%
associate-+l+99.8%
associate--l-99.8%
associate-/l*99.8%
associate-/l*99.8%
associate-/r*99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in sinTheta_i around 0 99.8%
associate--l+99.8%
associate-*r/99.8%
+-commutative99.8%
exp-sum99.8%
associate--l+99.8%
exp-sum99.8%
add-exp-log99.8%
associate-*r/99.8%
Applied egg-rr99.8%
add-cube-cbrt99.8%
exp-prod99.8%
cbrt-unprod99.8%
pow299.8%
sub-div99.8%
fma-neg99.8%
metadata-eval99.8%
sub-div99.8%
fma-neg99.8%
metadata-eval99.8%
Applied egg-rr99.8%
metadata-eval99.8%
fma-neg99.8%
*-commutative99.8%
fma-neg99.8%
metadata-eval99.8%
metadata-eval99.8%
fma-neg99.8%
*-commutative99.8%
fma-neg99.8%
metadata-eval99.8%
Simplified99.8%
add-cube-cbrt99.9%
pow399.9%
Applied egg-rr99.9%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (let* ((t_0 (/ (fma cosTheta_i cosTheta_O -1.0) v))) (* (* (/ 0.5 v) (pow (exp (cbrt (pow t_0 2.0))) (cbrt t_0))) (exp 0.6931))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
float t_0 = fmaf(cosTheta_i, cosTheta_O, -1.0f) / v;
return ((0.5f / v) * powf(expf(cbrtf(powf(t_0, 2.0f))), cbrtf(t_0))) * expf(0.6931f);
}
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) t_0 = Float32(fma(cosTheta_i, cosTheta_O, Float32(-1.0)) / v) return Float32(Float32(Float32(Float32(0.5) / v) * (exp(cbrt((t_0 ^ Float32(2.0)))) ^ cbrt(t_0))) * exp(Float32(0.6931))) end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{\mathsf{fma}\left(cosTheta\_i, cosTheta\_O, -1\right)}{v}\\
\left(\frac{0.5}{v} \cdot {\left(e^{\sqrt[3]{{t\_0}^{2}}}\right)}^{\left(\sqrt[3]{t\_0}\right)}\right) \cdot e^{0.6931}
\end{array}
\end{array}
Initial program 99.8%
associate-+l+99.8%
associate--l-99.8%
associate-/l*99.8%
associate-/l*99.8%
associate-/r*99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in sinTheta_i around 0 99.8%
associate--l+99.8%
associate-*r/99.8%
+-commutative99.8%
exp-sum99.8%
associate--l+99.8%
exp-sum99.8%
add-exp-log99.8%
associate-*r/99.8%
Applied egg-rr99.8%
add-cube-cbrt99.8%
exp-prod99.8%
cbrt-unprod99.8%
pow299.8%
sub-div99.8%
fma-neg99.8%
metadata-eval99.8%
sub-div99.8%
fma-neg99.8%
metadata-eval99.8%
Applied egg-rr99.8%
metadata-eval99.8%
fma-neg99.8%
*-commutative99.8%
fma-neg99.8%
metadata-eval99.8%
metadata-eval99.8%
fma-neg99.8%
*-commutative99.8%
fma-neg99.8%
metadata-eval99.8%
Simplified99.8%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (exp (* cosTheta_i (/ (+ (+ 0.6931 (log (/ 0.5 v))) (/ -1.0 v)) cosTheta_i))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return expf((cosTheta_i * (((0.6931f + logf((0.5f / v))) + (-1.0f / v)) / 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 = exp((costheta_i * (((0.6931e0 + log((0.5e0 / v))) + ((-1.0e0) / v)) / costheta_i)))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return exp(Float32(cosTheta_i * Float32(Float32(Float32(Float32(0.6931) + log(Float32(Float32(0.5) / v))) + Float32(Float32(-1.0) / v)) / cosTheta_i))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = exp((cosTheta_i * (((single(0.6931) + log((single(0.5) / v))) + (single(-1.0) / v)) / cosTheta_i))); end
\begin{array}{l}
\\
e^{cosTheta\_i \cdot \frac{\left(0.6931 + \log \left(\frac{0.5}{v}\right)\right) + \frac{-1}{v}}{cosTheta\_i}}
\end{array}
Initial program 99.8%
associate-+l+99.8%
associate--l-99.8%
associate-/l*99.8%
associate-/l*99.8%
associate-/r*99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in sinTheta_i around 0 99.8%
Taylor expanded in cosTheta_i around -inf 99.8%
associate-*r*99.8%
mul-1-neg99.8%
mul-1-neg99.8%
unsub-neg99.8%
associate-*r/99.8%
mul-1-neg99.8%
associate--l+99.8%
sub-neg99.8%
distribute-neg-frac99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in cosTheta_i around 0 99.8%
Final simplification99.8%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (exp (+ (+ 0.6931 (log (/ 0.5 v))) (/ -1.0 v))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return expf(((0.6931f + logf((0.5f / v))) + (-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 = exp(((0.6931e0 + log((0.5e0 / v))) + ((-1.0e0) / v)))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return exp(Float32(Float32(Float32(0.6931) + log(Float32(Float32(0.5) / v))) + Float32(Float32(-1.0) / v))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = exp(((single(0.6931) + log((single(0.5) / v))) + (single(-1.0) / v))); end
\begin{array}{l}
\\
e^{\left(0.6931 + \log \left(\frac{0.5}{v}\right)\right) + \frac{-1}{v}}
\end{array}
Initial program 99.8%
associate-+l+99.8%
associate--l-99.8%
associate-/l*99.8%
associate-/l*99.8%
associate-/r*99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in sinTheta_i around 0 99.8%
Taylor expanded in cosTheta_O around 0 99.8%
Final simplification99.8%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (if (<= sinTheta_O 8.000000025368615e-28) (exp (/ (* cosTheta_i cosTheta_O) v)) (exp (/ (* sinTheta_O sinTheta_i) v))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
float tmp;
if (sinTheta_O <= 8.000000025368615e-28f) {
tmp = expf(((cosTheta_i * cosTheta_O) / v));
} else {
tmp = expf(((sinTheta_O * sinTheta_i) / 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) :: tmp
if (sintheta_o <= 8.000000025368615e-28) then
tmp = exp(((costheta_i * costheta_o) / v))
else
tmp = exp(((sintheta_o * sintheta_i) / v))
end if
code = tmp
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = Float32(0.0) if (sinTheta_O <= Float32(8.000000025368615e-28)) tmp = exp(Float32(Float32(cosTheta_i * cosTheta_O) / v)); else tmp = exp(Float32(Float32(sinTheta_O * sinTheta_i) / v)); end return tmp end
function tmp_2 = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = single(0.0); if (sinTheta_O <= single(8.000000025368615e-28)) tmp = exp(((cosTheta_i * cosTheta_O) / v)); else tmp = exp(((sinTheta_O * sinTheta_i) / v)); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;sinTheta\_O \leq 8.000000025368615 \cdot 10^{-28}:\\
\;\;\;\;e^{\frac{cosTheta\_i \cdot cosTheta\_O}{v}}\\
\mathbf{else}:\\
\;\;\;\;e^{\frac{sinTheta\_O \cdot sinTheta\_i}{v}}\\
\end{array}
\end{array}
if sinTheta_O < 8.00000003e-28Initial program 99.8%
associate-+l+99.8%
associate--l-99.8%
associate-/l*99.8%
associate-/l*99.8%
associate-/r*99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in sinTheta_i around -inf 99.8%
mul-1-neg99.8%
cancel-sign-sub-inv99.8%
mul-1-neg99.8%
associate--l+99.8%
associate-/l*99.8%
metadata-eval99.8%
*-lft-identity99.8%
Simplified99.8%
Taylor expanded in cosTheta_O around inf 14.7%
if 8.00000003e-28 < sinTheta_O Initial program 99.8%
associate-+l+99.8%
associate--l-99.8%
associate-/l*99.8%
associate-/l*99.8%
associate-/r*99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in sinTheta_i around -inf 99.8%
mul-1-neg99.8%
cancel-sign-sub-inv99.8%
mul-1-neg99.8%
associate--l+99.8%
associate-/l*99.8%
metadata-eval99.8%
*-lft-identity99.8%
Simplified99.8%
Taylor expanded in sinTheta_i around inf 14.1%
associate-*r/14.1%
associate-*r*14.1%
neg-mul-114.1%
*-commutative14.1%
Simplified14.1%
Taylor expanded in sinTheta_i around inf 14.1%
mul-1-neg14.1%
associate-/l*14.1%
distribute-rgt-neg-in14.1%
distribute-neg-frac214.1%
Simplified14.1%
associate-*r/14.1%
add-sqr-sqrt-0.0%
sqrt-unprod23.8%
sqr-neg23.8%
sqrt-unprod17.2%
add-sqr-sqrt17.2%
Applied egg-rr17.2%
Final simplification15.6%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (if (<= sinTheta_O 8.000000025368615e-28) (exp (/ (* cosTheta_i cosTheta_O) v)) (exp (* sinTheta_O (/ sinTheta_i v)))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
float tmp;
if (sinTheta_O <= 8.000000025368615e-28f) {
tmp = expf(((cosTheta_i * cosTheta_O) / v));
} else {
tmp = expf((sinTheta_O * (sinTheta_i / 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) :: tmp
if (sintheta_o <= 8.000000025368615e-28) then
tmp = exp(((costheta_i * costheta_o) / v))
else
tmp = exp((sintheta_o * (sintheta_i / v)))
end if
code = tmp
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = Float32(0.0) if (sinTheta_O <= Float32(8.000000025368615e-28)) tmp = exp(Float32(Float32(cosTheta_i * cosTheta_O) / v)); else tmp = exp(Float32(sinTheta_O * Float32(sinTheta_i / v))); end return tmp end
function tmp_2 = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = single(0.0); if (sinTheta_O <= single(8.000000025368615e-28)) tmp = exp(((cosTheta_i * cosTheta_O) / v)); else tmp = exp((sinTheta_O * (sinTheta_i / v))); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;sinTheta\_O \leq 8.000000025368615 \cdot 10^{-28}:\\
\;\;\;\;e^{\frac{cosTheta\_i \cdot cosTheta\_O}{v}}\\
\mathbf{else}:\\
\;\;\;\;e^{sinTheta\_O \cdot \frac{sinTheta\_i}{v}}\\
\end{array}
\end{array}
if sinTheta_O < 8.00000003e-28Initial program 99.8%
associate-+l+99.8%
associate--l-99.8%
associate-/l*99.8%
associate-/l*99.8%
associate-/r*99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in sinTheta_i around -inf 99.8%
mul-1-neg99.8%
cancel-sign-sub-inv99.8%
mul-1-neg99.8%
associate--l+99.8%
associate-/l*99.8%
metadata-eval99.8%
*-lft-identity99.8%
Simplified99.8%
Taylor expanded in cosTheta_O around inf 14.7%
if 8.00000003e-28 < sinTheta_O Initial program 99.8%
associate-+l+99.8%
associate--l-99.8%
associate-/l*99.8%
associate-/l*99.8%
associate-/r*99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in sinTheta_i around -inf 99.8%
mul-1-neg99.8%
cancel-sign-sub-inv99.8%
mul-1-neg99.8%
associate--l+99.8%
associate-/l*99.8%
metadata-eval99.8%
*-lft-identity99.8%
Simplified99.8%
Taylor expanded in sinTheta_i around inf 14.1%
associate-*r/14.1%
associate-*r*14.1%
neg-mul-114.1%
*-commutative14.1%
Simplified14.1%
Taylor expanded in sinTheta_i around inf 14.1%
mul-1-neg14.1%
associate-/l*14.1%
distribute-rgt-neg-in14.1%
distribute-neg-frac214.1%
Simplified14.1%
add-log-exp14.1%
*-un-lft-identity14.1%
log-prod14.1%
metadata-eval14.1%
add-log-exp14.1%
add-sqr-sqrt-0.0%
sqrt-unprod26.1%
sqr-neg26.1%
sqrt-unprod17.2%
add-sqr-sqrt17.2%
Applied egg-rr17.2%
+-lft-identity17.2%
Simplified17.2%
Final simplification15.6%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (if (<= sinTheta_O 8.000000025368615e-28) (exp (* cosTheta_O (/ cosTheta_i v))) (exp (* sinTheta_O (/ sinTheta_i v)))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
float tmp;
if (sinTheta_O <= 8.000000025368615e-28f) {
tmp = expf((cosTheta_O * (cosTheta_i / v)));
} else {
tmp = expf((sinTheta_O * (sinTheta_i / 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) :: tmp
if (sintheta_o <= 8.000000025368615e-28) then
tmp = exp((costheta_o * (costheta_i / v)))
else
tmp = exp((sintheta_o * (sintheta_i / v)))
end if
code = tmp
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = Float32(0.0) if (sinTheta_O <= Float32(8.000000025368615e-28)) tmp = exp(Float32(cosTheta_O * Float32(cosTheta_i / v))); else tmp = exp(Float32(sinTheta_O * Float32(sinTheta_i / v))); end return tmp end
function tmp_2 = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = single(0.0); if (sinTheta_O <= single(8.000000025368615e-28)) tmp = exp((cosTheta_O * (cosTheta_i / v))); else tmp = exp((sinTheta_O * (sinTheta_i / v))); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;sinTheta\_O \leq 8.000000025368615 \cdot 10^{-28}:\\
\;\;\;\;e^{cosTheta\_O \cdot \frac{cosTheta\_i}{v}}\\
\mathbf{else}:\\
\;\;\;\;e^{sinTheta\_O \cdot \frac{sinTheta\_i}{v}}\\
\end{array}
\end{array}
if sinTheta_O < 8.00000003e-28Initial program 99.8%
associate-+l+99.8%
associate--l-99.8%
associate-/l*99.8%
associate-/l*99.8%
associate-/r*99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in sinTheta_i around -inf 99.8%
mul-1-neg99.8%
cancel-sign-sub-inv99.8%
mul-1-neg99.8%
associate--l+99.8%
associate-/l*99.8%
metadata-eval99.8%
*-lft-identity99.8%
Simplified99.8%
Taylor expanded in cosTheta_O around inf 14.7%
associate-*r/14.7%
Simplified14.7%
if 8.00000003e-28 < sinTheta_O Initial program 99.8%
associate-+l+99.8%
associate--l-99.8%
associate-/l*99.8%
associate-/l*99.8%
associate-/r*99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in sinTheta_i around -inf 99.8%
mul-1-neg99.8%
cancel-sign-sub-inv99.8%
mul-1-neg99.8%
associate--l+99.8%
associate-/l*99.8%
metadata-eval99.8%
*-lft-identity99.8%
Simplified99.8%
Taylor expanded in sinTheta_i around inf 14.1%
associate-*r/14.1%
associate-*r*14.1%
neg-mul-114.1%
*-commutative14.1%
Simplified14.1%
Taylor expanded in sinTheta_i around inf 14.1%
mul-1-neg14.1%
associate-/l*14.1%
distribute-rgt-neg-in14.1%
distribute-neg-frac214.1%
Simplified14.1%
add-log-exp14.1%
*-un-lft-identity14.1%
log-prod14.1%
metadata-eval14.1%
add-log-exp14.1%
add-sqr-sqrt-0.0%
sqrt-unprod26.1%
sqr-neg26.1%
sqrt-unprod17.2%
add-sqr-sqrt17.2%
Applied egg-rr17.2%
+-lft-identity17.2%
Simplified17.2%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* 0.5 (/ (exp (+ 0.6931 (/ -1.0 v))) v)))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return 0.5f * (expf((0.6931f + (-1.0f / v))) / 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 * (exp((0.6931e0 + ((-1.0e0) / v))) / v)
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(0.5) * Float32(exp(Float32(Float32(0.6931) + Float32(Float32(-1.0) / v))) / v)) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = single(0.5) * (exp((single(0.6931) + (single(-1.0) / v))) / v); end
\begin{array}{l}
\\
0.5 \cdot \frac{e^{0.6931 + \frac{-1}{v}}}{v}
\end{array}
Initial program 99.8%
associate-+l+99.8%
associate--l-99.8%
associate-/l*99.8%
associate-/l*99.8%
associate-/r*99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in sinTheta_i around 0 99.8%
associate--l+99.8%
associate-*r/99.8%
+-commutative99.8%
exp-sum99.8%
associate--l+99.8%
exp-sum99.8%
add-exp-log99.8%
associate-*r/99.8%
Applied egg-rr99.8%
Taylor expanded in cosTheta_O around 0 99.8%
prod-exp99.8%
distribute-neg-frac99.8%
metadata-eval99.8%
Simplified99.8%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (exp (/ (+ -1.0 (* cosTheta_i cosTheta_O)) v)))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return expf(((-1.0f + (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 = exp((((-1.0e0) + (costheta_i * costheta_o)) / v))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return exp(Float32(Float32(Float32(-1.0) + Float32(cosTheta_i * cosTheta_O)) / v)) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = exp(((single(-1.0) + (cosTheta_i * cosTheta_O)) / v)); end
\begin{array}{l}
\\
e^{\frac{-1 + cosTheta\_i \cdot cosTheta\_O}{v}}
\end{array}
Initial program 99.8%
associate-+l+99.8%
associate--l-99.8%
associate-/l*99.8%
associate-/l*99.8%
associate-/r*99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in sinTheta_i around 0 99.8%
Taylor expanded in cosTheta_i around -inf 99.8%
associate-*r*99.8%
mul-1-neg99.8%
mul-1-neg99.8%
unsub-neg99.8%
associate-*r/99.8%
mul-1-neg99.8%
associate--l+99.8%
sub-neg99.8%
distribute-neg-frac99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in v around 0 98.4%
mul-1-neg98.4%
distribute-neg-frac298.4%
neg-mul-198.4%
+-commutative98.4%
distribute-lft-in98.4%
rgt-mult-inverse98.4%
*-commutative98.4%
distribute-lft-neg-out98.4%
unsub-neg98.4%
*-commutative98.4%
Simplified98.4%
Final simplification98.4%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (exp (* cosTheta_O (/ cosTheta_i v))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return expf((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((costheta_o * (costheta_i / v)))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return exp(Float32(cosTheta_O * Float32(cosTheta_i / v))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = exp((cosTheta_O * (cosTheta_i / v))); end
\begin{array}{l}
\\
e^{cosTheta\_O \cdot \frac{cosTheta\_i}{v}}
\end{array}
Initial program 99.8%
associate-+l+99.8%
associate--l-99.8%
associate-/l*99.8%
associate-/l*99.8%
associate-/r*99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in sinTheta_i around -inf 99.8%
mul-1-neg99.8%
cancel-sign-sub-inv99.8%
mul-1-neg99.8%
associate--l+99.8%
associate-/l*99.8%
metadata-eval99.8%
*-lft-identity99.8%
Simplified99.8%
Taylor expanded in cosTheta_O around inf 13.5%
associate-*r/13.5%
Simplified13.5%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 1.0)
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return 1.0f;
}
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 = 1.0e0
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(1.0) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = single(1.0); end
\begin{array}{l}
\\
1
\end{array}
Initial program 99.8%
associate-+l+99.8%
associate--l-99.8%
associate-/l*99.8%
associate-/l*99.8%
associate-/r*99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in sinTheta_i around -inf 99.8%
mul-1-neg99.8%
cancel-sign-sub-inv99.8%
mul-1-neg99.8%
associate--l+99.8%
associate-/l*99.8%
metadata-eval99.8%
*-lft-identity99.8%
Simplified99.8%
Taylor expanded in cosTheta_O around inf 13.5%
associate-*r/13.5%
Simplified13.5%
Taylor expanded in cosTheta_O around 0 6.3%
herbie shell --seed 2024136
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v)
:name "HairBSDF, Mp, lower"
:precision binary32
:pre (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))) (and (<= -1.5707964 v) (<= v 0.1)))
(exp (+ (+ (- (- (/ (* cosTheta_i cosTheta_O) v) (/ (* sinTheta_i sinTheta_O) v)) (/ 1.0 v)) 0.6931) (log (/ 1.0 (* 2.0 v))))))