loadpath 0.1.1__py3-none-any.whl

loadpath/data/published_experiments.json

@@ -0,0 +1,71 @@
+{
+ "about": "Published compression tests of printed lattices with fully specified geometry, used to check loadpath's lattice models against measurements (bench/published_experiments.py, docs/VALIDATION.md). Values are copied from the papers' tables; mean +- standard deviation where the paper gives one. Each specimen is a cube of lattice with cut-open side faces, compressed between flat platens along the build direction.",
+ "sources": {
+  "shevchenko2023": {
+   "citation": "Shevchenko V, Balabanov S, Sychov M, Karimova L. Prediction of Cellular Structure Mechanical Properties with the Geometry of Triply Periodic Minimal Surfaces (TPMS). ACS Omega 8(30):26895, 2023.",
+   "doi": "10.1021/acsomega.3c01631",
+   "url": "https://pmc.ncbi.nlm.nih.gov/articles/PMC10398688/",
+   "tables": "Table 1 (solid compressive E and strength of each printed material, same lab), Tables 3-5 (lattices)",
+   "notes": "30 mm cubes, 3x3x3 sheet cells (10 mm), 0.8 mm wall; SLS PA2200 (EOS Formiga P110, 0.1 mm layers), FDM PLA (Anet A6, 0.2 mm layers), LCD urethane-acrylate resin (Anycubic Photon S, 0.05 mm). ISO 604, 5 mm/min, along the build axis; modulus from the stress-strain curve (crosshead strain, no extensometer). Solid modulus and strength measured in compression by the same lab, so the material card is the paper's own. Strength is the first peak (sigma_B). Lattice phase at the faces not stated."
+  },
+  "guerrasilva2021": {
+   "citation": "Guerra Silva R, Salinas Estay C, Morales Pavez G, Zahr Vinuela J, Torres MJ. Influence of Geometric and Manufacturing Parameters on the Compressive Behavior of 3D Printed Polymer Lattice Structures. Materials 14(6):1462, 2021.",
+   "doi": "10.3390/ma14061462",
+   "url": "https://pmc.ncbi.nlm.nih.gov/articles/PMC8002545/",
+   "tables": "Table 2 (feedstock properties, from datasheets), Tables 3 and 5 (configurations, measured modulus and plateau stress)",
+   "notes": "35 mm cubes, FDM (UP mini 2, 0.4 mm nozzle, raft, no supports) in ABS and PC; truncated-octahedron (Kelvin) and cubic diamond strut cells, 5 or 7 mm cells, 1.0 or 1.5 mm round struts, 0.15 or 0.25 mm layers; two specimens each, ASTM D1621, 2 mm/min, crosshead strain. Solid modulus and strength are filament datasheet values, not printed coupons. Strength reported is the plateau stress, not first yield."
+  },
+  "abouali2022": {
+   "citation": "Abou-Ali AM, Lee DW, Abu Al-Rub RK. On the Effect of Lattice Topology on Mechanical Properties of SLS Additively Manufactured Sheet-, Ligament-, and Strut-Based Polymeric Metamaterials. Polymers 14(21):4583, 2022.",
+   "doi": "10.3390/polym14214583",
+   "url": "https://pmc.ncbi.nlm.nih.gov/articles/PMC9658147/",
+   "tables": "Tables 2-4 (designed and measured relative density), Table 5 (power-law fits E = C RD^n, yield = C RD^n over measured RD)",
+   "notes": "40 mm cubes, 5x5x5 cells of 8 mm, SLS PA 1102 black (polyamide 11, EOS Formiga P110, 80 um layers), 0.001/s. Only power-law fits over the measured relative density are published, so each 'measured' value here is the paper's fit evaluated at its measured RD (two specimens per point); the solid modulus is the EOS datasheet's (1560 MPa x/y, ISO 527), strength 48 MPa. Sheet lattices weighed 1.4-2.9x their designed density (sintered powder stuck to the large surface): loadpath models the designed geometry."
+  }
+ },
+ "materials": {
+  "shevchenko-pa12-sls": {"E": 1260.0, "strength": 57.3, "nu": 0.38, "density": 1.01, "basis": "measured in compression by the same lab (Table 1)"},
+  "shevchenko-pla-fdm": {"E": 1832.0, "strength": 78.7, "nu": 0.35, "density": 1.21, "basis": "measured in compression by the same lab (Table 1)"},
+  "shevchenko-resin-lcd": {"E": 1287.0, "strength": 58.3, "nu": 0.4, "density": 1.18, "basis": "measured in compression by the same lab (Table 1)"},
+  "guerrasilva-abs": {"E": 2180.0, "strength": 30.0, "nu": 0.35, "density": 1.195, "basis": "filament datasheet (paper's Table 2)"},
+  "guerrasilva-pc": {"E": 1958.0, "strength": 64.0, "nu": 0.37, "density": 1.2, "basis": "filament datasheet (paper's Table 2)"},
+  "abouali-pa11-sls": {"E": 1560.0, "strength": 48.0, "nu": 0.4, "density": 0.99, "basis": "EOS PA 1102 black datasheet (x/y, ISO 527)"}
+ },
+ "specimens": [
+  {"id": "shev-sls-primitive", "source": "shevchenko2023", "material": "shevchenko-pa12-sls", "process": "SLS PA12", "lattice": {"type": "schwarz-p", "cell": 10.0, "rho": 0.18, "wall": 0.8}, "size": [30, 30, 30], "measured": {"E": [62.8, 3.2], "strength": [2.67, 0.12], "strengthKind": "first peak"}},
+  {"id": "shev-sls-diamond", "source": "shevchenko2023", "material": "shevchenko-pa12-sls", "process": "SLS PA12", "lattice": {"type": "diamond", "cell": 10.0, "rho": 0.33, "wall": 0.8}, "size": [30, 30, 30], "measured": {"E": [105.4, 4.1], "strength": [5.90, 0.19], "strengthKind": "first peak"}},
+  {"id": "shev-sls-gyroid", "source": "shevchenko2023", "material": "shevchenko-pa12-sls", "process": "SLS PA12", "lattice": {"type": "gyroid", "cell": 10.0, "rho": 0.24, "wall": 0.8}, "size": [30, 30, 30], "measured": {"E": [68.9, 3.7], "strength": [3.10, 0.15], "strengthKind": "first peak"}},
+  {"id": "shev-fdm-primitive", "source": "shevchenko2023", "material": "shevchenko-pla-fdm", "process": "FDM PLA", "lattice": {"type": "schwarz-p", "cell": 10.0, "rho": 0.18, "wall": 0.8}, "size": [30, 30, 30], "measured": {"E": [38.3, 2.1], "strength": [1.72, 0.15], "strengthKind": "first peak"}},
+  {"id": "shev-fdm-diamond", "source": "shevchenko2023", "material": "shevchenko-pla-fdm", "process": "FDM PLA", "lattice": {"type": "diamond", "cell": 10.0, "rho": 0.33, "wall": 0.8}, "size": [30, 30, 30], "measured": {"E": [134.1, 4.2], "strength": [8.13, 0.21], "strengthKind": "first peak"}},
+  {"id": "shev-fdm-gyroid", "source": "shevchenko2023", "material": "shevchenko-pla-fdm", "process": "FDM PLA", "lattice": {"type": "gyroid", "cell": 10.0, "rho": 0.24, "wall": 0.8}, "size": [30, 30, 30], "measured": {"E": [79.2, 2.2], "strength": [4.73, 0.16], "strengthKind": "first peak"}},
+  {"id": "shev-lcd-primitive", "source": "shevchenko2023", "material": "shevchenko-resin-lcd", "process": "LCD resin", "lattice": {"type": "schwarz-p", "cell": 10.0, "rho": 0.18, "wall": 0.8}, "size": [30, 30, 30], "measured": {"E": [31.3, 0.8], "strength": [1.86, 0.07], "strengthKind": "first peak"}},
+  {"id": "shev-lcd-diamond", "source": "shevchenko2023", "material": "shevchenko-resin-lcd", "process": "LCD resin", "lattice": {"type": "diamond", "cell": 10.0, "rho": 0.33, "wall": 0.8}, "size": [30, 30, 30], "measured": {"E": [113.9, 2.1], "strength": [6.40, 0.11], "strengthKind": "first peak"}},
+  {"id": "shev-lcd-gyroid", "source": "shevchenko2023", "material": "shevchenko-resin-lcd", "process": "LCD resin", "lattice": {"type": "gyroid", "cell": 10.0, "rho": 0.24, "wall": 0.8}, "size": [30, 30, 30], "measured": {"E": [68.4, 0.9], "strength": [3.88, 0.06], "strengthKind": "first peak"}},
+
+  {"id": "gs-001", "source": "guerrasilva2021", "material": "guerrasilva-abs", "process": "FDM ABS 0.25 mm", "lattice": {"type": "kelvin", "cell": 5.0, "diameter": 1.0}, "size": [35, 35, 35], "measured": {"E": [63.43, 4.72], "strength": [4.17, 0.03], "strengthKind": "plateau"}},
+  {"id": "gs-002", "source": "guerrasilva2021", "material": "guerrasilva-abs", "process": "FDM ABS 0.15 mm", "lattice": {"type": "kelvin", "cell": 7.0, "diameter": 1.0}, "size": [35, 35, 35], "measured": {"E": [15.48, 0.30], "strength": [1.00, 0.05], "strengthKind": "plateau"}},
+  {"id": "gs-003", "source": "guerrasilva2021", "material": "guerrasilva-abs", "process": "FDM ABS 0.15 mm", "lattice": {"type": "diamond-strut", "cell": 5.0, "diameter": 1.0}, "size": [35, 35, 35], "measured": {"E": [19.65, 2.33], "strength": [1.10, 0.01], "strengthKind": "plateau"}},
+  {"id": "gs-004", "source": "guerrasilva2021", "material": "guerrasilva-abs", "process": "FDM ABS 0.25 mm", "lattice": {"type": "diamond-strut", "cell": 7.0, "diameter": 1.0}, "size": [35, 35, 35], "measured": {"E": [4.49, 0.10], "strength": [0.28, 0.02], "strengthKind": "plateau"}},
+  {"id": "gs-005", "source": "guerrasilva2021", "material": "guerrasilva-abs", "process": "FDM ABS 0.15 mm", "lattice": {"type": "kelvin", "cell": 5.0, "diameter": 1.5}, "size": [35, 35, 35], "measured": {"E": [89.62, 0.66], "strength": [7.82, 0.21], "strengthKind": "plateau"}},
+  {"id": "gs-006", "source": "guerrasilva2021", "material": "guerrasilva-abs", "process": "FDM ABS 0.25 mm", "lattice": {"type": "kelvin", "cell": 7.0, "diameter": 1.5}, "size": [35, 35, 35], "measured": {"E": [42.27, 1.34], "strength": [2.68, 0.03], "strengthKind": "plateau"}},
+  {"id": "gs-007", "source": "guerrasilva2021", "material": "guerrasilva-abs", "process": "FDM ABS 0.25 mm", "lattice": {"type": "diamond-strut", "cell": 5.0, "diameter": 1.5}, "size": [35, 35, 35], "measured": {"E": [57.57, 3.00], "strength": [3.08, 0.02], "strengthKind": "plateau"}},
+  {"id": "gs-008", "source": "guerrasilva2021", "material": "guerrasilva-abs", "process": "FDM ABS 0.15 mm", "lattice": {"type": "diamond-strut", "cell": 7.0, "diameter": 1.5}, "size": [35, 35, 35], "measured": {"E": [18.21, 2.47], "strength": [1.11, 0.01], "strengthKind": "plateau"}},
+  {"id": "gs-009", "source": "guerrasilva2021", "material": "guerrasilva-pc", "process": "FDM PC 0.15 mm", "lattice": {"type": "kelvin", "cell": 5.0, "diameter": 1.0}, "size": [35, 35, 35], "measured": {"E": [60.42, 4.55], "strength": [4.79, 0.03], "strengthKind": "plateau"}},
+  {"id": "gs-010", "source": "guerrasilva2021", "material": "guerrasilva-pc", "process": "FDM PC 0.25 mm", "lattice": {"type": "kelvin", "cell": 7.0, "diameter": 1.0}, "size": [35, 35, 35], "measured": {"E": [20.39, 2.48], "strength": [1.52, 0.03], "strengthKind": "plateau"}},
+  {"id": "gs-011", "source": "guerrasilva2021", "material": "guerrasilva-pc", "process": "FDM PC 0.25 mm", "lattice": {"type": "diamond-strut", "cell": 5.0, "diameter": 1.0}, "size": [35, 35, 35], "measured": {"E": [21.19, 1.39], "strength": [1.55, 0.07], "strengthKind": "plateau"}},
+  {"id": "gs-012", "source": "guerrasilva2021", "material": "guerrasilva-pc", "process": "FDM PC 0.15 mm", "lattice": {"type": "diamond-strut", "cell": 7.0, "diameter": 1.0}, "size": [35, 35, 35], "measured": {"E": [3.05, 1.02], "strength": [0.29, 0.01], "strengthKind": "plateau"}},
+  {"id": "gs-013", "source": "guerrasilva2021", "material": "guerrasilva-pc", "process": "FDM PC 0.25 mm", "lattice": {"type": "kelvin", "cell": 5.0, "diameter": 1.5}, "size": [35, 35, 35], "measured": {"E": [130.66, 0.22], "strength": [12.84, 0.13], "strengthKind": "plateau"}},
+  {"id": "gs-014", "source": "guerrasilva2021", "material": "guerrasilva-pc", "process": "FDM PC 0.15 mm", "lattice": {"type": "kelvin", "cell": 7.0, "diameter": 1.5}, "size": [35, 35, 35], "measured": {"E": [33.72, 1.90], "strength": [2.48, 0.93], "strengthKind": "plateau"}},
+  {"id": "gs-015", "source": "guerrasilva2021", "material": "guerrasilva-pc", "process": "FDM PC 0.15 mm", "lattice": {"type": "diamond-strut", "cell": 5.0, "diameter": 1.5}, "size": [35, 35, 35], "measured": {"E": [59.83, 0.68], "strength": [3.88, 0.06], "strengthKind": "plateau"}},
+  {"id": "gs-016", "source": "guerrasilva2021", "material": "guerrasilva-pc", "process": "FDM PC 0.25 mm", "lattice": {"type": "diamond-strut", "cell": 7.0, "diameter": 1.5}, "size": [35, 35, 35], "measured": {"E": [25.60, 4.62], "strength": [1.54, 0.01], "strengthKind": "plateau"}},
+
+  {"id": "aa-sheet-gyroid-10", "source": "abouali2022", "material": "abouali-pa11-sls", "process": "SLS PA11", "lattice": {"type": "gyroid", "cell": 8.0, "rho": 0.104}, "size": [40, 40, 40], "measured": {"rho": 0.181, "fit": {"E": [329.7, 1.08], "strength": [8.9, 1.01]}, "strengthKind": "yield (fit)"}},
+  {"id": "aa-sheet-gyroid-18", "source": "abouali2022", "material": "abouali-pa11-sls", "process": "SLS PA11", "lattice": {"type": "gyroid", "cell": 8.0, "rho": 0.182}, "size": [40, 40, 40], "measured": {"rho": 0.250, "fit": {"E": [329.7, 1.08], "strength": [8.9, 1.01]}, "strengthKind": "yield (fit)"}},
+  {"id": "aa-sheet-diamond-13", "source": "abouali2022", "material": "abouali-pa11-sls", "process": "SLS PA11", "lattice": {"type": "diamond", "cell": 8.0, "rho": 0.126}, "size": [40, 40, 40], "measured": {"rho": 0.279, "fit": {"E": [225.5, 0.75], "strength": [17.4, 1.52]}, "strengthKind": "yield (fit)"}},
+  {"id": "aa-sheet-diamond-21", "source": "abouali2022", "material": "abouali-pa11-sls", "process": "SLS PA11", "lattice": {"type": "diamond", "cell": 8.0, "rho": 0.209}, "size": [40, 40, 40], "measured": {"rho": 0.342, "fit": {"E": [225.5, 0.75], "strength": [17.4, 1.52]}, "strengthKind": "yield (fit)"}},
+  {"id": "aa-octet-11", "source": "abouali2022", "material": "abouali-pa11-sls", "process": "SLS PA11", "lattice": {"type": "octet", "cell": 8.0, "rho": 0.114}, "size": [40, 40, 40], "measured": {"rho": 0.118, "fit": {"E": [79.0, 0.74], "strength": [4.5, 0.91]}, "strengthKind": "yield (fit)"}},
+  {"id": "aa-octet-25", "source": "abouali2022", "material": "abouali-pa11-sls", "process": "SLS PA11", "lattice": {"type": "octet", "cell": 8.0, "rho": 0.25}, "size": [40, 40, 40], "measured": {"rho": 0.284, "fit": {"E": [79.0, 0.74], "strength": [4.5, 0.91]}, "strengthKind": "yield (fit)"}},
+  {"id": "aa-kelvin-14", "source": "abouali2022", "material": "abouali-pa11-sls", "process": "SLS PA11", "lattice": {"type": "kelvin", "cell": 8.0, "rho": 0.14}, "size": [40, 40, 40], "measured": {"rho": 0.153, "fit": {"E": [808.5, 1.85], "strength": [32.7, 2.00]}, "strengthKind": "yield (fit)"}},
+  {"id": "aa-kelvin-25", "source": "abouali2022", "material": "abouali-pa11-sls", "process": "SLS PA11", "lattice": {"type": "kelvin", "cell": 8.0, "rho": 0.247}, "size": [40, 40, 40], "measured": {"rho": 0.276, "fit": {"E": [808.5, 1.85], "strength": [32.7, 2.00]}, "strengthKind": "yield (fit)"}}
+ ]
+}

