[Python] Add unstructured grid example and fix array casting

Introduced an unstructured grid example in example.py and updated several method calls to use the .cast() method for array arguments to ensure correct type handling. Also made minor updates in fesapi.ipynb for array casting and kernel metadata.

Author: philippeVerney

python/example/example.py

@@ -20,7 +20,7 @@ def serialize_grid(repo: fesapi.DataObjectRepository):
         resqml_points.setitem(i, value)
     single_cell_ijk_grid.setGeometryAsCoordinateLineNodes(
         fesapi.resqml20__PillarShape_vertical, fesapi.resqml20__KDirection_down, False,
-        resqml_points)
+        resqml_points.cast())
 
     # ONE CONTINUOUS PROPERTY ON A PARTIAL 3 cells I=1 J=1 K=3 GRID
     partial_ijk_grid = repo.createPartialIjkGridRepresentation("bc92a216-aa17-4a1f-9253-8b3ab094bf84", "partial grid")
@@ -33,7 +33,48 @@ def serialize_grid(repo: fesapi.DataObjectRepository):
     stats = fesapi.FloatArrayStatistics()
     stats.setMaximum(3.3)
     stats.setMinimum(1.1)
-    continuous_prop.pushBackFloatArray3dOfValuesPlusStatistics(resqml_values, 1, 1, 3, stats)
+    continuous_prop.pushBackFloatArray3dOfValuesPlusStatistics(resqml_values.cast(), 1, 1, 3, stats)
+
+    # unstructured grid example
+    unstructured_grid = repo.createUnstructuredGridRepresentation(
+        "9283cd33-5e52-4110-b7b1-616abde2b303",
+        "One tetrahedron + prism grid",
+        2
+    )
+    python_points = [0, 0, 300, 375, 0, 300, 0, 150, 300, # points for shared face between tetra and wedge
+		0, 0, 500, # point for tetra
+		0, 0, 0, 375, 0, 0, 0, 150, 0]
+    resqml_points = fesapi.DoubleArray(21)
+    for i, value in enumerate(python_points):
+        resqml_points.setitem(i, value)
+    python_node_indices_per_face = [0, 1, 2, # shared face
+		1, 2, 3, 0, 1, 3, 0, 2, 3, # faces for tetra
+		0, 2, 6, 4, 2, 1, 5, 6, 0, 1, 5, 4, 4, 5, 6] # faces for wedge
+    resqml_node_indices_per_face = fesapi.UInt64Array(27)
+    for i, value in enumerate(python_node_indices_per_face):
+        resqml_node_indices_per_face.setitem(i, value)
+    python_node_indices_cumulative_count_per_face = [3, # shared face
+		6, 9, 12, # faces for tetra
+		16, 20, 24, 27] # faces for wedge
+    resqml_node_indices_cumulative_count_per_face = fesapi.UInt64Array(8)
+    for i, value in enumerate(python_node_indices_cumulative_count_per_face):
+        resqml_node_indices_cumulative_count_per_face.setitem(i, value)
+    python_face_indices_per_cell = [0, 1, 2, 3, # tetra
+		0, 4, 5, 6, 7] # wedge
+    resqml_face_indices_per_cell = fesapi.UInt64Array(9)
+    for i, value in enumerate(python_face_indices_per_cell):
+        resqml_face_indices_per_cell.setitem(i, value)
+    python_face_indices_cumulative_count_per_cell = [ 4, 9 ];
+    resqml_face_indices_cumulative_count_per_cell = fesapi.UInt64Array(2)
+    for i, value in enumerate(python_face_indices_cumulative_count_per_cell):
+        resqml_face_indices_cumulative_count_per_cell.setitem(i, value)
+    python_face_right_handness = [ 1, 0, 0, 1, 0, 0, 0, 1, 1 ];
+    resqml_face_right_handness = fesapi.UInt8Array(9)
+    for i, value in enumerate(python_face_right_handness):
+        resqml_face_right_handness.setitem(i, value)
+    unstructured_grid.setGeometry(resqml_face_right_handness.cast(), resqml_points.cast(), 7, None, resqml_face_indices_per_cell.cast(), resqml_face_indices_cumulative_count_per_cell.cast(), 8,
+                                  resqml_node_indices_per_face.cast(), resqml_node_indices_cumulative_count_per_face.cast(), fesapi.resqml20__CellShape_prism);
+
 
 def serialize(file_name: str):
     """
@@ -122,15 +163,15 @@ def show_ijk_grid(ijk_grid: fesapi.Resqml2_AbstractIjkGridRepresentation):
         for patch_index in range(patch_count):
             nb_xyz_points_in_patch = ijk_grid.getXyzPointCountOfPatch(patch_index)
             xyz_points_in_patch = fesapi.DoubleArray(nb_xyz_points_in_patch * 3)
-            ijk_grid.getXyzPointsOfPatch(patch_index, xyz_points_in_patch)
+            ijk_grid.getXyzPointsOfPatch(patch_index, xyz_points_in_patch.cast())
             for vertex_index in range(nb_xyz_points_in_patch):
                 x = xyz_points_in_patch.getitem(vertex_index * 3)
                 y = xyz_points_in_patch.getitem(vertex_index * 3 + 1)
                 z = xyz_points_in_patch.getitem(vertex_index * 3 + 2)
                 print("Vertex ", vertex_index, " : ", x, " ", y, " ", z)
 
         xyz_points = fesapi.DoubleArray(nb_xyz_points * 3)
-        ijk_grid.getXyzPointsOfAllPatches(xyz_points)
+        ijk_grid.getXyzPointsOfAllPatches(xyz_points.cast())
 
         ijk_grid.loadSplitInformation()
         for cell_corner in range(8):

python/example/fesapi.ipynb

@@ -116,7 +116,7 @@
     "resqml_points = fesapi.DoubleArray(24)\n",
     "for i in range(6):\n",
     "    resqml_points.setitem(i, i*100)\n",
-    "horizon_grid_2d_representation.setGeometryAsArray2dOfExplicitZ(resqml_points, 2, 3, hdf_proxy,\n",
+    "horizon_grid_2d_representation.setGeometryAsArray2dOfExplicitZ(resqml_points.cast(), 2, 3, hdf_proxy,\n",
     "                                                                0.0, 0.0, 0.0,\n",
     "                                                               1.0, 0.0, 0.0, 25.0,\n",
     "                                                               0.0, 1.0, 0.0, 50.0)"
@@ -132,15 +132,12 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 6,
    "id": "b5d289e6-c699-4d88-a870-3dae2baf14ee",
    "metadata": {},
-   "outputs": [],
    "source": [
     "# Create an EPC Document in the current folder nad named TestingFesapiWithPython.epc\n",
     "epc_file = fesapi.EpcDocument(file_name_without_extension + \".epc\")\n",
-    "epc_file.serializeFrom(repo)\n",
-    "epc_file.close()"
+    "epc_file.serializeFrom(repo)"
    ]
   },
   {
@@ -164,7 +161,7 @@
  ],
  "metadata": {
   "kernelspec": {
-   "display_name": "Python 3 (ipykernel)",
+   "display_name": ".venv (3.11.9)",
    "language": "python",
    "name": "python3"
   },
@@ -178,7 +175,7 @@
    "name": "python",
    "nbconvert_exporter": "python",
    "pygments_lexer": "ipython3",
-   "version": "3.11.5"
+   "version": "3.11.9"
   }
  },
  "nbformat": 4,