sqrt(x*x+y*y) = hypot(x,y)
log(1+x) = log1p(x)
expm1(x) = exp x - 1