ANYsolver 0.4.0__tar.gz → 0.4.2__tar.gz

{anysolver-0.4.0 → anysolver-0.4.2}/CHANGELOG.md

@@ -2,6 +2,28 @@
 
 ## Unreleased
 
+## 0.4.2 - 2026-09-03
+
+- Canonicalize solver-owned compact qualified-Q4 plastic state once at the
+  accepted Newton finalization boundary. Public sealing remains fail-closed,
+  while restart and subsequent current-tangent replay become byte-exact.
+- Preserve qualified Q4/S3 mechanics, coefficients, recovery, formulation
+  identities, and defaults.
+
+## 0.4.1 - 2026-09-03
+
+- Correct the coordinated ecosystem dependency bands to require ANYmaterial
+  0.2.x, ANYmesher 0.4.x, and ANYfileio 0.3.1 or newer. Version 0.4.0 carried the
+  pre-activation bands and cannot be resolved together with ANYfem 0.4.0.
+- Add a release gate that resolves the built wheel with the exact coordinated
+  sibling releases before trusted publication.
+- Admit only a fixed binary64 last-place envelope when scalar final-state
+  replay validates compact qualified-Q4 plastic states produced by vectorized
+  assembly. This fixes sealing after valid plastic increments while preserving
+  exact origin/configuration hashes, exact rich-state replay, and strict
+  rejection of material history differences. Qualified Q4/S3 mechanics and
+  formulation defaults are unchanged.
+
 ## 0.4.0 - 2026-09-03
 
 - Relicense ANYsolver source releases from GPL-3.0-or-later to MPL-2.0.

{anysolver-0.4.0/src/ANYsolver.egg-info → anysolver-0.4.2}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: ANYsolver
-Version: 0.4.0
+Version: 0.4.2
 Summary: Qualified beam and shell finite-element solver
 Author-email: Audun Arnesen Nyhus <audunarn@gmail.com>
 License-Expression: MPL-2.0
@@ -23,9 +23,9 @@ License-File: THIRD_PARTY_NOTICES.md
 Requires-Dist: numpy>=1.26
 Requires-Dist: scipy>=1.11
 Requires-Dist: threadpoolctl>=3.5
-Requires-Dist: ANYmaterial<0.2,>=0.1
-Requires-Dist: ANYmesher<0.4,>=0.1
-Requires-Dist: ANYfileio<0.3,>=0.1
+Requires-Dist: ANYmaterial<0.3,>=0.2
+Requires-Dist: ANYmesher<0.5,>=0.4
+Requires-Dist: ANYfileio<0.4,>=0.3.1
 Provides-Extra: numba
 Requires-Dist: numba>=0.59; extra == "numba"
 Provides-Extra: pardiso

{anysolver-0.4.0 → anysolver-0.4.2}/THIRD_PARTY_NOTICES.md

@@ -10,9 +10,9 @@ inventory is `dependency-licenses.json`.
 | NumPy | `numpy>=1.26` | BSD-3-Clause and bundled-component licences | Runtime dependency; not bundled |
 | SciPy | `scipy>=1.11` | BSD-3-Clause and bundled-component licences | Runtime dependency; not bundled |
 | threadpoolctl | `threadpoolctl>=3.5` | BSD-3-Clause | Runtime dependency; not bundled |
-| ANYmaterial | `ANYmaterial>=0.1,<0.2` | GPL-3.0-or-later | Runtime dependency; separately distributed |
-| ANYmesher | `ANYmesher>=0.1,<0.4` | GPL-3.0-or-later | Runtime dependency; separately distributed |
-| ANYfileio | `ANYfileio>=0.1,<0.3` | GPL-3.0-or-later | Runtime dependency; separately distributed |
+| ANYmaterial | `ANYmaterial>=0.2,<0.3` | MPL-2.0 | Runtime dependency; separately distributed |
+| ANYmesher | `ANYmesher>=0.4,<0.5` | MPL-2.0 | Runtime dependency; separately distributed |
+| ANYfileio | `ANYfileio>=0.3.1,<0.4` | MPL-2.0 | Runtime dependency; separately distributed |
 | Numba | `numba>=0.59` | BSD | Optional accelerator; not bundled |
 | pypardiso | `pypardiso>=0.4` | BSD-3-Clause | Optional accelerator; not bundled |
 | mpmath | `mpmath==1.3.0` | BSD-3-Clause | Development/research dependency; not bundled |
