I'm going to present a work on numerical stability of the acoustic terms in #SPH, and as I'm finalizing the conference paper with some benchmarks, I've come across an interesting situation. The cornerstone of the work involves replacing a powf() function with a combination of exp1mf() and logp1f() for higher accuracy.
How much of a performance loss does this higher accuracy involve?
And here's where things get weird.