add accurate log(1 + x) via log1p()