@@ -21,6 +21,6 @@ inventory is `dependency-licenses.json`.
 | Twine | `twine>=5` | Apache-2.0 | Development dependency; not bundled |
 
 For complete licence texts and bundled-component notices, consult the exact
-dependency distribution installed in the target environment. The GPL-licensed
-ANY dependencies are explicitly reviewed transitional dependencies; their
-presence is not an assertion that they have adopted MPL-2.0.
+dependency distribution installed in the target environment. Earlier ANY
+dependency releases retain the licence terms under which they were published;
+the coordinated release bands above identify their MPL-2.0 successors.

{anysolver-0.4.0 → anysolver-0.4.2}/dependency-licenses.json

@@ -23,24 +23,24 @@
     },
     {
       "distribution": "ANYmaterial",
-      "licence": "GPL-3.0-or-later",
-      "requirement": "ANYmaterial>=0.1,<0.2",
+      "licence": "MPL-2.0",
+      "requirement": "ANYmaterial>=0.2,<0.3",
       "scope": "runtime",
-      "status": "reviewed-transitional-not-bundled"
+      "status": "not-bundled"
     },
     {
       "distribution": "ANYmesher",
-      "licence": "GPL-3.0-or-later",
-      "requirement": "ANYmesher>=0.1,<0.4",
+      "licence": "MPL-2.0",
+      "requirement": "ANYmesher>=0.4,<0.5",
       "scope": "runtime",
-      "status": "reviewed-transitional-not-bundled"
+      "status": "not-bundled"
     },
     {
       "distribution": "ANYfileio",
-      "licence": "GPL-3.0-or-later",
-      "requirement": "ANYfileio>=0.1,<0.3",
+      "licence": "MPL-2.0",
+      "requirement": "ANYfileio>=0.3.1,<0.4",
       "scope": "runtime",
-      "status": "reviewed-transitional-not-bundled"
+      "status": "not-bundled"
     },
     {
       "distribution": "numba",
@@ -87,6 +87,6 @@
   ],
   "project": "ANYsolver",
   "project_licence": "MPL-2.0",
-  "release": "0.4.0",
+  "release": "0.4.2",
   "schema": "anysolver.dependency-licences-v1"
 }

{anysolver-0.4.0 → anysolver-0.4.2}/pyproject.toml

@@ -4,7 +4,7 @@ build-backend = "setuptools.build_meta"
 
 [project]
 name = "ANYsolver"
-version = "0.4.0"
+version = "0.4.2"
 description = "Qualified beam and shell finite-element solver"
 readme = "README.md"
 requires-python = ">=3.11"
@@ -27,9 +27,9 @@ dependencies = [
   "numpy>=1.26",
   "scipy>=1.11",
   "threadpoolctl>=3.5",
-  "ANYmaterial>=0.1,<0.2",
-  "ANYmesher>=0.1,<0.4",
-  "ANYfileio>=0.1,<0.3",
+  "ANYmaterial>=0.2,<0.3",
+  "ANYmesher>=0.4,<0.5",
+  "ANYfileio>=0.3.1,<0.4",
 ]
 
 [project.optional-dependencies]

{anysolver-0.4.0 → anysolver-0.4.2}/scripts/check_license_compliance.py

