resqpy 4.13.1__tar.gz → 4.13.3__tar.gz

{resqpy-4.13.1 → resqpy-4.13.3}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: resqpy
-Version: 4.13.1
+Version: 4.13.3
 Summary: Python API for working with RESQML models
 Home-page: https://github.com/bp/resqpy
 License: MIT

{resqpy-4.13.1 → resqpy-4.13.3}/pyproject.toml

@@ -9,7 +9,7 @@ build-backend = "poetry.masonry.api"
 
 [tool.poetry]
 name = "resqpy"
-version = "4.13.1" # Set at build time
+version = "4.13.3" # Set at build time
 description = "Python API for working with RESQML models"
 authors = ["BP"]
 license = "MIT"

{resqpy-4.13.1 → resqpy-4.13.3}/resqpy/__init__.py

@@ -28,6 +28,6 @@
 
 import logging
 
-__version__ = "4.13.1"  # Set at build time
+__version__ = "4.13.3"  # Set at build time
 log = logging.getLogger(__name__)
 log.info(f"Imported resqpy version {__version__}")

{resqpy-4.13.1 → resqpy-4.13.3}/resqpy/lines/_polyline.py

@@ -427,6 +427,10 @@ class Polyline(rql_c._BasePolyline):
         """
 
         assert 0.0 <= fraction <= 1.0
+        if maths.isclose(fraction, 0.0):
+            return self.coordinates[0]
+        if maths.isclose(fraction, 1.0):
+            return (self.coordinates[0] if self.isclosed else self.coordinates[-1])
         target = fraction * self.full_length(in_xy = in_xy)
         seg_index = 0
         while True:

{resqpy-4.13.1 → resqpy-4.13.3}/resqpy/surface/_surface.py

@@ -539,6 +539,16 @@ class Surface(rqsb.BaseSurface):
         assert saucer_parameter is None or 0.0 <= saucer_parameter < 1.0
         crs = rqc.Crs(self.model, uuid = point_set.crs_uuid)
         p = point_set.full_array_ref()
+        assert p.ndim >= 2
+        assert p.shape[-1] == 3
+        p = p.reshape((-1, 3))
+        nan_mask = np.isnan(p)
+        if np.any(nan_mask):
+            row_mask = np.logical_not(np.any(nan_mask, axis = -1))
+            log.info(
+                f'removing {len(p) - np.count_nonzero(row_mask)} NaN points from point set {point_set.title} prior to surface triangulation'
+            )
+            p = p[row_mask, :]
         if crs.xy_units == crs.z_units or not reorient:
             unit_adjusted_p = p
         else: