z = geometricMean(r, 2*s);