log1pf
extern float log1pf(float x);
return log1pf(x-1 + sqrtf((x-1)*(x-1)+2*(x-1)));
x = log1pf(x + x*x/(sqrtf(x*x+1)+1));
y = 0.5f*log1pf(2*y + 2*y*y/(1-y));
y = 0.5f*log1pf(2*(y/(1-y)));