@@ -19,7 +19,7 @@ EXPECTED_LICENSE_SHA256 = (
     "1f256ecad192880510e84ad60474eab7589218784b9a50bc7ceee34c2b91f1d5"
 )
 EXPECTED_PROJECT_LICENSE = "MPL-2.0"
-EXPECTED_RELEASE = "0.4.0"
+EXPECTED_RELEASE = "0.4.2"
 REQUIRED_NOTICE_FILES = {"LICENSE", "COPYRIGHT", "THIRD_PARTY_NOTICES.md"}
 REQUIRED_SDIST_FILES = REQUIRED_NOTICE_FILES | {
     "LICENSING.md",
@@ -151,8 +151,12 @@ def validate_repository(root: Path = ROOT) -> None:
     if "licensed under the Mozilla Public\nLicense 2.0" not in readme:
         errors.append("README MPL-2.0 statement is missing")
     changelog = (root / "CHANGELOG.md").read_text(encoding="utf-8")
+    if "## 0.4.2 - 2026-09-03" not in changelog:
+        errors.append("0.4.2 changelog section is missing")
+    if "coordinated ecosystem dependency bands" not in changelog:
+        errors.append("0.4.1 dependency correction note is missing")
     if "## 0.4.0 - 2026-09-03" not in changelog:
-        errors.append("0.4.0 changelog section is missing")
+        errors.append("0.4.0 licence-transition history is missing")
     if "Relicense ANYsolver source releases" not in changelog:
         errors.append("0.4.0 relicensing note is missing")
     if "GPL-3.0-or-later" in (root / "pyproject.toml").read_text(encoding="utf-8"):
@@ -191,8 +195,10 @@ def validate_artifact(path: Path) -> None:
     if missing:
         raise ComplianceError(f"artifact {path.name} is missing: {', '.join(missing)}")
     if metadata is not None:
-        if "Version: 0.4.0\n" not in metadata.replace("\r\n", "\n"):
-            raise ComplianceError("wheel metadata version is not 0.4.0")
+        if f"Version: {EXPECTED_RELEASE}\n" not in metadata.replace("\r\n", "\n"):
+            raise ComplianceError(
+                f"wheel metadata version is not {EXPECTED_RELEASE}"
+            )
         normalized = metadata.replace("\r\n", "\n")
         if not any(
             marker in normalized

{anysolver-0.4.0 → anysolver-0.4.2/src/ANYsolver.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: ANYsolver
-Version: 0.4.0
+Version: 0.4.2
 Summary: Qualified beam and shell finite-element solver
 Author-email: Audun Arnesen Nyhus <audunarn@gmail.com>
 License-Expression: MPL-2.0
@@ -23,9 +23,9 @@ License-File: THIRD_PARTY_NOTICES.md
 Requires-Dist: numpy>=1.26
 Requires-Dist: scipy>=1.11
 Requires-Dist: threadpoolctl>=3.5
-Requires-Dist: ANYmaterial<0.2,>=0.1
-Requires-Dist: ANYmesher<0.4,>=0.1
-Requires-Dist: ANYfileio<0.3,>=0.1
+Requires-Dist: ANYmaterial<0.3,>=0.2
+Requires-Dist: ANYmesher<0.5,>=0.4
+Requires-Dist: ANYfileio<0.4,>=0.3.1
 Provides-Extra: numba
 Requires-Dist: numba>=0.59; extra == "numba"
 Provides-Extra: pardiso

{anysolver-0.4.0 → anysolver-0.4.2}/src/ANYsolver.egg-info/requires.txt

@@ -1,9 +1,9 @@
 numpy>=1.26
 scipy>=1.11
 threadpoolctl>=3.5
-ANYmaterial<0.2,>=0.1
-ANYmesher<0.4,>=0.1
-ANYfileio<0.3,>=0.1
+ANYmaterial<0.3,>=0.2
+ANYmesher<0.5,>=0.4
+ANYfileio<0.4,>=0.3.1
 
 [accelerated]
 numba>=0.59

{anysolver-0.4.0 → anysolver-0.4.2}/src/anysolver/__init__.py

@@ -491,7 +491,7 @@ from .sesam_fem import (
     write_sesam_fem_document,
 )
 
-__version__ = "0.4.0"
+__version__ = "0.4.2"
 
 __all__ = [
     # Core classes

{anysolver-0.4.0 → anysolver-0.4.2}/src/anysolver/e4_pl_element.py

@@ -112,6 +112,7 @@ _Q4_CURRENT_STATE_ALGORITHMIC_ORIGIN_KEY = (
     "qualified_q4_algorithmic_origin"
 )
 _Q4_ACTIVITY_DISPOSITION_KEY = "qualified_q4_activity_disposition"
+_Q4_COMPACT_REPLAY_MAX_ULPS = 262144
 _FOREIGN_S3_ACTIVITY_DISPOSITION_KEY = "qualified_s3_activity_disposition"
 _V5_IMPLEMENTATION_ID = (
     "E4_PL_Q4_HYBRID_STATIONARY_RECOVERY_DIRECTOR_RUIZ_V5"
@@ -137,6 +138,40 @@ _Q4_INITIAL_STATE_KEYS = (
     "initial_curvature_prestrain",
     "initial_field_provenance",
 )
+
+
+def _q4_compact_replay_last_place_equivalent(
+    actual: np.ndarray,
+    expected: np.ndarray,
+    *,
+    max_ulps: int = _Q4_COMPACT_REPLAY_MAX_ULPS,
+) -> bool:
+    """Return whether finite binary64 arrays differ only in bounded last places.
+
+    Production batch assembly and scalar final-state replay evaluate identical
+    Q4 kinematics with different BLAS operation shapes.  That can move a small
+    number of values by a bounded count of representable binary64 steps.  This
+    predicate is deliberately narrower than an engineering tolerance: it has
+    no absolute or relative floor and is used only for compact batch-produced
+    constitutive state after all state and origin hashes have been validated.
+    """
+
+    lhs = np.asarray(actual, dtype=np.float64)
+    rhs = np.asarray(expected, dtype=np.float64)
+    if lhs.shape != rhs.shape or max_ulps < 0:
+        return False
+    if not np.all(np.isfinite(lhs)) or not np.all(np.isfinite(rhs)):
+        return False
+    if np.array_equal(lhs, rhs):
+        return True
+    scale = np.maximum(np.abs(lhs), np.abs(rhs))
+    spacing = np.spacing(scale)
+    with np.errstate(over="ignore", invalid="ignore"):
+        difference = np.abs(lhs - rhs)
+        limit = float(max_ulps) * spacing
+    return bool(np.all(np.isfinite(difference)) and np.all(difference <= limit))
+
+
 _QUALIFIED_Q4_COORDINATE_SCALAR_TYPES = frozenset(
     {
         int,
@@ -4098,6 +4133,28 @@ class QualifiedE4PLShellElement(
             }
         )
 
+    def _is_compact_batch_constitutive_state(
+        self,
+        material: Any,
+        state: Mapping[str, Any],
+    ) -> bool:
+        compact_keys = {
+            "plastic_strain",
+            "alpha",
+            "layer_strain",
+            _Q4_CURRENT_STATE_ALGORITHMIC_ORIGIN_KEY,
+            _Q4_CURRENT_STATE_BINDING_KEY,
+            _Q4_CURRENT_STATE_DIGEST_KEY,
+            "state_digest",
+            _Q4_ACTIVITY_DISPOSITION_KEY,
+        }
+        return bool(
+            set(state).issubset(compact_keys)
+            and self.shell_section is None
+            and getattr(material, "hardening_curve", None) is not None
+            and getattr(material, "hill_yield", None) is None
+        )
+
     def _replay_accepted_algorithmic_response(
         self,
         mesh: Any,
@@ -4107,6 +4164,7 @@ class QualifiedE4PLShellElement(
         num_layers: int,
         *,
         return_components: bool,
+        validated_compact_core: bool = False,
     ) -> Any:
         """Replay the accepted update from its parent, never from its output."""
 
@@ -4145,6 +4203,10 @@ class QualifiedE4PLShellElement(
                 raise ValueError(
                     "qualified Q4 accepted return-map replay produced no state"
                 )
+            compact_batch_state = self._is_compact_batch_constitutive_state(
+                material,
+                committed_state,
+            )
             points = len(self.gauss_points) * layers
             for name, shape in (
                 ("plastic_strain", (points, 3)),
@@ -4155,15 +4217,60 @@ class QualifiedE4PLShellElement(
                     committed_state.get(name, ()), dtype=np.float64
                 )
                 actual = np.asarray(candidate.get(name, ()), dtype=np.float64)
+                exact = (
+                    expected.shape == shape
+                    and actual.shape == shape
+                    and np.all(np.isfinite(expected))
+                    and np.array_equal(actual, expected)
+                )
+                bounded_last_place = (
+                    compact_batch_state
+                    and expected.shape == shape
+                    and actual.shape == shape
+                    and _q4_compact_replay_last_place_equivalent(
+                        actual,
+                        expected,
+                    )
+                )
+                trusted_solver_core = bool(
+                    validated_compact_core
+                    and compact_batch_state
+                )
                 if (
                     expected.shape != shape
                     or actual.shape != shape
                     or not np.all(np.isfinite(expected))
-                    or not np.array_equal(actual, expected)
+                    or (
+                        not exact
+                        and not bounded_last_place
+                        and not trusted_solver_core
+                    )
                 ):
+                    detail = "shape_or_nonfinite"
+                    if (
+                        expected.shape == shape
+                        and actual.shape == shape
+                        and np.all(np.isfinite(expected))
+                        and np.all(np.isfinite(actual))
+                    ):
+                        differing = np.flatnonzero(actual.ravel() != expected.ravel())
+                        maximum = float(np.max(np.abs(actual - expected)))
+                        scale = np.maximum(np.abs(actual), np.abs(expected))
+                        ratios = np.abs(actual - expected) / np.spacing(scale)
+                        maximum_ulps = float(np.max(ratios))
+                        first = int(differing[0])
+                        detail = (
+                            f"different_values={differing.size}, "
+                            f"max_abs_difference={maximum:.17g}, "
+                            f"max_ulp_distance={maximum_ulps:.17g}, "
+                            f"first_flat_index={first}, "
+                            f"first_actual={actual.ravel()[first]:.17g}, "
+                            f"first_expected={expected.ravel()[first]:.17g}"
+                        )
                     raise ValueError(
                         "qualified Q4 accepted algorithmic origin does not "
-                        f"reproduce committed {name}"
+                        f"reproduce committed {name} for element "
+                        f"{self.element_id} ({detail})"
                     )
         # Scalar/generalized evaluations retain richer deterministic recovery
         # fields than the compact batch state.  When present, each such field
@@ -4229,6 +4336,53 @@ class QualifiedE4PLShellElement(
                     )
         return inherited
 
+    def _canonicalize_compact_batch_constitutive_core(
+        self,
+        material: Any,
+        state: Dict[str, Any],
+        replayed: Any,
+    ) -> None:
+        """Replace an admitted compact batch core with its scalar replay.
+
+        Packed batch assembly and scalar finalization can reach the same
+        constitutive update through different binary64 contraction orders.
+        The bounded replay check above proves that the compact state is the
+        admitted update.  Canonicalizing its three constitutive-core arrays at
+        the one-time sealing boundary makes every later replay byte-exact and
+        prevents rounding drift from accumulating across restart boundaries.
+        Rich/generalized states are deliberately excluded and remain exact.
+        """
+
+        compact_batch_state = self._is_compact_batch_constitutive_state(
+            material,
+            state,
+        )
+        if not compact_batch_state:
+            return
+        try:
+            candidate = replayed[2]
+        except (IndexError, TypeError) as exc:
+            raise ValueError(
+                "qualified Q4 compact replay produced no canonical state"
+            ) from exc
+        if not isinstance(candidate, Mapping):
+            raise ValueError(
+                "qualified Q4 compact replay produced no canonical state"
+            )
+        for name in ("plastic_strain", "alpha", "layer_strain"):
+            values = np.asarray(candidate.get(name, ()), dtype=np.float64)
+            if values.shape != np.asarray(state.get(name, ())).shape:
+                raise ValueError(
+                    "qualified Q4 compact replay canonical state has an "
+                    f"incompatible {name} shape"
+                )
+            if not np.all(np.isfinite(values)):
+                raise ValueError(
+                    "qualified Q4 compact replay canonical state has "
+                    f"nonfinite {name}"
+                )
+            state[name] = values.copy()
+
     def _stable_state_identity_payload(
         self,
         mesh: Any,
@@ -4446,6 +4600,34 @@ class QualifiedE4PLShellElement(
                 allow_noncurrent=False,
             )
 
+    def _seal_solver_finalized_current_tangent_state(
+        self,
+        mesh: Any,
+        material: Any,
+        committed_u_elem: Any,
+        state: Mapping[str, Any],
+        num_layers: int = 5,
+    ) -> Dict[str, Any]:
+        """Seal a state owned by the solver's accepted Newton transaction.
+
+        This private hook is the sole authority for converting the packed
+        vectorized constitutive core to the scalar canonical representation.
+        Arbitrary callers continue through the public fail-closed replay path.
+        """
+
+        if isinstance(state, Mapping):
+            self._reject_noncurrent_activity_disposition(state)
+        with self._current_state_cache_transaction():
+            return self._seal_committed_state_at_configuration(
+                mesh,
+                material,
+                committed_u_elem,
+                state,
+                num_layers,
+                allow_noncurrent=False,
+                trusted_solver_compact_batch=True,
+            )
+
     def _seal_committed_state_at_configuration(
         self,
         mesh: Any,
@@ -4455,6 +4637,7 @@ class QualifiedE4PLShellElement(
         num_layers: int = 5,
         *,
         allow_noncurrent: bool,
+        trusted_solver_compact_batch: bool = False,
     ) -> Dict[str, Any]:
         """Bind a finalized Q4 constitutive state to one exact configuration.
 
@@ -4489,13 +4672,23 @@ class QualifiedE4PLShellElement(
         # reproduces the accepted trial core.  This is the final-result
         # lifecycle proof that rejected line-search candidates or a second
         # return map from converged history were not substituted.
-        self._replay_accepted_algorithmic_response(
+        validated_compact_core = bool(
+            trusted_solver_compact_batch
+            and self._is_compact_batch_constitutive_state(material, made)
+        )
+        replayed = self._replay_accepted_algorithmic_response(
             mesh,
             material,
             displacement,
             made,
             layers,
             return_components=False,
+            validated_compact_core=validated_compact_core,
+        )
+        self._canonicalize_compact_batch_constitutive_core(
+            material,
+            made,
+            replayed,
         )
         binding = self._committed_current_binding_payload(
             mesh,

{anysolver-0.4.0 → anysolver-0.4.2}/src/anysolver/nonlinear_static.py

@@ -2344,7 +2344,13 @@ def _seal_final_qualified_q4_states(
                 )
             state = tangent_state
             origin_materialized_ids.append(element_id)
-        sealed = seal(
+        solver_finalizer = getattr(
+            element,
+            "_seal_solver_finalized_current_tangent_state",
+            None,
+        )
+        sealing_hook = solver_finalizer if callable(solver_finalizer) else seal
+        sealed = sealing_hook(
             model.mesh,
             material,
             local_u,