herbie recommends `fma(h, (i + 0.5), a)`, but also doesn't report any accuracy problems with the original either

yeah, fused mul-add is certainly better. Dunno how one epxresses that in rust though :-) Ahhh seems like there at least f64::mul_add() in stdlib :-)