deploy: a14f1f3

Deltares · Sep 2, 2024 · 5e97a1e · 5e97a1e
1 parent b463b19
commit 5e97a1e
Show file tree

Hide file tree

Showing 101 changed files with 1,773 additions and 127 deletions.
diff --git a/.doctrees/api/api/pandamesh.GmshMesher.doctree b/.doctrees/api/api/pandamesh.GmshMesher.doctree
diff --git a/.doctrees/api/api/pandamesh.TriangleMesher.doctree b/.doctrees/api/api/pandamesh.TriangleMesher.doctree
diff --git a/.doctrees/api/api/pandamesh.find_edge_intersections.doctree b/.doctrees/api/api/pandamesh.find_edge_intersections.doctree
diff --git a/.doctrees/api/index.doctree b/.doctrees/api/index.doctree
diff --git a/.doctrees/environment.pickle b/.doctrees/environment.pickle
diff --git a/.doctrees/examples/gmsh-basic.doctree b/.doctrees/examples/gmsh-basic.doctree
diff --git a/.doctrees/examples/gmsh-fields.doctree b/.doctrees/examples/gmsh-fields.doctree
diff --git a/.doctrees/examples/preprocessing.doctree b/.doctrees/examples/preprocessing.doctree
diff --git a/.doctrees/examples/sg_execution_times.doctree b/.doctrees/examples/sg_execution_times.doctree
diff --git a/.doctrees/examples/triangle-basic.doctree b/.doctrees/examples/triangle-basic.doctree
diff --git a/.doctrees/examples/triangle-geospatial.doctree b/.doctrees/examples/triangle-geospatial.doctree
diff --git a/.doctrees/sg_execution_times.doctree b/.doctrees/sg_execution_times.doctree
diff --git a/_downloads/5bc19a906f7aee7609a9fe736a8fcec5/preprocessing.py b/_downloads/5bc19a906f7aee7609a9fe736a8fcec5/preprocessing.py
@@ -221,6 +221,15 @@
 gdf["cellsize"] = [2.0, 1.0, 1.0]
 gdf.plot(edgecolor="k")
 
+# %%
+# We can identify these problematic intersections using
+# :func:`pandamesh.find_edge_intersections`:
+
+intersections = pm.find_edge_intersections(gdf.geometry)
+fig, ax = plt.subplots()
+gdf.plot(ax=ax, facecolor="none")
+intersections.plot(ax=ax)
+
 # %%
 # A first step is to remove line segments that do not fall in any polygon:
 
@@ -234,9 +243,10 @@
 # %%
 # However, this doesn't create suitable input for meshing. The ``GmshMesher``
 # appears to hang on this input, and Triangle generates a grid with very small
-# triangles:
+# triangles. Pandamesh errors on these intersections by default, but way may
+# proceed:
 
-vertices, faces = pm.TriangleMesher(resolved).generate()
+vertices, faces = pm.TriangleMesher(resolved, intersecting_edges="warn").generate()
 pm.plot(vertices, faces)
 
 # %%

diff --git a/_downloads/79a425db58e6180cef408093f1ab28a8/preprocessing.zip b/_downloads/79a425db58e6180cef408093f1ab28a8/preprocessing.zip
diff --git a/_downloads/7a7ce2204255f277838397c00eecb9b7/gmsh-basic.zip b/_downloads/7a7ce2204255f277838397c00eecb9b7/gmsh-basic.zip
diff --git a/_downloads/8d89d0cac7b99393a56c43d6bed625ed/triangle-basic.zip b/_downloads/8d89d0cac7b99393a56c43d6bed625ed/triangle-basic.zip
diff --git a/_downloads/a977adb4611357edd665dfb42e7a7790/preprocessing.ipynb b/_downloads/a977adb4611357edd665dfb42e7a7790/preprocessing.ipynb
@@ -198,6 +198,24 @@
         "donut = sg.Polygon(\n    [\n        [0.0, 0.0],\n        [10.0, 0.0],\n        [10.0, 10.0],\n        [0.0, 10.0],\n    ],\n    holes=[\n        [\n            [2.0, 5.0],\n            [5.0, 8.0],\n            [8.0, 5.0],\n            [5.0, 2.0],\n        ]\n    ],\n)\nline0 = shapely.LineString(\n    [\n        [-2.0, 0.0],\n        [12.0, 10.0],\n    ]\n)\nline1 = shapely.LineString(\n    [\n        [5.5, 9.0],\n        [9.0, 5.5],\n    ]\n)\n\ngdf = gpd.GeoDataFrame(geometry=[donut, line0, line1])\ngdf[\"cellsize\"] = [2.0, 1.0, 1.0]\ngdf.plot(edgecolor=\"k\")"
       ]
     },
