log1pl
extern long double log1pl(long double x);
return log1pl(x-1 + sqrtl((x-1)*(x-1)+2*(x-1)));
x = log1pl(x + x*x/(sqrtl(x*x+1)+1));
x = 0.5*log1pl(2*x + 2*x*x/(1-x));
x = 0.5*log1pl(2*(x/(1-x)));