cmath.rst 7.05 KB
Newer Older
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
: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.

16
.. note::
17

18 19 20 21 22 23 24
   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.


25 26
Conversions to and from polar coordinates
-----------------------------------------
27

28 29 30 31
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::
32

33
   z == z.real + z.imag*1j
34

35 36 37 38 39 40
*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
angle from the positive x-axis to the line segment that joins the
origin to *z*.
41

42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59
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
60 61 62 63


.. note::

64 65 66
   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.
67 68


69
.. function:: polar(x)
70

71 72 73 74
   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))``.
75 76


77
.. function:: rect(r, phi)
78

79 80
   Return the complex number *x* with polar coordinates *r* and *phi*.
   Equivalent to ``r * (math.cos(phi) + math.sin(phi)*1j)``.
81 82


83 84
Power and logarithmic functions
-------------------------------
85

86
.. function:: exp(x)
87

88
   Return the exponential value ``e**x``.
89 90


91
.. function:: log(x[, base])
92

93 94 95
   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.
96

97 98
   .. versionchanged:: 2.4
      *base* argument added.
99

100

101
.. function:: log10(x)
102

103 104
   Return the base-10 logarithm of *x*. This has the same branch cut as
   :func:`log`.
105 106


107
.. function:: sqrt(x)
108

109
   Return the square root of *x*. This has the same branch cut as :func:`log`.
110 111


112 113
Trigonometric functions
-----------------------
114

115
.. function:: acos(x)
116

117 118 119
   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.
120 121


122 123 124
.. function:: asin(x)

   Return the arc sine of *x*. This has the same branch cuts as :func:`acos`.
125

126 127 128 129

.. function:: atan(x)

   Return the arc tangent of *x*. There are two branch cuts: One extends from
130
   ``1j`` along the imaginary axis to ``∞j``, continuous from the right. The
131
   other extends from ``-1j`` along the imaginary axis to ``-∞j``, continuous
132 133
   from the left.

134 135 136 137 138 139

.. function:: cos(x)

   Return the cosine of *x*.


140
.. function:: sin(x)
141

142
   Return the sine of *x*.
143 144


145
.. function:: tan(x)
146

147
   Return the tangent of *x*.
148 149


150 151
Hyperbolic functions
--------------------
152

153
.. function:: acosh(x)
154

155 156
   Return the hyperbolic arc cosine of *x*. There is one branch cut, extending left
   from 1 along the real axis to -∞, continuous from above.
157 158


159
.. function:: asinh(x)
160

161 162 163 164
   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.
165 166


167
.. function:: atanh(x)
168

169 170 171 172
   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.
173 174


175
.. function:: cosh(x)
176

177
   Return the hyperbolic cosine of *x*.
178 179 180 181 182 183 184


.. function:: sinh(x)

   Return the hyperbolic sine of *x*.


185
.. function:: tanh(x)
186

187
   Return the hyperbolic tangent of *x*.
188 189


190 191
Classification functions
------------------------
192

193 194 195 196
.. function:: isinf(x)

   Return *True* if the real or the imaginary part of x is positive
   or negative infinity.
197 198


199 200 201
.. function:: isnan(x)

   Return *True* if the real or imaginary part of x is not a number (NaN).
202 203


204 205
Constants
---------
206 207 208 209


.. data:: pi

210
   The mathematical constant *π*, as a float.
211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240


.. 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.

241