polytope-python 2.1.5__tar.gz → 2.1.7__tar.gz

{polytope_python-2.1.5 → polytope_python-2.1.7}/.github/ci-config.yml

@@ -1,6 +1,6 @@
 dependencies: |
   ecmwf/ecbuild
-  MathisRosenhauer/libaec@refs/tags/v1.1.3
+  MathisRosenhauer/libaec@refs/tags/v1.1.4
   ecmwf/eccodes
   ecmwf/eckit@develop
   ecmwf/metkit

{polytope_python-2.1.5 → polytope_python-2.1.7}/.github/workflows/label-public-pr.yml

@@ -7,4 +7,4 @@ on:
 
 jobs:
   label:
-    uses: ecmwf-actions/reusable-workflows/.github/workflows/label-pr.yml@v2
+    uses: ecmwf/reusable-workflows/.github/workflows/label-pr.yml@v2

{polytope_python-2.1.5 → polytope_python-2.1.7}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: polytope-python
-Version: 2.1.5
+Version: 2.1.7
 Summary: Polytope datacube feature extraction library
 Author-email: "European Centre for Medium-Range Weather Forecasts (ECMWF)" <software.support@ecmwf.int>
 Maintainer-email: James Hawkes <James.Hawkes@ecmwf.int>, Mathilde Leuridan <Mathilde.Leuridan@ecmwf.int>

{polytope_python-2.1.5 → polytope_python-2.1.7}/examples/3D_shipping_route.py

@@ -39,13 +39,8 @@ class Test:
                 lines.append(point_list)
 
         speed_km_hr = 100
-        initial_distance = math.sqrt(
-            lines[0][0] ** 2 + lines[0][1] ** 2
-        )  # This is technically only the degrees...
-        time_hours = [
-            (math.sqrt(pt[0] ** 2 + pt[1] ** 2) - initial_distance) / speed_km_hr
-            for pt in lines
-        ]
+        initial_distance = math.sqrt(lines[0][0] ** 2 + lines[0][1] ** 2)  # This is technically only the degrees...
+        time_hours = [(math.sqrt(pt[0] ** 2 + pt[1] ** 2) - initial_distance) / speed_km_hr for pt in lines]
         step_list_seconds = np.array(time_hours)
         step_list_seconds = step_list_seconds * 3600 * 1000 * 1000 * 1000
         step_list_seconds = [int(a) for a in step_list_seconds]
@@ -55,23 +50,17 @@ class Test:
 
         new_points = []
         for point, step in zip(lines[:7], step_list):
-            new_points.append(
-                [point[1], point[0], step + np.timedelta64(3600000000000, "ns")]
-            )
+            new_points.append([point[1], point[0], step + np.timedelta64(3600000000000, "ns")])
 
         # Pad the shipping route with an initial shape
 
         padded_point_upper = [0.24, 0.24, np.timedelta64(3590, "s")]
         padded_point_lower = [-0.24, -0.24, np.timedelta64(-3590, "s")]
-        initial_shape = Ellipsoid(
-            ["latitude", "longitude", "step"], padded_point_lower, padded_point_upper
-        )
+        initial_shape = Ellipsoid(["latitude", "longitude", "step"], padded_point_lower, padded_point_upper)
 
         # Then somehow make this list of points into just a sequence of points
 
-        ship_route_polytope = Path(
-            ["latitude", "longitude", "step"], initial_shape, *new_points
-        )
+        ship_route_polytope = Path(["latitude", "longitude", "step"], initial_shape, *new_points)
         request = Request(
             ship_route_polytope,
             Select("number", [0]),
@@ -101,9 +90,7 @@ class Test:
         fig, ax = plt.subplots(figsize=(12, 6))
         worldmap.plot(color="darkgrey", ax=ax)
         ax.scatter(longs, lats, s=8, c=parameter_values, cmap="viridis")
-        norm = mpl.colors.Normalize(
-            vmin=min(parameter_values), vmax=max(parameter_values)
-        )
+        norm = mpl.colors.Normalize(vmin=min(parameter_values), vmax=max(parameter_values))
 
         sm = plt.cm.ScalarMappable(cmap="viridis", norm=norm)
         sm.set_array([])

{polytope_python-2.1.5 → polytope_python-2.1.7}/examples/3D_shipping_route_wave_model.py

@@ -25,9 +25,7 @@ class TestReducedLatLonGrid:
                     }
                 }
             },
