Open
Description
When trying to pass nest=True to healpy.query_strip, the Python interpreter crashes:
>>> import healpy
>>> healpy.query_strip(2, 1, 2, nest=True)
Error encountered at /Users/runner/work/healpy/healpy/healpixsubmodule/src/cxx/Healpix_cxx/healpix_base.cc, line 1058
(function void T_Healpix_Base<I>::query_strip_internal(double, double, bool, rangeset<T>&) const [with I = long long int])
query_strip not yet implemented for NESTED
terminate called after throwing an instance of 'PlanckError'It is not hard to work around this issue, e.g.:
>>> import healpy
>>> healpy.ring2nest(2, healpy.query_strip(2, 1, 2))
array([19, 0, 23, 4, 27, 8, 31, 12, 17, 22, 21, 26, 25, 30, 29, 18, 16,
35, 20, 39, 24, 43, 28, 47])but it would be great if healpy would do this workaround internally to make users less confused :-)
Metadata
Metadata
Assignees
Labels
No labels