+    {
+      "cell_type": "markdown",
+      "metadata": {},
+      "source": [
+        "We can identify these problematic intersections using\n:func:`pandamesh.find_edge_intersections`:\n\n"
+      ]
+    },
+    {
+      "cell_type": "code",
+      "execution_count": null,
+      "metadata": {
+        "collapsed": false
+      },
+      "outputs": [],
+      "source": [
+        "intersections = pm.find_edge_intersections(gdf.geometry)\nfig, ax = plt.subplots()\ngdf.plot(ax=ax, facecolor=\"none\")\nintersections.plot(ax=ax)"
+      ]
+    },
     {
       "cell_type": "markdown",
       "metadata": {},
@@ -220,7 +238,7 @@
       "cell_type": "markdown",
       "metadata": {},
       "source": [
-        "However, this doesn't create suitable input for meshing. The ``GmshMesher``\nappears to hang on this input, and Triangle generates a grid with very small\ntriangles:\n\n"
+        "However, this doesn't create suitable input for meshing. The ``GmshMesher``\nappears to hang on this input, and Triangle generates a grid with very small\ntriangles. Pandamesh errors on these intersections by default, but way may\nproceed:\n\n"
       ]
     },
     {
@@ -231,7 +249,7 @@
       },
       "outputs": [],
       "source": [
-        "vertices, faces = pm.TriangleMesher(resolved).generate()\npm.plot(vertices, faces)"
+        "vertices, faces = pm.TriangleMesher(resolved, intersecting_edges=\"warn\").generate()\npm.plot(vertices, faces)"
       ]
     },
     {

diff --git a/_downloads/ae350e4fcfd5c3cf82fb43e401fc1d34/triangle-geospatial.zip b/_downloads/ae350e4fcfd5c3cf82fb43e401fc1d34/triangle-geospatial.zip
diff --git a/_downloads/d85ad210e21a91643987461e0d983111/gmsh-fields.zip b/_downloads/d85ad210e21a91643987461e0d983111/gmsh-fields.zip
diff --git a/_images/sphx_glr_preprocessing_010.png b/_images/sphx_glr_preprocessing_010.png
diff --git a/_images/sphx_glr_preprocessing_011.png b/_images/sphx_glr_preprocessing_011.png
diff --git a/_images/sphx_glr_preprocessing_012.png b/_images/sphx_glr_preprocessing_012.png
diff --git a/_images/sphx_glr_preprocessing_013.png b/_images/sphx_glr_preprocessing_013.png
diff --git a/_images/sphx_glr_preprocessing_014.png b/_images/sphx_glr_preprocessing_014.png
diff --git a/_images/sphx_glr_preprocessing_015.png b/_images/sphx_glr_preprocessing_015.png
diff --git a/_images/sphx_glr_preprocessing_016.png b/_images/sphx_glr_preprocessing_016.png
diff --git a/_images/sphx_glr_preprocessing_017.png b/_images/sphx_glr_preprocessing_017.png
diff --git a/_images/sphx_glr_preprocessing_018.png b/_images/sphx_glr_preprocessing_018.png
diff --git a/_images/sphx_glr_preprocessing_019.png b/_images/sphx_glr_preprocessing_019.png
diff --git a/_modules/index.html b/_modules/index.html
@@ -354,7 +354,8 @@
                 <article class="bd-article">
 
   <h1>All modules for which code is available</h1>
-<ul><li><a href="pandamesh/gmsh_enums.html">pandamesh.gmsh_enums</a></li>
+<ul><li><a href="pandamesh/common.html">pandamesh.common</a></li>
+<li><a href="pandamesh/gmsh_enums.html">pandamesh.gmsh_enums</a></li>
 <li><a href="pandamesh/gmsh_mesher.html">pandamesh.gmsh_mesher</a></li>
 <li><a href="pandamesh/preprocessor.html">pandamesh.preprocessor</a></li>
 <li><a href="pandamesh/triangle_enums.html">pandamesh.triangle_enums</a></li>