-            "date": {
-                "transformation": {"merge": {"with": "time", "linkers": ["T", "00"]}}
-            },
+            "date": {"transformation": {"merge": {"with": "time", "linkers": ["T", "00"]}}},
             "step": {"transformation": {"type_change": "int"}},
             "number": {"transformation": {"type_change": "int"}},
             "longitude": {"transformation": {"cyclic": [0, 360]}},
@@ -35,9 +33,7 @@ class TestReducedLatLonGrid:
         self.config = {"class": "od", "stream": "wave"}
         self.fdbdatacube = FDBDatacube(self.config, axis_options=self.options)
         self.slicer = HullSlicer()
-        self.API = Polytope(
-            datacube=self.fdbdatacube, engine=self.slicer, axis_options=self.options
-        )
+        self.API = Polytope(datacube=self.fdbdatacube, engine=self.slicer, axis_options=self.options)
 
     def find_nearest_latlon(self, grib_file, target_lat, target_lon):
         messages = grib_file
@@ -77,15 +73,11 @@ class TestReducedLatLonGrid:
 
         padded_point_upper = [0.24, 0.24, 1]
         padded_point_lower = [-0.24, -0.24, 1]
-        initial_shape = Ellipsoid(
-            ["latitude", "longitude", "step"], padded_point_lower, padded_point_upper
-        )
+        initial_shape = Ellipsoid(["latitude", "longitude", "step"], padded_point_lower, padded_point_upper)
 
         # Then somehow make this list of points into just a sequence of points
 
-        ship_route_polytope = Path(
-            ["latitude", "longitude", "step"], initial_shape, *new_points
-        )
+        ship_route_polytope = Path(["latitude", "longitude", "step"], initial_shape, *new_points)
 
         request = Request(
             ship_route_polytope,

{polytope_python-2.1.5 → polytope_python-2.1.7}/examples/4D_flight_path.py

@@ -21,9 +21,7 @@ class Test:
             array = array.sortby(dim)
         self.array = array
         self.slicer = HullSlicer()
-        self.API = Polytope(
-            datacube=array, engine=self.slicer, axis_options=axis_options
-        )
+        self.API = Polytope(datacube=array, engine=self.slicer, axis_options=axis_options)
 
     def test_slice_shipping_route(self):
         colorscale = [
@@ -56,9 +54,7 @@ class Test:
         texture = np.asarray(Image.open("./examples/data/earth_image.jpg")).T
         radius = 1
         x, y, z = sphere(radius, texture)
-        surf = go.Surface(
-            x=x, y=y, z=z, surfacecolor=texture, colorscale=colorscale, hoverinfo="none"
-        )
+        surf = go.Surface(x=x, y=y, z=z, surfacecolor=texture, colorscale=colorscale, hoverinfo="none")
         layout = go.Layout(scene=dict(aspectratio=dict(x=1, y=1, z=1)), spikedistance=0)
         fig = go.Figure(data=[surf], layout=layout)
         fig.update_layout(scene=dict(xaxis_showspikes=False, yaxis_showspikes=False))
@@ -96,9 +92,7 @@ class Test:
             *route_point_CDG_LHR,
         )
 
-        request = Request(
-            flight_route_polytope, Select("time", ["2022-12-02T12:00:00"])
-        )
+        request = Request(flight_route_polytope, Select("time", ["2022-12-02T12:00:00"]))
         result = self.API.retrieve(request)
 
         lats = []
@@ -128,9 +122,7 @@ class Test:
         for i in range(len(lats)):
             x1 = radius * np.cos(longs[i]) * np.cos(lats[i])
             y1 = radius * np.sin(longs[i]) * np.cos(lats[i])
-            z1 = (
-                radius * np.sin(lats[i]) + alts[i] / 300
-            )  # the 300 here should really be the earth of the radius
+            z1 = radius * np.sin(lats[i]) + alts[i] / 300  # the 300 here should really be the earth of the radius
             x.append(x1)
             y.append(y1)
             z.append(z1)
@@ -157,15 +149,9 @@ class Test:
                 landcolor="rgb(230, 145, 56)",
                 lakecolor="rgb(0, 255, 255)",
                 oceancolor="rgb(0, 255, 255)",
-                projection=dict(
-                    type="orthographic", rotation=dict(lon=-100, lat=40, roll=0)
-                ),
-                lonaxis=dict(
-                    showgrid=True, gridcolor="rgb(102, 102, 102)", gridwidth=0.5
-                ),
-                lataxis=dict(
-                    showgrid=True, gridcolor="rgb(102, 102, 102)", gridwidth=0.5
-                ),
+                projection=dict(type="orthographic", rotation=dict(lon=-100, lat=40, roll=0)),
+                lonaxis=dict(showgrid=True, gridcolor="rgb(102, 102, 102)", gridwidth=0.5),
+                lataxis=dict(showgrid=True, gridcolor="rgb(102, 102, 102)", gridwidth=0.5),
             ),
         )
 

