\left(0.5 \cdot \sin re\right) \cdot \left(e^{0.0 - im} + e^{im}\right)\mathsf{fma}\left(0.5 \cdot \sin re, e^{0.0 - im}, \left(0.5 \cdot \sin re\right) \cdot e^{im}\right)double code(double re, double im) {
return ((0.5 * sin(re)) * (exp((0.0 - im)) + exp(im)));
}
double code(double re, double im) {
return fma((0.5 * sin(re)), exp((0.0 - im)), ((0.5 * sin(re)) * exp(im)));
}



Bits error versus re



Bits error versus im
Results
Initial program 0.1
rmApplied distribute-lft-in0.1
rmApplied fma-def0.1
Final simplification0.1
herbie shell --seed 2020078 +o rules:numerics
(FPCore (re im)
:name "math.sin on complex, real part"
:precision binary64
(* (* 0.5 (sin re)) (+ (exp (- 0.0 im)) (exp im))))