loadpath/decimate.py

@@ -0,0 +1,177 @@
+"""Adaptive decimation of a field mesh against the field itself.
+
+Marching cubes puts triangles a pitch apart everywhere: a flat face gets
+as many as a fillet. `decimate` collapses edges, shortest first, wherever
+the result stays within `tol` mm of the true surface, so flat and gently
+curved regions thin out and creases, small fillets and thin lattice necks
+keep their triangles.
+
+Each pass works on many edges at once (numpy, no priority queue):
+
+1. candidates: edges whose two ends share exactly two neighbours (the
+   link condition, so the collapse keeps the mesh a closed manifold) and
+   touch no locked vertex;
+2. an independent set: an edge goes ahead when it is the shortest
+   candidate around both its ends and their neighbours, so no two
+   collapses in a pass share a triangle;
+3. each collapse moves the pair to its midpoint projected onto the field's
+   zero level (two Newton steps) and is kept only if every triangle it
+   reshapes stays within `tol` of the surface at its centroid and edge
+   midpoints (|field| there, which is the distance for the distance-like
+   fields in fields.py) and does not turn over (normal within ~60 degrees
+   of before).
+
+Passes repeat until one removes under 1% of the triangles. The mesh
+stays watertight, and every vertex stays on the surface.
+"""
+from __future__ import annotations
+
+import numpy as np
+
+__all__ = ["decimate"]
+
+
+def _normals(T):
+    n = np.cross(T[:, 1] - T[:, 0], T[:, 2] - T[:, 0])
+    return n / np.maximum(np.linalg.norm(n, axis=1, keepdims=True), 1e-300)
+
+
+def decimate(V, F, field, tol, pitch, locked=None, max_passes=80, min_gain=0.01, progress=None,
+             target=None, place=None, turn=0.5, labels=None, allow=None, max_reject=3, quality=None):
+    """(V, F, stats, kept): the mesh with edges collapsed wherever it stays within
+    `tol` mm of field = 0. `locked` vertices (a mask) never move or go:
+    the streamed export's seams. Unused vertices are dropped (locked ones
+    keep their order relative to each other); `kept` is the mask of the
+    input vertices that remain.
+
+    With `field` None there is no surface to hold to: a collapse of edges
+    (a, b) puts the new vertex at `place(V, a, b)` (the midpoints if not
+    given), and passes stop once the mesh is down to `target` faces.
+    Shortest edges first makes that an even coarsening. `turn` is the
+    least cosine between a reshaped face's normal and its old one.
+    `labels` (one per face) ride along with the faces and come back in
+    stats["labels"]; `allow(E, F, labels)` filters each pass's candidate
+    edges (a mask), say to keep region seams in place. A vertex whose
+    collapses failed `max_reject` times is left alone. `quality(Told, Tnew,
+    labels)` (the reshaped triangles' corners before and after, (k, 3, 3),
+    and their labels) can veto collapses too."""
+    from .fields import _project
+    from .narrowband import peval
+
+    V = np.array(V, float)
+    F = np.array(F, np.int64)
+    n = len(V)
+    lock = np.zeros(n, bool) if locked is None else np.array(locked, bool)
+    alive = np.ones(n, bool)
+    rej = np.zeros(n, np.int64)
+    f0 = len(F)
+    lab = None if labels is None else np.asarray(labels).copy()
+    rng = np.random.default_rng(0)
+    passes = 0
+    for passes in range(1, max_passes + 1):
+        nf = len(F)
+        a = np.r_[F[:, 0], F[:, 1], F[:, 2]]
+        b = np.r_[F[:, 1], F[:, 2], F[:, 0]]
+        # closed manifold: every edge is a half-edge each way, so the edges
+        # are the half-edges with a < b
+        h = a < b
+        E = np.stack([a[h], b[h]], axis=1)
+        # an edge that failed before yields to its neighbours, or it would
+        # block them pass after pass; a vertex that failed three times is done
+        frozen = lock | (rej >= max_reject)
+        E = E[~(frozen[E[:, 0]] | frozen[E[:, 1]])]
+        if allow is not None and len(E):
+            E = E[allow(E, F, lab)]
+        if not len(E):
+            break
+        d = V[E[:, 0]] - V[E[:, 1]]
+        prio = (np.einsum("ij,ij->i", d, d) * np.exp2(2 * (rej[E[:, 0]] + rej[E[:, 1]]))
+                * (1.0 + 1e-6 * rng.random(len(E))))
+        # independent set: shortest around both ends and their neighbours
+        m = np.full(n, np.inf)
+        np.minimum.at(m, E[:, 0], prio)
+        np.minimum.at(m, E[:, 1], prio)
+        M = m.copy()
+        np.minimum.at(M, a, m[b])
+        win = (prio == M[E[:, 0]]) & (prio == M[E[:, 1]])
+        W = E[win]
+        if not len(W):
+            break
+        # the triangles around each collapse: no triangle holds ends of two
+        # collapses (their ends are not neighbours), so each belongs to at
+        # most one
+        wid = np.full(n, -1)
+        wid[W[:, 0]] = np.arange(len(W))
+        wid[W[:, 1]] = np.arange(len(W))
+        fw = wid[F].max(axis=1)
+        fi = np.flatnonzero(fw >= 0)
+        tri, w = F[fi], fw[fi]
+        wa, wb = W[w, 0], W[w, 1]
+        in_a = (tri == wa[:, None]).any(axis=1)
+        in_b = (tri == wb[:, None]).any(axis=1)
+        # link condition: the two ends share exactly the two opposite vertices
+        nk = []
+        for side in (in_a, in_b):
+            t, ww = tri[side], w[side]
+            v = t.ravel()
+            k = np.repeat(ww, 3) * n + v
+            k = k[(v != np.repeat(W[ww, 0], 3)) & (v != np.repeat(W[ww, 1], 3))]
+            nk.append(np.unique(k))
+        common = np.bincount(np.intersect1d(nk[0], nk[1], assume_unique=True) // n, minlength=len(W))
+        link = common == 2
+        rej[W[~link].ravel()] += 1
+        sel = link[w] & ~(in_a & in_b)            # the two triangles on the edge go
+        tri, w, wa, wb = tri[sel], w[sel], wa[sel], wb[sel]
+        pr = np.stack([w, fi[sel]], axis=1)
+        # the new vertex: midpoint, onto the surface
+        P = np.zeros((len(W), 3))
+        li = np.flatnonzero(link)
+        if len(li) and field is None:
+            P[li] = (place(V, W[li, 0], W[li, 1]) if place is not None
+                     else 0.5 * (V[W[li, 0]] + V[W[li, 1]]))
+        elif len(li):
+            q = _project(field, 0.5 * (V[W[li, 0]] + V[W[li, 1]]), pitch)
+            P[li] = _project(field, q, pitch)
+        Told = V[tri]
+        Tnew = Told.copy()
+        hit = (tri == wa[:, None]) | (tri == wb[:, None])
+        Tnew[hit] = np.repeat(P[pr[:, 0]], hit.sum(axis=1), axis=0)
+        good = (_normals(Told) * _normals(Tnew)).sum(axis=1) > turn
+        if quality is not None and len(Tnew):
+            good &= quality(Told, Tnew, None if lab is None else lab[pr[:, 1]])
+        if field is not None:
+            S = np.concatenate([Tnew.mean(axis=1), 0.5 * (Tnew[:, 0] + Tnew[:, 1]),
+                                0.5 * (Tnew[:, 1] + Tnew[:, 2]), 0.5 * (Tnew[:, 2] + Tnew[:, 0])])
+            dev = np.abs(peval(field, S)).reshape(4, -1).max(axis=0)
+            good &= dev <= tol
+        bad = ~link
+        np.logical_or.at(bad, pr[:, 0], ~good)
+        acc = ~bad
+        rej[W[bad & link].ravel()] += 1
+        if target is not None:              # each collapse takes two triangles: do not overshoot
+            room = max(0, (len(F) - int(target)) // 2)
+            acc[np.flatnonzero(acc)[room:]] = False
+        if not acc.any():
+            break
+        Wa, Wb = W[acc, 0], W[acc, 1]
+        V[Wa] = P[acc]
+        remap = np.arange(n)
+        remap[Wb] = Wa
+        alive[Wb] = False
+        F = remap[F]
+        keepf = (F[:, 0] != F[:, 1]) & (F[:, 1] != F[:, 2]) & (F[:, 0] != F[:, 2])
+        F = F[keepf]
+        if lab is not None:
+            lab = lab[keepf]
+        if progress:
+            progress({"stage": "decimate", "pass": passes, "faces": int(len(F))})
+        if nf - len(F) < min_gain * nf or (target is not None and len(F) <= target):
+            break
+    used = np.zeros(n, bool)
+    used[F] = True
+    used |= lock & alive
+    new = np.cumsum(used) - 1
+    stats = {"tolerance": float(tol), "facesBefore": int(f0), "faces": int(len(F)), "passes": int(passes)}
+    if lab is not None:
+        stats["labels"] = lab
+    return V[used], new[F], stats, used

loadpath/delivery.py

@@ -0,0 +1,190 @@
+"""
+The part to hand over, in one place: report.json["deliver"].
+
+A run's report lists every outcome it tried; an agent wants the one it
+should give the user, whether it passes, which number to quote as its safety
+factor (and from which analysis), what set its size, and every warning that
+matters for it, gathered from the outcome, its print and process checks, its
+durability and the run. Agents testing the brief evals kept missing print and
+fatigue warnings that sat three levels down, and quoting the voxel safety
+factor where the report had replaced it with a finer one.
+"""
+from __future__ import annotations
+
+_VF_CEILING = 0.7   # study._VF_MAX: the search never goes past it
+_VF_FLOOR = 0.02    # study._VF_MIN
+
+
+def attention(o: dict, study: dict) -> list[str]:
+    """What SKILL.md tells an agent to look for, already looked for."""
+    out = []
+    if o.get("passes") is False:
+        # which limit it misses, so an agent can tell a weak part from a floppy one
+        goal = study.get("goal") or {}
+        why = []
+        sf, want = o.get("safetyFactor"), goal.get("safetyFactor")
+        if want and sf is not None and sf < want:
+            why.append(f"safety factor {sf:.2f} under {want:g}")
+        dmax = goal.get("maxDeflection")
+        defl = max((c.get("deflection") or 0.0 for c in o.get("cases") or []), default=0.0)
+        if dmax and defl > dmax:
+            why.append(f"deflection {defl:.2f} mm over {dmax:g}")
+        out.append("does not meet the goal" + (f" ({', '.join(why)})" if why else "") +
+                   ": raise the volume fraction or the material, or check the loads")
+    if (o.get("bodies") or 1) > 1:
+        out.append(f"{o['bodies']} separate bodies: something is not connected to a support")
+    if (o.get("keepOutIntrusion") or 0) > 1.0:
+        out.append(f"keepOutIntrusion {o['keepOutIntrusion']:.1f} mm³: material in a keep-out (report a bug)")
+    if (o.get("preserveMissing") or 0) > 1.0:
+        out.append(f"preserveMissing {o['preserveMissing']:.1f} mm³: preserved geometry lost (report a bug)")
+    mm = study.get("minMember")
+    p5 = _p5(o)
+    if mm and p5 is not None and p5 < 0.8 * mm:
+        out.append(f"p5Member {p5:.2f} mm is under minMember {mm} mm: raise the resolution")
+    if o.get("closedBores"):
+        out.append(f"closed bores: {o['closedBores']}")
+    dur = o.get("durability") or {}
+    if dur.get("safetyFactor") is not None and dur["safetyFactor"] < dur.get("target", 1.5):
+        out.append(f"durability safety factor {dur['safetyFactor']:.2f} (fatigue/creep) is under "
+                   f"{dur.get('target', 1.5):g}: set durability.constrain to size for it")
+    out.extend(dur.get("warnings") or [])
+    return out
+
+
+def _p5(o: dict):
+    if o.get("p5Member") is not None:
+        return o["p5Member"]
+    return (o.get("confirm") or {}).get("p5Member")
+
+
+def _sf_method(o: dict) -> str:
+    c = o.get("confirm")
+    if not c:
+        return "voxel"
+    return c.get("method") or "voxel-fine"
+
+
+def _governing(o: dict, goal: dict) -> str | None:
+    """Which of the goal's limits the delivered part is closest to."""
+    ms = {}
+    if goal.get("safetyFactor") and o.get("safetyFactor"):
+        ms["safetyFactor"] = o["safetyFactor"] / goal["safetyFactor"]
+    if goal.get("maxDeflection"):
+        d = max((c.get("deflection") or 0.0) for c in o.get("cases") or [{}])
+        if d > 0:
+            ms["maxDeflection"] = goal["maxDeflection"] / d
+    dur = o.get("durability") or {}
+    if dur.get("safetyFactor") is not None and dur.get("constrain"):
+        ms["durability"] = dur["safetyFactor"] / dur.get("target", 1.5)
+    return min(ms, key=ms.get) if ms else None
+
+
+def delivery(summary: dict, study: dict) -> dict | None:
+    """report.json["deliver"]: the best outcome, flattened, with every
+    warning that applies to it."""
+    outs = summary.get("outcomes") or []
+    best = summary.get("best")
+    if best is None:
+        if not outs:
+            return None
+        best = max(range(len(outs)), key=lambda i: outs[i].get("safetyFactor") or 0.0)
+    o = outs[best]
+    goal = study.get("goal") or {}
+    warn = []
+    for w in (list(summary.get("warnings") or []) + list(o.get("warnings") or [])
+              + list((o.get("print") or {}).get("warnings") or [])
+              + list((o.get("process") or {}).get("warnings") or []) + attention(o, study)):
+        if w not in warn:
+            warn.append(w)
+    peaks = {c.get("name"): c.get("peakStress") for c in (o.get("confirm") or {}).get("cases") or []}
+    cases = [{"name": c.get("name"), "safetyFactor": c.get("safetyFactor"), "deflection": c.get("deflection"),
+              "peakStress": c.get("peakStress") if c.get("peakStress") is not None else peaks.get(c.get("name")),
+              **({"peakAt": c["peakAt"]} if c.get("peakAt") is not None else {}),
+              **({"worstDirection": c["worstDirection"]} if c.get("worstDirection") else {})}
+             for c in o.get("cases") or []]
+    d = {
+        "outcome": best,
+        "passes": o.get("passes") if summary.get("best") is not None else False,
+        "stl": o.get("stl"),
+        **({"3mf": o["3mf"]} if o.get("3mf") else {}),
+        **({"render": o["render"]} if o.get("render") else {}),
+        "material": o.get("material"),
+        "mass": o.get("mass"),
+        "volumeFraction": o.get("volumeFraction"),
+        "safetyFactor": o.get("safetyFactor"),
+        # which analysis that number comes from: "tet" (body-fitted cross-check),
+        # "voxel-fine" (a finer voxel re-analysis) or "voxel" (the design grid)
+        "safetyFactorMethod": _sf_method(o),
+        **({"target": {k: goal[k] for k in ("safetyFactor", "maxDeflection", "maxMass") if k in goal}}
+           if any(k in goal for k in ("safetyFactor", "maxDeflection", "maxMass")) else {}),
+        "governedBy": _governing(o, goal),
+        "maxDeflection": max((c.get("deflection") or 0.0) for c in o.get("cases") or [{}]) or None,
+        "cases": cases,
+        "p5Member": _p5(o),
+        **({"durabilitySafetyFactor": o["durability"].get("safetyFactor")} if o.get("durability") else {}),
+        **{k: o[k] for k in ("bodies", "keepOutIntrusion", "preserveMissing") if o.get(k) is not None},
+        **({"process": o["process"]} if o.get("process") else {}),
+        **({"print": {k: v for k, v in o["print"].items() if k != "warnings"}} if o.get("print") else {}),
+        "warnings": warn,
+    }
+    cap = goal.get("maxMass")
+    if cap and d.get("mass") is not None and d["mass"] > cap:
+        d["passes"] = False
+        d["warnings"].insert(0, f"over the mass cap: {d['mass']:.1f} g against goal.maxMass {cap:g} g; this is "
+                                "already the lightest that meets the other limits, so use a lighter material, "
+                                "trim the preserves or lower the safety factor")
+    dur = o.get("durability") or {}
+    if d["passes"] and dur.get("safetyFactor") is not None and dur["safetyFactor"] < dur.get("target", 1.5):
+        # a permanently loaded PLA bracket read static SF 2.6 and passes: true
+        # at creep SF 1.09; "passes" means the whole part, durability included
+        d["passes"] = False
+        d["warnings"].insert(0, f"fails on durability: fatigue/creep safety factor {dur['safetyFactor']:.2f} is "
+                                f"under {dur.get('target', 1.5):g} (the static goal passes); set "
+                                "durability.constrain: true to size the part for it")
+    proc = study.get("process")
+    kind = proc.get("type") if isinstance(proc, dict) else proc
+    if str(kind or "").lower() == "fdm" and "slicer" not in o:
+        d["slicer"] = "not run: no slicer found or none asked for (install PrusaSlicer or set LOADPATH_SLICER " \
+                      "for print time and filament)"
+    # the one number to quote: the factor of whichever limit sets the part
+    d["governingSafetyFactor"] = d.get("durabilitySafetyFactor") if d["governedBy"] == "durability" \
+        else d.get("safetyFactor")
+    if d["governedBy"] == "durability" and d.get("durabilitySafetyFactor") is not None:
+        d["warnings"].insert(0, f"durability (fatigue/creep) sets this part: quote its safety factor "
+                                f"{d['durabilitySafetyFactor']:.2f}, not the static {d['safetyFactor']:.2f}")
+    if summary.get("best") is not None and (o.get("volumeFraction") or 1.0) <= _VF_FLOOR + 1e-6:
+        d["governedBy"] = "volumeFraction floor"
+        d["warnings"].insert(0, f"the search stopped at its smallest volume fraction ({_VF_FLOOR:.0%}) with the "
+                                "safety factor still over target: the preserved interfaces, not the loads, set "
+                                "the mass; trim the preserves or the design space to go lighter")
+    if summary.get("best") is None:
+        d["warnings"].insert(0, "no outcome meets the goal: this is the strongest one tried, not a passing part")
+        sat = _saturated(outs)
+        if sat:
+            d["warnings"].insert(1, sat)
+        top = max((x.get("volumeFraction") or 0.0) for x in outs)
+        if top >= _VF_CEILING - 1e-6:
+            d["warnings"].insert(1, f"the search stopped at the volume-fraction ceiling ({_VF_CEILING:.0%} of the design "
+                                    "space): more material is not allowed, so change the material, the design space or the loads")
+    return d
+
+
+def _saturated(outs: list) -> str | None:
+    """When more material hardly helps, the shape is not the problem: a notch
+    (a keep-out's corner, a thin interface) caps the safety factor."""
+    pts = sorted((o["volumeFraction"], o["safetyFactor"], o) for o in outs
+                 if o.get("volumeFraction") and o.get("safetyFactor"))
+    import math
+
+    # a sound layout's safety factor grows at least as fast as material
+    # (roughly vf^1 to vf^2); under vf^0.6 the material is not where it fails
+    if len(pts) < 2 or pts[-1][0] < 1.3 * pts[0][0] or pts[-1][1] <= 0 or pts[0][1] <= 0 or \
+            math.log(pts[-1][1] / pts[0][1]) > 0.6 * math.log(pts[-1][0] / pts[0][0]):
+        return None
+    worst = min(pts[-1][2].get("cases") or [{}], key=lambda c: c.get("safetyFactor") or float("inf"))
+    at = worst.get("peakAt")
+    where = f" at {[round(v, 1) for v in at]} (case {worst.get('name')!r})" if at else ""
+    return (f"the safety factor barely rises with material ({pts[0][1]:.2f} at {pts[0][0]:.0%} of the design "
+            f"space, {pts[-1][1]:.2f} at {pts[-1][0]:.0%}): a local hot spot{where} sets it, so look at the "
+            "geometry there (a keep-out notching an interface, a thin preserved ring, a fastener's head path) "
+            "rather than adding material")