updated version, added xarray export test

.github/workflows/tests.yml

@@ -21,7 +21,7 @@ jobs:
     - name: Install dependencies
       run: |
         python -m pip install --upgrade pip
-        pip install -e .
+        pip install -e .[raster]
 
     - name: Build C extensions
       run: |

CHANGELOG.md

@@ -13,4 +13,7 @@
 - fixed history tracking bug
 - fixed a bug in the rng algorithm that impacted reproducibility in some cases
 - odd-sized grids are now supported
-- fixed a bug preventing early exit in some situations during optimization when the convergence criterion is met
+- fixed a bug preventing early exit in some situations during optimization when the convergence criterion is met
+
+# 1.3.20251010
+- fixed critical bug preventing export to xarray datasets

docs/conf.py

@@ -14,7 +14,7 @@
 project = 'PyOCN'
 copyright = '2025, Alexander S. Fox'
 author = 'Alexander S. Fox'
-release = '1.2.20251010'
+release = '1.3.20251010'
 
 # -- General configuration ---------------------------------------------------
 # https://www.sphinx-doc.org/en/master/usage/configuration.html#general-configuration

pyproject.toml

@@ -4,7 +4,7 @@ build-backend = "setuptools.build_meta"
 
 [project]
 name = "PyOCN"
-version = "1.2.20251010"
+version = "1.3.20251010"
 description = "Optimal Channel Networks (OCN) in Python with a C core"
 readme = "readme.md"
 requires-python = ">=3.10"

tests/test_basic.py

@@ -11,6 +11,23 @@
 class TestBasicOCN(unittest.TestCase):
     """Basic tests for OCN creation and energy computation."""
 
+    def test_ocn_to_xarray(self):
+        """Test that OCN can be converted to xarray and back."""
+        ocn = po.OCN.from_net_type(
+            net_type="V",
+            dims=(32, 32),
+            random_state=1234,
+        )
+        ocn.to_xarray()
+
+        ocn = po.OCN.from_net_type(
+            net_type="V",
+            dims=(32, 32),
+            random_state=1234,
+            wrap=True,
+        )
+        ocn.to_xarray()
+
     def test_ocn_energy(self):
         """Test that OCN creation produces expected energy value."""
         ocn = po.OCN.from_net_type(