Displaying 1 result from an estimated 1 matches for "logspace_sup".
Did you mean:
logspace_sub
2009 Apr 27
0
Patch proposal for logspace_sub
...om the logarithms of the individual terms; so the
functions logspace_add and logspace_sub which are part of R's API come
in handy.
However, I noticed that logspace_sub can have problems if both
arguments are (very) small or the difference between the arguments are
vary small. The logic behind logspace_sup, when called with arguments
lx and ly, is:
log( exp(lx) - exp(ly) ) = log( exp(lx) * ( 1 - exp(ly - lx) ) )
= lx + log( 1 - exp(ly - lx) )
= lx + log1p( - exp(ly - lx) )
and it is the last expression that is evaluated by logspace_sub.
However the...