log1pf
float log1pf(float);
extern float log1pf(float);
F_IM(ans) = half * log1pf((ay + ay) / t);
F_IM(ans) = half * log1pf(two / (-t));