{polytope_python-2.1.5 → polytope_python-2.1.7}/examples/country_slicing.py

@@ -18,9 +18,7 @@ class Test:
         axis_options = {"longitude": {"cyclic": [0, 360.0]}}
         self.xarraydatacube = XArrayDatacube(array)
         self.slicer = HullSlicer()
-        self.API = Polytope(
-            datacube=array, engine=self.slicer, axis_options=axis_options
-        )
+        self.API = Polytope(datacube=array, engine=self.slicer, axis_options=axis_options)
 
     def test_slice_country(self):
         # Read a shapefile for a given country and extract the geometry polygons

{polytope_python-2.1.5 → polytope_python-2.1.7}/examples/cyclic_route_around_earth.py

@@ -17,15 +17,11 @@ class Test:
         axis_options = {"longitude": {"cyclic": [0, 360.0]}}
         self.xarraydatacube = XArrayDatacube(array)
         self.slicer = HullSlicer()
-        self.API = Polytope(
-            datacube=array, engine=self.slicer, axis_options=axis_options
-        )
+        self.API = Polytope(datacube=array, engine=self.slicer, axis_options=axis_options)
 
     def test_slice_country(self):
         bounding_box = Box(["latitude", "longitude"], [-0.1, -0.1], [0.1, 0.1])
