logsum

Functions for computing log(a + b) and log(a - b) when direct computation is numerically unstable