
(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 15 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}
NOTE: cosTheta_i and cosTheta_O should be sorted in increasing order before calling this function. (FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (/ cosTheta_i (* (- (exp (/ 1.0 v)) (exp (/ -1.0 v))) (/ v (/ cosTheta_O v)))))
assert(cosTheta_i < cosTheta_O);
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return cosTheta_i / ((expf((1.0f / v)) - expf((-1.0f / v))) * (v / (cosTheta_O / v)));
}
NOTE: cosTheta_i and cosTheta_O should be sorted in increasing order before calling this function.
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 / ((exp((1.0e0 / v)) - exp(((-1.0e0) / v))) * (v / (costheta_o / v)))
end function
cosTheta_i, cosTheta_O = sort([cosTheta_i, cosTheta_O]) function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(cosTheta_i / Float32(Float32(exp(Float32(Float32(1.0) / v)) - exp(Float32(Float32(-1.0) / v))) * Float32(v / Float32(cosTheta_O / v)))) end
cosTheta_i, cosTheta_O = num2cell(sort([cosTheta_i, cosTheta_O])){:}
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v)
tmp = cosTheta_i / ((exp((single(1.0) / v)) - exp((single(-1.0) / v))) * (v / (cosTheta_O / v)));
end
\begin{array}{l}
[cosTheta_i, cosTheta_O] = \mathsf{sort}([cosTheta_i, cosTheta_O])\\
\\
\frac{cosTheta_i}{\left(e^{\frac{1}{v}} - e^{\frac{-1}{v}}\right) \cdot \frac{v}{\frac{cosTheta_O}{v}}}
\end{array}
Initial program 98.5%
*-commutative98.5%
associate-*r/98.5%
*-commutative98.5%
associate-*l/98.6%
*-commutative98.6%
*-commutative98.6%
associate-*r*98.6%
associate-/l/98.6%
exp-neg98.6%
associate-/l/98.6%
associate-/r*98.6%
metadata-eval98.6%
associate-*l/98.6%
*-commutative98.6%
exp-prod98.6%
Simplified98.6%
Taylor expanded in cosTheta_i around 0 98.6%
associate-/l*98.7%
*-commutative98.7%
*-commutative98.7%
rec-exp98.7%
unpow298.7%
associate-*r/98.7%
Simplified98.7%
Taylor expanded in sinTheta_i around 0 98.7%
associate-/l*98.6%
rec-exp98.6%
associate-/r/98.7%
unpow298.7%
rec-exp98.7%
distribute-neg-frac98.7%
metadata-eval98.7%
Simplified98.7%
pow198.7%
associate-/l*98.7%
Applied egg-rr98.7%
unpow198.7%
*-commutative98.7%
Simplified98.7%
Final simplification98.7%
NOTE: cosTheta_i and cosTheta_O should be sorted in increasing order before calling this function. (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))))))
assert(cosTheta_i < cosTheta_O);
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))));
}
NOTE: cosTheta_i and cosTheta_O should be sorted in increasing order before calling this function.
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
cosTheta_i, cosTheta_O = sort([cosTheta_i, cosTheta_O]) 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
cosTheta_i, cosTheta_O = num2cell(sort([cosTheta_i, cosTheta_O])){:}
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}
[cosTheta_i, cosTheta_O] = \mathsf{sort}([cosTheta_i, cosTheta_O])\\
\\
\frac{cosTheta_i}{v \cdot v} \cdot \frac{cosTheta_O}{e^{\frac{1}{v}} - e^{\frac{-1}{v}}}
\end{array}
Initial program 98.5%
*-commutative98.5%
associate-*r/98.5%
*-commutative98.5%
associate-*l/98.6%
*-commutative98.6%
*-commutative98.6%
associate-*r*98.6%
associate-/l/98.6%
exp-neg98.6%
associate-/l/98.6%
associate-/r*98.6%
metadata-eval98.6%
associate-*l/98.6%
*-commutative98.6%
exp-prod98.6%
Simplified98.6%
Taylor expanded in sinTheta_O around 0 98.6%
times-frac98.7%
unpow298.7%
rec-exp98.7%
distribute-neg-frac98.7%
metadata-eval98.7%
Simplified98.7%
Final simplification98.7%
NOTE: cosTheta_i and cosTheta_O should be sorted in increasing order before calling this function. (FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (/ cosTheta_i (* (- (exp (/ 1.0 v)) (exp (/ -1.0 v))) (* v (/ v cosTheta_O)))))
assert(cosTheta_i < cosTheta_O);
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return cosTheta_i / ((expf((1.0f / v)) - expf((-1.0f / v))) * (v * (v / cosTheta_O)));
}
NOTE: cosTheta_i and cosTheta_O should be sorted in increasing order before calling this function.
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 / ((exp((1.0e0 / v)) - exp(((-1.0e0) / v))) * (v * (v / costheta_o)))
end function
cosTheta_i, cosTheta_O = sort([cosTheta_i, cosTheta_O]) function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(cosTheta_i / Float32(Float32(exp(Float32(Float32(1.0) / v)) - exp(Float32(Float32(-1.0) / v))) * Float32(v * Float32(v / cosTheta_O)))) end
cosTheta_i, cosTheta_O = num2cell(sort([cosTheta_i, cosTheta_O])){:}
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v)
tmp = cosTheta_i / ((exp((single(1.0) / v)) - exp((single(-1.0) / v))) * (v * (v / cosTheta_O)));
end
\begin{array}{l}
[cosTheta_i, cosTheta_O] = \mathsf{sort}([cosTheta_i, cosTheta_O])\\
\\
\frac{cosTheta_i}{\left(e^{\frac{1}{v}} - e^{\frac{-1}{v}}\right) \cdot \left(v \cdot \frac{v}{cosTheta_O}\right)}
\end{array}
Initial program 98.5%
*-commutative98.5%
associate-*r/98.5%
*-commutative98.5%
associate-*l/98.6%
*-commutative98.6%
*-commutative98.6%
associate-*r*98.6%
associate-/l/98.6%
exp-neg98.6%
associate-/l/98.6%
associate-/r*98.6%
metadata-eval98.6%
associate-*l/98.6%
*-commutative98.6%
exp-prod98.6%
Simplified98.6%
Taylor expanded in cosTheta_i around 0 98.6%
associate-/l*98.7%
*-commutative98.7%
*-commutative98.7%
rec-exp98.7%
unpow298.7%
associate-*r/98.7%
Simplified98.7%
Taylor expanded in sinTheta_i around 0 98.7%
associate-/l*98.6%
rec-exp98.6%
associate-/r/98.7%
unpow298.7%
rec-exp98.7%
distribute-neg-frac98.7%
metadata-eval98.7%
Simplified98.7%
pow198.7%
associate-/l*98.7%
Applied egg-rr98.7%
unpow198.7%
*-commutative98.7%
associate-/r/98.6%
Simplified98.6%
Final simplification98.6%
NOTE: cosTheta_i and cosTheta_O should be sorted in increasing order before calling this function.
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v)
:precision binary32
(if (<= v 0.33000001311302185)
(/ cosTheta_i (* (/ (* v v) cosTheta_O) (+ (exp (/ 1.0 v)) -1.0)))
(/
cosTheta_i
(+
(* (/ v cosTheta_O) 2.0)
(+
(* 0.3333333333333333 (/ 1.0 (* v cosTheta_O)))
(* 0.016666666666666666 (/ 1.0 (* cosTheta_O (pow v 3.0)))))))))assert(cosTheta_i < cosTheta_O);
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
float tmp;
if (v <= 0.33000001311302185f) {
tmp = cosTheta_i / (((v * v) / cosTheta_O) * (expf((1.0f / v)) + -1.0f));
} else {
tmp = cosTheta_i / (((v / cosTheta_O) * 2.0f) + ((0.3333333333333333f * (1.0f / (v * cosTheta_O))) + (0.016666666666666666f * (1.0f / (cosTheta_O * powf(v, 3.0f))))));
}
return tmp;
}
NOTE: cosTheta_i and cosTheta_O should be sorted in increasing order before calling this function.
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.33000001311302185e0) then
tmp = costheta_i / (((v * v) / costheta_o) * (exp((1.0e0 / v)) + (-1.0e0)))
else
tmp = costheta_i / (((v / costheta_o) * 2.0e0) + ((0.3333333333333333e0 * (1.0e0 / (v * costheta_o))) + (0.016666666666666666e0 * (1.0e0 / (costheta_o * (v ** 3.0e0))))))
end if
code = tmp
end function
cosTheta_i, cosTheta_O = sort([cosTheta_i, cosTheta_O]) function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = Float32(0.0) if (v <= Float32(0.33000001311302185)) tmp = Float32(cosTheta_i / Float32(Float32(Float32(v * v) / cosTheta_O) * Float32(exp(Float32(Float32(1.0) / v)) + Float32(-1.0)))); else tmp = Float32(cosTheta_i / Float32(Float32(Float32(v / cosTheta_O) * Float32(2.0)) + Float32(Float32(Float32(0.3333333333333333) * Float32(Float32(1.0) / Float32(v * cosTheta_O))) + Float32(Float32(0.016666666666666666) * Float32(Float32(1.0) / Float32(cosTheta_O * (v ^ Float32(3.0)))))))); end return tmp end
cosTheta_i, cosTheta_O = num2cell(sort([cosTheta_i, cosTheta_O])){:}
function tmp_2 = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v)
tmp = single(0.0);
if (v <= single(0.33000001311302185))
tmp = cosTheta_i / (((v * v) / cosTheta_O) * (exp((single(1.0) / v)) + single(-1.0)));
else
tmp = cosTheta_i / (((v / cosTheta_O) * single(2.0)) + ((single(0.3333333333333333) * (single(1.0) / (v * cosTheta_O))) + (single(0.016666666666666666) * (single(1.0) / (cosTheta_O * (v ^ single(3.0)))))));
end
tmp_2 = tmp;
end
\begin{array}{l}
[cosTheta_i, cosTheta_O] = \mathsf{sort}([cosTheta_i, cosTheta_O])\\
\\
\begin{array}{l}
\mathbf{if}\;v \leq 0.33000001311302185:\\
\;\;\;\;\frac{cosTheta_i}{\frac{v \cdot v}{cosTheta_O} \cdot \left(e^{\frac{1}{v}} + -1\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{cosTheta_i}{\frac{v}{cosTheta_O} \cdot 2 + \left(0.3333333333333333 \cdot \frac{1}{v \cdot cosTheta_O} + 0.016666666666666666 \cdot \frac{1}{cosTheta_O \cdot {v}^{3}}\right)}\\
\end{array}
\end{array}
if v < 0.330000013Initial program 98.3%
*-commutative98.3%
associate-*r/98.2%
*-commutative98.2%
associate-*l/98.3%
*-commutative98.3%
*-commutative98.3%
associate-*r*98.3%
associate-/l/98.3%
exp-neg98.3%
associate-/l/98.3%
associate-/r*98.3%
metadata-eval98.3%
associate-*l/98.3%
*-commutative98.3%
exp-prod98.3%
Simplified98.3%
Taylor expanded in cosTheta_i around 0 98.4%
associate-/l*98.4%
*-commutative98.4%
*-commutative98.4%
rec-exp98.4%
unpow298.4%
associate-*r/98.4%
Simplified98.4%
Taylor expanded in sinTheta_i around 0 98.4%
associate-/l*98.2%
rec-exp98.2%
associate-/r/98.4%
unpow298.4%
rec-exp98.4%
distribute-neg-frac98.4%
metadata-eval98.4%
Simplified98.4%
Taylor expanded in v around inf 76.0%
if 0.330000013 < v Initial program 98.7%
*-commutative98.7%
associate-*r/98.7%
*-commutative98.7%
associate-*l/98.9%
*-commutative98.9%
*-commutative98.9%
associate-*r*98.9%
associate-/l/98.9%
exp-neg98.9%
associate-/l/98.9%
associate-/r*98.9%
metadata-eval98.9%
associate-*l/98.9%
*-commutative98.9%
exp-prod98.9%
Simplified98.9%
Taylor expanded in cosTheta_i around 0 98.8%
associate-/l*98.9%
*-commutative98.9%
*-commutative98.9%
rec-exp99.0%
unpow299.0%
associate-*r/99.0%
Simplified99.0%
Taylor expanded in sinTheta_i around 0 98.9%
associate-/l*99.0%
rec-exp98.9%
associate-/r/98.9%
unpow298.9%
rec-exp99.0%
distribute-neg-frac99.0%
metadata-eval99.0%
Simplified99.0%
Taylor expanded in v around inf 76.8%
Final simplification76.5%
NOTE: cosTheta_i and cosTheta_O should be sorted in increasing order before calling this function.
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v)
:precision binary32
(if (<= v 0.5199999809265137)
(/ cosTheta_i (* (/ (* v v) cosTheta_O) (+ (exp (/ 1.0 v)) -1.0)))
(*
(* cosTheta_i (/ cosTheta_O v))
(+
0.5
(-
(/ 0.009722222222222222 (pow v 4.0))
(/ 0.08333333333333333 (* v v)))))))assert(cosTheta_i < cosTheta_O);
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
float tmp;
if (v <= 0.5199999809265137f) {
tmp = cosTheta_i / (((v * v) / cosTheta_O) * (expf((1.0f / v)) + -1.0f));
} else {
tmp = (cosTheta_i * (cosTheta_O / v)) * (0.5f + ((0.009722222222222222f / powf(v, 4.0f)) - (0.08333333333333333f / (v * v))));
}
return tmp;
}
NOTE: cosTheta_i and cosTheta_O should be sorted in increasing order before calling this function.
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.5199999809265137e0) then
tmp = costheta_i / (((v * v) / costheta_o) * (exp((1.0e0 / v)) + (-1.0e0)))
else
tmp = (costheta_i * (costheta_o / v)) * (0.5e0 + ((0.009722222222222222e0 / (v ** 4.0e0)) - (0.08333333333333333e0 / (v * v))))
end if
code = tmp
end function
cosTheta_i, cosTheta_O = sort([cosTheta_i, cosTheta_O]) function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = Float32(0.0) if (v <= Float32(0.5199999809265137)) tmp = Float32(cosTheta_i / Float32(Float32(Float32(v * v) / cosTheta_O) * Float32(exp(Float32(Float32(1.0) / v)) + Float32(-1.0)))); else tmp = Float32(Float32(cosTheta_i * Float32(cosTheta_O / v)) * Float32(Float32(0.5) + Float32(Float32(Float32(0.009722222222222222) / (v ^ Float32(4.0))) - Float32(Float32(0.08333333333333333) / Float32(v * v))))); end return tmp end
cosTheta_i, cosTheta_O = num2cell(sort([cosTheta_i, cosTheta_O])){:}
function tmp_2 = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v)
tmp = single(0.0);
if (v <= single(0.5199999809265137))
tmp = cosTheta_i / (((v * v) / cosTheta_O) * (exp((single(1.0) / v)) + single(-1.0)));
else
tmp = (cosTheta_i * (cosTheta_O / v)) * (single(0.5) + ((single(0.009722222222222222) / (v ^ single(4.0))) - (single(0.08333333333333333) / (v * v))));
end
tmp_2 = tmp;
end
\begin{array}{l}
[cosTheta_i, cosTheta_O] = \mathsf{sort}([cosTheta_i, cosTheta_O])\\
\\
\begin{array}{l}
\mathbf{if}\;v \leq 0.5199999809265137:\\
\;\;\;\;\frac{cosTheta_i}{\frac{v \cdot v}{cosTheta_O} \cdot \left(e^{\frac{1}{v}} + -1\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(cosTheta_i \cdot \frac{cosTheta_O}{v}\right) \cdot \left(0.5 + \left(\frac{0.009722222222222222}{{v}^{4}} - \frac{0.08333333333333333}{v \cdot v}\right)\right)\\
\end{array}
\end{array}
if v < 0.519999981Initial program 98.4%
*-commutative98.4%
associate-*r/98.3%
*-commutative98.3%
associate-*l/98.4%
*-commutative98.4%
*-commutative98.4%
associate-*r*98.4%
associate-/l/98.4%
exp-neg98.4%
associate-/l/98.4%
associate-/r*98.4%
metadata-eval98.4%
associate-*l/98.4%
*-commutative98.4%
exp-prod98.4%
Simplified98.4%
Taylor expanded in cosTheta_i around 0 98.4%
associate-/l*98.5%
*-commutative98.5%
*-commutative98.5%
rec-exp98.5%
unpow298.5%
associate-*r/98.5%
Simplified98.5%
Taylor expanded in sinTheta_i around 0 98.5%
associate-/l*98.3%
rec-exp98.3%
associate-/r/98.5%
unpow298.5%
rec-exp98.5%
distribute-neg-frac98.5%
metadata-eval98.5%
Simplified98.5%
Taylor expanded in v around inf 71.7%
if 0.519999981 < v Initial program 98.8%
*-commutative98.8%
associate-*r/98.8%
*-commutative98.8%
associate-*l/99.0%
*-commutative99.0%
*-commutative99.0%
associate-*r*99.0%
associate-/l/99.0%
exp-neg99.0%
associate-/l/99.0%
associate-/r*99.0%
metadata-eval99.0%
associate-*l/99.0%
*-commutative99.0%
exp-prod99.0%
Simplified99.0%
expm1-log1p-u98.9%
Applied egg-rr98.9%
Taylor expanded in sinTheta_O around 0 98.9%
associate-/r*99.1%
rec-exp99.1%
distribute-neg-frac99.1%
metadata-eval99.1%
Simplified99.1%
Taylor expanded in v around inf 74.9%
associate--l+74.9%
associate-*r/74.9%
metadata-eval74.9%
unpow274.9%
associate-*r/74.9%
metadata-eval74.9%
Simplified74.9%
Final simplification73.0%
NOTE: cosTheta_i and cosTheta_O should be sorted in increasing order before calling this function.
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v)
:precision binary32
(if (<= v 0.5199999809265137)
(/ cosTheta_i (* (/ (* v v) cosTheta_O) (+ (exp (/ 1.0 v)) -1.0)))
(*
(* cosTheta_i (/ cosTheta_O v))
(+
0.5
(-
(/ 0.009722222222222222 (pow v 4.0))
(/ (/ 0.08333333333333333 v) v))))))assert(cosTheta_i < cosTheta_O);
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
float tmp;
if (v <= 0.5199999809265137f) {
tmp = cosTheta_i / (((v * v) / cosTheta_O) * (expf((1.0f / v)) + -1.0f));
} else {
tmp = (cosTheta_i * (cosTheta_O / v)) * (0.5f + ((0.009722222222222222f / powf(v, 4.0f)) - ((0.08333333333333333f / v) / v)));
}
return tmp;
}
NOTE: cosTheta_i and cosTheta_O should be sorted in increasing order before calling this function.
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.5199999809265137e0) then
tmp = costheta_i / (((v * v) / costheta_o) * (exp((1.0e0 / v)) + (-1.0e0)))
else
tmp = (costheta_i * (costheta_o / v)) * (0.5e0 + ((0.009722222222222222e0 / (v ** 4.0e0)) - ((0.08333333333333333e0 / v) / v)))
end if
code = tmp
end function
cosTheta_i, cosTheta_O = sort([cosTheta_i, cosTheta_O]) function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = Float32(0.0) if (v <= Float32(0.5199999809265137)) tmp = Float32(cosTheta_i / Float32(Float32(Float32(v * v) / cosTheta_O) * Float32(exp(Float32(Float32(1.0) / v)) + Float32(-1.0)))); else tmp = Float32(Float32(cosTheta_i * Float32(cosTheta_O / v)) * Float32(Float32(0.5) + Float32(Float32(Float32(0.009722222222222222) / (v ^ Float32(4.0))) - Float32(Float32(Float32(0.08333333333333333) / v) / v)))); end return tmp end
cosTheta_i, cosTheta_O = num2cell(sort([cosTheta_i, cosTheta_O])){:}
function tmp_2 = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v)
tmp = single(0.0);
if (v <= single(0.5199999809265137))
tmp = cosTheta_i / (((v * v) / cosTheta_O) * (exp((single(1.0) / v)) + single(-1.0)));
else
tmp = (cosTheta_i * (cosTheta_O / v)) * (single(0.5) + ((single(0.009722222222222222) / (v ^ single(4.0))) - ((single(0.08333333333333333) / v) / v)));
end
tmp_2 = tmp;
end
\begin{array}{l}
[cosTheta_i, cosTheta_O] = \mathsf{sort}([cosTheta_i, cosTheta_O])\\
\\
\begin{array}{l}
\mathbf{if}\;v \leq 0.5199999809265137:\\
\;\;\;\;\frac{cosTheta_i}{\frac{v \cdot v}{cosTheta_O} \cdot \left(e^{\frac{1}{v}} + -1\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(cosTheta_i \cdot \frac{cosTheta_O}{v}\right) \cdot \left(0.5 + \left(\frac{0.009722222222222222}{{v}^{4}} - \frac{\frac{0.08333333333333333}{v}}{v}\right)\right)\\
\end{array}
\end{array}
if v < 0.519999981Initial program 98.4%
*-commutative98.4%
associate-*r/98.3%
*-commutative98.3%
associate-*l/98.4%
*-commutative98.4%
*-commutative98.4%
associate-*r*98.4%
associate-/l/98.4%
exp-neg98.4%
associate-/l/98.4%
associate-/r*98.4%
metadata-eval98.4%
associate-*l/98.4%
*-commutative98.4%
exp-prod98.4%
Simplified98.4%
Taylor expanded in cosTheta_i around 0 98.4%
associate-/l*98.5%
*-commutative98.5%
*-commutative98.5%
rec-exp98.5%
unpow298.5%
associate-*r/98.5%
Simplified98.5%
Taylor expanded in sinTheta_i around 0 98.5%
associate-/l*98.3%
rec-exp98.3%
associate-/r/98.5%
unpow298.5%
rec-exp98.5%
distribute-neg-frac98.5%
metadata-eval98.5%
Simplified98.5%
Taylor expanded in v around inf 71.7%
if 0.519999981 < v Initial program 98.8%
*-commutative98.8%
associate-*r/98.8%
*-commutative98.8%
associate-*l/99.0%
*-commutative99.0%
*-commutative99.0%
associate-*r*99.0%
associate-/l/99.0%
exp-neg99.0%
associate-/l/99.0%
associate-/r*99.0%
metadata-eval99.0%
associate-*l/99.0%
*-commutative99.0%
exp-prod99.0%
Simplified99.0%
expm1-log1p-u98.9%
Applied egg-rr98.9%
Taylor expanded in sinTheta_O around 0 98.9%
associate-/r*99.1%
rec-exp99.1%
distribute-neg-frac99.1%
metadata-eval99.1%
Simplified99.1%
Taylor expanded in v around inf 74.9%
associate--l+74.9%
associate-*r/74.9%
metadata-eval74.9%
unpow274.9%
associate-*r/74.9%
metadata-eval74.9%
associate-/r*74.9%
Simplified74.9%
Final simplification73.0%
NOTE: cosTheta_i and cosTheta_O should be sorted in increasing order before calling this function.
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v)
:precision binary32
(if (<= v 0.46000000834465027)
(/ cosTheta_i (* (/ (* v v) cosTheta_O) (+ (exp (/ 1.0 v)) -1.0)))
(/
cosTheta_i
(*
(* v (/ v cosTheta_O))
(+ (/ 2.0 v) (/ 0.3333333333333333 (pow v 3.0)))))))assert(cosTheta_i < cosTheta_O);
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
float tmp;
if (v <= 0.46000000834465027f) {
tmp = cosTheta_i / (((v * v) / cosTheta_O) * (expf((1.0f / v)) + -1.0f));
} else {
tmp = cosTheta_i / ((v * (v / cosTheta_O)) * ((2.0f / v) + (0.3333333333333333f / powf(v, 3.0f))));
}
return tmp;
}
NOTE: cosTheta_i and cosTheta_O should be sorted in increasing order before calling this function.
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.46000000834465027e0) then
tmp = costheta_i / (((v * v) / costheta_o) * (exp((1.0e0 / v)) + (-1.0e0)))
else
tmp = costheta_i / ((v * (v / costheta_o)) * ((2.0e0 / v) + (0.3333333333333333e0 / (v ** 3.0e0))))
end if
code = tmp
end function
cosTheta_i, cosTheta_O = sort([cosTheta_i, cosTheta_O]) function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = Float32(0.0) if (v <= Float32(0.46000000834465027)) tmp = Float32(cosTheta_i / Float32(Float32(Float32(v * v) / cosTheta_O) * Float32(exp(Float32(Float32(1.0) / v)) + Float32(-1.0)))); else tmp = Float32(cosTheta_i / Float32(Float32(v * Float32(v / cosTheta_O)) * Float32(Float32(Float32(2.0) / v) + Float32(Float32(0.3333333333333333) / (v ^ Float32(3.0)))))); end return tmp end
cosTheta_i, cosTheta_O = num2cell(sort([cosTheta_i, cosTheta_O])){:}
function tmp_2 = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v)
tmp = single(0.0);
if (v <= single(0.46000000834465027))
tmp = cosTheta_i / (((v * v) / cosTheta_O) * (exp((single(1.0) / v)) + single(-1.0)));
else
tmp = cosTheta_i / ((v * (v / cosTheta_O)) * ((single(2.0) / v) + (single(0.3333333333333333) / (v ^ single(3.0)))));
end
tmp_2 = tmp;
end
\begin{array}{l}
[cosTheta_i, cosTheta_O] = \mathsf{sort}([cosTheta_i, cosTheta_O])\\
\\
\begin{array}{l}
\mathbf{if}\;v \leq 0.46000000834465027:\\
\;\;\;\;\frac{cosTheta_i}{\frac{v \cdot v}{cosTheta_O} \cdot \left(e^{\frac{1}{v}} + -1\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{cosTheta_i}{\left(v \cdot \frac{v}{cosTheta_O}\right) \cdot \left(\frac{2}{v} + \frac{0.3333333333333333}{{v}^{3}}\right)}\\
\end{array}
\end{array}
if v < 0.460000008Initial program 98.5%
*-commutative98.5%
associate-*r/98.4%
*-commutative98.4%
associate-*l/98.4%
*-commutative98.4%
*-commutative98.4%
associate-*r*98.4%
associate-/l/98.4%
exp-neg98.4%
associate-/l/98.4%
associate-/r*98.4%
metadata-eval98.4%
associate-*l/98.4%
*-commutative98.4%
exp-prod98.4%
Simplified98.4%
Taylor expanded in cosTheta_i around 0 98.5%
associate-/l*98.6%
*-commutative98.6%
*-commutative98.6%
rec-exp98.6%
unpow298.6%
associate-*r/98.6%
Simplified98.6%
Taylor expanded in sinTheta_i around 0 98.5%
associate-/l*98.4%
rec-exp98.4%
associate-/r/98.6%
unpow298.6%
rec-exp98.6%
distribute-neg-frac98.6%
metadata-eval98.6%
Simplified98.6%
Taylor expanded in v around inf 73.6%
if 0.460000008 < v Initial program 98.7%
*-commutative98.7%
associate-*r/98.7%
*-commutative98.7%
associate-*l/98.9%
*-commutative98.9%
*-commutative98.9%
associate-*r*98.9%
associate-/l/98.9%
exp-neg98.9%
associate-/l/98.9%
associate-/r*98.9%
metadata-eval98.9%
associate-*l/98.9%
*-commutative98.9%
exp-prod98.9%
Simplified98.9%
Taylor expanded in cosTheta_i around 0 98.7%
associate-/l*98.8%
*-commutative98.8%
*-commutative98.8%
rec-exp98.9%
unpow298.9%
associate-*r/98.9%
Simplified98.9%
Taylor expanded in sinTheta_i around 0 98.9%
associate-/l*99.0%
rec-exp99.0%
associate-/r/98.9%
unpow298.9%
rec-exp98.9%
distribute-neg-frac98.9%
metadata-eval98.9%
Simplified98.9%
Taylor expanded in v around inf 69.7%
associate-*r/69.7%
metadata-eval69.7%
associate-*r/69.7%
metadata-eval69.7%
Simplified69.7%
associate-/l*69.7%
associate-/r/69.7%
Applied egg-rr69.7%
Final simplification71.9%
NOTE: cosTheta_i and cosTheta_O should be sorted in increasing order before calling this function.
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v)
:precision binary32
(let* ((t_0 (/ (* v v) cosTheta_O)))
(if (<= v 0.46000000834465027)
(/ cosTheta_i (* t_0 (+ (exp (/ 1.0 v)) -1.0)))
(/ cosTheta_i (* t_0 (+ (/ 2.0 v) (/ 0.3333333333333333 (pow v 3.0))))))))assert(cosTheta_i < cosTheta_O);
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
float t_0 = (v * v) / cosTheta_O;
float tmp;
if (v <= 0.46000000834465027f) {
tmp = cosTheta_i / (t_0 * (expf((1.0f / v)) + -1.0f));
} else {
tmp = cosTheta_i / (t_0 * ((2.0f / v) + (0.3333333333333333f / powf(v, 3.0f))));
}
return tmp;
}
NOTE: cosTheta_i and cosTheta_O should be sorted in increasing order before calling this function.
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 * v) / costheta_o
if (v <= 0.46000000834465027e0) then
tmp = costheta_i / (t_0 * (exp((1.0e0 / v)) + (-1.0e0)))
else
tmp = costheta_i / (t_0 * ((2.0e0 / v) + (0.3333333333333333e0 / (v ** 3.0e0))))
end if
code = tmp
end function
cosTheta_i, cosTheta_O = sort([cosTheta_i, cosTheta_O]) function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) t_0 = Float32(Float32(v * v) / cosTheta_O) tmp = Float32(0.0) if (v <= Float32(0.46000000834465027)) tmp = Float32(cosTheta_i / Float32(t_0 * Float32(exp(Float32(Float32(1.0) / v)) + Float32(-1.0)))); else tmp = Float32(cosTheta_i / Float32(t_0 * Float32(Float32(Float32(2.0) / v) + Float32(Float32(0.3333333333333333) / (v ^ Float32(3.0)))))); end return tmp end
cosTheta_i, cosTheta_O = num2cell(sort([cosTheta_i, cosTheta_O])){:}
function tmp_2 = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v)
t_0 = (v * v) / cosTheta_O;
tmp = single(0.0);
if (v <= single(0.46000000834465027))
tmp = cosTheta_i / (t_0 * (exp((single(1.0) / v)) + single(-1.0)));
else
tmp = cosTheta_i / (t_0 * ((single(2.0) / v) + (single(0.3333333333333333) / (v ^ single(3.0)))));
end
tmp_2 = tmp;
end
\begin{array}{l}
[cosTheta_i, cosTheta_O] = \mathsf{sort}([cosTheta_i, cosTheta_O])\\
\\
\begin{array}{l}
t_0 := \frac{v \cdot v}{cosTheta_O}\\
\mathbf{if}\;v \leq 0.46000000834465027:\\
\;\;\;\;\frac{cosTheta_i}{t_0 \cdot \left(e^{\frac{1}{v}} + -1\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{cosTheta_i}{t_0 \cdot \left(\frac{2}{v} + \frac{0.3333333333333333}{{v}^{3}}\right)}\\
\end{array}
\end{array}
if v < 0.460000008Initial program 98.5%
*-commutative98.5%
associate-*r/98.4%
*-commutative98.4%
associate-*l/98.4%
*-commutative98.4%
*-commutative98.4%
associate-*r*98.4%
associate-/l/98.4%
exp-neg98.4%
associate-/l/98.4%
associate-/r*98.4%
metadata-eval98.4%
associate-*l/98.4%
*-commutative98.4%
exp-prod98.4%
Simplified98.4%
Taylor expanded in cosTheta_i around 0 98.5%
associate-/l*98.6%
*-commutative98.6%
*-commutative98.6%
rec-exp98.6%
unpow298.6%
associate-*r/98.6%
Simplified98.6%
Taylor expanded in sinTheta_i around 0 98.5%
associate-/l*98.4%
rec-exp98.4%
associate-/r/98.6%
unpow298.6%
rec-exp98.6%
distribute-neg-frac98.6%
metadata-eval98.6%
Simplified98.6%
Taylor expanded in v around inf 73.6%
if 0.460000008 < v Initial program 98.7%
*-commutative98.7%
associate-*r/98.7%
*-commutative98.7%
associate-*l/98.9%
*-commutative98.9%
*-commutative98.9%
associate-*r*98.9%
associate-/l/98.9%
exp-neg98.9%
associate-/l/98.9%
associate-/r*98.9%
metadata-eval98.9%
associate-*l/98.9%
*-commutative98.9%
exp-prod98.9%
Simplified98.9%
Taylor expanded in cosTheta_i around 0 98.7%
associate-/l*98.8%
*-commutative98.8%
*-commutative98.8%
rec-exp98.9%
unpow298.9%
associate-*r/98.9%
Simplified98.9%
Taylor expanded in sinTheta_i around 0 98.9%
associate-/l*99.0%
rec-exp99.0%
associate-/r/98.9%
unpow298.9%
rec-exp98.9%
distribute-neg-frac98.9%
metadata-eval98.9%
Simplified98.9%
Taylor expanded in v around inf 69.7%
associate-*r/69.7%
metadata-eval69.7%
associate-*r/69.7%
metadata-eval69.7%
Simplified69.7%
Final simplification71.9%
NOTE: cosTheta_i and cosTheta_O should be sorted in increasing order before calling this function.
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v)
:precision binary32
(let* ((t_0 (+ (/ 0.3333333333333333 v) (* v 2.0))))
(if (<= v 0.46000000834465027)
(/ cosTheta_i (* (/ v (/ cosTheta_O v)) (+ (exp (/ 1.0 v)) -1.0)))
(-
(/ cosTheta_i (/ t_0 cosTheta_O))
(* (/ sinTheta_i v) (/ (* cosTheta_i (* cosTheta_O sinTheta_O)) t_0))))))assert(cosTheta_i < cosTheta_O);
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
float t_0 = (0.3333333333333333f / v) + (v * 2.0f);
float tmp;
if (v <= 0.46000000834465027f) {
tmp = cosTheta_i / ((v / (cosTheta_O / v)) * (expf((1.0f / v)) + -1.0f));
} else {
tmp = (cosTheta_i / (t_0 / cosTheta_O)) - ((sinTheta_i / v) * ((cosTheta_i * (cosTheta_O * sinTheta_O)) / t_0));
}
return tmp;
}
NOTE: cosTheta_i and cosTheta_O should be sorted in increasing order before calling this function.
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 = (0.3333333333333333e0 / v) + (v * 2.0e0)
if (v <= 0.46000000834465027e0) then
tmp = costheta_i / ((v / (costheta_o / v)) * (exp((1.0e0 / v)) + (-1.0e0)))
else
tmp = (costheta_i / (t_0 / costheta_o)) - ((sintheta_i / v) * ((costheta_i * (costheta_o * sintheta_o)) / t_0))
end if
code = tmp
end function
cosTheta_i, cosTheta_O = sort([cosTheta_i, cosTheta_O]) function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) t_0 = Float32(Float32(Float32(0.3333333333333333) / v) + Float32(v * Float32(2.0))) tmp = Float32(0.0) if (v <= Float32(0.46000000834465027)) tmp = Float32(cosTheta_i / Float32(Float32(v / Float32(cosTheta_O / 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 / v) * Float32(Float32(cosTheta_i * Float32(cosTheta_O * sinTheta_O)) / t_0))); end return tmp end
cosTheta_i, cosTheta_O = num2cell(sort([cosTheta_i, cosTheta_O])){:}
function tmp_2 = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v)
t_0 = (single(0.3333333333333333) / v) + (v * single(2.0));
tmp = single(0.0);
if (v <= single(0.46000000834465027))
tmp = cosTheta_i / ((v / (cosTheta_O / v)) * (exp((single(1.0) / v)) + single(-1.0)));
else
tmp = (cosTheta_i / (t_0 / cosTheta_O)) - ((sinTheta_i / v) * ((cosTheta_i * (cosTheta_O * sinTheta_O)) / t_0));
end
tmp_2 = tmp;
end
\begin{array}{l}
[cosTheta_i, cosTheta_O] = \mathsf{sort}([cosTheta_i, cosTheta_O])\\
\\
\begin{array}{l}
t_0 := \frac{0.3333333333333333}{v} + v \cdot 2\\
\mathbf{if}\;v \leq 0.46000000834465027:\\
\;\;\;\;\frac{cosTheta_i}{\frac{v}{\frac{cosTheta_O}{v}} \cdot \left(e^{\frac{1}{v}} + -1\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{cosTheta_i}{\frac{t_0}{cosTheta_O}} - \frac{sinTheta_i}{v} \cdot \frac{cosTheta_i \cdot \left(cosTheta_O \cdot sinTheta_O\right)}{t_0}\\
\end{array}
\end{array}
if v < 0.460000008Initial program 98.5%
*-commutative98.5%
associate-*r/98.4%
*-commutative98.4%
associate-*l/98.4%
*-commutative98.4%
*-commutative98.4%
associate-*r*98.4%
associate-/l/98.4%
exp-neg98.4%
associate-/l/98.4%
associate-/r*98.4%
metadata-eval98.4%
associate-*l/98.4%
*-commutative98.4%
exp-prod98.4%
Simplified98.4%
Taylor expanded in cosTheta_i around 0 98.5%
associate-/l*98.6%
*-commutative98.6%
*-commutative98.6%
rec-exp98.6%
unpow298.6%
associate-*r/98.6%
Simplified98.6%
Taylor expanded in sinTheta_i around 0 98.5%
associate-/l*98.4%
rec-exp98.4%
associate-/r/98.6%
unpow298.6%
rec-exp98.6%
distribute-neg-frac98.6%
metadata-eval98.6%
Simplified98.6%
pow198.6%
associate-/l*98.6%
Applied egg-rr98.6%
unpow198.6%
*-commutative98.6%
Simplified98.6%
Taylor expanded in v around inf 73.6%
if 0.460000008 < v Initial program 98.7%
associate-*l/98.7%
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 69.7%
Taylor expanded in sinTheta_O around 0 69.7%
associate-/l*69.7%
associate-*r/69.7%
metadata-eval69.7%
*-commutative69.7%
mul-1-neg69.7%
times-frac69.7%
*-commutative69.7%
associate-*r/69.7%
metadata-eval69.7%
*-commutative69.7%
Simplified69.7%
Final simplification71.9%
NOTE: cosTheta_i and cosTheta_O should be sorted in increasing order before calling this function.
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v)
:precision binary32
(let* ((t_0 (+ (/ 0.3333333333333333 v) (* v 2.0))))
(if (<= v 0.46000000834465027)
(/ cosTheta_i (* (/ (* v v) cosTheta_O) (+ (exp (/ 1.0 v)) -1.0)))
(-
(/ cosTheta_i (/ t_0 cosTheta_O))
(* (/ sinTheta_i v) (/ (* cosTheta_i (* cosTheta_O sinTheta_O)) t_0))))))assert(cosTheta_i < cosTheta_O);
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
float t_0 = (0.3333333333333333f / v) + (v * 2.0f);
float tmp;
if (v <= 0.46000000834465027f) {
tmp = cosTheta_i / (((v * v) / cosTheta_O) * (expf((1.0f / v)) + -1.0f));
} else {
tmp = (cosTheta_i / (t_0 / cosTheta_O)) - ((sinTheta_i / v) * ((cosTheta_i * (cosTheta_O * sinTheta_O)) / t_0));
}
return tmp;
}
NOTE: cosTheta_i and cosTheta_O should be sorted in increasing order before calling this function.
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 = (0.3333333333333333e0 / v) + (v * 2.0e0)
if (v <= 0.46000000834465027e0) then
tmp = costheta_i / (((v * v) / costheta_o) * (exp((1.0e0 / v)) + (-1.0e0)))
else
tmp = (costheta_i / (t_0 / costheta_o)) - ((sintheta_i / v) * ((costheta_i * (costheta_o * sintheta_o)) / t_0))
end if
code = tmp
end function
cosTheta_i, cosTheta_O = sort([cosTheta_i, cosTheta_O]) function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) t_0 = Float32(Float32(Float32(0.3333333333333333) / v) + Float32(v * Float32(2.0))) tmp = Float32(0.0) if (v <= Float32(0.46000000834465027)) tmp = Float32(cosTheta_i / Float32(Float32(Float32(v * v) / cosTheta_O) * 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 / v) * Float32(Float32(cosTheta_i * Float32(cosTheta_O * sinTheta_O)) / t_0))); end return tmp end
cosTheta_i, cosTheta_O = num2cell(sort([cosTheta_i, cosTheta_O])){:}
function tmp_2 = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v)
t_0 = (single(0.3333333333333333) / v) + (v * single(2.0));
tmp = single(0.0);
if (v <= single(0.46000000834465027))
tmp = cosTheta_i / (((v * v) / cosTheta_O) * (exp((single(1.0) / v)) + single(-1.0)));
else
tmp = (cosTheta_i / (t_0 / cosTheta_O)) - ((sinTheta_i / v) * ((cosTheta_i * (cosTheta_O * sinTheta_O)) / t_0));
end
tmp_2 = tmp;
end
\begin{array}{l}
[cosTheta_i, cosTheta_O] = \mathsf{sort}([cosTheta_i, cosTheta_O])\\
\\
\begin{array}{l}
t_0 := \frac{0.3333333333333333}{v} + v \cdot 2\\
\mathbf{if}\;v \leq 0.46000000834465027:\\
\;\;\;\;\frac{cosTheta_i}{\frac{v \cdot v}{cosTheta_O} \cdot \left(e^{\frac{1}{v}} + -1\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{cosTheta_i}{\frac{t_0}{cosTheta_O}} - \frac{sinTheta_i}{v} \cdot \frac{cosTheta_i \cdot \left(cosTheta_O \cdot sinTheta_O\right)}{t_0}\\
\end{array}
\end{array}
if v < 0.460000008Initial program 98.5%
*-commutative98.5%
associate-*r/98.4%
*-commutative98.4%
associate-*l/98.4%
*-commutative98.4%
*-commutative98.4%
associate-*r*98.4%
associate-/l/98.4%
exp-neg98.4%
associate-/l/98.4%
associate-/r*98.4%
metadata-eval98.4%
associate-*l/98.4%
*-commutative98.4%
exp-prod98.4%
Simplified98.4%
Taylor expanded in cosTheta_i around 0 98.5%
associate-/l*98.6%
*-commutative98.6%
*-commutative98.6%
rec-exp98.6%
unpow298.6%
associate-*r/98.6%
Simplified98.6%
Taylor expanded in sinTheta_i around 0 98.5%
associate-/l*98.4%
rec-exp98.4%
associate-/r/98.6%
unpow298.6%
rec-exp98.6%
distribute-neg-frac98.6%
metadata-eval98.6%
Simplified98.6%
Taylor expanded in v around inf 73.6%
if 0.460000008 < v Initial program 98.7%
associate-*l/98.7%
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 69.7%
Taylor expanded in sinTheta_O around 0 69.7%
associate-/l*69.7%
associate-*r/69.7%
metadata-eval69.7%
*-commutative69.7%
mul-1-neg69.7%
times-frac69.7%
*-commutative69.7%
associate-*r/69.7%
metadata-eval69.7%
*-commutative69.7%
Simplified69.7%
Final simplification71.9%
NOTE: cosTheta_i and cosTheta_O should be sorted in increasing order before calling this function.
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v)
:precision binary32
(let* ((t_0 (+ (/ 0.3333333333333333 v) (* v 2.0))))
(-
(/ cosTheta_i (/ t_0 cosTheta_O))
(* (/ sinTheta_i v) (/ (* cosTheta_i (* cosTheta_O sinTheta_O)) t_0)))))assert(cosTheta_i < cosTheta_O);
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
float t_0 = (0.3333333333333333f / v) + (v * 2.0f);
return (cosTheta_i / (t_0 / cosTheta_O)) - ((sinTheta_i / v) * ((cosTheta_i * (cosTheta_O * sinTheta_O)) / t_0));
}
NOTE: cosTheta_i and cosTheta_O should be sorted in increasing order before calling this function.
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 = (0.3333333333333333e0 / v) + (v * 2.0e0)
code = (costheta_i / (t_0 / costheta_o)) - ((sintheta_i / v) * ((costheta_i * (costheta_o * sintheta_o)) / t_0))
end function
cosTheta_i, cosTheta_O = sort([cosTheta_i, cosTheta_O]) function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) t_0 = Float32(Float32(Float32(0.3333333333333333) / v) + Float32(v * Float32(2.0))) return Float32(Float32(cosTheta_i / Float32(t_0 / cosTheta_O)) - Float32(Float32(sinTheta_i / v) * Float32(Float32(cosTheta_i * Float32(cosTheta_O * sinTheta_O)) / t_0))) end
cosTheta_i, cosTheta_O = num2cell(sort([cosTheta_i, cosTheta_O])){:}
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v)
t_0 = (single(0.3333333333333333) / v) + (v * single(2.0));
tmp = (cosTheta_i / (t_0 / cosTheta_O)) - ((sinTheta_i / v) * ((cosTheta_i * (cosTheta_O * sinTheta_O)) / t_0));
end
\begin{array}{l}
[cosTheta_i, cosTheta_O] = \mathsf{sort}([cosTheta_i, cosTheta_O])\\
\\
\begin{array}{l}
t_0 := \frac{0.3333333333333333}{v} + v \cdot 2\\
\frac{cosTheta_i}{\frac{t_0}{cosTheta_O}} - \frac{sinTheta_i}{v} \cdot \frac{cosTheta_i \cdot \left(cosTheta_O \cdot sinTheta_O\right)}{t_0}
\end{array}
\end{array}
Initial program 98.5%
associate-*l/98.5%
times-frac98.7%
exp-neg98.7%
associate-*l/98.7%
*-lft-identity98.7%
associate-/l/98.7%
associate-*l*98.7%
associate-*l*98.6%
*-commutative98.6%
*-commutative98.6%
associate-*l/98.6%
Simplified98.6%
Taylor expanded in v around inf 64.8%
Taylor expanded in sinTheta_O around 0 64.8%
associate-/l*64.8%
associate-*r/64.8%
metadata-eval64.8%
*-commutative64.8%
mul-1-neg64.8%
times-frac64.8%
*-commutative64.8%
associate-*r/64.8%
metadata-eval64.8%
*-commutative64.8%
Simplified64.8%
Final simplification64.8%
NOTE: cosTheta_i and cosTheta_O should be sorted in increasing order before calling this function. (FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (/ cosTheta_i (/ (+ (/ 0.3333333333333333 v) (* v 2.0)) cosTheta_O)))
assert(cosTheta_i < cosTheta_O);
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return cosTheta_i / (((0.3333333333333333f / v) + (v * 2.0f)) / cosTheta_O);
}
NOTE: cosTheta_i and cosTheta_O should be sorted in increasing order before calling this function.
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 / (((0.3333333333333333e0 / v) + (v * 2.0e0)) / costheta_o)
end function
cosTheta_i, cosTheta_O = sort([cosTheta_i, cosTheta_O]) function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(cosTheta_i / Float32(Float32(Float32(Float32(0.3333333333333333) / v) + Float32(v * Float32(2.0))) / cosTheta_O)) end
cosTheta_i, cosTheta_O = num2cell(sort([cosTheta_i, cosTheta_O])){:}
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v)
tmp = cosTheta_i / (((single(0.3333333333333333) / v) + (v * single(2.0))) / cosTheta_O);
end
\begin{array}{l}
[cosTheta_i, cosTheta_O] = \mathsf{sort}([cosTheta_i, cosTheta_O])\\
\\
\frac{cosTheta_i}{\frac{\frac{0.3333333333333333}{v} + v \cdot 2}{cosTheta_O}}
\end{array}
Initial program 98.5%
associate-*l/98.5%
times-frac98.7%
exp-neg98.7%
associate-*l/98.7%
*-lft-identity98.7%
associate-/l/98.7%
associate-*l*98.7%
associate-*l*98.6%
*-commutative98.6%
*-commutative98.6%
associate-*l/98.6%
Simplified98.6%
Taylor expanded in v around inf 64.8%
Taylor expanded in sinTheta_O around 0 64.8%
associate-/l*64.8%
associate-*r/64.8%
metadata-eval64.8%
*-commutative64.8%
Simplified64.8%
Final simplification64.8%
NOTE: cosTheta_i and cosTheta_O should be sorted in increasing order before calling this function. (FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* 0.5 (* cosTheta_O (/ cosTheta_i v))))
assert(cosTheta_i < cosTheta_O);
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return 0.5f * (cosTheta_O * (cosTheta_i / v));
}
NOTE: cosTheta_i and cosTheta_O should be sorted in increasing order before calling this function.
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
cosTheta_i, cosTheta_O = sort([cosTheta_i, cosTheta_O]) function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(0.5) * Float32(cosTheta_O * Float32(cosTheta_i / v))) end
cosTheta_i, cosTheta_O = num2cell(sort([cosTheta_i, cosTheta_O])){:}
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}
[cosTheta_i, cosTheta_O] = \mathsf{sort}([cosTheta_i, cosTheta_O])\\
\\
0.5 \cdot \left(cosTheta_O \cdot \frac{cosTheta_i}{v}\right)
\end{array}
Initial program 98.5%
*-commutative98.5%
associate-*r/98.5%
*-commutative98.5%
associate-*l/98.6%
*-commutative98.6%
*-commutative98.6%
associate-*r*98.6%
associate-/l/98.6%
exp-neg98.6%
associate-/l/98.6%
associate-/r*98.6%
metadata-eval98.6%
associate-*l/98.6%
*-commutative98.6%
exp-prod98.6%
Simplified98.6%
Taylor expanded in v around inf 59.2%
associate-*l/59.2%
Applied egg-rr59.2%
Final simplification59.2%
NOTE: cosTheta_i and cosTheta_O should be sorted in increasing order before calling this function. (FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (/ (* 0.5 (* cosTheta_i cosTheta_O)) v))
assert(cosTheta_i < cosTheta_O);
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return (0.5f * (cosTheta_i * cosTheta_O)) / v;
}
NOTE: cosTheta_i and cosTheta_O should be sorted in increasing order before calling this function.
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
cosTheta_i, cosTheta_O = sort([cosTheta_i, cosTheta_O]) function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(Float32(0.5) * Float32(cosTheta_i * cosTheta_O)) / v) end
cosTheta_i, cosTheta_O = num2cell(sort([cosTheta_i, cosTheta_O])){:}
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}
[cosTheta_i, cosTheta_O] = \mathsf{sort}([cosTheta_i, cosTheta_O])\\
\\
\frac{0.5 \cdot \left(cosTheta_i \cdot cosTheta_O\right)}{v}
\end{array}
Initial program 98.5%
*-commutative98.5%
associate-*r/98.5%
*-commutative98.5%
associate-*l/98.6%
*-commutative98.6%
*-commutative98.6%
associate-*r*98.6%
associate-/l/98.6%
exp-neg98.6%
associate-/l/98.6%
associate-/r*98.6%
metadata-eval98.6%
associate-*l/98.6%
*-commutative98.6%
exp-prod98.6%
Simplified98.6%
Taylor expanded in v around inf 59.2%
associate-*r/59.2%
*-commutative59.2%
Applied egg-rr59.2%
Final simplification59.2%
NOTE: cosTheta_i and cosTheta_O should be sorted in increasing order before calling this function. (FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (/ cosTheta_i 0.0))
assert(cosTheta_i < cosTheta_O);
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return cosTheta_i / 0.0f;
}
NOTE: cosTheta_i and cosTheta_O should be sorted in increasing order before calling this function.
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 / 0.0e0
end function
cosTheta_i, cosTheta_O = sort([cosTheta_i, cosTheta_O]) function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(cosTheta_i / Float32(0.0)) end
cosTheta_i, cosTheta_O = num2cell(sort([cosTheta_i, cosTheta_O])){:}
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v)
tmp = cosTheta_i / single(0.0);
end
\begin{array}{l}
[cosTheta_i, cosTheta_O] = \mathsf{sort}([cosTheta_i, cosTheta_O])\\
\\
\frac{cosTheta_i}{0}
\end{array}
Initial program 98.5%
*-commutative98.5%
associate-*r/98.5%
*-commutative98.5%
associate-*l/98.6%
*-commutative98.6%
*-commutative98.6%
associate-*r*98.6%
associate-/l/98.6%
exp-neg98.6%
associate-/l/98.6%
associate-/r*98.6%
metadata-eval98.6%
associate-*l/98.6%
*-commutative98.6%
exp-prod98.6%
Simplified98.6%
Taylor expanded in cosTheta_i around 0 98.6%
associate-/l*98.7%
*-commutative98.7%
*-commutative98.7%
rec-exp98.7%
unpow298.7%
associate-*r/98.7%
Simplified98.7%
Taylor expanded in sinTheta_i around 0 98.7%
associate-/l*98.6%
rec-exp98.6%
associate-/r/98.7%
unpow298.7%
rec-exp98.7%
distribute-neg-frac98.7%
metadata-eval98.7%
Simplified98.7%
pow198.7%
associate-/l*98.7%
Applied egg-rr98.7%
unpow198.7%
*-commutative98.7%
Simplified98.7%
expm1-log1p-u45.8%
expm1-udef45.8%
add-sqr-sqrt45.8%
sqrt-unprod45.8%
frac-times45.8%
metadata-eval45.8%
metadata-eval45.8%
frac-times45.8%
sqrt-unprod-0.0%
add-sqr-sqrt3.2%
associate-/r/3.2%
Applied egg-rr3.2%
expm1-def3.5%
expm1-log1p3.5%
+-inverses3.5%
mul0-lft3.2%
Simplified3.2%
Final simplification3.2%
herbie shell --seed 2023240
(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)))