2007-08-15 11:28:22 -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 21:31:39 -03:00
|
|
|
.. note::
|
2007-08-15 11:28:22 -03:00
|
|
|
|
2008-04-18 21:31:39 -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:31:03 -03:00
|
|
|
Conversions to and from polar coordinates
|
|
|
|
-----------------------------------------
|
2008-04-18 21:31:39 -03:00
|
|
|
|
2009-07-28 13:31:03 -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 21:31:39 -03:00
|
|
|
|
2009-07-28 13:31:03 -03:00
|
|
|
z == z.real + z.imag*1j
|
2008-04-18 21:31:39 -03:00
|
|
|
|
2009-07-28 13:31:03 -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:33:10 -04:00
|
|
|
angle, measured in radians, from the positive x-axis to the line
|
|
|
|
segment that joins the origin to *z*.
|
2008-04-18 21:31:39 -03:00
|
|
|
|
2009-07-28 13:31:03 -03:00
|
|
|
The following functions can be used to convert from the native
|
|
|
|
rectangular coordinates to polar coordinates and back.
|
|
|
|
|
|
|
|
.. function:: phase(x)
|
|
|
|
|
|
|
|
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::
|
|
|
|
|
|
|
|
>>> phase(complex(-1.0, 0.0))
|
|
|
|
3.141592653589793
|
|
|
|
>>> phase(complex(-1.0, -0.0))
|
|
|
|
-3.141592653589793
|
2008-04-18 21:31:39 -03:00
|
|
|
|
|
|
|
|
|
|
|
.. note::
|
|
|
|
|
2009-07-28 13:31:03 -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 21:31:39 -03:00
|
|
|
|
|
|
|
|
2009-07-28 13:31:03 -03:00
|
|
|
.. function:: polar(x)
|
2008-04-18 21:31:39 -03:00
|
|
|
|
2009-07-28 13:31:03 -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 21:31:39 -03:00
|
|
|
|
|
|
|
|
2009-07-28 13:31:03 -03:00
|
|
|
.. function:: rect(r, phi)
|
2008-04-18 21:31:39 -03:00
|
|
|
|
2009-07-28 13:31:03 -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 21:31:39 -03:00
|
|
|
|
|
|
|
|
2009-07-28 13:31:03 -03:00
|
|
|
Power and logarithmic functions
|
|
|
|
-------------------------------
|
2008-04-18 21:31:39 -03:00
|
|
|
|
2009-07-28 13:31:03 -03:00
|
|
|
.. function:: exp(x)
|
2008-04-18 21:31:39 -03:00
|
|
|
|
2009-07-28 13:31:03 -03:00
|
|
|
Return the exponential value ``e**x``.
|
2008-04-18 21:31:39 -03:00
|
|
|
|
|
|
|
|
2009-07-28 13:31:03 -03:00
|
|
|
.. function:: log(x[, base])
|
2008-04-18 21:31:39 -03:00
|
|
|
|
2009-07-28 13:31:03 -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.
|
2008-04-18 21:31:39 -03:00
|
|
|
|
2007-08-15 11:28:22 -03:00
|
|
|
|
2009-07-28 13:31:03 -03:00
|
|
|
.. function:: log10(x)
|
2007-08-15 11:28:22 -03:00
|
|
|
|
2009-07-28 13:31:03 -03:00
|
|
|
Return the base-10 logarithm of *x*. This has the same branch cut as
|
|
|
|
:func:`log`.
|
2007-08-15 11:28:22 -03:00
|
|
|
|
|
|
|
|
2009-07-28 13:31:03 -03:00
|
|
|
.. function:: sqrt(x)
|
2007-08-15 11:28:22 -03:00
|
|
|
|
2009-07-28 13:31:03 -03:00
|
|
|
Return the square root of *x*. This has the same branch cut as :func:`log`.
|
2007-08-15 11:28:22 -03:00
|
|
|
|
|
|
|
|
2009-07-28 13:31:03 -03:00
|
|
|
Trigonometric functions
|
|
|
|
-----------------------
|
2007-08-15 11:28:22 -03:00
|
|
|
|
2009-07-28 13:31:03 -03:00
|
|
|
.. function:: acos(x)
|
2007-08-15 11:28:22 -03:00
|
|
|
|
2009-07-28 13:31:03 -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.
|
2007-08-15 11:28:22 -03:00
|
|
|
|
|
|
|
|
2009-07-28 13:31:03 -03:00
|
|
|
.. function:: asin(x)
|
|
|
|
|
|
|
|
Return the arc sine of *x*. This has the same branch cuts as :func:`acos`.
|
2008-04-18 21:31:39 -03:00
|
|
|
|
2007-08-15 11:28:22 -03:00
|
|
|
|
|
|
|
.. function:: atan(x)
|
|
|
|
|
|
|
|
Return the arc tangent of *x*. There are two branch cuts: One extends from
|
2008-04-18 21:31:39 -03:00
|
|
|
``1j`` along the imaginary axis to ``∞j``, continuous from the right. The
|
2007-08-15 11:28:22 -03:00
|
|
|
other extends from ``-1j`` along the imaginary axis to ``-∞j``, continuous
|
2008-04-18 21:31:39 -03:00
|
|
|
from the left.
|
|
|
|
|
2007-08-15 11:28:22 -03:00
|
|
|
|
|
|
|
.. function:: cos(x)
|
|
|
|
|
|
|
|
Return the cosine of *x*.
|
|
|
|
|
|
|
|
|
2009-07-28 13:31:03 -03:00
|
|
|
.. function:: sin(x)
|
2007-08-15 11:28:22 -03:00
|
|
|
|
2009-07-28 13:31:03 -03:00
|
|
|
Return the sine of *x*.
|
2007-08-15 11:28:22 -03:00
|
|
|
|
|
|
|
|
2009-07-28 13:31:03 -03:00
|
|
|
.. function:: tan(x)
|
2007-08-15 11:28:22 -03:00
|
|
|
|
2009-07-28 13:31:03 -03:00
|
|
|
Return the tangent of *x*.
|
2008-04-18 21:31:39 -03:00
|
|
|
|
|
|
|
|
2009-07-28 13:31:03 -03:00
|
|
|
Hyperbolic functions
|
|
|
|
--------------------
|
2008-04-18 21:31:39 -03:00
|
|
|
|
2009-07-28 13:31:03 -03:00
|
|
|
.. function:: acosh(x)
|
2008-04-18 21:31:39 -03:00
|
|
|
|
2015-01-11 09:08:05 -04:00
|
|
|
Return the inverse hyperbolic cosine of *x*. There is one branch cut,
|
|
|
|
extending left from 1 along the real axis to -∞, continuous from above.
|
2008-04-18 21:31:39 -03:00
|
|
|
|
|
|
|
|
2009-07-28 13:31:03 -03:00
|
|
|
.. function:: asinh(x)
|
2007-08-15 11:28:22 -03:00
|
|
|
|
2015-01-11 09:08:05 -04:00
|
|
|
Return the inverse hyperbolic sine of *x*. There are two branch cuts:
|
2009-07-28 13:31:03 -03:00
|
|
|
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.
|
2007-08-15 11:28:22 -03:00
|
|
|
|
|
|
|
|
2009-07-28 13:31:03 -03:00
|
|
|
.. function:: atanh(x)
|
2007-08-15 11:28:22 -03:00
|
|
|
|
2015-01-11 09:08:05 -04:00
|
|
|
Return the inverse hyperbolic tangent of *x*. There are two branch cuts: One
|
2009-07-28 13:31:03 -03:00
|
|
|
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:22 -03:00
|
|
|
|
|
|
|
|
2009-07-28 13:31:03 -03:00
|
|
|
.. function:: cosh(x)
|
2007-08-15 11:28:22 -03:00
|
|
|
|
2009-07-28 13:31:03 -03:00
|
|
|
Return the hyperbolic cosine of *x*.
|
2007-08-15 11:28:22 -03:00
|
|
|
|
|
|
|
|
|
|
|
.. function:: sinh(x)
|
|
|
|
|
|
|
|
Return the hyperbolic sine of *x*.
|
|
|
|
|
|
|
|
|
2009-07-28 13:31:03 -03:00
|
|
|
.. function:: tanh(x)
|
2007-08-15 11:28:22 -03:00
|
|
|
|
2009-07-28 13:31:03 -03:00
|
|
|
Return the hyperbolic tangent of *x*.
|
2007-08-15 11:28:22 -03:00
|
|
|
|
|
|
|
|
2009-07-28 13:31:03 -03:00
|
|
|
Classification functions
|
|
|
|
------------------------
|
2007-08-15 11:28:22 -03:00
|
|
|
|
2010-07-11 14:38:24 -03:00
|
|
|
.. function:: isfinite(x)
|
|
|
|
|
2010-07-11 16:47:37 -03:00
|
|
|
Return ``True`` if both the real and imaginary parts of *x* are finite, and
|
|
|
|
``False`` otherwise.
|
|
|
|
|
|
|
|
.. versionadded:: 3.2
|
2010-07-11 14:38:24 -03:00
|
|
|
|
|
|
|
|
2009-07-28 13:31:03 -03:00
|
|
|
.. function:: isinf(x)
|
|
|
|
|
2010-07-11 16:47:37 -03:00
|
|
|
Return ``True`` if either the real or the imaginary part of *x* is an
|
|
|
|
infinity, and ``False`` otherwise.
|
2007-08-15 11:28:22 -03:00
|
|
|
|
|
|
|
|
2009-07-28 13:31:03 -03:00
|
|
|
.. function:: isnan(x)
|
|
|
|
|
2010-07-11 16:47:37 -03:00
|
|
|
Return ``True`` if either the real or the imaginary part of *x* is a NaN,
|
|
|
|
and ``False`` otherwise.
|
2007-08-15 11:28:22 -03:00
|
|
|
|
|
|
|
|
2015-05-31 16:05:00 -03:00
|
|
|
.. function:: isclose(a, b, *, rel_tol=1e-09, abs_tol=0.0)
|
|
|
|
|
|
|
|
Return ``True`` if the values *a* and *b* are close to each other and
|
|
|
|
``False`` otherwise.
|
|
|
|
|
|
|
|
Whether or not two values are considered close is determined according to
|
|
|
|
given absolute and relative tolerances.
|
|
|
|
|
|
|
|
*rel_tol* is the relative tolerance -- it is the maximum allowed difference
|
|
|
|
between *a* and *b*, relative to the larger absolute value of *a* or *b*.
|
|
|
|
For example, to set a tolerance of 5%, pass ``rel_tol=0.05``. The default
|
|
|
|
tolerance is ``1e-09``, which assures that the two values are the same
|
|
|
|
within about 9 decimal digits. *rel_tol* must be greater than zero.
|
|
|
|
|
|
|
|
*abs_tol* is the minimum absolute tolerance -- useful for comparisons near
|
|
|
|
zero. *abs_tol* must be at least zero.
|
|
|
|
|
|
|
|
If no errors occur, the result will be:
|
|
|
|
``abs(a-b) <= max(rel_tol * max(abs(a), abs(b)), abs_tol)``.
|
|
|
|
|
|
|
|
The IEEE 754 special values of ``NaN``, ``inf``, and ``-inf`` will be
|
|
|
|
handled according to IEEE rules. Specifically, ``NaN`` is not considered
|
|
|
|
close to any other value, including ``NaN``. ``inf`` and ``-inf`` are only
|
|
|
|
considered close to themselves.
|
|
|
|
|
|
|
|
.. versionadded:: 3.5
|
|
|
|
|
|
|
|
.. seealso::
|
|
|
|
|
|
|
|
:pep:`485` -- A function for testing approximate equality
|
|
|
|
|
|
|
|
|
2009-07-28 13:31:03 -03:00
|
|
|
Constants
|
|
|
|
---------
|
2007-08-15 11:28:22 -03:00
|
|
|
|
|
|
|
|
|
|
|
.. data:: pi
|
|
|
|
|
2009-07-28 13:31:03 -03:00
|
|
|
The mathematical constant *π*, as a float.
|
2007-08-15 11:28:22 -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 21:31:39 -03:00
|
|
|
|