
(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 8 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 (* (/ (exp (/ (- (* 0.6931 v) 1.0) v)) v) 0.5))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return (expf((((0.6931f * v) - 1.0f) / v)) / v) * 0.5f;
}
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 * v) - 1.0e0) / v)) / v) * 0.5e0
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(exp(Float32(Float32(Float32(Float32(0.6931) * v) - Float32(1.0)) / v)) / v) * Float32(0.5)) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = (exp((((single(0.6931) * v) - single(1.0)) / v)) / v) * single(0.5); end
\begin{array}{l}
\\
\frac{e^{\frac{0.6931 \cdot v - 1}{v}}}{v} \cdot 0.5
\end{array}
Initial program 99.8%
Taylor expanded in sinTheta_i around 0
associate-+r+N/A
associate--l+N/A
exp-sumN/A
lower-*.f32N/A
+-commutativeN/A
exp-sumN/A
lower-*.f32N/A
rem-exp-logN/A
lower-/.f32N/A
lower-exp.f32N/A
lower-exp.f32N/A
div-subN/A
lower-/.f32N/A
lower--.f32N/A
lower-*.f3299.8
Applied rewrites99.8%
Taylor expanded in cosTheta_i around 0
Applied rewrites99.8%
Taylor expanded in v around 0
Applied rewrites99.8%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* (exp (+ (/ -1.0 v) 0.6931)) (/ 0.5 v)))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return expf(((-1.0f / v) + 0.6931f)) * (0.5f / 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) + 0.6931e0)) * (0.5e0 / v)
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(exp(Float32(Float32(Float32(-1.0) / v) + Float32(0.6931))) * Float32(Float32(0.5) / v)) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = exp(((single(-1.0) / v) + single(0.6931))) * (single(0.5) / v); end
\begin{array}{l}
\\
e^{\frac{-1}{v} + 0.6931} \cdot \frac{0.5}{v}
\end{array}
Initial program 99.8%
Taylor expanded in sinTheta_i around 0
associate-+r+N/A
associate--l+N/A
exp-sumN/A
lower-*.f32N/A
+-commutativeN/A
exp-sumN/A
lower-*.f32N/A
rem-exp-logN/A
lower-/.f32N/A
lower-exp.f32N/A
lower-exp.f32N/A
div-subN/A
lower-/.f32N/A
lower--.f32N/A
lower-*.f3299.8
Applied rewrites99.8%
Taylor expanded in cosTheta_i around 0
Applied rewrites99.8%
Applied rewrites99.8%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (if (<= (* sinTheta_i sinTheta_O) 5.200000222756716e-37) (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_i * sinTheta_O) <= 5.200000222756716e-37f) {
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_i * sintheta_o) <= 5.200000222756716e-37) 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 (Float32(sinTheta_i * sinTheta_O) <= Float32(5.200000222756716e-37)) tmp = exp(Float32(Float32(cosTheta_O * cosTheta_i) / v)); else tmp = exp(Float32(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_i * sinTheta_O) <= single(5.200000222756716e-37)) 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\_i \cdot sinTheta\_O \leq 5.200000222756716 \cdot 10^{-37}:\\
\;\;\;\;e^{\frac{cosTheta\_O \cdot cosTheta\_i}{v}}\\
\mathbf{else}:\\
\;\;\;\;e^{\left(-sinTheta\_O\right) \cdot \frac{sinTheta\_i}{v}}\\
\end{array}
\end{array}
if (*.f32 sinTheta_i sinTheta_O) < 5.20000022e-37Initial program 99.7%
Taylor expanded in sinTheta_i around inf
*-commutativeN/A
times-fracN/A
associate-*r/N/A
metadata-evalN/A
associate-*l/N/A
associate-/l*N/A
div-addN/A
*-commutativeN/A
associate-/r*N/A
associate--l-N/A
div-subN/A
lower-*.f32N/A
Applied rewrites99.7%
Taylor expanded in cosTheta_i around inf
lower-/.f32N/A
lower-*.f3213.7
Applied rewrites13.7%
if 5.20000022e-37 < (*.f32 sinTheta_i sinTheta_O) Initial program 99.8%
Taylor expanded in sinTheta_i around inf
*-commutativeN/A
times-fracN/A
associate-*r/N/A
metadata-evalN/A
associate-*l/N/A
associate-/l*N/A
div-addN/A
*-commutativeN/A
associate-/r*N/A
associate--l-N/A
div-subN/A
lower-*.f32N/A
Applied rewrites99.8%
Taylor expanded in cosTheta_i around inf
lower-/.f32N/A
lower-*.f3215.1
Applied rewrites15.1%
Taylor expanded in sinTheta_i around inf
mul-1-negN/A
associate-/l*N/A
distribute-lft-neg-inN/A
lower-*.f32N/A
lower-neg.f32N/A
lower-/.f3241.3
Applied rewrites41.3%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* (/ (exp (/ -1.0 v)) v) 0.5))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return (expf((-1.0f / v)) / v) * 0.5f;
}
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)) / v) * 0.5e0
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(exp(Float32(Float32(-1.0) / v)) / v) * Float32(0.5)) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = (exp((single(-1.0) / v)) / v) * single(0.5); end
\begin{array}{l}
\\
\frac{e^{\frac{-1}{v}}}{v} \cdot 0.5
\end{array}
Initial program 99.8%
Taylor expanded in sinTheta_i around 0
associate-+r+N/A
associate--l+N/A
exp-sumN/A
lower-*.f32N/A
+-commutativeN/A
exp-sumN/A
lower-*.f32N/A
rem-exp-logN/A
lower-/.f32N/A
lower-exp.f32N/A
lower-exp.f32N/A
div-subN/A
lower-/.f32N/A
lower--.f32N/A
lower-*.f3299.8
Applied rewrites99.8%
Taylor expanded in cosTheta_i around 0
Applied rewrites99.8%
Taylor expanded in v around 0
Applied rewrites98.0%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (exp (/ (- (* cosTheta_O cosTheta_i) 1.0) v)))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return expf((((cosTheta_O * cosTheta_i) - 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((((costheta_o * costheta_i) - 1.0e0) / v))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return exp(Float32(Float32(Float32(cosTheta_O * cosTheta_i) - Float32(1.0)) / v)) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = exp((((cosTheta_O * cosTheta_i) - single(1.0)) / v)); end
\begin{array}{l}
\\
e^{\frac{cosTheta\_O \cdot cosTheta\_i - 1}{v}}
\end{array}
Initial program 99.8%
Taylor expanded in sinTheta_i around 0
associate-+r+N/A
associate--l+N/A
exp-sumN/A
lower-*.f32N/A
+-commutativeN/A
exp-sumN/A
lower-*.f32N/A
rem-exp-logN/A
lower-/.f32N/A
lower-exp.f32N/A
lower-exp.f32N/A
div-subN/A
lower-/.f32N/A
lower--.f32N/A
lower-*.f3299.8
Applied rewrites99.8%
Applied rewrites99.8%
Taylor expanded in v around 0
Applied rewrites97.8%
(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(Float32(cosTheta_O * 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^{\frac{cosTheta\_O \cdot cosTheta\_i}{v}}
\end{array}
Initial program 99.8%
Taylor expanded in sinTheta_i around inf
*-commutativeN/A
times-fracN/A
associate-*r/N/A
metadata-evalN/A
associate-*l/N/A
associate-/l*N/A
div-addN/A
*-commutativeN/A
associate-/r*N/A
associate--l-N/A
div-subN/A
lower-*.f32N/A
Applied rewrites99.8%
Taylor expanded in cosTheta_i around inf
lower-/.f32N/A
lower-*.f3214.0
Applied rewrites14.0%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (exp (* (/ cosTheta_i v) cosTheta_O)))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return expf(((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 = exp(((costheta_i / v) * costheta_o))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return exp(Float32(Float32(cosTheta_i / v) * cosTheta_O)) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = exp(((cosTheta_i / v) * cosTheta_O)); end
\begin{array}{l}
\\
e^{\frac{cosTheta\_i}{v} \cdot cosTheta\_O}
\end{array}
Initial program 99.8%
Taylor expanded in sinTheta_i around inf
*-commutativeN/A
times-fracN/A
associate-*r/N/A
metadata-evalN/A
associate-*l/N/A
associate-/l*N/A
div-addN/A
*-commutativeN/A
associate-/r*N/A
associate--l-N/A
div-subN/A
lower-*.f32N/A
Applied rewrites99.8%
Taylor expanded in cosTheta_i around inf
lower-/.f32N/A
lower-*.f3214.0
Applied rewrites14.0%
Applied rewrites14.0%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (/ (* 0.5 (exp 0.6931)) v))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return (0.5f * expf(0.6931f)) / 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)) / v
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(Float32(0.5) * exp(Float32(0.6931))) / v) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = (single(0.5) * exp(single(0.6931))) / v; end
\begin{array}{l}
\\
\frac{0.5 \cdot e^{0.6931}}{v}
\end{array}
Initial program 99.8%
Taylor expanded in v around inf
associate-+r+N/A
exp-sumN/A
rem-exp-logN/A
lower-*.f32N/A
exp-sumN/A
rem-exp-logN/A
lower-*.f32N/A
lower-exp.f32N/A
lower-/.f324.6
Applied rewrites4.6%
Applied rewrites4.6%
herbie shell --seed 2024326
(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))))))