Disclosed are methods to compute the precise value of the logarithm of a sum and the logarithm of a subtraction.
A number of practical problems can result in having too big or small values in intermediate values of a calculation. Then one tries to take logarithm of these values and operate on logarithms instead.
In the case: log (p • q) = log p + log q is very easy to compute, but the problem is then to compute (or approximate) the value of log (a + b) from the value of log a and log b.
Let us assume that log a and log b are known, and that we want to approximate log (a + b).
Most basic solution would be calculating
sum_log = log(exp(a_log) + exp(b_log),
where a_log = log a and b_log = log b and therefore sum_log = log(a + b).
But the method I propose requires calling only one exp() and one log(), instead of two exp() and one log() in the basic solution.
Additionally, the proposed method has the critical advantage of not overflowing in case of large numbers of a and b.
The method is based on the notion that
ln(a + b) = ln{exp[ln(a)  ln(b)] + 1} + ln(b).
In a similar manner we can compute the logarithm of subtraction, based on the notion that
ln(a  b) = ln{exp[ln(a)  ln(b)]  1} + ln(b), if a > b.
Roland Pihlakas (2020). Methods for calculating precise logarithm of a sum and subtraction (https://www.mathworks.com/matlabcentral/fileexchange/25273methodsforcalculatingpreciselogarithmofasumandsubtraction), MATLAB Central File Exchange. Retrieved .
1.3.0.0  Seems like the updated zip file did not upload with my previous post. I will now resend it. 

1.2.0.0  1) Included method for calculating logarithm of subtraction.

Create scripts with code, output, and formatted text in a single executable document.
Darren Rowland (view profile)
I can see this being a useful utility for some applications. However the code could be improved. I would like to see you add
 a H1 line
 an example (or two) showing usage, particularly comparing your method with a method which will produce overflow
 Comments within the code explaining the procedure so a user does not need to consult the pdf file
Lastly I think the actual computation should be on just one line
i.e. R = log(exp(a_ln  b_ln) + 1) + b_ln;
rather than spread out over many lines. Let comments explain the shifting.
Darren