tmp2 = 0.5 / g * (w1 * w1 + w2 * w2 - w1 * w2)  # OK*/