spacetelescope / spherical_geometry

A Python package for handling spherical polygons that represent arbitrary regions of the sky
http://spherical-geometry.readthedocs.io/
62 stars 31 forks source link

TST: 27 test failures with numpy 2.1.dev #261

Closed pllim closed 5 months ago

pllim commented 5 months ago

Example log: https://github.com/spacetelescope/spherical_geometry/actions/runs/8595276994/job/23549707417

Example traceback (one of many):

>       A = polygon.SphericalPolygon(a)

spherical_geometry/tests/test_union.py:329: 
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ 
spherical_geometry/polygon.py:772: in __init__
    polygons.extend(g.disjoint_polygons())
spherical_geometry/graph.py:458: in disjoint_polygons
    changed = self._find_all_intersections() or changed
spherical_geometry/graph.py:643: in _find_all_intersections
    changed = self._find_point_to_arc_intersections() or changed
spherical_geometry/graph.py:542: in _find_point_to_arc_intersections
    intersects = gca.intersects_point(
spherical_geometry/great_circle_arc.py:286: in intersects_point
    total_length = length(A, B)
spherical_geometry/great_circle_arc.py:227: in length
    dot = inner1d(A, B)
spherical_geometry/great_circle_arc.py:31: in _inner1d_np
    return np.multiply(x, y).sum(axis=1)
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ 

a = array([0.09374197, 0.00091684, 0.90534119]), axis = 1, dtype = None
out = None, keepdims = False, initial = <no value>, where = True

    def _sum(a, axis=None, dtype=None, out=None, keepdims=False,
             initial=_NoValue, where=True):
>       return umr_sum(a, axis, dtype, out, keepdims, initial, where)
E       numpy.exceptions.AxisError: axis 1 is out of bounds for array of dimension 1