\left(0.5 \cdot \sin re\right) \cdot \left(e^{0 - im} + e^{im}\right)
\left(0.5 \cdot e^{im}\right) \cdot \sin re + \sin re \cdot \frac{0.5}{e^{im}}
(FPCore (re im) :precision binary64 (* (* 0.5 (sin re)) (+ (exp (- 0.0 im)) (exp im))))
(FPCore (re im) :precision binary64 (+ (* (* 0.5 (exp im)) (sin re)) (* (sin re) (/ 0.5 (exp im)))))
double code(double re, double im) {
return (0.5 * sin(re)) * (exp(0.0 - im) + exp(im));
}
double code(double re, double im) {
return ((0.5 * exp(im)) * sin(re)) + (sin(re) * (0.5 / exp(im)));
}



Bits error versus re



Bits error versus im
Results
Initial program 0.0
Simplified0.0
Applied fma-udef_binary640.0
Applied distribute-rgt-in_binary640.0
Final simplification0.0
herbie shell --seed 2022081
(FPCore (re im)
:name "math.sin on complex, real part"
:precision binary64
(* (* 0.5 (sin re)) (+ (exp (- 0.0 im)) (exp im))))