Logarithm of complex sine and cosine avoiding overflow

version (2.27 KB) by David Young
Functions to compute log(sin(Z)) and log(cos(Z)) avoiding overflow for large abs(imag(Z)).


Updated 19 Sep 2011

View License

log(cos(Z)) and log(sin(Z)) return infinities for abs(imag(Z)) > 711 in double precision, or abs(imag(Z)) > 90 in single precision. The overflow occurs in the trigonometric function; the correct final result has about the same magnitude as Z, and can be readily computed.

Using logcos(Z) or logsin(Z) instead of log(cos(Z)) or log(sin(Z)) respectively avoids the overflow, greatly extending the range of arguments which give a useful result.

The functions exploit a simple approximation; this is only applied when the approximation error is less than the precision to which numbers can be represented, so there is no loss of accuracy. Both single and double precision computation are supported.

Note that there is no point in using these functions when computing with real numbers, or with complex numbers if it is known that the imaginary part is close to zero.

Cite As

David Young (2022). Logarithm of complex sine and cosine avoiding overflow (https://www.mathworks.com/matlabcentral/fileexchange/32947-logarithm-of-complex-sine-and-cosine-avoiding-overflow), MATLAB Central File Exchange. Retrieved .

MATLAB Release Compatibility
Created with R2011a
Compatible with any release
Platform Compatibility
Windows macOS Linux

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!