1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 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 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265
|
************
Known Issues
************
.. contents::
:local:
:depth: 2
While most bugs and issues are managed using the `astropy issue
tracker <https://github.com/astropy/astropy/issues>`_, this document
lists issues that are too difficult to fix, may require some
intervention from the user to work around, or are caused by bugs in other
projects or packages.
Issues listed on this page are grouped into two categories: The first is known
issues and shortcomings in actual algorithms and interfaces that currently do
not have fixes or workarounds, and that users should be aware of when writing
code that uses ``astropy``. Some of those issues are still platform-specific,
while others are very general. The second category is of common issues that come
up when configuring, building, or installing ``astropy``. This also includes
cases where the test suite can report false negatives depending on the context/
platform on which it was run.
Known Deficiencies
==================
.. _quantity_issues:
Quantities Lose Their Units with Some Operations
------------------------------------------------
Quantities are subclassed from ``numpy``'s `~numpy.ndarray` and while we have
ensured that ``numpy`` functions will work well with them, they do not always
work in functions from ``scipy`` or other packages that use ``numpy``
internally, but ignore the subclass. Furthermore, at a few places in ``numpy``
itself we cannot control the behaviour. For instance, care must be taken when
setting array slices using Quantities::
>>> import astropy.units as u
>>> import numpy as np
>>> a = np.ones(4)
>>> a[2:3] = 2*u.kg
>>> a # doctest: +FLOAT_CMP
array([1., 1., 2., 1.])
::
>>> a = np.ones(4)
>>> a[2:3] = 1*u.cm/u.m
>>> a # doctest: +FLOAT_CMP
array([1., 1., 1., 1.])
Either set single array entries or use lists of Quantities::
>>> a = np.ones(4)
>>> a[2] = 1*u.cm/u.m
>>> a # doctest: +FLOAT_CMP
array([1. , 1. , 0.01, 1. ])
::
>>> a = np.ones(4)
>>> a[2:3] = [1*u.cm/u.m]
>>> a # doctest: +FLOAT_CMP
array([1. , 1. , 0.01, 1. ])
Both will throw an exception if units do not cancel, e.g.::
>>> a = np.ones(4)
>>> a[2] = 1*u.cm
Traceback (most recent call last):
...
TypeError: only dimensionless scalar quantities can be converted to Python scalars
See: https://github.com/astropy/astropy/issues/7582
Numpy array creation functions cannot be used to initialize Quantity
--------------------------------------------------------------------
Trying the following example will ignore the unit:
>>> np.full(10, 1 * u.m)
array([1., 1., 1., 1., 1., 1., 1., 1., 1., 1.])
A workaround for this at the moment would be to do::
>>> np.full(10, 1) << u.m
<Quantity [1., 1., 1., 1., 1., 1., 1., 1., 1., 1.] m>
As well as with `~numpy.full` one cannot do `~numpy.zeros`, `~numpy.ones`, and `~numpy.empty`.
The `~numpy.arange` function does not work either::
>>> np.arange(0 * u.m, 10 * u.m, 1 * u.m)
Traceback (most recent call last):
...
TypeError: only dimensionless scalar quantities can be converted to Python scalars
Workarounds include moving the units outside of the call to
`~numpy.arange`::
>>> np.arange(0, 10, 1) * u.m
<Quantity [0., 1., 2., 3., 4., 5., 6., 7., 8., 9.] m>
Also, `~numpy.linspace` does work:
>>> np.linspace(0 * u.m, 9 * u.m, 10)
<Quantity [0., 1., 2., 3., 4., 5., 6., 7., 8., 9.] m>
Quantities Lose Their Units When Broadcasted
--------------------------------------------
When broadcasting Quantities, it is necessary to pass ``subok=True`` to
`~numpy.broadcast_to`, or else a bare `~numpy.ndarray` will be returned::
>>> q = u.Quantity(np.arange(10.), u.m)
>>> b = np.broadcast_to(q, (2, len(q)))
>>> b # doctest: +FLOAT_CMP
array([[0., 1., 2., 3., 4., 5., 6., 7., 8., 9.],
[0., 1., 2., 3., 4., 5., 6., 7., 8., 9.]])
>>> b2 = np.broadcast_to(q, (2, len(q)), subok=True)
>>> b2 # doctest: +FLOAT_CMP
<Quantity [[0., 1., 2., 3., 4., 5., 6., 7., 8., 9.],
[0., 1., 2., 3., 4., 5., 6., 7., 8., 9.]] m>
This is analogous to the case of passing a Quantity to `~numpy.array`::
>>> a = np.array(q)
>>> a # doctest: +FLOAT_CMP
array([0., 1., 2., 3., 4., 5., 6., 7., 8., 9.])
>>> a2 = np.array(q, subok=True)
>>> a2 # doctest: +FLOAT_CMP
<Quantity [0., 1., 2., 3., 4., 5., 6., 7., 8., 9.] m>
See: https://github.com/astropy/astropy/issues/7832
mmap Support for ``astropy.io.fits`` on GNU Hurd
------------------------------------------------
On Hurd and possibly other platforms, ``flush()`` on memory-mapped files are not
implemented, so writing changes to a mmap'd FITS file may not be reliable and is
thus disabled. Attempting to open a FITS file in writeable mode with mmap will
result in a warning (and mmap will be disabled on the file automatically).
See: https://github.com/astropy/astropy/issues/968
Color Printing on Windows
-------------------------
Colored printing of log messages and other colored text does work in Windows,
but only when running in the IPython console. Colors are not currently
supported in the basic Python command-line interpreter on Windows.
``numpy.int64`` does not decompose input ``Quantity`` objects
-------------------------------------------------------------
Python's ``int()`` goes through ``__index__``
while ``numpy.int64`` or ``numpy.int_`` do not go through ``__index__``. This
means that an upstream fix in NumPy is required in order for
``astropy.units`` to control decomposing the input in these functions::
>>> np.int64((15 * u.km) / (15 * u.imperial.foot))
1
>>> np.int_((15 * u.km) / (15 * u.imperial.foot))
1
>>> int((15 * u.km) / (15 * u.imperial.foot))
3280
To convert a dimensionless `~astropy.units.Quantity` to an integer, it is
therefore recommended to use ``int(...)``.
Inconsistent behavior when converting complex numbers to floats
---------------------------------------------------------------
Attempting to use `float` or NumPy's ``numpy.float`` on a standard
complex number (e.g., ``5 + 6j``) results in a `TypeError`. In
contrast, using `float` or ``numpy.float`` on a complex number from
NumPy (e.g., ``numpy.complex128``) drops the imaginary component and
issues a ``numpy.ComplexWarning``. This inconsistency persists between
`~astropy.units.Quantity` instances based on standard and NumPy
complex numbers. To get the real part of a complex number, it is
recommended to use ``numpy.real``.
.. _structured_unit_deserialization_segfault:
Structured units deserialization segfaults in big-endian
--------------------------------------------------------
Structured units deserialization with ``pickle`` may cause segmentation
fault in big-endian machine with ``numpy<1.21.1``.
Build/Installation/Test Issues
==============================
Anaconda Users Should Upgrade with ``conda``, Not ``pip``
---------------------------------------------------------
Upgrading ``astropy`` in the Anaconda Python distribution using ``pip`` can result
in a corrupted install with a mix of files from the old version and the new
version. Anaconda users should update with ``conda update astropy``. There
may be a brief delay between the release of ``astropy`` on PyPI and its release
via the ``conda`` package manager; users can check the availability of new
versions with ``conda search astropy``.
Locale Errors in MacOS X and Linux
----------------------------------
On MacOS X, you may see the following error when running ``pip``::
...
ValueError: unknown locale: UTF-8
This is due to the ``LC_CTYPE`` environment variable being incorrectly set to
``UTF-8`` by default, which is not a valid locale setting.
On MacOS X or Linux (or other platforms) you may also encounter the following
error::
...
stderr = stderr.decode(stdio_encoding)
TypeError: decode() argument 1 must be str, not None
This also indicates that your locale is not set correctly.
To fix either of these issues, set this environment variable, as well as the
``LANG`` and ``LC_ALL`` environment variables to e.g. ``en_US.UTF-8`` using, in
the case of ``bash``::
export LANG="en_US.UTF-8"
export LC_ALL="en_US.UTF-8"
export LC_CTYPE="en_US.UTF-8"
To avoid any issues in future, you should add this line to your e.g.
``~/.bash_profile`` or ``.bashrc`` file.
To test these changes, open a new terminal and type ``locale``, and you should
see something like::
$ locale
LANG="en_US.UTF-8"
LC_COLLATE="en_US.UTF-8"
LC_CTYPE="en_US.UTF-8"
LC_MESSAGES="en_US.UTF-8"
LC_MONETARY="en_US.UTF-8"
LC_NUMERIC="en_US.UTF-8"
LC_TIME="en_US.UTF-8"
LC_ALL="en_US.UTF-8"
If so, you can go ahead and try running ``pip`` again (in the new
terminal).
Failing Logging Tests When Running the Tests in IPython
-------------------------------------------------------
When running the Astropy tests using ``astropy.test()`` in an IPython
interpreter, some of the tests in the ``astropy/tests/test_logger.py`` *might*
fail depending on the version of IPython or other factors.
This is due to mutually incompatible behaviors in IPython and pytest, and is
not due to a problem with the test itself or the feature being tested.
See: https://github.com/astropy/astropy/issues/717
|