
(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 19 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}
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v)
:precision binary32
(*
(* cosTheta_i cosTheta_O)
(/
(/ -1.0 v)
(/
(-
-1.0
(*
sinTheta_i
(+
(/ sinTheta_O v)
(*
sinTheta_i
(+
(/
(*
0.16666666666666666
(* sinTheta_i (* sinTheta_O (* sinTheta_O sinTheta_O))))
(* v (* v v)))
(* 0.5 (/ (/ (* sinTheta_O sinTheta_O) v) v)))))))
(/ (/ 0.5 v) (sinh (/ 1.0 v)))))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return (cosTheta_i * cosTheta_O) * ((-1.0f / v) / ((-1.0f - (sinTheta_i * ((sinTheta_O / v) + (sinTheta_i * (((0.16666666666666666f * (sinTheta_i * (sinTheta_O * (sinTheta_O * sinTheta_O)))) / (v * (v * v))) + (0.5f * (((sinTheta_O * sinTheta_O) / v) / v))))))) / ((0.5f / v) / sinhf((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 = (costheta_i * costheta_o) * (((-1.0e0) / v) / (((-1.0e0) - (sintheta_i * ((sintheta_o / v) + (sintheta_i * (((0.16666666666666666e0 * (sintheta_i * (sintheta_o * (sintheta_o * sintheta_o)))) / (v * (v * v))) + (0.5e0 * (((sintheta_o * sintheta_o) / v) / v))))))) / ((0.5e0 / v) / sinh((1.0e0 / v)))))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(cosTheta_i * cosTheta_O) * Float32(Float32(Float32(-1.0) / v) / Float32(Float32(Float32(-1.0) - Float32(sinTheta_i * Float32(Float32(sinTheta_O / v) + Float32(sinTheta_i * Float32(Float32(Float32(Float32(0.16666666666666666) * Float32(sinTheta_i * Float32(sinTheta_O * Float32(sinTheta_O * sinTheta_O)))) / Float32(v * Float32(v * v))) + Float32(Float32(0.5) * Float32(Float32(Float32(sinTheta_O * sinTheta_O) / v) / v))))))) / Float32(Float32(Float32(0.5) / v) / sinh(Float32(Float32(1.0) / v)))))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = (cosTheta_i * cosTheta_O) * ((single(-1.0) / v) / ((single(-1.0) - (sinTheta_i * ((sinTheta_O / v) + (sinTheta_i * (((single(0.16666666666666666) * (sinTheta_i * (sinTheta_O * (sinTheta_O * sinTheta_O)))) / (v * (v * v))) + (single(0.5) * (((sinTheta_O * sinTheta_O) / v) / v))))))) / ((single(0.5) / v) / sinh((single(1.0) / v))))); end
\begin{array}{l}
\\
\left(cosTheta\_i \cdot cosTheta\_O\right) \cdot \frac{\frac{-1}{v}}{\frac{-1 - sinTheta\_i \cdot \left(\frac{sinTheta\_O}{v} + sinTheta\_i \cdot \left(\frac{0.16666666666666666 \cdot \left(sinTheta\_i \cdot \left(sinTheta\_O \cdot \left(sinTheta\_O \cdot sinTheta\_O\right)\right)\right)}{v \cdot \left(v \cdot v\right)} + 0.5 \cdot \frac{\frac{sinTheta\_O \cdot sinTheta\_O}{v}}{v}\right)\right)}{\frac{\frac{0.5}{v}}{\sinh \left(\frac{1}{v}\right)}}}
\end{array}
Initial program 98.6%
/-lowering-/.f32N/A
exp-negN/A
associate-*l/N/A
*-lft-identityN/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
exp-lowering-exp.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
associate-*l*N/A
*-lowering-*.f32N/A
sinh-lowering-sinh.f32N/A
/-lowering-/.f32N/A
*-commutativeN/A
*-lowering-*.f3298.6%
Simplified98.6%
div-invN/A
associate-/l*N/A
remove-double-negN/A
div-invN/A
distribute-rgt-neg-inN/A
distribute-lft-neg-inN/A
associate-*l*N/A
*-lowering-*.f32N/A
neg-sub0N/A
--lowering--.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
Applied egg-rr98.9%
Taylor expanded in sinTheta_i around 0
+-lowering-+.f32N/A
*-lowering-*.f32N/A
+-commutativeN/A
+-lowering-+.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
Simplified98.9%
clear-numN/A
un-div-invN/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
Applied egg-rr98.9%
Final simplification98.9%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v)
:precision binary32
(*
(* cosTheta_i cosTheta_O)
(*
(/ -1.0 v)
(/
(/ (/ 0.5 v) (sinh (/ 1.0 v)))
(-
-1.0
(*
sinTheta_i
(+
(/ sinTheta_O v)
(*
sinTheta_i
(+
(/
(*
0.16666666666666666
(* sinTheta_i (* sinTheta_O (* sinTheta_O sinTheta_O))))
(* v (* v v)))
(/ (* (* sinTheta_O sinTheta_O) 0.5) (* v v)))))))))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return (cosTheta_i * cosTheta_O) * ((-1.0f / v) * (((0.5f / v) / sinhf((1.0f / v))) / (-1.0f - (sinTheta_i * ((sinTheta_O / v) + (sinTheta_i * (((0.16666666666666666f * (sinTheta_i * (sinTheta_O * (sinTheta_O * sinTheta_O)))) / (v * (v * v))) + (((sinTheta_O * sinTheta_O) * 0.5f) / (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 = (costheta_i * costheta_o) * (((-1.0e0) / v) * (((0.5e0 / v) / sinh((1.0e0 / v))) / ((-1.0e0) - (sintheta_i * ((sintheta_o / v) + (sintheta_i * (((0.16666666666666666e0 * (sintheta_i * (sintheta_o * (sintheta_o * sintheta_o)))) / (v * (v * v))) + (((sintheta_o * sintheta_o) * 0.5e0) / (v * v)))))))))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(cosTheta_i * cosTheta_O) * Float32(Float32(Float32(-1.0) / v) * Float32(Float32(Float32(Float32(0.5) / v) / sinh(Float32(Float32(1.0) / v))) / Float32(Float32(-1.0) - Float32(sinTheta_i * Float32(Float32(sinTheta_O / v) + Float32(sinTheta_i * Float32(Float32(Float32(Float32(0.16666666666666666) * Float32(sinTheta_i * Float32(sinTheta_O * Float32(sinTheta_O * sinTheta_O)))) / Float32(v * Float32(v * v))) + Float32(Float32(Float32(sinTheta_O * sinTheta_O) * Float32(0.5)) / Float32(v * v)))))))))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = (cosTheta_i * cosTheta_O) * ((single(-1.0) / v) * (((single(0.5) / v) / sinh((single(1.0) / v))) / (single(-1.0) - (sinTheta_i * ((sinTheta_O / v) + (sinTheta_i * (((single(0.16666666666666666) * (sinTheta_i * (sinTheta_O * (sinTheta_O * sinTheta_O)))) / (v * (v * v))) + (((sinTheta_O * sinTheta_O) * single(0.5)) / (v * v))))))))); end
\begin{array}{l}
\\
\left(cosTheta\_i \cdot cosTheta\_O\right) \cdot \left(\frac{-1}{v} \cdot \frac{\frac{\frac{0.5}{v}}{\sinh \left(\frac{1}{v}\right)}}{-1 - sinTheta\_i \cdot \left(\frac{sinTheta\_O}{v} + sinTheta\_i \cdot \left(\frac{0.16666666666666666 \cdot \left(sinTheta\_i \cdot \left(sinTheta\_O \cdot \left(sinTheta\_O \cdot sinTheta\_O\right)\right)\right)}{v \cdot \left(v \cdot v\right)} + \frac{\left(sinTheta\_O \cdot sinTheta\_O\right) \cdot 0.5}{v \cdot v}\right)\right)}\right)
\end{array}
Initial program 98.6%
/-lowering-/.f32N/A
exp-negN/A
associate-*l/N/A
*-lft-identityN/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
exp-lowering-exp.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
associate-*l*N/A
*-lowering-*.f32N/A
sinh-lowering-sinh.f32N/A
/-lowering-/.f32N/A
*-commutativeN/A
*-lowering-*.f3298.6%
Simplified98.6%
div-invN/A
associate-/l*N/A
remove-double-negN/A
div-invN/A
distribute-rgt-neg-inN/A
distribute-lft-neg-inN/A
associate-*l*N/A
*-lowering-*.f32N/A
neg-sub0N/A
--lowering--.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
Applied egg-rr98.9%
Taylor expanded in sinTheta_i around 0
+-lowering-+.f32N/A
*-lowering-*.f32N/A
+-commutativeN/A
+-lowering-+.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
Simplified98.9%
Final simplification98.9%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v)
:precision binary32
(*
(* cosTheta_i cosTheta_O)
(/
(/ -1.0 v)
(/
(-
-1.0
(*
sinTheta_i
(+
(/ sinTheta_O v)
(* sinTheta_i (/ (* (* sinTheta_O sinTheta_O) 0.5) (* v v))))))
(/ (/ 0.5 v) (sinh (/ 1.0 v)))))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return (cosTheta_i * cosTheta_O) * ((-1.0f / v) / ((-1.0f - (sinTheta_i * ((sinTheta_O / v) + (sinTheta_i * (((sinTheta_O * sinTheta_O) * 0.5f) / (v * v)))))) / ((0.5f / v) / sinhf((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 = (costheta_i * costheta_o) * (((-1.0e0) / v) / (((-1.0e0) - (sintheta_i * ((sintheta_o / v) + (sintheta_i * (((sintheta_o * sintheta_o) * 0.5e0) / (v * v)))))) / ((0.5e0 / v) / sinh((1.0e0 / v)))))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(cosTheta_i * cosTheta_O) * Float32(Float32(Float32(-1.0) / v) / Float32(Float32(Float32(-1.0) - Float32(sinTheta_i * Float32(Float32(sinTheta_O / v) + Float32(sinTheta_i * Float32(Float32(Float32(sinTheta_O * sinTheta_O) * Float32(0.5)) / Float32(v * v)))))) / Float32(Float32(Float32(0.5) / v) / sinh(Float32(Float32(1.0) / v)))))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = (cosTheta_i * cosTheta_O) * ((single(-1.0) / v) / ((single(-1.0) - (sinTheta_i * ((sinTheta_O / v) + (sinTheta_i * (((sinTheta_O * sinTheta_O) * single(0.5)) / (v * v)))))) / ((single(0.5) / v) / sinh((single(1.0) / v))))); end
\begin{array}{l}
\\
\left(cosTheta\_i \cdot cosTheta\_O\right) \cdot \frac{\frac{-1}{v}}{\frac{-1 - sinTheta\_i \cdot \left(\frac{sinTheta\_O}{v} + sinTheta\_i \cdot \frac{\left(sinTheta\_O \cdot sinTheta\_O\right) \cdot 0.5}{v \cdot v}\right)}{\frac{\frac{0.5}{v}}{\sinh \left(\frac{1}{v}\right)}}}
\end{array}
Initial program 98.6%
/-lowering-/.f32N/A
exp-negN/A
associate-*l/N/A
*-lft-identityN/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
exp-lowering-exp.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
associate-*l*N/A
*-lowering-*.f32N/A
sinh-lowering-sinh.f32N/A
/-lowering-/.f32N/A
*-commutativeN/A
*-lowering-*.f3298.6%
Simplified98.6%
div-invN/A
associate-/l*N/A
remove-double-negN/A
div-invN/A
distribute-rgt-neg-inN/A
distribute-lft-neg-inN/A
associate-*l*N/A
*-lowering-*.f32N/A
neg-sub0N/A
--lowering--.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
Applied egg-rr98.9%
Taylor expanded in sinTheta_i around 0
+-lowering-+.f32N/A
*-lowering-*.f32N/A
+-commutativeN/A
+-lowering-+.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
Simplified98.9%
clear-numN/A
un-div-invN/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
Applied egg-rr98.9%
Taylor expanded in sinTheta_i around 0
associate-*r/N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f32N/A
unpow2N/A
*-lowering-*.f3298.9%
Simplified98.9%
Final simplification98.9%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v)
:precision binary32
(*
(/
(-
(/
(*
(* cosTheta_i cosTheta_O)
(+
(* sinTheta_i sinTheta_O)
(/ (* sinTheta_O (* sinTheta_i (* sinTheta_i sinTheta_O))) (/ v -0.5))))
v)
(* cosTheta_i cosTheta_O))
v)
(/ (/ -0.5 v) (sinh (/ 1.0 v)))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return (((((cosTheta_i * cosTheta_O) * ((sinTheta_i * sinTheta_O) + ((sinTheta_O * (sinTheta_i * (sinTheta_i * sinTheta_O))) / (v / -0.5f)))) / v) - (cosTheta_i * cosTheta_O)) / v) * ((-0.5f / v) / sinhf((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 = (((((costheta_i * costheta_o) * ((sintheta_i * sintheta_o) + ((sintheta_o * (sintheta_i * (sintheta_i * sintheta_o))) / (v / (-0.5e0))))) / v) - (costheta_i * costheta_o)) / v) * (((-0.5e0) / v) / sinh((1.0e0 / v)))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(Float32(Float32(Float32(Float32(cosTheta_i * cosTheta_O) * Float32(Float32(sinTheta_i * sinTheta_O) + Float32(Float32(sinTheta_O * Float32(sinTheta_i * Float32(sinTheta_i * sinTheta_O))) / Float32(v / Float32(-0.5))))) / v) - Float32(cosTheta_i * cosTheta_O)) / v) * Float32(Float32(Float32(-0.5) / v) / sinh(Float32(Float32(1.0) / v)))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = (((((cosTheta_i * cosTheta_O) * ((sinTheta_i * sinTheta_O) + ((sinTheta_O * (sinTheta_i * (sinTheta_i * sinTheta_O))) / (v / single(-0.5))))) / v) - (cosTheta_i * cosTheta_O)) / v) * ((single(-0.5) / v) / sinh((single(1.0) / v))); end
\begin{array}{l}
\\
\frac{\frac{\left(cosTheta\_i \cdot cosTheta\_O\right) \cdot \left(sinTheta\_i \cdot sinTheta\_O + \frac{sinTheta\_O \cdot \left(sinTheta\_i \cdot \left(sinTheta\_i \cdot sinTheta\_O\right)\right)}{\frac{v}{-0.5}}\right)}{v} - cosTheta\_i \cdot cosTheta\_O}{v} \cdot \frac{\frac{-0.5}{v}}{\sinh \left(\frac{1}{v}\right)}
\end{array}
Initial program 98.6%
/-lowering-/.f32N/A
exp-negN/A
associate-*l/N/A
*-lft-identityN/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
exp-lowering-exp.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
associate-*l*N/A
*-lowering-*.f32N/A
sinh-lowering-sinh.f32N/A
/-lowering-/.f32N/A
*-commutativeN/A
*-lowering-*.f3298.6%
Simplified98.6%
Taylor expanded in v around -inf
mul-1-negN/A
distribute-neg-frac2N/A
/-lowering-/.f32N/A
Simplified98.5%
Applied egg-rr98.7%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* (* cosTheta_i cosTheta_O) (* (/ (/ 0.5 v) (sinh (/ 1.0 v))) (/ 1.0 v))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return (cosTheta_i * cosTheta_O) * (((0.5f / v) / sinhf((1.0f / 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 = (costheta_i * costheta_o) * (((0.5e0 / v) / sinh((1.0e0 / v))) * (1.0e0 / v))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(cosTheta_i * cosTheta_O) * Float32(Float32(Float32(Float32(0.5) / v) / sinh(Float32(Float32(1.0) / v))) * Float32(Float32(1.0) / v))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = (cosTheta_i * cosTheta_O) * (((single(0.5) / v) / sinh((single(1.0) / v))) * (single(1.0) / v)); end
\begin{array}{l}
\\
\left(cosTheta\_i \cdot cosTheta\_O\right) \cdot \left(\frac{\frac{0.5}{v}}{\sinh \left(\frac{1}{v}\right)} \cdot \frac{1}{v}\right)
\end{array}
Initial program 98.6%
/-lowering-/.f32N/A
exp-negN/A
associate-*l/N/A
*-lft-identityN/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
exp-lowering-exp.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
associate-*l*N/A
*-lowering-*.f32N/A
sinh-lowering-sinh.f32N/A
/-lowering-/.f32N/A
*-commutativeN/A
*-lowering-*.f3298.6%
Simplified98.6%
div-invN/A
associate-/l*N/A
remove-double-negN/A
div-invN/A
distribute-rgt-neg-inN/A
distribute-lft-neg-inN/A
associate-*l*N/A
*-lowering-*.f32N/A
neg-sub0N/A
--lowering--.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f32N/A
Applied egg-rr98.9%
Taylor expanded in sinTheta_i around 0
Simplified98.7%
Final simplification98.7%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* (/ (/ 0.5 v) (sinh (/ 1.0 v))) (/ (* cosTheta_i cosTheta_O) v)))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return ((0.5f / v) / sinhf((1.0f / v))) * ((cosTheta_i * cosTheta_O) / v);
}
real(4) function code(costheta_i, costheta_o, sintheta_i, sintheta_o, v)
real(4), intent (in) :: costheta_i
real(4), intent (in) :: costheta_o
real(4), intent (in) :: sintheta_i
real(4), intent (in) :: sintheta_o
real(4), intent (in) :: v
code = ((0.5e0 / v) / sinh((1.0e0 / v))) * ((costheta_i * costheta_o) / v)
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(Float32(Float32(0.5) / v) / sinh(Float32(Float32(1.0) / v))) * Float32(Float32(cosTheta_i * cosTheta_O) / v)) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = ((single(0.5) / v) / sinh((single(1.0) / v))) * ((cosTheta_i * cosTheta_O) / v); end
\begin{array}{l}
\\
\frac{\frac{0.5}{v}}{\sinh \left(\frac{1}{v}\right)} \cdot \frac{cosTheta\_i \cdot cosTheta\_O}{v}
\end{array}
Initial program 98.6%
/-lowering-/.f32N/A
exp-negN/A
associate-*l/N/A
*-lft-identityN/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
exp-lowering-exp.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
associate-*l*N/A
*-lowering-*.f32N/A
sinh-lowering-sinh.f32N/A
/-lowering-/.f32N/A
*-commutativeN/A
*-lowering-*.f3298.6%
Simplified98.6%
associate-/l/N/A
*-commutativeN/A
times-fracN/A
*-commutativeN/A
associate-*r*N/A
times-fracN/A
*-commutativeN/A
sinh-undefN/A
*-lowering-*.f32N/A
Applied egg-rr98.5%
Taylor expanded in sinTheta_i around 0
Simplified98.3%
frac-timesN/A
div-invN/A
associate-*l*N/A
metadata-evalN/A
associate-/r/N/A
div-invN/A
clear-numN/A
*-lowering-*.f32N/A
associate-*r/N/A
*-commutativeN/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
sinh-lowering-sinh.f32N/A
/-lowering-/.f3298.6%
Applied egg-rr98.6%
Final simplification98.6%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* (/ (/ 0.5 v) (sinh (/ 1.0 v))) (/ cosTheta_i (/ v cosTheta_O))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return ((0.5f / v) / sinhf((1.0f / v))) * (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 = ((0.5e0 / v) / sinh((1.0e0 / v))) * (costheta_i / (v / costheta_o))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(Float32(Float32(0.5) / v) / sinh(Float32(Float32(1.0) / v))) * Float32(cosTheta_i / Float32(v / cosTheta_O))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = ((single(0.5) / v) / sinh((single(1.0) / v))) * (cosTheta_i / (v / cosTheta_O)); end
\begin{array}{l}
\\
\frac{\frac{0.5}{v}}{\sinh \left(\frac{1}{v}\right)} \cdot \frac{cosTheta\_i}{\frac{v}{cosTheta\_O}}
\end{array}
Initial program 98.6%
/-lowering-/.f32N/A
exp-negN/A
associate-*l/N/A
*-lft-identityN/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
exp-lowering-exp.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
associate-*l*N/A
*-lowering-*.f32N/A
sinh-lowering-sinh.f32N/A
/-lowering-/.f32N/A
*-commutativeN/A
*-lowering-*.f3298.6%
Simplified98.6%
associate-/l/N/A
*-commutativeN/A
times-fracN/A
*-commutativeN/A
associate-*r*N/A
times-fracN/A
*-commutativeN/A
sinh-undefN/A
*-lowering-*.f32N/A
Applied egg-rr98.5%
Taylor expanded in sinTheta_i around 0
Simplified98.3%
*-commutativeN/A
associate-/l/N/A
associate-/r*N/A
frac-timesN/A
associate-/l*N/A
metadata-evalN/A
div-invN/A
associate-*r/N/A
clear-numN/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
associate-*r/N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
sinh-lowering-sinh.f32N/A
/-lowering-/.f3293.5%
Applied egg-rr93.5%
clear-numN/A
div-invN/A
associate-/l*N/A
times-fracN/A
frac-2negN/A
un-div-invN/A
metadata-evalN/A
frac-2negN/A
un-div-invN/A
clear-numN/A
*-lowering-*.f32N/A
frac-2negN/A
metadata-evalN/A
un-div-invN/A
frac-2negN/A
associate-*r/N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
*-commutativeN/A
Applied egg-rr98.2%
Final simplification98.2%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* cosTheta_O (/ (/ cosTheta_i v) (/ (sinh (/ 1.0 v)) (/ 0.5 v)))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return cosTheta_O * ((cosTheta_i / v) / (sinhf((1.0f / v)) / (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 = costheta_o * ((costheta_i / v) / (sinh((1.0e0 / v)) / (0.5e0 / v)))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(cosTheta_O * Float32(Float32(cosTheta_i / v) / Float32(sinh(Float32(Float32(1.0) / v)) / Float32(Float32(0.5) / v)))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = cosTheta_O * ((cosTheta_i / v) / (sinh((single(1.0) / v)) / (single(0.5) / v))); end
\begin{array}{l}
\\
cosTheta\_O \cdot \frac{\frac{cosTheta\_i}{v}}{\frac{\sinh \left(\frac{1}{v}\right)}{\frac{0.5}{v}}}
\end{array}
Initial program 98.6%
/-lowering-/.f32N/A
exp-negN/A
associate-*l/N/A
*-lft-identityN/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
exp-lowering-exp.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
associate-*l*N/A
*-lowering-*.f32N/A
sinh-lowering-sinh.f32N/A
/-lowering-/.f32N/A
*-commutativeN/A
*-lowering-*.f3298.6%
Simplified98.6%
associate-/l/N/A
*-commutativeN/A
times-fracN/A
*-commutativeN/A
associate-*r*N/A
times-fracN/A
*-commutativeN/A
sinh-undefN/A
*-lowering-*.f32N/A
Applied egg-rr98.5%
Taylor expanded in sinTheta_i around 0
Simplified98.3%
frac-timesN/A
associate-/l*N/A
*-lowering-*.f32N/A
associate-*l*N/A
metadata-evalN/A
associate-/r/N/A
div-invN/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
sinh-lowering-sinh.f32N/A
/-lowering-/.f32N/A
/-lowering-/.f3298.5%
Applied egg-rr98.5%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v)
:precision binary32
(*
(-
(/
(*
(* cosTheta_i cosTheta_O)
(+
(* sinTheta_i sinTheta_O)
(/ (* sinTheta_O (* sinTheta_i (* sinTheta_i sinTheta_O))) (/ v -0.5))))
v)
(* cosTheta_i cosTheta_O))
(/
-1.0
(*
v
(-
(/ (+ 0.3333333333333333 (/ 0.016666666666666666 (* v v))) (* v v))
-2.0)))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return ((((cosTheta_i * cosTheta_O) * ((sinTheta_i * sinTheta_O) + ((sinTheta_O * (sinTheta_i * (sinTheta_i * sinTheta_O))) / (v / -0.5f)))) / v) - (cosTheta_i * cosTheta_O)) * (-1.0f / (v * (((0.3333333333333333f + (0.016666666666666666f / (v * v))) / (v * 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 = ((((costheta_i * costheta_o) * ((sintheta_i * sintheta_o) + ((sintheta_o * (sintheta_i * (sintheta_i * sintheta_o))) / (v / (-0.5e0))))) / v) - (costheta_i * costheta_o)) * ((-1.0e0) / (v * (((0.3333333333333333e0 + (0.016666666666666666e0 / (v * v))) / (v * v)) - (-2.0e0))))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(Float32(Float32(Float32(cosTheta_i * cosTheta_O) * Float32(Float32(sinTheta_i * sinTheta_O) + Float32(Float32(sinTheta_O * Float32(sinTheta_i * Float32(sinTheta_i * sinTheta_O))) / Float32(v / Float32(-0.5))))) / v) - Float32(cosTheta_i * cosTheta_O)) * Float32(Float32(-1.0) / Float32(v * Float32(Float32(Float32(Float32(0.3333333333333333) + Float32(Float32(0.016666666666666666) / Float32(v * v))) / Float32(v * v)) - Float32(-2.0))))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = ((((cosTheta_i * cosTheta_O) * ((sinTheta_i * sinTheta_O) + ((sinTheta_O * (sinTheta_i * (sinTheta_i * sinTheta_O))) / (v / single(-0.5))))) / v) - (cosTheta_i * cosTheta_O)) * (single(-1.0) / (v * (((single(0.3333333333333333) + (single(0.016666666666666666) / (v * v))) / (v * v)) - single(-2.0)))); end
\begin{array}{l}
\\
\left(\frac{\left(cosTheta\_i \cdot cosTheta\_O\right) \cdot \left(sinTheta\_i \cdot sinTheta\_O + \frac{sinTheta\_O \cdot \left(sinTheta\_i \cdot \left(sinTheta\_i \cdot sinTheta\_O\right)\right)}{\frac{v}{-0.5}}\right)}{v} - cosTheta\_i \cdot cosTheta\_O\right) \cdot \frac{-1}{v \cdot \left(\frac{0.3333333333333333 + \frac{0.016666666666666666}{v \cdot v}}{v \cdot v} - -2\right)}
\end{array}
Initial program 98.6%
/-lowering-/.f32N/A
exp-negN/A
associate-*l/N/A
*-lft-identityN/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
exp-lowering-exp.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
associate-*l*N/A
*-lowering-*.f32N/A
sinh-lowering-sinh.f32N/A
/-lowering-/.f32N/A
*-commutativeN/A
*-lowering-*.f3298.6%
Simplified98.6%
Taylor expanded in v around -inf
mul-1-negN/A
distribute-neg-frac2N/A
/-lowering-/.f32N/A
Simplified98.5%
Applied egg-rr98.4%
Taylor expanded in v around -inf
mul-1-negN/A
*-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f32N/A
Simplified67.3%
Final simplification67.3%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* (- (* cosTheta_i cosTheta_O) (/ (* cosTheta_O (* cosTheta_i (* sinTheta_i sinTheta_O))) v)) (/ 1.0 (* v (+ 2.0 (/ 0.3333333333333333 (* v v)))))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return ((cosTheta_i * cosTheta_O) - ((cosTheta_O * (cosTheta_i * (sinTheta_i * sinTheta_O))) / v)) * (1.0f / (v * (2.0f + (0.3333333333333333f / (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 = ((costheta_i * costheta_o) - ((costheta_o * (costheta_i * (sintheta_i * sintheta_o))) / v)) * (1.0e0 / (v * (2.0e0 + (0.3333333333333333e0 / (v * v)))))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(Float32(cosTheta_i * cosTheta_O) - Float32(Float32(cosTheta_O * Float32(cosTheta_i * Float32(sinTheta_i * sinTheta_O))) / v)) * Float32(Float32(1.0) / Float32(v * Float32(Float32(2.0) + Float32(Float32(0.3333333333333333) / Float32(v * v)))))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = ((cosTheta_i * cosTheta_O) - ((cosTheta_O * (cosTheta_i * (sinTheta_i * sinTheta_O))) / v)) * (single(1.0) / (v * (single(2.0) + (single(0.3333333333333333) / (v * v))))); end
\begin{array}{l}
\\
\left(cosTheta\_i \cdot cosTheta\_O - \frac{cosTheta\_O \cdot \left(cosTheta\_i \cdot \left(sinTheta\_i \cdot sinTheta\_O\right)\right)}{v}\right) \cdot \frac{1}{v \cdot \left(2 + \frac{0.3333333333333333}{v \cdot v}\right)}
\end{array}
Initial program 98.6%
/-lowering-/.f32N/A
exp-negN/A
associate-*l/N/A
*-lft-identityN/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
exp-lowering-exp.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
associate-*l*N/A
*-lowering-*.f32N/A
sinh-lowering-sinh.f32N/A
/-lowering-/.f32N/A
*-commutativeN/A
*-lowering-*.f3298.6%
Simplified98.6%
Taylor expanded in v around -inf
mul-1-negN/A
distribute-neg-frac2N/A
/-lowering-/.f32N/A
Simplified98.5%
Applied egg-rr98.4%
Taylor expanded in v around inf
*-lowering-*.f32N/A
+-lowering-+.f32N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f32N/A
unpow2N/A
*-lowering-*.f3261.0%
Simplified61.0%
Taylor expanded in sinTheta_i around 0
+-commutativeN/A
mul-1-negN/A
sub-negN/A
--lowering--.f32N/A
*-lowering-*.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
*-lowering-*.f3261.0%
Simplified61.0%
Final simplification61.0%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v)
:precision binary32
(*
(/
cosTheta_O
(/
(-
(/ (+ 0.3333333333333333 (/ 0.016666666666666666 (* v v))) (* v v))
-2.0)
v))
(/ (/ cosTheta_i v) v)))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return (cosTheta_O / ((((0.3333333333333333f + (0.016666666666666666f / (v * v))) / (v * v)) - -2.0f) / v)) * ((cosTheta_i / 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 = (costheta_o / ((((0.3333333333333333e0 + (0.016666666666666666e0 / (v * v))) / (v * v)) - (-2.0e0)) / v)) * ((costheta_i / v) / v)
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(cosTheta_O / Float32(Float32(Float32(Float32(Float32(0.3333333333333333) + Float32(Float32(0.016666666666666666) / Float32(v * v))) / Float32(v * v)) - Float32(-2.0)) / v)) * Float32(Float32(cosTheta_i / v) / v)) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = (cosTheta_O / ((((single(0.3333333333333333) + (single(0.016666666666666666) / (v * v))) / (v * v)) - single(-2.0)) / v)) * ((cosTheta_i / v) / v); end
\begin{array}{l}
\\
\frac{cosTheta\_O}{\frac{\frac{0.3333333333333333 + \frac{0.016666666666666666}{v \cdot v}}{v \cdot v} - -2}{v}} \cdot \frac{\frac{cosTheta\_i}{v}}{v}
\end{array}
Initial program 98.6%
/-lowering-/.f32N/A
exp-negN/A
associate-*l/N/A
*-lft-identityN/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
exp-lowering-exp.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
associate-*l*N/A
*-lowering-*.f32N/A
sinh-lowering-sinh.f32N/A
/-lowering-/.f32N/A
*-commutativeN/A
*-lowering-*.f3298.6%
Simplified98.6%
associate-/l/N/A
*-commutativeN/A
times-fracN/A
*-commutativeN/A
associate-*r*N/A
times-fracN/A
*-commutativeN/A
sinh-undefN/A
*-lowering-*.f32N/A
Applied egg-rr98.5%
Taylor expanded in sinTheta_i around 0
Simplified98.3%
Taylor expanded in v around -inf
mul-1-negN/A
distribute-neg-frac2N/A
mul-1-negN/A
/-lowering-/.f32N/A
Simplified67.2%
Final simplification67.2%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* (/ (/ cosTheta_i v) v) (/ cosTheta_O (/ (+ 2.0 (/ 0.3333333333333333 (* v v))) v))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return ((cosTheta_i / v) / v) * (cosTheta_O / ((2.0f + (0.3333333333333333f / (v * 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 = ((costheta_i / v) / v) * (costheta_o / ((2.0e0 + (0.3333333333333333e0 / (v * v))) / v))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(Float32(cosTheta_i / v) / v) * Float32(cosTheta_O / Float32(Float32(Float32(2.0) + Float32(Float32(0.3333333333333333) / Float32(v * v))) / v))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = ((cosTheta_i / v) / v) * (cosTheta_O / ((single(2.0) + (single(0.3333333333333333) / (v * v))) / v)); end
\begin{array}{l}
\\
\frac{\frac{cosTheta\_i}{v}}{v} \cdot \frac{cosTheta\_O}{\frac{2 + \frac{0.3333333333333333}{v \cdot v}}{v}}
\end{array}
Initial program 98.6%
/-lowering-/.f32N/A
exp-negN/A
associate-*l/N/A
*-lft-identityN/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
exp-lowering-exp.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
associate-*l*N/A
*-lowering-*.f32N/A
sinh-lowering-sinh.f32N/A
/-lowering-/.f32N/A
*-commutativeN/A
*-lowering-*.f3298.6%
Simplified98.6%
associate-/l/N/A
*-commutativeN/A
times-fracN/A
*-commutativeN/A
associate-*r*N/A
times-fracN/A
*-commutativeN/A
sinh-undefN/A
*-lowering-*.f32N/A
Applied egg-rr98.5%
Taylor expanded in sinTheta_i around 0
Simplified98.3%
Taylor expanded in v around inf
/-lowering-/.f32N/A
+-lowering-+.f32N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f32N/A
unpow2N/A
*-lowering-*.f3261.0%
Simplified61.0%
Final simplification61.0%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* cosTheta_O (/ cosTheta_i (* v (+ 2.0 (/ 0.3333333333333333 (* v v)))))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return cosTheta_O * (cosTheta_i / (v * (2.0f + (0.3333333333333333f / (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 = costheta_o * (costheta_i / (v * (2.0e0 + (0.3333333333333333e0 / (v * v)))))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(cosTheta_O * Float32(cosTheta_i / Float32(v * Float32(Float32(2.0) + Float32(Float32(0.3333333333333333) / Float32(v * v)))))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = cosTheta_O * (cosTheta_i / (v * (single(2.0) + (single(0.3333333333333333) / (v * v))))); end
\begin{array}{l}
\\
cosTheta\_O \cdot \frac{cosTheta\_i}{v \cdot \left(2 + \frac{0.3333333333333333}{v \cdot v}\right)}
\end{array}
Initial program 98.6%
/-lowering-/.f32N/A
exp-negN/A
associate-*l/N/A
*-lft-identityN/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
exp-lowering-exp.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
associate-*l*N/A
*-lowering-*.f32N/A
sinh-lowering-sinh.f32N/A
/-lowering-/.f32N/A
*-commutativeN/A
*-lowering-*.f3298.6%
Simplified98.6%
Taylor expanded in v around -inf
mul-1-negN/A
distribute-neg-frac2N/A
/-lowering-/.f32N/A
Simplified98.5%
Applied egg-rr98.4%
Taylor expanded in v around inf
*-lowering-*.f32N/A
+-lowering-+.f32N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f32N/A
unpow2N/A
*-lowering-*.f3261.0%
Simplified61.0%
Taylor expanded in sinTheta_i around 0
associate-/l*N/A
*-lowering-*.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
+-lowering-+.f32N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f32N/A
unpow2N/A
*-lowering-*.f3261.0%
Simplified61.0%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (/ 1.0 (/ (/ v cosTheta_i) (* cosTheta_O 0.5))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return 1.0f / ((v / cosTheta_i) / (cosTheta_O * 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 = 1.0e0 / ((v / costheta_i) / (costheta_o * 0.5e0))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(1.0) / Float32(Float32(v / cosTheta_i) / Float32(cosTheta_O * Float32(0.5)))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = single(1.0) / ((v / cosTheta_i) / (cosTheta_O * single(0.5))); end
\begin{array}{l}
\\
\frac{1}{\frac{\frac{v}{cosTheta\_i}}{cosTheta\_O \cdot 0.5}}
\end{array}
Initial program 98.6%
/-lowering-/.f32N/A
exp-negN/A
associate-*l/N/A
*-lft-identityN/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
exp-lowering-exp.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
associate-*l*N/A
*-lowering-*.f32N/A
sinh-lowering-sinh.f32N/A
/-lowering-/.f32N/A
*-commutativeN/A
*-lowering-*.f3298.6%
Simplified98.6%
Taylor expanded in v around inf
associate-*r/N/A
/-lowering-/.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
*-lowering-*.f3255.4%
Simplified55.4%
*-commutativeN/A
associate-/l*N/A
frac-2negN/A
*-commutativeN/A
un-div-invN/A
metadata-evalN/A
frac-2negN/A
*-lowering-*.f32N/A
*-commutativeN/A
distribute-lft-neg-outN/A
frac-2negN/A
metadata-evalN/A
div-invN/A
associate-/l*N/A
distribute-neg-frac2N/A
distribute-rgt-neg-inN/A
remove-double-negN/A
*-lowering-*.f32N/A
/-lowering-/.f3255.3%
Applied egg-rr55.3%
associate-*r*N/A
clear-numN/A
un-div-invN/A
clear-numN/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f3255.5%
Applied egg-rr55.5%
Final simplification55.5%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (/ 1.0 (/ v (* (* cosTheta_i cosTheta_O) 0.5))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return 1.0f / (v / ((cosTheta_i * cosTheta_O) * 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 = 1.0e0 / (v / ((costheta_i * costheta_o) * 0.5e0))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(1.0) / Float32(v / Float32(Float32(cosTheta_i * cosTheta_O) * Float32(0.5)))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = single(1.0) / (v / ((cosTheta_i * cosTheta_O) * single(0.5))); end
\begin{array}{l}
\\
\frac{1}{\frac{v}{\left(cosTheta\_i \cdot cosTheta\_O\right) \cdot 0.5}}
\end{array}
Initial program 98.6%
/-lowering-/.f32N/A
exp-negN/A
associate-*l/N/A
*-lft-identityN/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
exp-lowering-exp.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
associate-*l*N/A
*-lowering-*.f32N/A
sinh-lowering-sinh.f32N/A
/-lowering-/.f32N/A
*-commutativeN/A
*-lowering-*.f3298.6%
Simplified98.6%
Taylor expanded in v around inf
associate-*r/N/A
/-lowering-/.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
*-lowering-*.f3255.4%
Simplified55.4%
*-commutativeN/A
associate-/l*N/A
frac-2negN/A
*-commutativeN/A
un-div-invN/A
metadata-evalN/A
frac-2negN/A
*-lowering-*.f32N/A
*-commutativeN/A
distribute-lft-neg-outN/A
frac-2negN/A
metadata-evalN/A
div-invN/A
associate-/l*N/A
distribute-neg-frac2N/A
distribute-rgt-neg-inN/A
remove-double-negN/A
*-lowering-*.f32N/A
/-lowering-/.f3255.3%
Applied egg-rr55.3%
*-commutativeN/A
associate-*r/N/A
*-commutativeN/A
associate-*l/N/A
clear-numN/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
*-lowering-*.f3255.5%
Applied egg-rr55.5%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (/ 0.5 (/ v (* cosTheta_i cosTheta_O))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return 0.5f / (v / (cosTheta_i * 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 = 0.5e0 / (v / (costheta_i * costheta_o))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(0.5) / Float32(v / Float32(cosTheta_i * cosTheta_O))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = single(0.5) / (v / (cosTheta_i * cosTheta_O)); end
\begin{array}{l}
\\
\frac{0.5}{\frac{v}{cosTheta\_i \cdot cosTheta\_O}}
\end{array}
Initial program 98.6%
/-lowering-/.f32N/A
exp-negN/A
associate-*l/N/A
*-lft-identityN/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
exp-lowering-exp.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
associate-*l*N/A
*-lowering-*.f32N/A
sinh-lowering-sinh.f32N/A
/-lowering-/.f32N/A
*-commutativeN/A
*-lowering-*.f3298.6%
Simplified98.6%
Taylor expanded in v around inf
associate-*r/N/A
/-lowering-/.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
*-lowering-*.f3255.4%
Simplified55.4%
*-commutativeN/A
associate-/l*N/A
frac-2negN/A
*-commutativeN/A
un-div-invN/A
metadata-evalN/A
frac-2negN/A
*-lowering-*.f32N/A
*-commutativeN/A
distribute-lft-neg-outN/A
frac-2negN/A
metadata-evalN/A
div-invN/A
associate-/l*N/A
distribute-neg-frac2N/A
distribute-rgt-neg-inN/A
remove-double-negN/A
*-lowering-*.f32N/A
/-lowering-/.f3255.3%
Applied egg-rr55.3%
associate-*r/N/A
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f32N/A
/-lowering-/.f3255.3%
Applied egg-rr55.3%
associate-*r/N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f3255.5%
Applied egg-rr55.5%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* (* cosTheta_i cosTheta_O) (/ 0.5 v)))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return (cosTheta_i * cosTheta_O) * (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 = (costheta_i * costheta_o) * (0.5e0 / v)
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(cosTheta_i * cosTheta_O) * Float32(Float32(0.5) / v)) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = (cosTheta_i * cosTheta_O) * (single(0.5) / v); end
\begin{array}{l}
\\
\left(cosTheta\_i \cdot cosTheta\_O\right) \cdot \frac{0.5}{v}
\end{array}
Initial program 98.6%
/-lowering-/.f32N/A
exp-negN/A
associate-*l/N/A
*-lft-identityN/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
exp-lowering-exp.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
associate-*l*N/A
*-lowering-*.f32N/A
sinh-lowering-sinh.f32N/A
/-lowering-/.f32N/A
*-commutativeN/A
*-lowering-*.f3298.6%
Simplified98.6%
Taylor expanded in v around inf
associate-*r/N/A
/-lowering-/.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
*-lowering-*.f3255.4%
Simplified55.4%
associate-/l*N/A
*-commutativeN/A
*-lowering-*.f32N/A
/-lowering-/.f32N/A
*-commutativeN/A
*-lowering-*.f3255.4%
Applied egg-rr55.4%
Final simplification55.4%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* cosTheta_i (/ cosTheta_O (/ v 0.5))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return cosTheta_i * (cosTheta_O / (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 = costheta_i * (costheta_o / (v / 0.5e0))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(cosTheta_i * Float32(cosTheta_O / Float32(v / Float32(0.5)))) end
function tmp = code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) tmp = cosTheta_i * (cosTheta_O / (v / single(0.5))); end
\begin{array}{l}
\\
cosTheta\_i \cdot \frac{cosTheta\_O}{\frac{v}{0.5}}
\end{array}
Initial program 98.6%
/-lowering-/.f32N/A
exp-negN/A
associate-*l/N/A
*-lft-identityN/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
exp-lowering-exp.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
associate-*l*N/A
*-lowering-*.f32N/A
sinh-lowering-sinh.f32N/A
/-lowering-/.f32N/A
*-commutativeN/A
*-lowering-*.f3298.6%
Simplified98.6%
Taylor expanded in v around inf
associate-*r/N/A
/-lowering-/.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
*-lowering-*.f3255.4%
Simplified55.4%
associate-/l*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f32N/A
metadata-evalN/A
associate-/r*N/A
*-commutativeN/A
un-div-invN/A
/-lowering-/.f32N/A
metadata-evalN/A
div-invN/A
/-lowering-/.f3255.4%
Applied egg-rr55.4%
(FPCore (cosTheta_i cosTheta_O sinTheta_i sinTheta_O v) :precision binary32 (* 0.5 (* cosTheta_O (/ cosTheta_i v))))
float code(float cosTheta_i, float cosTheta_O, float sinTheta_i, float sinTheta_O, float v) {
return 0.5f * (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 = 0.5e0 * (costheta_o * (costheta_i / v))
end function
function code(cosTheta_i, cosTheta_O, sinTheta_i, sinTheta_O, v) return Float32(Float32(0.5) * Float32(cosTheta_O * Float32(cosTheta_i / v))) end
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}
\\
0.5 \cdot \left(cosTheta\_O \cdot \frac{cosTheta\_i}{v}\right)
\end{array}
Initial program 98.6%
/-lowering-/.f32N/A
exp-negN/A
associate-*l/N/A
*-lft-identityN/A
/-lowering-/.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
exp-lowering-exp.f32N/A
/-lowering-/.f32N/A
*-lowering-*.f32N/A
associate-*l*N/A
*-lowering-*.f32N/A
sinh-lowering-sinh.f32N/A
/-lowering-/.f32N/A
*-commutativeN/A
*-lowering-*.f3298.6%
Simplified98.6%
Taylor expanded in v around inf
associate-*r/N/A
/-lowering-/.f32N/A
*-commutativeN/A
*-lowering-*.f32N/A
*-lowering-*.f3255.4%
Simplified55.4%
*-commutativeN/A
associate-/l*N/A
frac-2negN/A
*-commutativeN/A
un-div-invN/A
metadata-evalN/A
frac-2negN/A
*-lowering-*.f32N/A
*-commutativeN/A
distribute-lft-neg-outN/A
frac-2negN/A
metadata-evalN/A
div-invN/A
associate-/l*N/A
distribute-neg-frac2N/A
distribute-rgt-neg-inN/A
remove-double-negN/A
*-lowering-*.f32N/A
/-lowering-/.f3255.3%
Applied egg-rr55.3%
herbie shell --seed 2024164
(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)))