sum_squares
rx = log1p(4 * ax / sum_squares(ax - 1, ay)) / 4;
rx = log1pf(4 * ax / sum_squares(ax - 1, ay)) / 4;
rx = log1pl(4 * ax / sum_squares(ax - 1, ay)) / 4;