
(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 12 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 (cbrt (+ 0.6931 (log (/ 0.5 v))))))
(exp
(fma
(pow t_0 2.0)
t_0
(-
(/ cosTheta_i (/ v cosTheta_O))
(fma sinTheta_O (/ sinTheta_i v) (/ 1.0 v)))))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
float t_0 = cbrtf((0.6931f + logf((0.5f / v))));
return expf(fmaf(powf(t_0, 2.0f), t_0, ((cosTheta_i / (v / cosTheta_O)) - fmaf(sinTheta_O, (sinTheta_i / v), (1.0f / v)))));
}
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) t_0 = cbrt(Float32(Float32(0.6931) + log(Float32(Float32(0.5) / v)))) return exp(fma((t_0 ^ Float32(2.0)), t_0, Float32(Float32(cosTheta_i / Float32(v / cosTheta_O)) - fma(sinTheta_O, Float32(sinTheta_i / v), Float32(Float32(1.0) / v))))) end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt[3]{0.6931 + \log \left(\frac{0.5}{v}\right)}\\
e^{\mathsf{fma}\left({t_0}^{2}, t_0, \frac{cosTheta_i}{\frac{v}{cosTheta_O}} - \mathsf{fma}\left(sinTheta_O, \frac{sinTheta_i}{v}, \frac{1}{v}\right)\right)}
\end{array}
\end{array}
Initial program 99.8%
associate-+l+99.8%
sub-neg99.8%
associate-+l-99.8%
associate-+l-99.8%
sub-neg99.8%
associate--l-99.8%
associate-/l*99.8%
associate-/r*99.8%
metadata-eval99.8%
Simplified99.8%
+-commutative99.8%
add-cube-cbrt99.8%
fma-def99.8%
pow299.8%
associate-*l/99.8%
*-commutative99.8%
fma-def99.8%
Applied egg-rr99.8%
Final simplification99.8%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v)
:precision binary32
(exp
(+
(+
0.6931
(+
(- (/ (* cosTheta_i cosTheta_O) v) (/ (* sinTheta_O sinTheta_i) v))
(/ -1.0 v)))
(log (/ 1.0 (* v 2.0))))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return expf(((0.6931f + ((((cosTheta_i * cosTheta_O) / v) - ((sinTheta_O * sinTheta_i) / v)) + (-1.0f / v))) + logf((1.0f / (v * 2.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 = exp(((0.6931e0 + ((((costheta_i * costheta_o) / v) - ((sintheta_o * sintheta_i) / v)) + ((-1.0e0) / v))) + log((1.0e0 / (v * 2.0e0)))))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return exp(Float32(Float32(Float32(0.6931) + Float32(Float32(Float32(Float32(cosTheta_i * cosTheta_O) / v) - Float32(Float32(sinTheta_O * sinTheta_i) / v)) + Float32(Float32(-1.0) / v))) + log(Float32(Float32(1.0) / Float32(v * Float32(2.0)))))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = exp(((single(0.6931) + ((((cosTheta_i * cosTheta_O) / v) - ((sinTheta_O * sinTheta_i) / v)) + (single(-1.0) / v))) + log((single(1.0) / (v * single(2.0)))))); end
\begin{array}{l}
\\
e^{\left(0.6931 + \left(\left(\frac{cosTheta_i \cdot cosTheta_O}{v} - \frac{sinTheta_O \cdot sinTheta_i}{v}\right) + \frac{-1}{v}\right)\right) + \log \left(\frac{1}{v \cdot 2}\right)}
\end{array}
Initial program 99.8%
Final simplification99.8%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v)
:precision binary32
(exp
(+
(+ 0.6931 (log (/ 0.5 v)))
(-
(/ cosTheta_i (/ v cosTheta_O))
(+ (/ 1.0 v) (/ (* sinTheta_O sinTheta_i) 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))) + ((cosTheta_i / (v / cosTheta_O)) - ((1.0f / v) + ((sinTheta_O * sinTheta_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(((0.6931e0 + log((0.5e0 / v))) + ((costheta_i / (v / costheta_o)) - ((1.0e0 / v) + ((sintheta_o * sintheta_i) / 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(cosTheta_i / Float32(v / cosTheta_O)) - Float32(Float32(Float32(1.0) / v) + Float32(Float32(sinTheta_O * sinTheta_i) / v))))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = exp(((single(0.6931) + log((single(0.5) / v))) + ((cosTheta_i / (v / cosTheta_O)) - ((single(1.0) / v) + ((sinTheta_O * sinTheta_i) / v))))); end
\begin{array}{l}
\\
e^{\left(0.6931 + \log \left(\frac{0.5}{v}\right)\right) + \left(\frac{cosTheta_i}{\frac{v}{cosTheta_O}} - \left(\frac{1}{v} + \frac{sinTheta_O \cdot sinTheta_i}{v}\right)\right)}
\end{array}
Initial program 99.8%
associate-+l+99.8%
sub-neg99.8%
associate-+l-99.8%
associate-+l-99.8%
sub-neg99.8%
associate--l-99.8%
associate-/l*99.8%
associate-/r*99.8%
metadata-eval99.8%
Simplified99.8%
Final simplification99.8%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* (/ 0.5 v) (pow E (+ 0.6931 (/ (- -1.0 (* sinTheta_O sinTheta_i)) v)))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return (0.5f / v) * powf(((float) M_E), (0.6931f + ((-1.0f - (sinTheta_O * sinTheta_i)) / v)));
}
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(Float32(0.5) / v) * (Float32(exp(1)) ^ Float32(Float32(0.6931) + Float32(Float32(Float32(-1.0) - Float32(sinTheta_O * sinTheta_i)) / v)))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = (single(0.5) / v) * (single(2.71828182845904523536) ^ (single(0.6931) + ((single(-1.0) - (sinTheta_O * sinTheta_i)) / v))); end
\begin{array}{l}
\\
\frac{0.5}{v} \cdot {e}^{\left(0.6931 + \frac{-1 - sinTheta_O \cdot sinTheta_i}{v}\right)}
\end{array}
Initial program 99.8%
exp-sum99.7%
Simplified99.7%
*-un-lft-identity99.7%
exp-prod99.7%
+-commutative99.7%
+-commutative99.7%
associate-/r/99.7%
*-commutative99.7%
Applied egg-rr99.7%
exp-1-e99.7%
associate-+l+99.7%
associate-/r/99.7%
*-commutative99.7%
*-commutative99.7%
associate-*l/99.7%
Simplified99.7%
Taylor expanded in cosTheta_O around 0 99.7%
sub-neg99.7%
distribute-neg-in99.7%
sub-neg99.7%
distribute-neg-frac99.7%
metadata-eval99.7%
*-commutative99.7%
div-sub99.7%
Simplified99.7%
Final simplification99.7%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (/ (* (/ 0.5 v) (exp 0.6931)) (exp (/ (+ 1.0 (* sinTheta_O sinTheta_i)) v))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return ((0.5f / v) * expf(0.6931f)) / expf(((1.0f + (sinTheta_O * sinTheta_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 / v) * exp(0.6931e0)) / exp(((1.0e0 + (sintheta_o * sintheta_i)) / v))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(Float32(Float32(0.5) / v) * exp(Float32(0.6931))) / exp(Float32(Float32(Float32(1.0) + Float32(sinTheta_O * sinTheta_i)) / v))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = ((single(0.5) / v) * exp(single(0.6931))) / exp(((single(1.0) + (sinTheta_O * sinTheta_i)) / v)); end
\begin{array}{l}
\\
\frac{\frac{0.5}{v} \cdot e^{0.6931}}{e^{\frac{1 + sinTheta_O \cdot sinTheta_i}{v}}}
\end{array}
Initial program 99.8%
associate-+l+99.8%
sub-neg99.8%
associate-+l-99.8%
associate-+l-99.8%
sub-neg99.8%
associate--l-99.8%
associate-/l*99.8%
associate-/r*99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in cosTheta_i around 0 99.8%
exp-diff99.7%
exp-sum99.7%
rem-exp-log99.7%
+-commutative99.7%
associate-*l/99.7%
*-commutative99.7%
fma-udef99.7%
Simplified99.7%
Taylor expanded in v around 0 99.7%
Final simplification99.7%
(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%
sub-neg99.8%
associate-+l-99.8%
associate-+l-99.8%
sub-neg99.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.7%
Taylor expanded in cosTheta_i around 0 99.7%
Final simplification99.7%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* (/ 0.5 v) (exp (+ 0.6931 (/ (- -1.0 (* sinTheta_O sinTheta_i)) v)))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return (0.5f / v) * expf((0.6931f + ((-1.0f - (sinTheta_O * sinTheta_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 / v) * exp((0.6931e0 + (((-1.0e0) - (sintheta_o * sintheta_i)) / v)))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(Float32(0.5) / v) * exp(Float32(Float32(0.6931) + Float32(Float32(Float32(-1.0) - Float32(sinTheta_O * sinTheta_i)) / v)))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = (single(0.5) / v) * exp((single(0.6931) + ((single(-1.0) - (sinTheta_O * sinTheta_i)) / v))); end
\begin{array}{l}
\\
\frac{0.5}{v} \cdot e^{0.6931 + \frac{-1 - sinTheta_O \cdot sinTheta_i}{v}}
\end{array}
Initial program 99.8%
exp-sum99.7%
Simplified99.7%
*-un-lft-identity99.7%
exp-prod99.7%
+-commutative99.7%
+-commutative99.7%
associate-/r/99.7%
*-commutative99.7%
Applied egg-rr99.7%
exp-1-e99.7%
associate-+l+99.7%
associate-/r/99.7%
*-commutative99.7%
*-commutative99.7%
associate-*l/99.7%
Simplified99.7%
Taylor expanded in cosTheta_O around 0 99.7%
sub-neg99.7%
distribute-neg-in99.7%
sub-neg99.7%
distribute-neg-frac99.7%
metadata-eval99.7%
*-commutative99.7%
div-sub99.7%
Simplified99.7%
pow-to-exp99.6%
e-exp-199.6%
add-log-exp99.7%
*-un-lft-identity99.7%
Applied egg-rr99.7%
Final simplification99.7%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* (/ 0.5 v) (exp (+ 0.6931 (/ -1.0 v)))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return (0.5f / v) * expf((0.6931f + (-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 = (0.5e0 / v) * exp((0.6931e0 + ((-1.0e0) / v)))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(Float32(0.5) / v) * exp(Float32(Float32(0.6931) + Float32(Float32(-1.0) / v)))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = (single(0.5) / v) * exp((single(0.6931) + (single(-1.0) / v))); end
\begin{array}{l}
\\
\frac{0.5}{v} \cdot e^{0.6931 + \frac{-1}{v}}
\end{array}
Initial program 99.8%
associate-+l+99.8%
sub-neg99.8%
associate-+l-99.8%
associate-+l-99.8%
sub-neg99.8%
associate--l-99.8%
associate-/l*99.8%
associate-/r*99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in cosTheta_i around 0 99.8%
exp-diff99.7%
exp-sum99.7%
rem-exp-log99.7%
+-commutative99.7%
associate-*l/99.7%
*-commutative99.7%
fma-udef99.7%
Simplified99.7%
Taylor expanded in sinTheta_O around 0 99.7%
associate-*r/99.7%
times-frac99.7%
div-exp99.7%
Simplified99.7%
Final simplification99.7%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (exp (/ (- -1.0 (* sinTheta_O sinTheta_i)) v)))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return expf(((-1.0f - (sinTheta_O * sinTheta_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((((-1.0e0) - (sintheta_o * sintheta_i)) / v))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return exp(Float32(Float32(Float32(-1.0) - Float32(sinTheta_O * sinTheta_i)) / v)) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = exp(((single(-1.0) - (sinTheta_O * sinTheta_i)) / v)); end
\begin{array}{l}
\\
e^{\frac{-1 - sinTheta_O \cdot sinTheta_i}{v}}
\end{array}
Initial program 99.8%
associate-+l+99.8%
sub-neg99.8%
associate-+l-99.8%
associate-+l-99.8%
sub-neg99.8%
associate--l-99.8%
associate-/l*99.8%
associate-/r*99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in v around 0 98.1%
Taylor expanded in cosTheta_i around 0 98.1%
associate-*r/98.1%
mul-1-neg98.1%
*-commutative98.1%
distribute-neg-in98.1%
metadata-eval98.1%
*-commutative98.1%
Simplified98.1%
Taylor expanded in sinTheta_i around inf 98.1%
associate-*r/98.1%
neg-mul-198.1%
distribute-neg-in98.1%
metadata-eval98.1%
*-commutative98.1%
sub-neg98.1%
Simplified98.1%
Final simplification98.1%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* (/ 0.5 v) (exp (/ -1.0 v))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return (0.5f / 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 = (0.5e0 / v) * exp(((-1.0e0) / v))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(Float32(0.5) / v) * exp(Float32(Float32(-1.0) / v))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = (single(0.5) / v) * exp((single(-1.0) / v)); end
\begin{array}{l}
\\
\frac{0.5}{v} \cdot e^{\frac{-1}{v}}
\end{array}
Initial program 99.8%
associate-+l+99.8%
sub-neg99.8%
associate-+l-99.8%
associate-+l-99.8%
sub-neg99.8%
associate--l-99.8%
associate-/l*99.8%
associate-/r*99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in cosTheta_i around 0 99.8%
exp-diff99.7%
exp-sum99.7%
rem-exp-log99.7%
+-commutative99.7%
associate-*l/99.7%
*-commutative99.7%
fma-udef99.7%
Simplified99.7%
Taylor expanded in sinTheta_O around 0 99.7%
associate-*r/99.7%
times-frac99.7%
div-exp99.7%
Simplified99.7%
Taylor expanded in v around 0 98.3%
Final simplification98.3%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (exp (/ -1.0 v)))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return 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 = exp(((-1.0e0) / v))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return exp(Float32(Float32(-1.0) / v)) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = exp((single(-1.0) / v)); end
\begin{array}{l}
\\
e^{\frac{-1}{v}}
\end{array}
Initial program 99.8%
associate-+l+99.8%
sub-neg99.8%
associate-+l-99.8%
associate-+l-99.8%
sub-neg99.8%
associate--l-99.8%
associate-/l*99.8%
associate-/r*99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in v around 0 98.1%
Taylor expanded in cosTheta_i around 0 98.1%
associate-*r/98.1%
mul-1-neg98.1%
*-commutative98.1%
distribute-neg-in98.1%
metadata-eval98.1%
*-commutative98.1%
Simplified98.1%
Taylor expanded in sinTheta_i around 0 98.1%
Final simplification98.1%
(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%
sub-neg99.8%
associate-+l-99.8%
associate-+l-99.8%
sub-neg99.8%
associate--l-99.8%
associate-/l*99.8%
associate-/r*99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in v around 0 98.1%
Taylor expanded in v around inf 6.4%
Final simplification6.4%
herbie shell --seed 2023174
(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))))))