Describe the bug
Healpix map spatial queries result in bad writes when one or more ranks did a query near its boundary. The end result is that the set of valid pixels is not the same in the written dataset, even after accounting for re-partitoning.
To Reproduce
There is a known failing case in test/parallel/test_healpixmap_mpi.py::test_healpix_write
Expected behavior
Healpix maps should warn users if they are doing a query that does not fully overlap with the map, but writing should still proceed correctly.
Desktop (please complete the following information):
- OS: Any
- Python Version: Any
- Using MPI? yes
Describe the bug
Healpix map spatial queries result in bad writes when one or more ranks did a query near its boundary. The end result is that the set of valid pixels is not the same in the written dataset, even after accounting for re-partitoning.
To Reproduce
There is a known failing case in
test/parallel/test_healpixmap_mpi.py::test_healpix_writeExpected behavior
Healpix maps should warn users if they are doing a query that does not fully overlap with the map, but writing should still proceed correctly.
Desktop (please complete the following information):