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;