[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.
See #367

Author: philippeVerney

.github/workflows/github-actions.yml

@@ -71,11 +71,10 @@ jobs:
       fail-fast: false  # 'false' means Don't stop matrix workflows even if some matrix failed.
       matrix:
         include: [
-          { xcc_name: 'gcc 9',      xcc_pkg: gcc-9,       cc: gcc-9,      cxx: g++-9 },
           { xcc_name: 'gcc 10',     xcc_pkg: gcc-10,      cc: gcc-10,     cxx: g++-10 },
-          { xcc_name: 'clang 11',   xcc_pkg: clang-11,    cc: clang-11,   cxx: clang++-11 },
-          { xcc_name: 'clang 12',   xcc_pkg: clang-12,    cc: clang-12,   cxx: clang++-12 },
+          { xcc_name: 'gcc 12',     xcc_pkg: gcc-12,      cc: gcc-12,     cxx: g++-12 },
           { xcc_name: 'clang 13',   xcc_pkg: clang-13,    cc: clang-13,   cxx: clang++-13 },
+          { xcc_name: 'clang 15',   xcc_pkg: clang-15,    cc: clang-15,   cxx: clang++-15 },
         ]
     env:
       XCC: ${{ matrix.xcc_name }}

example/example.cpp

@@ -1229,7 +1229,7 @@ void serializeGrid(COMMON_NS::DataObjectRepository * repo, EML2_NS::AbstractHdfP
 	// The cumulative count of points per face i.e. first face contains 3 points, second face contains 6-3=3 points, third face contains 9-6=3 points etc...
 	uint64_t nodeIndicesCumulativeCountPerFace[8] = { 3, // shared face
 		6, 9, 12, // faces for tetra
-		16, 20, 24, 27 //  faces for wedge
+		16, 20, 24, 27 // faces for wedge
 	};
 	// The face indices of each cell. 
 	uint64_t faceIndicesPerCell[9] = { 0, 1, 2, 3, // tetra

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,