wehub-resource-sync
8 lines
209 B
Matlab
Executable File
8 lines
209 B
Matlab
Executable File
function y = log1mexp(x)
|
|
% Accurately compute y = log(1-exp(x))
|
|
% reference: Accurately Computing log(1-exp(-|a|)) Martin Machler
|
|
y = x;
|
|
i = x < -log(2);
|
|
y(i) = log1p(-exp(x(i)));
|
|
y(~i) = log(-expm1(x(~i)));
|