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.
|
|
|
|
|
|
|
|
|
|
|
|
Complex coordinates
|
|
|
|
-------------------
|
|
|
|
|
|
|
|
Complex numbers can be expressed by two important coordinate systems.
|
|
|
|
Python's :class:`complex` type uses rectangular coordinates where a number
|
|
|
|
on the complex plain is defined by two floats, the real part and the imaginary
|
|
|
|
part.
|
|
|
|
|
|
|
|
Definition::
|
|
|
|
|
|
|
|
z = x + 1j * y
|
|
|
|
|
|
|
|
x := real(z)
|
|
|
|
y := imag(z)
|
|
|
|
|
|
|
|
In engineering the polar coordinate system is popular for complex numbers. In
|
|
|
|
polar coordinates a complex number is defined by the radius *r* and the phase
|
2008-06-22 16:07:59 -03:00
|
|
|
angle *phi*. The radius *r* is the absolute value of the complex, which can be
|
2008-04-18 20:13:07 -03:00
|
|
|
viewed as distance from (0, 0). The radius *r* is always 0 or a positive float.
|
2008-06-22 16:07:59 -03:00
|
|
|
The phase angle *phi* is the counter clockwise angle from the positive x axis,
|
2008-04-18 20:13:07 -03:00
|
|
|
e.g. *1* has the angle *0*, *1j* has the angle *π/2* and *-1* the angle *-π*.
|
|
|
|
|
|
|
|
.. note::
|
|
|
|
While :func:`phase` and func:`polar` return *+π* for a negative real they
|
|
|
|
may return *-π* for a complex with a very small negative imaginary
|
|
|
|
part, e.g. *-1-1E-300j*.
|
|
|
|
|
|
|
|
|
|
|
|
Definition::
|
|
|
|
|
2008-06-22 16:07:59 -03:00
|
|
|
z = r * exp(1j * phi)
|
|
|
|
z = r * cis(phi)
|
2008-04-18 20:13:07 -03:00
|
|
|
|
|
|
|
r := abs(z) := sqrt(real(z)**2 + imag(z)**2)
|
|
|
|
phi := phase(z) := atan2(imag(z), real(z))
|
2008-06-22 16:07:59 -03:00
|
|
|
cis(phi) := cos(phi) + 1j * sin(phi)
|
2008-04-18 20:13:07 -03:00
|
|
|
|
|
|
|
|
|
|
|
.. function:: phase(x)
|
|
|
|
|
|
|
|
Return phase, also known as the argument, of a complex.
|
|
|
|
|
|
|
|
.. versionadded:: 2.6
|
|
|
|
|
|
|
|
|
|
|
|
.. function:: polar(x)
|
|
|
|
|
Merged revisions 68133-68134,68141-68142,68145-68146,68148-68149,68159-68162,68166,68171-68174,68179,68195-68196,68210,68214-68215,68217-68222 via svnmerge from
svn+ssh://pythondev@svn.python.org/python/trunk
........
r68133 | antoine.pitrou | 2009-01-01 16:38:03 +0100 (Thu, 01 Jan 2009) | 1 line
fill in actual issue number in tests
........
r68134 | hirokazu.yamamoto | 2009-01-01 16:45:39 +0100 (Thu, 01 Jan 2009) | 2 lines
Issue #4797: IOError.filename was not set when _fileio.FileIO failed to open
file with `str' filename on Windows.
........
r68141 | benjamin.peterson | 2009-01-01 17:43:12 +0100 (Thu, 01 Jan 2009) | 1 line
fix highlighting
........
r68142 | benjamin.peterson | 2009-01-01 18:29:49 +0100 (Thu, 01 Jan 2009) | 2 lines
welcome to 2009, Python!
........
r68145 | amaury.forgeotdarc | 2009-01-02 01:03:54 +0100 (Fri, 02 Jan 2009) | 5 lines
#4801 _collections module fails to build on cygwin.
_PyObject_GC_TRACK is the macro version of PyObject_GC_Track,
and according to documentation it should not be used for extension modules.
........
r68146 | ronald.oussoren | 2009-01-02 11:44:46 +0100 (Fri, 02 Jan 2009) | 2 lines
Fix for issue4472: "configure --enable-shared doesn't work on OSX"
........
r68148 | ronald.oussoren | 2009-01-02 11:48:31 +0100 (Fri, 02 Jan 2009) | 2 lines
Forgot to add a NEWS item in my previous checkin
........
r68149 | ronald.oussoren | 2009-01-02 11:50:48 +0100 (Fri, 02 Jan 2009) | 2 lines
Fix for issue4780
........
r68159 | ronald.oussoren | 2009-01-02 15:48:17 +0100 (Fri, 02 Jan 2009) | 2 lines
Fix for issue 1627952
........
r68160 | ronald.oussoren | 2009-01-02 15:52:09 +0100 (Fri, 02 Jan 2009) | 2 lines
Fix for issue r1737832
........
r68161 | ronald.oussoren | 2009-01-02 16:00:05 +0100 (Fri, 02 Jan 2009) | 3 lines
Fix for issue 1149804
........
r68162 | ronald.oussoren | 2009-01-02 16:06:00 +0100 (Fri, 02 Jan 2009) | 3 lines
Fix for issue 4472 is incompatible with Cygwin, this patch
should fix that.
........
r68166 | benjamin.peterson | 2009-01-02 19:26:23 +0100 (Fri, 02 Jan 2009) | 1 line
document PyMemberDef
........
r68171 | georg.brandl | 2009-01-02 21:25:14 +0100 (Fri, 02 Jan 2009) | 3 lines
#4811: fix markup glitches (mostly remains of the conversion),
found by Gabriel Genellina.
........
r68172 | martin.v.loewis | 2009-01-02 21:32:55 +0100 (Fri, 02 Jan 2009) | 2 lines
Issue #4075: Use OutputDebugStringW in Py_FatalError.
........
r68173 | martin.v.loewis | 2009-01-02 21:40:14 +0100 (Fri, 02 Jan 2009) | 2 lines
Issue #4051: Prevent conflict of UNICODE macros in cPickle.
........
r68174 | benjamin.peterson | 2009-01-02 21:47:27 +0100 (Fri, 02 Jan 2009) | 1 line
fix compilation on non-Windows platforms
........
r68179 | raymond.hettinger | 2009-01-02 22:26:45 +0100 (Fri, 02 Jan 2009) | 1 line
Issue #4615. Document how to use itertools for de-duping.
........
r68195 | georg.brandl | 2009-01-03 14:45:15 +0100 (Sat, 03 Jan 2009) | 2 lines
Remove useless string literal.
........
r68196 | georg.brandl | 2009-01-03 15:29:53 +0100 (Sat, 03 Jan 2009) | 2 lines
Fix indentation.
........
r68210 | georg.brandl | 2009-01-03 20:10:12 +0100 (Sat, 03 Jan 2009) | 2 lines
Set eol-style correctly for mp_distributing.py.
........
r68214 | georg.brandl | 2009-01-03 20:44:48 +0100 (Sat, 03 Jan 2009) | 2 lines
Make indentation consistent.
........
r68215 | georg.brandl | 2009-01-03 21:15:14 +0100 (Sat, 03 Jan 2009) | 2 lines
Fix role name.
........
r68217 | georg.brandl | 2009-01-03 21:30:15 +0100 (Sat, 03 Jan 2009) | 2 lines
Add rstlint, a little tool to find subtle markup problems and inconsistencies in the Doc sources.
........
r68218 | georg.brandl | 2009-01-03 21:38:59 +0100 (Sat, 03 Jan 2009) | 2 lines
Recognize usage of the default role.
........
r68219 | georg.brandl | 2009-01-03 21:47:01 +0100 (Sat, 03 Jan 2009) | 2 lines
Fix uses of the default role.
........
r68220 | georg.brandl | 2009-01-03 21:55:06 +0100 (Sat, 03 Jan 2009) | 2 lines
Remove trailing whitespace.
........
r68221 | georg.brandl | 2009-01-03 22:04:55 +0100 (Sat, 03 Jan 2009) | 2 lines
Remove tabs from the documentation.
........
r68222 | georg.brandl | 2009-01-03 22:11:58 +0100 (Sat, 03 Jan 2009) | 2 lines
Disable the line length checker by default.
........
2009-01-03 17:55:17 -04:00
|
|
|
Convert a :class:`complex` from rectangular coordinates to polar
|
2008-04-18 20:13:07 -03:00
|
|
|
coordinates. The function returns a tuple with the two elements
|
Merged revisions 68133-68134,68141-68142,68145-68146,68148-68149,68159-68162,68166,68171-68174,68179,68195-68196,68210,68214-68215,68217-68222 via svnmerge from
svn+ssh://pythondev@svn.python.org/python/trunk
........
r68133 | antoine.pitrou | 2009-01-01 16:38:03 +0100 (Thu, 01 Jan 2009) | 1 line
fill in actual issue number in tests
........
r68134 | hirokazu.yamamoto | 2009-01-01 16:45:39 +0100 (Thu, 01 Jan 2009) | 2 lines
Issue #4797: IOError.filename was not set when _fileio.FileIO failed to open
file with `str' filename on Windows.
........
r68141 | benjamin.peterson | 2009-01-01 17:43:12 +0100 (Thu, 01 Jan 2009) | 1 line
fix highlighting
........
r68142 | benjamin.peterson | 2009-01-01 18:29:49 +0100 (Thu, 01 Jan 2009) | 2 lines
welcome to 2009, Python!
........
r68145 | amaury.forgeotdarc | 2009-01-02 01:03:54 +0100 (Fri, 02 Jan 2009) | 5 lines
#4801 _collections module fails to build on cygwin.
_PyObject_GC_TRACK is the macro version of PyObject_GC_Track,
and according to documentation it should not be used for extension modules.
........
r68146 | ronald.oussoren | 2009-01-02 11:44:46 +0100 (Fri, 02 Jan 2009) | 2 lines
Fix for issue4472: "configure --enable-shared doesn't work on OSX"
........
r68148 | ronald.oussoren | 2009-01-02 11:48:31 +0100 (Fri, 02 Jan 2009) | 2 lines
Forgot to add a NEWS item in my previous checkin
........
r68149 | ronald.oussoren | 2009-01-02 11:50:48 +0100 (Fri, 02 Jan 2009) | 2 lines
Fix for issue4780
........
r68159 | ronald.oussoren | 2009-01-02 15:48:17 +0100 (Fri, 02 Jan 2009) | 2 lines
Fix for issue 1627952
........
r68160 | ronald.oussoren | 2009-01-02 15:52:09 +0100 (Fri, 02 Jan 2009) | 2 lines
Fix for issue r1737832
........
r68161 | ronald.oussoren | 2009-01-02 16:00:05 +0100 (Fri, 02 Jan 2009) | 3 lines
Fix for issue 1149804
........
r68162 | ronald.oussoren | 2009-01-02 16:06:00 +0100 (Fri, 02 Jan 2009) | 3 lines
Fix for issue 4472 is incompatible with Cygwin, this patch
should fix that.
........
r68166 | benjamin.peterson | 2009-01-02 19:26:23 +0100 (Fri, 02 Jan 2009) | 1 line
document PyMemberDef
........
r68171 | georg.brandl | 2009-01-02 21:25:14 +0100 (Fri, 02 Jan 2009) | 3 lines
#4811: fix markup glitches (mostly remains of the conversion),
found by Gabriel Genellina.
........
r68172 | martin.v.loewis | 2009-01-02 21:32:55 +0100 (Fri, 02 Jan 2009) | 2 lines
Issue #4075: Use OutputDebugStringW in Py_FatalError.
........
r68173 | martin.v.loewis | 2009-01-02 21:40:14 +0100 (Fri, 02 Jan 2009) | 2 lines
Issue #4051: Prevent conflict of UNICODE macros in cPickle.
........
r68174 | benjamin.peterson | 2009-01-02 21:47:27 +0100 (Fri, 02 Jan 2009) | 1 line
fix compilation on non-Windows platforms
........
r68179 | raymond.hettinger | 2009-01-02 22:26:45 +0100 (Fri, 02 Jan 2009) | 1 line
Issue #4615. Document how to use itertools for de-duping.
........
r68195 | georg.brandl | 2009-01-03 14:45:15 +0100 (Sat, 03 Jan 2009) | 2 lines
Remove useless string literal.
........
r68196 | georg.brandl | 2009-01-03 15:29:53 +0100 (Sat, 03 Jan 2009) | 2 lines
Fix indentation.
........
r68210 | georg.brandl | 2009-01-03 20:10:12 +0100 (Sat, 03 Jan 2009) | 2 lines
Set eol-style correctly for mp_distributing.py.
........
r68214 | georg.brandl | 2009-01-03 20:44:48 +0100 (Sat, 03 Jan 2009) | 2 lines
Make indentation consistent.
........
r68215 | georg.brandl | 2009-01-03 21:15:14 +0100 (Sat, 03 Jan 2009) | 2 lines
Fix role name.
........
r68217 | georg.brandl | 2009-01-03 21:30:15 +0100 (Sat, 03 Jan 2009) | 2 lines
Add rstlint, a little tool to find subtle markup problems and inconsistencies in the Doc sources.
........
r68218 | georg.brandl | 2009-01-03 21:38:59 +0100 (Sat, 03 Jan 2009) | 2 lines
Recognize usage of the default role.
........
r68219 | georg.brandl | 2009-01-03 21:47:01 +0100 (Sat, 03 Jan 2009) | 2 lines
Fix uses of the default role.
........
r68220 | georg.brandl | 2009-01-03 21:55:06 +0100 (Sat, 03 Jan 2009) | 2 lines
Remove trailing whitespace.
........
r68221 | georg.brandl | 2009-01-03 22:04:55 +0100 (Sat, 03 Jan 2009) | 2 lines
Remove tabs from the documentation.
........
r68222 | georg.brandl | 2009-01-03 22:11:58 +0100 (Sat, 03 Jan 2009) | 2 lines
Disable the line length checker by default.
........
2009-01-03 17:55:17 -04:00
|
|
|
*r* and *phi*. *r* is the distance from 0 and *phi* the phase
|
2008-04-18 20:13:07 -03:00
|
|
|
angle.
|
|
|
|
|
|
|
|
.. versionadded:: 2.6
|
|
|
|
|
|
|
|
|
|
|
|
.. function:: rect(r, phi)
|
|
|
|
|
|
|
|
Convert from polar coordinates to rectangular coordinates and return
|
|
|
|
a :class:`complex`.
|
|
|
|
|
|
|
|
.. versionadded:: 2.6
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
cmath functions
|
|
|
|
---------------
|
2007-08-15 11:28:01 -03:00
|
|
|
|
|
|
|
.. function:: acos(x)
|
|
|
|
|
|
|
|
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.
|
|
|
|
|
|
|
|
|
|
|
|
.. function:: acosh(x)
|
|
|
|
|
|
|
|
Return the hyperbolic arc cosine of *x*. There is one branch cut, extending left
|
|
|
|
from 1 along the real axis to -∞, continuous from above.
|
|
|
|
|
|
|
|
|
|
|
|
.. function:: asin(x)
|
|
|
|
|
|
|
|
Return the arc sine of *x*. This has the same branch cuts as :func:`acos`.
|
|
|
|
|
|
|
|
|
|
|
|
.. function:: asinh(x)
|
|
|
|
|
2008-04-18 20:13:07 -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.
|
|
|
|
|
|
|
|
.. versionchanged:: 2.6
|
|
|
|
branch cuts moved to match those recommended by the C99 standard
|
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:: atanh(x)
|
|
|
|
|
|
|
|
Return the hyperbolic arc tangent of *x*. There are two branch cuts: One
|
2008-04-18 20:13:07 -03:00
|
|
|
extends from ``1`` along the real axis to ``∞``, continuous from below. The
|
2007-08-15 11:28:01 -03:00
|
|
|
other extends from ``-1`` along the real axis to ``-∞``, continuous from
|
2008-04-18 20:13:07 -03:00
|
|
|
above.
|
|
|
|
|
|
|
|
.. versionchanged:: 2.6
|
|
|
|
direction of continuity of right cut reversed
|
2007-08-15 11:28:01 -03:00
|
|
|
|
|
|
|
|
|
|
|
.. function:: cos(x)
|
|
|
|
|
|
|
|
Return the cosine of *x*.
|
|
|
|
|
|
|
|
|
|
|
|
.. function:: cosh(x)
|
|
|
|
|
|
|
|
Return the hyperbolic cosine of *x*.
|
|
|
|
|
|
|
|
|
|
|
|
.. function:: exp(x)
|
|
|
|
|
|
|
|
Return the exponential value ``e**x``.
|
|
|
|
|
|
|
|
|
2008-04-18 20:13:07 -03:00
|
|
|
.. function:: isinf(x)
|
|
|
|
|
|
|
|
Return *True* if the real or the imaginary part of x is positive
|
|
|
|
or negative infinity.
|
|
|
|
|
|
|
|
.. versionadded:: 2.6
|
|
|
|
|
|
|
|
|
|
|
|
.. function:: isnan(x)
|
|
|
|
|
|
|
|
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
|
|
|
.. function:: log(x[, base])
|
|
|
|
|
|
|
|
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.
|
|
|
|
|
|
|
|
.. versionchanged:: 2.4
|
|
|
|
*base* argument added.
|
|
|
|
|
|
|
|
|
|
|
|
.. function:: log10(x)
|
|
|
|
|
|
|
|
Return the base-10 logarithm of *x*. This has the same branch cut as
|
|
|
|
:func:`log`.
|
|
|
|
|
|
|
|
|
|
|
|
.. function:: sin(x)
|
|
|
|
|
|
|
|
Return the sine of *x*.
|
|
|
|
|
|
|
|
|
|
|
|
.. function:: sinh(x)
|
|
|
|
|
|
|
|
Return the hyperbolic sine of *x*.
|
|
|
|
|
|
|
|
|
|
|
|
.. function:: sqrt(x)
|
|
|
|
|
|
|
|
Return the square root of *x*. This has the same branch cut as :func:`log`.
|
|
|
|
|
|
|
|
|
|
|
|
.. function:: tan(x)
|
|
|
|
|
|
|
|
Return the tangent of *x*.
|
|
|
|
|
|
|
|
|
|
|
|
.. function:: tanh(x)
|
|
|
|
|
|
|
|
Return the hyperbolic tangent of *x*.
|
|
|
|
|
|
|
|
The module also defines two mathematical constants:
|
|
|
|
|
|
|
|
|
|
|
|
.. data:: pi
|
|
|
|
|
|
|
|
The mathematical constant *pi*, as a float.
|
|
|
|
|
|
|
|
|
|
|
|
.. 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
|
|
|
|