| 
									
										
										
										
											2007-08-15 14:28:22 +00: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-19 00:31:39 +00:00
										 |  |  | .. note::
 | 
					
						
							| 
									
										
										
										
											2007-08-15 14:28:22 +00:00
										 |  |  | 
 | 
					
						
							| 
									
										
										
										
											2008-04-19 00:31:39 +00: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 16:31:03 +00:00
										 |  |  | Conversions to and from polar coordinates
 | 
					
						
							|  |  |  | -----------------------------------------
 | 
					
						
							| 
									
										
										
										
											2008-04-19 00:31:39 +00:00
										 |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00: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-19 00:31:39 +00:00
										 |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00:00
										 |  |  |    z == z.real + z.imag*1j
 | 
					
						
							| 
									
										
										
										
											2008-04-19 00:31:39 +00:00
										 |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00: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 14:33:10 +00:00
										 |  |  | angle, measured in radians, from the positive x-axis to the line
 | 
					
						
							|  |  |  | segment that joins the origin to *z*.
 | 
					
						
							| 
									
										
										
										
											2008-04-19 00:31:39 +00:00
										 |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00: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-19 00:31:39 +00:00
										 |  |  | 
 | 
					
						
							|  |  |  | 
 | 
					
						
							|  |  |  | .. note::
 | 
					
						
							|  |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00: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-19 00:31:39 +00:00
										 |  |  | 
 | 
					
						
							|  |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00:00
										 |  |  | .. function:: polar(x)
 | 
					
						
							| 
									
										
										
										
											2008-04-19 00:31:39 +00:00
										 |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00: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-19 00:31:39 +00:00
										 |  |  | 
 | 
					
						
							|  |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00:00
										 |  |  | .. function:: rect(r, phi)
 | 
					
						
							| 
									
										
										
										
											2008-04-19 00:31:39 +00:00
										 |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00:00
										 |  |  |    Return the complex number *x* with polar coordinates *r* and *phi*.
 | 
					
						
							|  |  |  |    Equivalent to ``r * (math.cos(phi) + math.sin(phi)*1j)``.
 | 
					
						
							| 
									
										
										
										
											2008-04-19 00:31:39 +00:00
										 |  |  | 
 | 
					
						
							|  |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00:00
										 |  |  | Power and logarithmic functions
 | 
					
						
							|  |  |  | -------------------------------
 | 
					
						
							| 
									
										
										
										
											2008-04-19 00:31:39 +00:00
										 |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00:00
										 |  |  | .. function:: exp(x)
 | 
					
						
							| 
									
										
										
										
											2008-04-19 00:31:39 +00:00
										 |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00:00
										 |  |  |    Return the exponential value ``e**x``.
 | 
					
						
							| 
									
										
										
										
											2008-04-19 00:31:39 +00:00
										 |  |  | 
 | 
					
						
							|  |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00:00
										 |  |  | .. function:: log(x[, base])
 | 
					
						
							| 
									
										
										
										
											2008-04-19 00:31:39 +00:00
										 |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00: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-19 00:31:39 +00:00
										 |  |  | 
 | 
					
						
							| 
									
										
										
										
											2007-08-15 14:28:22 +00:00
										 |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00:00
										 |  |  | .. function:: log10(x)
 | 
					
						
							| 
									
										
										
										
											2007-08-15 14:28:22 +00:00
										 |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00:00
										 |  |  |    Return the base-10 logarithm of *x*. This has the same branch cut as
 | 
					
						
							|  |  |  |    :func:`log`.
 | 
					
						
							| 
									
										
										
										
											2007-08-15 14:28:22 +00:00
										 |  |  | 
 | 
					
						
							|  |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00:00
										 |  |  | .. function:: sqrt(x)
 | 
					
						
							| 
									
										
										
										
											2007-08-15 14:28:22 +00:00
										 |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00:00
										 |  |  |    Return the square root of *x*. This has the same branch cut as :func:`log`.
 | 
					
						
							| 
									
										
										
										
											2007-08-15 14:28:22 +00:00
										 |  |  | 
 | 
					
						
							|  |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00:00
										 |  |  | Trigonometric functions
 | 
					
						
							|  |  |  | -----------------------
 | 
					
						
							| 
									
										
										
										
											2007-08-15 14:28:22 +00:00
										 |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00:00
										 |  |  | .. function:: acos(x)
 | 
					
						
							| 
									
										
										
										
											2007-08-15 14:28:22 +00:00
										 |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00: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 14:28:22 +00:00
										 |  |  | 
 | 
					
						
							|  |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00:00
										 |  |  | .. function:: asin(x)
 | 
					
						
							|  |  |  | 
 | 
					
						
							|  |  |  |    Return the arc sine of *x*. This has the same branch cuts as :func:`acos`.
 | 
					
						
							| 
									
										
										
										
											2008-04-19 00:31:39 +00:00
										 |  |  | 
 | 
					
						
							| 
									
										
										
										
											2007-08-15 14:28:22 +00:00
										 |  |  | 
 | 
					
						
							|  |  |  | .. function:: atan(x)
 | 
					
						
							|  |  |  | 
 | 
					
						
							|  |  |  |    Return the arc tangent of *x*. There are two branch cuts: One extends from
 | 
					
						
							| 
									
										
										
										
											2008-04-19 00:31:39 +00:00
										 |  |  |    ``1j`` along the imaginary axis to ``∞j``, continuous from the right. The
 | 
					
						
							| 
									
										
										
										
											2007-08-15 14:28:22 +00:00
										 |  |  |    other extends from ``-1j`` along the imaginary axis to ``-∞j``, continuous
 | 
					
						
							| 
									
										
										
										
											2008-04-19 00:31:39 +00:00
										 |  |  |    from the left.
 | 
					
						
							|  |  |  | 
 | 
					
						
							| 
									
										
										
										
											2007-08-15 14:28:22 +00:00
										 |  |  | 
 | 
					
						
							|  |  |  | .. function:: cos(x)
 | 
					
						
							|  |  |  | 
 | 
					
						
							|  |  |  |    Return the cosine of *x*.
 | 
					
						
							|  |  |  | 
 | 
					
						
							|  |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00:00
										 |  |  | .. function:: sin(x)
 | 
					
						
							| 
									
										
										
										
											2007-08-15 14:28:22 +00:00
										 |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00:00
										 |  |  |    Return the sine of *x*.
 | 
					
						
							| 
									
										
										
										
											2007-08-15 14:28:22 +00:00
										 |  |  | 
 | 
					
						
							|  |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00:00
										 |  |  | .. function:: tan(x)
 | 
					
						
							| 
									
										
										
										
											2007-08-15 14:28:22 +00:00
										 |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00:00
										 |  |  |    Return the tangent of *x*.
 | 
					
						
							| 
									
										
										
										
											2008-04-19 00:31:39 +00:00
										 |  |  | 
 | 
					
						
							|  |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00:00
										 |  |  | Hyperbolic functions
 | 
					
						
							|  |  |  | --------------------
 | 
					
						
							| 
									
										
										
										
											2008-04-19 00:31:39 +00:00
										 |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00:00
										 |  |  | .. function:: acosh(x)
 | 
					
						
							| 
									
										
										
										
											2008-04-19 00:31:39 +00:00
										 |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00: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-19 00:31:39 +00:00
										 |  |  | 
 | 
					
						
							|  |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00:00
										 |  |  | .. function:: asinh(x)
 | 
					
						
							| 
									
										
										
										
											2007-08-15 14:28:22 +00:00
										 |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00: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.
 | 
					
						
							| 
									
										
										
										
											2007-08-15 14:28:22 +00:00
										 |  |  | 
 | 
					
						
							|  |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00:00
										 |  |  | .. function:: atanh(x)
 | 
					
						
							| 
									
										
										
										
											2007-08-15 14:28:22 +00:00
										 |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00: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 14:28:22 +00:00
										 |  |  | 
 | 
					
						
							|  |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00:00
										 |  |  | .. function:: cosh(x)
 | 
					
						
							| 
									
										
										
										
											2007-08-15 14:28:22 +00:00
										 |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00:00
										 |  |  |    Return the hyperbolic cosine of *x*.
 | 
					
						
							| 
									
										
										
										
											2007-08-15 14:28:22 +00:00
										 |  |  | 
 | 
					
						
							|  |  |  | 
 | 
					
						
							|  |  |  | .. function:: sinh(x)
 | 
					
						
							|  |  |  | 
 | 
					
						
							|  |  |  |    Return the hyperbolic sine of *x*.
 | 
					
						
							|  |  |  | 
 | 
					
						
							|  |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00:00
										 |  |  | .. function:: tanh(x)
 | 
					
						
							| 
									
										
										
										
											2007-08-15 14:28:22 +00:00
										 |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00:00
										 |  |  |    Return the hyperbolic tangent of *x*.
 | 
					
						
							| 
									
										
										
										
											2007-08-15 14:28:22 +00:00
										 |  |  | 
 | 
					
						
							|  |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00:00
										 |  |  | Classification functions
 | 
					
						
							|  |  |  | ------------------------
 | 
					
						
							| 
									
										
										
										
											2007-08-15 14:28:22 +00:00
										 |  |  | 
 | 
					
						
							| 
									
										
										
										
											2010-07-11 17:38:24 +00:00
										 |  |  | .. function:: isfinite(x)
 | 
					
						
							|  |  |  | 
 | 
					
						
							| 
									
										
										
										
											2010-07-11 19:47:37 +00:00
										 |  |  |    Return ``True`` if both the real and imaginary parts of *x* are finite, and
 | 
					
						
							|  |  |  |    ``False`` otherwise.
 | 
					
						
							|  |  |  | 
 | 
					
						
							|  |  |  |    .. versionadded:: 3.2
 | 
					
						
							| 
									
										
										
										
											2010-07-11 17:38:24 +00:00
										 |  |  | 
 | 
					
						
							|  |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00:00
										 |  |  | .. function:: isinf(x)
 | 
					
						
							|  |  |  | 
 | 
					
						
							| 
									
										
										
										
											2010-07-11 19:47:37 +00:00
										 |  |  |    Return ``True`` if either the real or the imaginary part of *x* is an
 | 
					
						
							|  |  |  |    infinity, and ``False`` otherwise.
 | 
					
						
							| 
									
										
										
										
											2007-08-15 14:28:22 +00:00
										 |  |  | 
 | 
					
						
							|  |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00:00
										 |  |  | .. function:: isnan(x)
 | 
					
						
							|  |  |  | 
 | 
					
						
							| 
									
										
										
										
											2010-07-11 19:47:37 +00:00
										 |  |  |    Return ``True`` if either the real or the imaginary part of *x* is a NaN,
 | 
					
						
							|  |  |  |    and ``False`` otherwise.
 | 
					
						
							| 
									
										
										
										
											2007-08-15 14:28:22 +00:00
										 |  |  | 
 | 
					
						
							|  |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00:00
										 |  |  | Constants
 | 
					
						
							|  |  |  | ---------
 | 
					
						
							| 
									
										
										
										
											2007-08-15 14:28:22 +00:00
										 |  |  | 
 | 
					
						
							|  |  |  | 
 | 
					
						
							|  |  |  | .. data:: pi
 | 
					
						
							|  |  |  | 
 | 
					
						
							| 
									
										
										
										
											2009-07-28 16:31:03 +00:00
										 |  |  |    The mathematical constant *π*, as a float.
 | 
					
						
							| 
									
										
										
										
											2007-08-15 14:28:22 +00: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-19 00:31:39 +00:00
										 |  |  | 
 |