-        request_obj = PathSegment(
-            ["latitude", "longitude"], bounding_box, [-88, -67], [68, 170]
-        )
+        request_obj = PathSegment(["latitude", "longitude"], bounding_box, [-88, -67], [68, 170])
         request = Request(
             request_obj,
             Select("number", [0]),

{polytope_python-2.1.5 → polytope_python-2.1.7}/examples/healpix_grid_box_example.py

@@ -12,9 +12,7 @@ from polytope_feature.shapes import Box, Select
 class TestOctahedralGrid:
     def setup_method(self, method):
         ds = data.from_source("file", "./tests/data/healpix.grib")
-        self.latlon_array = (
-            ds.to_xarray().isel(step=0).isel(time=0).isel(isobaricInhPa=0).z
-        )
+        self.latlon_array = ds.to_xarray().isel(step=0).isel(time=0).isel(isobaricInhPa=0).z
         self.xarraydatacube = XArrayDatacube(self.latlon_array)
         self.options = {
             "values": {
@@ -26,9 +24,7 @@ class TestOctahedralGrid:
             }
         }
         self.slicer = HullSlicer()
-        self.API = Polytope(
-            datacube=self.latlon_array, engine=self.slicer, axis_options=self.options
-        )
+        self.API = Polytope(datacube=self.latlon_array, engine=self.slicer, axis_options=self.options)
 
     def find_nearest_latlon(self, grib_file, target_lat, target_lon):
         # Open the GRIB file
@@ -75,9 +71,7 @@ class TestOctahedralGrid:
             lon = cubepath["longitude"]
             lats.append(lat)
             lons.append(lon)
-            nearest_points = self.find_nearest_latlon(
-                "./tests/data/healpix.grib", lat, lon
-            )
+            nearest_points = self.find_nearest_latlon("./tests/data/healpix.grib", lat, lon)
             eccodes_lat = nearest_points[0][0]["lat"]
             eccodes_lon = nearest_points[0][0]["lon"]
             eccodes_lats.append(eccodes_lat)

{polytope_python-2.1.5 → polytope_python-2.1.7}/examples/plotting_country_data.py

@@ -85,8 +85,7 @@ improved_polytope_data = [improved_polytope_data_original[i] for i in mask]
 country_list = [country_list_original[i] for i in mask]
 
 normalised_improvement = [
-    polytope / rectangle
-    for polytope, rectangle in zip(improved_polytope_data, naive_rectangle_data)
+    polytope / rectangle for polytope, rectangle in zip(improved_polytope_data, naive_rectangle_data)
 ]
 
 normalised_rectangle_data = list(np.ones(len(naive_rectangle_data)))

{polytope_python-2.1.5 → polytope_python-2.1.7}/examples/slicing_all_ecmwf_countries.py

@@ -18,9 +18,7 @@ class Test:
         axis_options = {"longitude": {"cyclic": [0, 360.0]}}
         self.xarraydatacube = XArrayDatacube(array)
         self.slicer = HullSlicer()
-        self.API = Polytope(
-            datacube=array, engine=self.slicer, axis_options=axis_options
-        )
+        self.API = Polytope(datacube=array, engine=self.slicer, axis_options=axis_options)
 
     def test_slice_country(self):
         # Read a shapefile for a given country and extract the geometry polygons
@@ -153,9 +151,7 @@ class Test:
             for geom in multi_polygon.geoms:
                 plt.plot(*geom.exterior.xy, color="black", linewidth=0.7)
         countries_temps = np.array(countries_temps)
-        plt.scatter(
-            countries_longs, countries_lats, s=8, c=countries_temps, cmap="YlOrRd"
-        )
+        plt.scatter(countries_longs, countries_lats, s=8, c=countries_temps, cmap="YlOrRd")
         plt.colorbar(label="Temperature")
 
         plt.show()

{polytope_python-2.1.5 → polytope_python-2.1.7}/examples/wind_farms.py

@@ -21,16 +21,12 @@ class Test:
         axis_options = {"longitude": {"cyclic": [0, 360.0]}}
         self.xarraydatacube = XArrayDatacube(array)
         self.slicer = HullSlicer()
-        self.API = Polytope(
-            datacube=array, engine=self.slicer, axis_options=axis_options
-        )
+        self.API = Polytope(datacube=array, engine=self.slicer, axis_options=axis_options)
 
     def test_slice_wind_farms(self):
         gdal.SetConfigOption("SHAPE_RESTORE_SHX", "YES")
 
-        shapefile = gpd.read_file(
-            "./examples/data/EMODnet_HA_WindFarms_pg_20220324.shp"
-        )
+        shapefile = gpd.read_file("./examples/data/EMODnet_HA_WindFarms_pg_20220324.shp")
         polygons = []
         for i in range(306):
             country = shapefile.iloc[i]
@@ -89,9 +85,7 @@ class Test:
         fig, ax = plt.subplots(figsize=(12, 6))
         worldmap.plot(color="darkgrey", ax=ax)
         ax.scatter(longs, lats, s=8, c=parameter_values, cmap="viridis")
-        norm = mpl.colors.Normalize(
-            vmin=min(parameter_values), vmax=max(parameter_values)
-        )
+        norm = mpl.colors.Normalize(vmin=min(parameter_values), vmax=max(parameter_values))
 
         sm = plt.cm.ScalarMappable(cmap="viridis", norm=norm)
         sm.set_array([])

{polytope_python-2.1.5 → polytope_python-2.1.7}/performance/fdb_performance.py

@@ -21,17 +21,13 @@ class TestSlicingFDBDatacube:
                     }
                 }
             },
-            "date": {
-                "transformation": {"merge": {"with": "time", "linkers": [" ", "00"]}}
-            },
+            "date": {"transformation": {"merge": {"with": "time", "linkers": [" ", "00"]}}},
             "step": {"transformation": {"type_change": "int"}},
         }
         self.config = {"class": "od", "expver": "0001", "levtype": "sfc", "step": 0}
         self.fdbdatacube = FDBDatacube(self.config, axis_options=self.options)
         self.slicer = HullSlicer()
