2007-08-15 11:28:01 -03:00
|
|
|
:mod:`cmath` --- Mathematical functions for complex numbers
|
|
|
|
===========================================================
|
|
|
|
|
|
|
|
.. module:: cmath
|
|
|
|
:synopsis: Mathematical functions for complex numbers.
|
|
|
|
|
|
|
|
|
|
|
|
This module is always available. It provides access to mathematical functions
|
|
|
|
for complex numbers. The functions in this module accept integers,
|
|
|
|
floating-point numbers or complex numbers as arguments. They will also accept
|
|
|
|
any Python object that has either a :meth:`__complex__` or a :meth:`__float__`
|
|
|
|
method: these methods are used to convert the object to a complex or
|
|
|
|
floating-point number, respectively, and the function is then applied to the
|
|
|
|
result of the conversion.
|
|
|
|
|
2008-04-18 20:13:07 -03:00
|
|
|
.. note::
|
2007-08-15 11:28:01 -03:00
|
|
|
|
2008-04-18 20:13:07 -03:00
|
|
|
On platforms with hardware and system-level support for signed
|
|
|
|
zeros, functions involving branch cuts are continuous on *both*
|
|
|
|
sides of the branch cut: the sign of the zero distinguishes one
|
|
|
|
side of the branch cut from the other. On platforms that do not
|
|
|
|
support signed zeros the continuity is as specified below.
|
|
|
|
|
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
Conversions to and from polar coordinates
|
|
|
|
-----------------------------------------
|
2008-04-18 20:13:07 -03:00
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
A Python complex number ``z`` is stored internally using *rectangular*
|
|
|
|
or *Cartesian* coordinates. It is completely determined by its *real
|
|
|
|
part* ``z.real`` and its *imaginary part* ``z.imag``. In other
|
|
|
|
words::
|
2008-04-18 20:13:07 -03:00
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
z == z.real + z.imag*1j
|
2008-04-18 20:13:07 -03:00
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
*Polar coordinates* give an alternative way to represent a complex
|
|
|
|
number. In polar coordinates, a complex number *z* is defined by the
|
|
|
|
modulus *r* and the phase angle *phi*. The modulus *r* is the distance
|
|
|
|
from *z* to the origin, while the phase *phi* is the counterclockwise
|
2010-01-02 10:31:13 -04:00
|
|
|
angle, measured in radians, from the positive x-axis to the line
|
|
|
|
segment that joins the origin to *z*.
|
2008-04-18 20:13:07 -03:00
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
The following functions can be used to convert from the native
|
|
|
|
rectangular coordinates to polar coordinates and back.
|
2008-04-18 20:13:07 -03:00
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
.. function:: phase(x)
|
2008-04-18 20:13:07 -03:00
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
Return the phase of *x* (also known as the *argument* of *x*), as a
|
|
|
|
float. ``phase(x)`` is equivalent to ``math.atan2(x.imag,
|
|
|
|
x.real)``. The result lies in the range [-π, π], and the branch
|
|
|
|
cut for this operation lies along the negative real axis,
|
|
|
|
continuous from above. On systems with support for signed zeros
|
|
|
|
(which includes most systems in current use), this means that the
|
|
|
|
sign of the result is the same as the sign of ``x.imag``, even when
|
|
|
|
``x.imag`` is zero::
|
2008-04-18 20:13:07 -03:00
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
>>> phase(complex(-1.0, 0.0))
|
|
|
|
3.1415926535897931
|
|
|
|
>>> phase(complex(-1.0, -0.0))
|
|
|
|
-3.1415926535897931
|
|
|
|
|
|
|
|
.. versionadded:: 2.6
|
2008-04-18 20:13:07 -03:00
|
|
|
|
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
.. note::
|
2008-04-18 20:13:07 -03:00
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
The modulus (absolute value) of a complex number *x* can be
|
|
|
|
computed using the built-in :func:`abs` function. There is no
|
|
|
|
separate :mod:`cmath` module function for this operation.
|
2008-04-18 20:13:07 -03:00
|
|
|
|
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
.. function:: polar(x)
|
2008-04-18 20:13:07 -03:00
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
Return the representation of *x* in polar coordinates. Returns a
|
|
|
|
pair ``(r, phi)`` where *r* is the modulus of *x* and phi is the
|
|
|
|
phase of *x*. ``polar(x)`` is equivalent to ``(abs(x),
|
|
|
|
phase(x))``.
|
2008-04-18 20:13:07 -03:00
|
|
|
|
|
|
|
.. versionadded:: 2.6
|
|
|
|
|
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
.. function:: rect(r, phi)
|
2008-04-18 20:13:07 -03:00
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
Return the complex number *x* with polar coordinates *r* and *phi*.
|
|
|
|
Equivalent to ``r * (math.cos(phi) + math.sin(phi)*1j)``.
|
2008-04-18 20:13:07 -03:00
|
|
|
|
|
|
|
.. versionadded:: 2.6
|
|
|
|
|
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
Power and logarithmic functions
|
|
|
|
-------------------------------
|
2008-04-18 20:13:07 -03:00
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
.. function:: exp(x)
|
2008-04-18 20:13:07 -03:00
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
Return the exponential value ``e**x``.
|
2008-04-18 20:13:07 -03:00
|
|
|
|
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
.. function:: log(x[, base])
|
2008-04-18 20:13:07 -03:00
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
Returns the logarithm of *x* to the given *base*. If the *base* is not
|
|
|
|
specified, returns the natural logarithm of *x*. There is one branch cut, from 0
|
|
|
|
along the negative real axis to -∞, continuous from above.
|
2007-08-15 11:28:01 -03:00
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
.. versionchanged:: 2.4
|
|
|
|
*base* argument added.
|
2007-08-15 11:28:01 -03:00
|
|
|
|
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
.. function:: log10(x)
|
2007-08-15 11:28:01 -03:00
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
Return the base-10 logarithm of *x*. This has the same branch cut as
|
|
|
|
:func:`log`.
|
2007-08-15 11:28:01 -03:00
|
|
|
|
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
.. function:: sqrt(x)
|
2007-08-15 11:28:01 -03:00
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
Return the square root of *x*. This has the same branch cut as :func:`log`.
|
2007-08-15 11:28:01 -03:00
|
|
|
|
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
Trigonometric functions
|
|
|
|
-----------------------
|
2007-08-15 11:28:01 -03:00
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
.. function:: acos(x)
|
2007-08-15 11:28:01 -03:00
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
Return the arc cosine of *x*. There are two branch cuts: One extends right from
|
|
|
|
1 along the real axis to ∞, continuous from below. The other extends left from
|
|
|
|
-1 along the real axis to -∞, continuous from above.
|
2008-04-18 20:13:07 -03:00
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
|
|
|
|
.. function:: asin(x)
|
|
|
|
|
|
|
|
Return the arc sine of *x*. This has the same branch cuts as :func:`acos`.
|
2007-08-15 11:28:01 -03:00
|
|
|
|
|
|
|
|
|
|
|
.. function:: atan(x)
|
|
|
|
|
|
|
|
Return the arc tangent of *x*. There are two branch cuts: One extends from
|
2008-04-18 20:13:07 -03:00
|
|
|
``1j`` along the imaginary axis to ``∞j``, continuous from the right. The
|
2007-08-15 11:28:01 -03:00
|
|
|
other extends from ``-1j`` along the imaginary axis to ``-∞j``, continuous
|
2008-04-18 20:13:07 -03:00
|
|
|
from the left.
|
|
|
|
|
|
|
|
.. versionchanged:: 2.6
|
|
|
|
direction of continuity of upper cut reversed
|
2007-08-15 11:28:01 -03:00
|
|
|
|
|
|
|
|
|
|
|
.. function:: cos(x)
|
|
|
|
|
|
|
|
Return the cosine of *x*.
|
|
|
|
|
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
.. function:: sin(x)
|
2007-08-15 11:28:01 -03:00
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
Return the sine of *x*.
|
2007-08-15 11:28:01 -03:00
|
|
|
|
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
.. function:: tan(x)
|
2007-08-15 11:28:01 -03:00
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
Return the tangent of *x*.
|
2007-08-15 11:28:01 -03:00
|
|
|
|
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
Hyperbolic functions
|
|
|
|
--------------------
|
2008-04-18 20:13:07 -03:00
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
.. function:: acosh(x)
|
2008-04-18 20:13:07 -03:00
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
Return the hyperbolic arc cosine of *x*. There is one branch cut, extending left
|
|
|
|
from 1 along the real axis to -∞, continuous from above.
|
2008-04-18 20:13:07 -03:00
|
|
|
|
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
.. function:: asinh(x)
|
2008-04-18 20:13:07 -03:00
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
Return the hyperbolic arc sine of *x*. There are two branch cuts:
|
|
|
|
One extends from ``1j`` along the imaginary axis to ``∞j``,
|
|
|
|
continuous from the right. The other extends from ``-1j`` along
|
|
|
|
the imaginary axis to ``-∞j``, continuous from the left.
|
2008-04-18 20:13:07 -03:00
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
.. versionchanged:: 2.6
|
|
|
|
branch cuts moved to match those recommended by the C99 standard
|
2008-04-18 20:13:07 -03:00
|
|
|
|
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
.. function:: atanh(x)
|
2007-08-15 11:28:01 -03:00
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
Return the hyperbolic arc tangent of *x*. There are two branch cuts: One
|
|
|
|
extends from ``1`` along the real axis to ``∞``, continuous from below. The
|
|
|
|
other extends from ``-1`` along the real axis to ``-∞``, continuous from
|
|
|
|
above.
|
2007-08-15 11:28:01 -03:00
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
.. versionchanged:: 2.6
|
|
|
|
direction of continuity of right cut reversed
|
2007-08-15 11:28:01 -03:00
|
|
|
|
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
.. function:: cosh(x)
|
2007-08-15 11:28:01 -03:00
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
Return the hyperbolic cosine of *x*.
|
2007-08-15 11:28:01 -03:00
|
|
|
|
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
.. function:: sinh(x)
|
2007-08-15 11:28:01 -03:00
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
Return the hyperbolic sine of *x*.
|
2007-08-15 11:28:01 -03:00
|
|
|
|
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
.. function:: tanh(x)
|
2007-08-15 11:28:01 -03:00
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
Return the hyperbolic tangent of *x*.
|
2007-08-15 11:28:01 -03:00
|
|
|
|
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
Classification functions
|
|
|
|
------------------------
|
2007-08-15 11:28:01 -03:00
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
.. function:: isinf(x)
|
2007-08-15 11:28:01 -03:00
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
Return *True* if the real or the imaginary part of x is positive
|
|
|
|
or negative infinity.
|
2007-08-15 11:28:01 -03:00
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
.. versionadded:: 2.6
|
2007-08-15 11:28:01 -03:00
|
|
|
|
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
.. function:: isnan(x)
|
2007-08-15 11:28:01 -03:00
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
Return *True* if the real or imaginary part of x is not a number (NaN).
|
|
|
|
|
|
|
|
.. versionadded:: 2.6
|
2007-08-15 11:28:01 -03:00
|
|
|
|
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
Constants
|
|
|
|
---------
|
2007-08-15 11:28:01 -03:00
|
|
|
|
|
|
|
|
|
|
|
.. data:: pi
|
|
|
|
|
2009-07-28 13:15:25 -03:00
|
|
|
The mathematical constant *π*, as a float.
|
2007-08-15 11:28:01 -03:00
|
|
|
|
|
|
|
|
|
|
|
.. data:: e
|
|
|
|
|
|
|
|
The mathematical constant *e*, as a float.
|
|
|
|
|
|
|
|
.. index:: module: math
|
|
|
|
|
|
|
|
Note that the selection of functions is similar, but not identical, to that in
|
|
|
|
module :mod:`math`. The reason for having two modules is that some users aren't
|
|
|
|
interested in complex numbers, and perhaps don't even know what they are. They
|
|
|
|
would rather have ``math.sqrt(-1)`` raise an exception than return a complex
|
|
|
|
number. Also note that the functions defined in :mod:`cmath` always return a
|
|
|
|
complex number, even if the answer can be expressed as a real number (in which
|
|
|
|
case the complex number has an imaginary part of zero).
|
|
|
|
|
|
|
|
A note on branch cuts: They are curves along which the given function fails to
|
|
|
|
be continuous. They are a necessary feature of many complex functions. It is
|
|
|
|
assumed that if you need to compute with complex functions, you will understand
|
|
|
|
about branch cuts. Consult almost any (not too elementary) book on complex
|
|
|
|
variables for enlightenment. For information of the proper choice of branch
|
|
|
|
cuts for numerical purposes, a good reference should be the following:
|
|
|
|
|
|
|
|
|
|
|
|
.. seealso::
|
|
|
|
|
|
|
|
Kahan, W: Branch cuts for complex elementary functions; or, Much ado about
|
|
|
|
nothing's sign bit. In Iserles, A., and Powell, M. (eds.), The state of the art
|
|
|
|
in numerical analysis. Clarendon Press (1987) pp165-211.
|
|
|
|
|
2008-04-18 20:13:07 -03:00
|
|
|
|