ANYsolver 0.4.1__tar.gz → 0.4.2__tar.gz

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

@@ -2,6 +2,14 @@
 
 ## 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

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

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: ANYsolver
-Version: 0.4.1
+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

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

@@ -87,6 +87,6 @@
   ],
   "project": "ANYsolver",
   "project_licence": "MPL-2.0",
-  "release": "0.4.1",
+  "release": "0.4.2",
   "schema": "anysolver.dependency-licences-v1"
 }

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

@@ -4,7 +4,7 @@ build-backend = "setuptools.build_meta"
 
 [project]
 name = "ANYsolver"
-version = "0.4.1"
+version = "0.4.2"
 description = "Qualified beam and shell finite-element solver"
 readme = "README.md"
 requires-python = ">=3.11"

{anysolver-0.4.1 → 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.1"
+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,8 @@ 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.1 - 2026-09-03" not in changelog:
-        errors.append("0.4.1 changelog section is missing")
+    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:

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

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: ANYsolver
-Version: 0.4.1
+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

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

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

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

@@ -4133,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,
@@ -4142,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."""
 
@@ -4180,21 +4203,9 @@ class QualifiedE4PLShellElement(
                 raise ValueError(
                     "qualified Q4 accepted return-map replay produced no state"
                 )
-            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,
-            }
-            compact_batch_state = (
-                set(committed_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
+            compact_batch_state = self._is_compact_batch_constitutive_state(
+                material,
+                committed_state,
             )
             points = len(self.gauss_points) * layers
             for name, shape in (
@@ -4221,11 +4232,19 @@ class QualifiedE4PLShellElement(
                         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 exact and not bounded_last_place)
+                    or (
+                        not exact
+                        and not bounded_last_place
+                        and not trusted_solver_core
+                    )
                 ):
                     detail = "shape_or_nonfinite"
                     if (
@@ -4317,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,
@@ -4534,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,
@@ -4543,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.
 
@@ -4577,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.1 → 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,

{anysolver-0.4.1 → anysolver-0.4.2}/tests/test_anyfileio_version_compatibility.py

@@ -187,7 +187,7 @@ def _assert_complete_source_graphs(roots: dict[str, str]) -> None:
                 "anyfileio": archives / "fileio-current",
             },
             {
-                "anysolver": "0.4.1",
+                "anysolver": "0.4.2",
                 "anymaterial": "0.2.0",
                 "anygeometry": "0.4.2",
                 "anymesher": "0.4.0",
@@ -632,7 +632,7 @@ def test_workflows_pin_compatibility_graph_and_actions() -> None:
 
     assert probe_environment(mesh_job) == "\n".join(
         (
-            '          EXPECTED_ANYSOLVER_VERSION: "0.4.1"',
+            '          EXPECTED_ANYSOLVER_VERSION: "0.4.2"',
             '          EXPECTED_ANYMATERIAL_VERSION: "0.2.0"',
             "          EXPECTED_ANYMESHER_VERSION: ${{ matrix.anymesher-version }}",
             "          EXPECTED_ANYGEOMETRY_VERSION: ${{ matrix.anygeometry-version }}",
@@ -642,7 +642,7 @@ def test_workflows_pin_compatibility_graph_and_actions() -> None:
     )
     assert probe_environment(fileio_job) == "\n".join(
         (
-            '          EXPECTED_ANYSOLVER_VERSION: "0.4.1"',
+            '          EXPECTED_ANYSOLVER_VERSION: "0.4.2"',
             '          EXPECTED_ANYMATERIAL_VERSION: "0.2.0"',
             "          EXPECTED_ANYFILEIO_VERSION: ${{ matrix.anyfileio-version }}",
             "          EXPECTED_ANYMESHER_VERSION: ${{ matrix.anymesher-version }}",

{anysolver-0.4.1 → anysolver-0.4.2}/tests/test_e4_pl_functional_constructor_policy.py

@@ -259,6 +259,10 @@ def test_functional_lane_direct_shell_calls_are_provably_non_q4() -> None:
                 "tests/test_e4_pl_q4_current_tangent.py",
                 "test_vectorized_q4_accepted_tangent_matches_sealed_scalar_replay_exactly",
             ): 1,
+            (
+                "tests/test_e4_pl_q4_current_tangent.py",
+                "test_distorted_vectorized_q4_plastic_history_replays_after_multiple_increments",
+            ): 1,
             ("tests/test_e4_pl_q4_state_lifecycle.py", "_loaded_q4_model"): 1,
             (
                 "tests/test_e4_pl_q4_state_lifecycle.py",
@@ -309,7 +313,7 @@ def test_functional_lane_direct_shell_calls_are_provably_non_q4() -> None:
         }
     )
     assert calls["QualifiedE4PLShellElement"] == expected_qualified_q4
-    assert sum(calls["QualifiedE4PLShellElement"].values()) == 53
+    assert sum(calls["QualifiedE4PLShellElement"].values()) == 54
 
 
 def test_functional_lane_legacy_q4_calls_are_explicit_and_closed_world() -> None:

{anysolver-0.4.1 → anysolver-0.4.2}/tests/test_e4_pl_q4_current_tangent.py

@@ -983,9 +983,20 @@ def test_vectorized_q4_accepted_tangent_matches_sealed_scalar_replay_exactly() -
         3,
         tangent_evaluated=True,
     )
-    sealed = element.seal_committed_current_tangent_state(
-        mesh, material, displacement, trial, 3
+    canonical_plastic = np.asarray(plastic[0]).copy()
+    perturbed = copy.deepcopy(trial)
+    positive = int(np.flatnonzero(canonical_plastic > 0.0)[0])
+    changed = canonical_plastic.ravel().copy()
+    changed[positive] += 568728.0 * np.spacing(abs(changed[positive]))
+    perturbed["plastic_strain"] = changed.reshape(canonical_plastic.shape)
+    with pytest.raises(ValueError, match="does not reproduce committed"):
+        element.seal_committed_current_tangent_state(
+            mesh, material, displacement, perturbed, 3
+        )
+    sealed = element._seal_solver_finalized_current_tangent_state(
+        mesh, material, displacement, perturbed, 3
     )
+    np.testing.assert_array_equal(sealed["plastic_strain"], canonical_plastic)
     components = element.compute_committed_current_tangent_components(
         mesh, material, displacement, sealed, 3
     )

{anysolver-0.4.1 → anysolver-0.4.2}/tests/test_e4_pl_s3_mixed_mesh_manifest.py

@@ -36,6 +36,12 @@ Q4_REPLAY_HOTFIX = (
     / "reference_cases"
     / "e4_pl_q4_0_4_1_replay_hotfix.json"
 )
+Q4_REPLAY_CANONICALIZATION = (
+    ROOT
+    / "docs"
+    / "reference_cases"
+    / "e4_pl_q4_0_4_2_replay_canonicalization.json"
+)
 Q4_MECHANICS_AUTHORITY_SCHEMA = "anysolver.python-ast-mechanics-closure-v1"
 Q4_MECHANICS_CANONICALIZATION = (
     "AST_NODE_ORDERED_FIELDS_COMPACT_JSON_WITHOUT_LOCATIONS_DOCSTRINGS_"
@@ -479,18 +485,28 @@ def test_qualification_contract_and_replay_hotfix_bind_q4_mechanics_symbols() ->
     )
     hotfix = _strict_json(Q4_REPLAY_HOTFIX.read_bytes())
     correction = hotfix["source_correction"]
+    successor = _strict_json(Q4_REPLAY_CANONICALIZATION.read_bytes())
+    successor_correction = successor["correction"]
     changed = sorted(
         name
         for name, digest in expected.items()
         if actual["symbols"][name] != digest
     )
     added = sorted(set(actual["symbols"]) - set(expected))
-    assert changed == correction["updated_existing_symbols"]
-    assert added == correction["added_symbols"]
+    assert changed == successor_correction["cumulative_updated_existing_symbols"]
+    assert added == successor_correction["cumulative_added_symbols"]
+    assert correction["updated_existing_symbols"] == changed
+    assert set(correction["added_symbols"]) < set(added)
     assert authority["aggregate_sha256"] == correction[
         "prior_mechanics_aggregate_sha256"
     ]
-    assert _q4_mechanics_aggregate(actual) == correction[
+    assert correction["current_mechanics_aggregate_sha256"] == successor_correction[
+        "prior_mechanics_aggregate_sha256"
+    ]
+    assert hashlib.sha256(Q4_REPLAY_HOTFIX.read_bytes()).hexdigest().upper() == successor[
+        "prior_evidence"
+    ]["sha256"]
+    assert _q4_mechanics_aggregate(actual) == successor_correction[
         "current_mechanics_aggregate_sha256"
     ]
 

{anysolver-0.4.1 → anysolver-0.4.2}/tests/test_extracted_package_wiring.py

@@ -37,4 +37,4 @@ def test_file_document_and_result_api_is_canonical_anyfileio() -> None:
 
 
 def test_distribution_version_marks_the_extraction_boundary() -> None:
-    assert anysolver.__version__ == "0.4.1"
+    assert anysolver.__version__ == "0.4.2"

{anysolver-0.4.1 → anysolver-0.4.2}/tests/test_s3_v2d_default_activation.py

@@ -42,6 +42,12 @@ Q4_REPLAY_HOTFIX = (
     / "reference_cases"
     / "e4_pl_q4_0_4_1_replay_hotfix.json"
 )
+Q4_REPLAY_CANONICALIZATION = (
+    ROOT
+    / "docs"
+    / "reference_cases"
+    / "e4_pl_q4_0_4_2_replay_canonicalization.json"
+)
 
 
 def _strict_object(pairs: list[tuple[str, object]]) -> dict[str, object]:
@@ -67,6 +73,20 @@ def _q4_replay_hotfix() -> dict[str, object]:
     return payload
 
 
+def _q4_replay_canonicalization() -> dict[str, object]:
+    raw = Q4_REPLAY_CANONICALIZATION.read_bytes()
+    payload = json.loads(raw, object_pairs_hook=_strict_object)
+    assert raw == (
+        json.dumps(payload, sort_keys=True, separators=(",", ":")) + "\n"
+    ).encode("utf-8")
+    prior = payload["prior_evidence"]
+    assert prior["path"] == Q4_REPLAY_HOTFIX.relative_to(ROOT).as_posix()
+    assert hashlib.sha256(Q4_REPLAY_HOTFIX.read_bytes()).hexdigest().upper() == prior[
+        "sha256"
+    ]
+    return payload
+
+
 def test_activation_candidate_binds_accepted_v6w_and_unchanged_q4() -> None:
     raw = ACTIVATION.read_bytes()
     payload = json.loads(raw, object_pairs_hook=_strict_object)
@@ -87,10 +107,14 @@ def test_activation_candidate_binds_accepted_v6w_and_unchanged_q4() -> None:
         f"blob {len(q4_data)}\0".encode("ascii") + q4_data
     ).hexdigest()
     hotfix = _q4_replay_hotfix()
+    successor = _q4_replay_canonicalization()
     assert payload["production_boundary"]["q4_mechanics_blob"] == hotfix[
         "source_correction"
     ]["prior_git_blob"]
-    assert q4_blob == hotfix["source_correction"]["current_git_blob"]
+    assert hotfix["source_correction"]["current_git_blob"] == successor["correction"][
+        "prior_e4_pl_element_blob"
+    ]
+    assert q4_blob == successor["correction"]["e4_pl_element_blob"]
     assert hotfix["qualification_boundary"]["q4_coefficients_or_operators_changed"] is False
     assert payload["production_boundary"]["q4_mechanics_unchanged"] is True
     assert payload["production_boundary"]["publication_authorized"] is False
@@ -124,10 +148,14 @@ def test_integration_repair_preserves_qualification_and_q4_boundary() -> None:
         f"blob {len(q4_data)}\0".encode("ascii") + q4_data
     ).hexdigest()
     hotfix = _q4_replay_hotfix()
+    successor = _q4_replay_canonicalization()
     assert payload["production_boundary"]["q4_mechanics_blob"] == hotfix[
         "source_correction"
     ]["prior_git_blob"]
-    assert q4_blob == hotfix["source_correction"]["current_git_blob"]
+    assert hotfix["source_correction"]["current_git_blob"] == successor["correction"][
+        "prior_e4_pl_element_blob"
+    ]
+    assert q4_blob == successor["correction"]["e4_pl_element_blob"]
     assert hotfix["qualification_boundary"]["s3_evidence_changed"] is False
     assert payload["production_boundary"]["v6w_evidence_changed"] is False
     assert payload["production_boundary"]["s3_linear_stiffness_operator_changed"] is False