-        self.API = Polytope(
-            datacube=self.fdbdatacube, engine=self.slicer, axis_options=self.options
-        )
+        self.API = Polytope(datacube=self.fdbdatacube, engine=self.slicer, axis_options=self.options)
 
     # Testing different shapes
     # @pytest.mark.skip(reason="can't install fdb branch on CI")

{polytope_python-2.1.5 → polytope_python-2.1.7}/performance/fdb_performance_3D.py

@@ -21,9 +21,7 @@ class TestSlicingFDBDatacube:
                     }
                 }
             },
-            "date": {
-                "transformation": {"merge": {"with": "time", "linkers": [" ", "00"]}}
-            },
+            "date": {"transformation": {"merge": {"with": "time", "linkers": [" ", "00"]}}},
             "step": {"transformation": {"type_change": "int"}},
             "levelist": {"transformation": {"type_change": "int"}},
             "longitude": {"transformation": {"cyclic": [0, 360]}},
@@ -31,9 +29,7 @@ class TestSlicingFDBDatacube:
         self.config = {"class": "od", "expver": "0001", "levtype": "sfc"}
         self.fdbdatacube = FDBDatacube(self.config, axis_options=self.options)
         self.slicer = HullSlicer()
-        self.API = Polytope(
-            datacube=self.fdbdatacube, engine=self.slicer, axis_options=self.options
-        )
+        self.API = Polytope(datacube=self.fdbdatacube, engine=self.slicer, axis_options=self.options)
 
     # Testing different shapes
     # @pytest.mark.skip(reason="can't install fdb branch on CI")

{polytope_python-2.1.5 → polytope_python-2.1.7}/performance/fdb_slice_many_numbers_timeseries.py

