
(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 (+ (log 0.5) (/ (* cosTheta_O cosTheta_i) v))) (/ -1.0 v))) v))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return expf(((0.6931f + (logf(0.5f) + ((cosTheta_O * cosTheta_i) / v))) + (-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 = exp(((0.6931e0 + (log(0.5e0) + ((costheta_o * costheta_i) / v))) + ((-1.0e0) / v))) / v
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(exp(Float32(Float32(Float32(0.6931) + Float32(log(Float32(0.5)) + Float32(Float32(cosTheta_O * cosTheta_i) / v))) + Float32(Float32(-1.0) / v))) / v) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = exp(((single(0.6931) + (log(single(0.5)) + ((cosTheta_O * cosTheta_i) / v))) + (single(-1.0) / v))) / v; end
\begin{array}{l}
\\
\frac{e^{\left(0.6931 + \left(\log 0.5 + \frac{cosTheta\_O \cdot cosTheta\_i}{v}\right)\right) + \frac{-1}{v}}}{v}
\end{array}
Initial program 99.4%
+-commutative99.4%
exp-sum99.4%
rem-exp-log99.4%
associate-/r*99.4%
metadata-eval99.4%
associate--l-99.4%
associate-/l*99.4%
associate-/l*99.4%
fma-define99.4%
Simplified99.4%
associate-*l/99.8%
clear-num99.8%
*-un-lft-identity99.8%
*-un-lft-identity99.8%
fma-define99.8%
associate--r+99.8%
associate-*r/99.8%
associate-*r/99.8%
sub-div99.8%
Applied egg-rr99.8%
add-exp-log99.9%
*-commutative99.9%
log-prod99.9%
add-log-exp99.9%
+-commutative99.9%
sub-div99.8%
*-commutative99.8%
Applied egg-rr99.8%
Taylor expanded in sinTheta_i around 0 99.9%
Final simplification99.9%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (/ (exp (+ (+ 0.6931 (log 0.5)) (/ -1.0 v))) v))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return expf(((0.6931f + logf(0.5f)) + (-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 = exp(((0.6931e0 + log(0.5e0)) + ((-1.0e0) / v))) / v
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(exp(Float32(Float32(Float32(0.6931) + log(Float32(0.5))) + Float32(Float32(-1.0) / v))) / v) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = exp(((single(0.6931) + log(single(0.5))) + (single(-1.0) / v))) / v; end
\begin{array}{l}
\\
\frac{e^{\left(0.6931 + \log 0.5\right) + \frac{-1}{v}}}{v}
\end{array}
Initial program 99.4%
+-commutative99.4%
exp-sum99.4%
rem-exp-log99.4%
associate-/r*99.4%
metadata-eval99.4%
associate--l-99.4%
associate-/l*99.4%
associate-/l*99.4%
fma-define99.4%
Simplified99.4%
associate-*l/99.8%
clear-num99.8%
*-un-lft-identity99.8%
*-un-lft-identity99.8%
fma-define99.8%
associate--r+99.8%
associate-*r/99.8%
associate-*r/99.8%
sub-div99.8%
Applied egg-rr99.8%
add-exp-log99.9%
*-commutative99.9%
log-prod99.9%
add-log-exp99.9%
+-commutative99.9%
sub-div99.8%
*-commutative99.8%
Applied egg-rr99.8%
Taylor expanded in sinTheta_i around 0 99.9%
Taylor expanded in cosTheta_O around 0 99.8%
Final simplification99.8%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* 0.5 (/ (exp (+ (+ 0.6931 (/ (* cosTheta_O cosTheta_i) v)) (/ -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 + ((cosTheta_O * cosTheta_i) / v)) + (-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 + ((costheta_o * costheta_i) / v)) + ((-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(Float32(0.6931) + Float32(Float32(cosTheta_O * cosTheta_i) / v)) + 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) + ((cosTheta_O * cosTheta_i) / v)) + (single(-1.0) / v))) / v); end
\begin{array}{l}
\\
0.5 \cdot \frac{e^{\left(0.6931 + \frac{cosTheta\_O \cdot cosTheta\_i}{v}\right) + \frac{-1}{v}}}{v}
\end{array}
Initial program 99.4%
+-commutative99.4%
exp-sum99.4%
rem-exp-log99.4%
associate-/r*99.4%
metadata-eval99.4%
associate--l-99.4%
associate-/l*99.4%
associate-/l*99.4%
fma-define99.4%
Simplified99.4%
Taylor expanded in sinTheta_i around 0 99.8%
Final simplification99.8%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (/ 1.0 (/ v (* 0.5 (exp (+ 0.6931 (/ -1.0 v)))))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return 1.0f / (v / (0.5f * 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 = 1.0e0 / (v / (0.5e0 * exp((0.6931e0 + ((-1.0e0) / v)))))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(1.0) / Float32(v / Float32(Float32(0.5) * 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(1.0) / (v / (single(0.5) * exp((single(0.6931) + (single(-1.0) / v))))); end
\begin{array}{l}
\\
\frac{1}{\frac{v}{0.5 \cdot e^{0.6931 + \frac{-1}{v}}}}
\end{array}
Initial program 99.4%
+-commutative99.4%
exp-sum99.4%
rem-exp-log99.4%
associate-/r*99.4%
metadata-eval99.4%
associate--l-99.4%
associate-/l*99.4%
associate-/l*99.4%
fma-define99.4%
Simplified99.4%
associate-*l/99.8%
clear-num99.8%
*-un-lft-identity99.8%
*-un-lft-identity99.8%
fma-define99.8%
associate--r+99.8%
associate-*r/99.8%
associate-*r/99.8%
sub-div99.8%
Applied egg-rr99.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 (* 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.4%
+-commutative99.4%
exp-sum99.4%
rem-exp-log99.4%
associate-/r*99.4%
metadata-eval99.4%
associate--l-99.4%
associate-/l*99.4%
associate-/l*99.4%
fma-define99.4%
Simplified99.4%
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 (pow E (/ (+ (* cosTheta_O cosTheta_i) -1.0) v)))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return powf(((float) M_E), (((cosTheta_O * cosTheta_i) + -1.0f) / v));
}
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(exp(1)) ^ 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 = single(2.71828182845904523536) ^ (((cosTheta_O * cosTheta_i) + single(-1.0)) / v); end
\begin{array}{l}
\\
{e}^{\left(\frac{cosTheta\_O \cdot cosTheta\_i + -1}{v}\right)}
\end{array}
Initial program 99.4%
+-commutative99.4%
exp-sum99.4%
rem-exp-log99.4%
associate-/r*99.4%
metadata-eval99.4%
associate--l-99.4%
associate-/l*99.4%
associate-/l*99.4%
fma-define99.4%
Simplified99.4%
*-commutative99.4%
add-exp-log99.4%
prod-exp99.4%
fma-define99.4%
associate-+r+99.4%
*-un-lft-identity99.4%
exp-prod99.4%
associate--r+99.4%
associate-+l-99.4%
Applied egg-rr99.4%
exp-1-e99.4%
associate--r-99.4%
sub-neg99.4%
fmm-def99.4%
distribute-rgt-neg-in99.4%
distribute-neg-frac99.4%
metadata-eval99.4%
+-commutative99.4%
Simplified99.4%
Taylor expanded in sinTheta_i around 0 99.4%
Taylor expanded in v around 0 97.8%
Final simplification97.8%
(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.4%
associate-+l+99.4%
associate--l-99.4%
associate-/l*99.4%
associate-/l*99.4%
associate-/r*99.4%
metadata-eval99.4%
Simplified99.4%
Taylor expanded in sinTheta_i around 0 99.4%
Taylor expanded in v around 0 97.8%
Taylor expanded in cosTheta_O around 0 97.8%
(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.4%
associate-+l+99.4%
associate--l-99.4%
associate-/l*99.4%
associate-/l*99.4%
associate-/r*99.4%
metadata-eval99.4%
Simplified99.4%
Taylor expanded in sinTheta_i around 0 99.4%
Taylor expanded in v around 0 97.8%
Taylor expanded in v around inf 6.4%
herbie shell --seed 2024155
(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))))))