@@ -21,9 +21,7 @@ options = {
         },
         {
             "axis_name": "date",
-            "transformations": [
-                {"name": "merge", "other_axis": "time", "linkers": ["T", "00"]}
-            ],
+            "transformations": [{"name": "merge", "other_axis": "time", "linkers": ["T", "00"]}],
         },
         {
             "axis_name": "values",

{polytope_python-2.1.5 → polytope_python-2.1.7}/performance/scalability_test.py

@@ -11,9 +11,7 @@ from polytope_feature.shapes import Box, Disk, Ellipsoid, Select
 
 class Test:
     def setup_method(self):
-        array = xr.open_dataset(
-            "../examples/data/temp_model_levels.grib", engine="cfgrib"
-        ).t
+        array = xr.open_dataset("../examples/data/temp_model_levels.grib", engine="cfgrib").t
         options = {"longitude": {"Cyclic": [0, 360.0]}}
         self.xarraydatacube = XArrayDatacube(array)
         for dim in array.dims:
@@ -26,9 +24,7 @@ class Test:
         time_start = time.time()
         print(time_start)
         box = Box(["latitude", "longitude"], [0, 0], [50, 360])
-        request = Request(
-            box, Select("step", [np.timedelta64(0, "ns")]), Select("hybrid", [1])
-        )
+        request = Request(box, Select("step", [np.timedelta64(0, "ns")]), Select("hybrid", [1]))
         result = self.API.retrieve(request)
         # result.pprint()
         print(len(result.leaves))
@@ -38,9 +34,7 @@ class Test:
         time_start = time.time()
         print(time_start)
         box = Box(["latitude", "longitude"], [0, 0], [100, 360])
-        request = Request(
-            box, Select("step", [np.timedelta64(0, "ns")]), Select("hybrid", [1])
-        )
+        request = Request(box, Select("step", [np.timedelta64(0, "ns")]), Select("hybrid", [1]))
         result = self.API.retrieve(request)
         print(len(result.leaves))
         print(time.time() - time_start)
@@ -49,9 +43,7 @@ class Test:
         time_start = time.time()
         print(time_start)
         box = Box(["latitude", "longitude"], [-50, 0], [100, 360])
-        request = Request(
-            box, Select("step", [np.timedelta64(0, "ns")]), Select("hybrid", [1])
-        )
+        request = Request(box, Select("step", [np.timedelta64(0, "ns")]), Select("hybrid", [1]))
         result = self.API.retrieve(request)
         print(len(result.leaves))
         print(time.time() - time_start)
@@ -60,9 +52,7 @@ class Test:
         time_start = time.time()
         print(time_start)
         box = Box(["latitude", "longitude"], [-100, 0], [100, 360])
-        request = Request(
-            box, Select("step", [np.timedelta64(0, "ns")]), Select("hybrid", [1])
-        )
+        request = Request(box, Select("step", [np.timedelta64(0, "ns")]), Select("hybrid", [1]))
         result = self.API.retrieve(request)
         print(len(result.leaves))
         print(time.time() - time_start)
@@ -71,9 +61,7 @@ class Test:
         time_start = time.time()
         print(time_start)
         box = Box(["latitude", "longitude"], [0, 0], [50, 180])
-        request = Request(
-            box, Select("step", [np.timedelta64(0, "ns")]), Select("hybrid", [1])
-        )
+        request = Request(box, Select("step", [np.timedelta64(0, "ns")]), Select("hybrid", [1]))
         result = self.API.retrieve(request)
         print(len(result.leaves))
         print(time.time() - time_start)
@@ -82,9 +70,7 @@ class Test:
         time_start = time.time()
         print(time_start)
         box = Box(["latitude", "longitude"], [-100, -180], [100, 360])
-        request = Request(
-            box, Select("step", [np.timedelta64(0, "ns")]), Select("hybrid", [1])
-        )
+        request = Request(box, Select("step", [np.timedelta64(0, "ns")]), Select("hybrid", [1]))
         result = self.API.retrieve(request)
         print(len(result.leaves))
         print(time.time() - time_start)
@@ -93,9 +79,7 @@ class Test:
         time_start = time.time()
         print(time_start)
         box = Box(["latitude", "longitude"], [-100, -360], [100, 360])
-        request = Request(
-            box, Select("step", [np.timedelta64(0, "ns")]), Select("hybrid", [1])
-        )
+        request = Request(box, Select("step", [np.timedelta64(0, "ns")]), Select("hybrid", [1]))
         result = self.API.retrieve(request)
         print(len(result.leaves))
         print(time.time() - time_start)
@@ -156,9 +140,7 @@ class Test:
         time_start = time.time()
         print(time_start)
         box = Disk(["latitude", "longitude"], [0, 0], [25, 180])
-        request = Request(
-            box, Select("step", [np.timedelta64(0, "ns")]), Select("hybrid", [1])
-        )
+        request = Request(box, Select("step", [np.timedelta64(0, "ns")]), Select("hybrid", [1]))
         result = self.API.retrieve(request)
         # result.pprint()
         print(len(result.leaves))
@@ -168,9 +150,7 @@ class Test:
         time_start = time.time()
         print(time_start)
         box = Disk(["latitude", "longitude"], [0, 0], [50, 180])
-        request = Request(
-            box, Select("step", [np.timedelta64(0, "ns")]), Select("hybrid", [1])
-        )
+        request = Request(box, Select("step", [np.timedelta64(0, "ns")]), Select("hybrid", [1]))
         result = self.API.retrieve(request)
         # result.pprint()
         print(len(result.leaves))
@@ -180,9 +160,7 @@ class Test:
         time_start = time.time()
         print(time_start)
         box = Disk(["latitude", "longitude"], [0, 0], [25, 90])
-        request = Request(
-            box, Select("step", [np.timedelta64(0, "ns")]), Select("hybrid", [1])
-        )
+        request = Request(box, Select("step", [np.timedelta64(0, "ns")]), Select("hybrid", [1]))
         result = self.API.retrieve(request)
         # result.pprint()
         print(len(result.leaves))

{polytope_python-2.1.5 → polytope_python-2.1.7}/performance/scalability_test_2.py

@@ -11,9 +11,7 @@ from polytope_feature.shapes import Box, Select, Union
 
 class Test:
     def setup_method(self):
-        array = xr.open_dataset(
-            "./examples/data/temp_model_levels.grib", engine="cfgrib"
-        ).t
+        array = xr.open_dataset("./examples/data/temp_model_levels.grib", engine="cfgrib").t
         options = {"longitude": {"Cyclic": [0, 360.0]}}
         self.xarraydatacube = XArrayDatacube(array)
         for dim in array.dims:
@@ -88,9 +86,7 @@ class Test:
                 union = Union(["latitude", "longitude"], union, box)
         time_start = time.time()
         print(time_start)
-        request = Request(
-            union, Select("step", [np.timedelta64(0, "ns")]), Select("hybrid", [1])
-        )
+        request = Request(union, Select("step", [np.timedelta64(0, "ns")]), Select("hybrid", [1]))
         result = self.API.retrieve(request)
         # result.pprint()
         print(len(result.leaves))

{polytope_python-2.1.5 → polytope_python-2.1.7}/performance_unstructured/octahedral_vs_unstructured_slicing.py

@@ -96,9 +96,7 @@ def set_up_slicing_unstructured(file_name):
             },
             {
                 "axis_name": "date",
-                "transformations": [
-                    {"name": "merge", "other_axis": "time", "linkers": ["T", "00"]}
-                ],
+                "transformations": [{"name": "merge", "other_axis": "time", "linkers": ["T", "00"]}],
             },
             {
                 "axis_name": "values",
@@ -181,9 +179,7 @@ def set_up_slicing_structured():
             },
             {
                 "axis_name": "date",
-                "transformations": [
-                    {"name": "merge", "other_axis": "time", "linkers": ["T", "00"]}
-                ],
+                "transformations": [{"name": "merge", "other_axis": "time", "linkers": ["T", "00"]}],
             },
             {
                 "axis_name": "values",

{polytope_python-2.1.5 → polytope_python-2.1.7}/performance_unstructured/plot_structured_vs_unstructured_slicing.py

@@ -30,13 +30,8 @@ tot_retrieve_time_unstructured = [
 ]
 
 
-pure_poly_time_structured = [
-    tot_retrieve_time_structured[i] - gj_retrieve_structured_time[i] for i in range(4)
-]
-pure_poly_time_unstructured = [
-    tot_retrieve_time_unstructured[i] - gj_retrieve_unstructured_time[i]
-    for i in range(4)
-]
+pure_poly_time_structured = [tot_retrieve_time_structured[i] - gj_retrieve_structured_time[i] for i in range(4)]
+pure_poly_time_unstructured = [tot_retrieve_time_unstructured[i] - gj_retrieve_unstructured_time[i] for i in range(4)]
 
 fig, ax = plt.subplots(1, 1)
 

{polytope_python-2.1.5 → polytope_python-2.1.7}/polytope_feature/_version.py

@@ -1,2 +1,2 @@
 # Do not change! Do not track in version control!
-__version__ = "2.1.5"
+__version__ = "2.1.7"

{polytope_python-2.1.5 → polytope_python-2.1.7}/polytope_feature/datacube/backends/datacube.py

@@ -49,9 +49,7 @@ class Datacube(ABC):
         """returns true if the input axes can be resolved against the datacube axes"""
         return validate_axes(list(self.axes.keys()), axes)
 
-    def _create_axes(
-        self, name, values, transformation_type_key, transformation_options
-    ):
+    def _create_axes(self, name, values, transformation_type_key, transformation_options):
         # first check what the final axes are for this axis name given transformations
         transformation_options = transformation_type_key
         final_axis_names, transformation = DatacubeAxisTransformation.get_final_axes(
@@ -98,9 +96,7 @@ class Datacube(ABC):
                 # etc axis properties add the specific transformation handled here to the relevant axes
                 # Modify the axis to update with the tag
 
-                if (
-                    transformation not in self._axes[axis_name].transformations
-                ):  # Avoids duplicates being stored
+                if transformation not in self._axes[axis_name].transformations:  # Avoids duplicates being stored
                     self._axes[axis_name].transformations.append(transformation)
             else:
                 # Means we have an unsliceable axis since we couln't transform values to desired type
@@ -145,9 +141,7 @@ class Datacube(ABC):
 
         idx_between = axis.find_indices_between(indexes, lower, upper, self, method)
 
-        logging.debug(
-            f"For axis {axis.name} between {lower} and {upper}, found indices {idx_between}"
-        )
+        logging.debug(f"For axis {axis.name} between {lower} and {upper}, found indices {idx_between}")
 
         return idx_between