weac 3.1.2__py3-none-any.whl → 3.1.3__py3-none-any.whl

weac/__init__.py

@@ -2,4 +2,4 @@
 WEAC - Weak Layer Anticrack Nucleation Model
 """
 
-__version__ = "3.1.2"
+__version__ = "3.1.3"

weac/analysis/criteria_evaluator.py

@@ -24,7 +24,7 @@ from weac.components import (
     WeakLayer,
 )
 from weac.constants import RHO_ICE
-from weac.core.slab_touchdown import TouchdownMode
+from weac.components.scenario_config import TouchdownMode
 from weac.core.system_model import SystemModel
 
 logger = logging.getLogger(__name__)
@@ -741,6 +741,9 @@ class CriteriaEvaluator:
             cut_length=cut_distance,
         )
         system_copy.update_scenario(segments=segments, scenario_config=scenario_config)
+        # Force the requested mode to avoid floating-point precision issues
+        # when l_AB/l_BC are recalculated with different scenario parameters
+        system_copy.set_forced_touchdown_mode(mode)
         touchdown_distance = system_copy.slab_touchdown.touchdown_distance
         analyzer = Analyzer(system_copy, printing_enabled=print_call_stats)
         energy_release_rate, _, _ = analyzer.differential_ERR(unit="J/m^2")

weac/components/__init__.py

@@ -7,7 +7,7 @@ from .criteria_config import CriteriaConfig
 from .layer import Layer, WeakLayer
 from .model_input import ModelInput
 from .segment import Segment
-from .scenario_config import ScenarioConfig, SystemType
+from .scenario_config import ScenarioConfig, SystemType, TouchdownMode
 
 __all__ = [
     "Config",
@@ -18,4 +18,5 @@ __all__ = [
     "ScenarioConfig",
     "ModelInput",
     "SystemType",
+    "TouchdownMode",
 ]

weac/components/config.py

@@ -12,6 +12,8 @@ field_name: type = Field(..., gt=0, description="Description")
 
 from pydantic import BaseModel, Field
 
+from weac.components.scenario_config import TouchdownMode
+
 
 class Config(BaseModel):
     """
@@ -21,11 +23,19 @@ class Config(BaseModel):
     ----------
     touchdown : bool
         Whether slab touchdown on the collapsed weak layer is considered.
+    forced_touchdown_mode : TouchdownMode | None
+        If set, forces the touchdown mode instead of calculating it from l_AB/l_BC.
+        This avoids floating-point precision issues when the mode boundary values
+        are recalculated with different scenario parameters.
     """
 
     touchdown: bool = Field(
         default=False, description="Whether to include slab touchdown in the analysis"
     )
+    forced_touchdown_mode: TouchdownMode | None = Field(
+        default=None,
+        description="Force a specific touchdown mode instead of auto-calculating",
+    )
 
 
 if __name__ == "__main__":

weac/components/scenario_config.py

@@ -11,6 +11,8 @@ SystemType = Literal[
     "skier", "skiers", "pst-", "-pst", "rot", "trans", "vpst-", "-vpst"
 ]
 
+TouchdownMode = Literal["A_free_hanging", "B_point_contact", "C_in_contact"]
+
 
 class ScenarioConfig(BaseModel):
     """

weac/core/slab_touchdown.py

@@ -9,7 +9,7 @@ from typing import Literal
 from scipy.optimize import brentq
 
 from weac.components.layer import WeakLayer
-from weac.components.scenario_config import ScenarioConfig
+from weac.components.scenario_config import ScenarioConfig, TouchdownMode
 from weac.components.segment import Segment
 from weac.constants import STIFFNESS_COLLAPSE_FACTOR
 from weac.core.eigensystem import Eigensystem
@@ -20,9 +20,6 @@ from weac.core.unknown_constants_solver import UnknownConstantsSolver
 logger = logging.getLogger(__name__)
 
 
-TouchdownMode = Literal["A_free_hanging", "B_point_contact", "C_in_contact"]
-
-
 class SlabTouchdown:  # pylint: disable=too-many-instance-attributes,too-few-public-methods
     """
     Handling the touchdown situation in a PST.
@@ -52,6 +49,9 @@ class SlabTouchdown:  # pylint: disable=too-many-instance-attributes,too-few-pub
     -----------
     scenario: `Scenario`
     eigensystem: `Eigensystem`
+    forced_mode: `TouchdownMode | None`
+        If provided, forces this touchdown mode instead of calculating from l_AB/l_BC.
+        This avoids floating-point precision issues when scenario parameters change.
 
     Attributes:
     -----------
@@ -81,9 +81,15 @@ class SlabTouchdown:  # pylint: disable=too-many-instance-attributes,too-few-pub
     touchdown_distance: float
     collapsed_weak_layer_kR: float | None = None
 
-    def __init__(self, scenario: Scenario, eigensystem: Eigensystem):
+    def __init__(
+        self,
+        scenario: Scenario,
+        eigensystem: Eigensystem,
+        forced_mode: TouchdownMode | None = None,
+    ):
         self.scenario = scenario
         self.eigensystem = eigensystem
+        self._forced_mode = forced_mode
 
         # Create a new scenario config with phi=0 (flat slab) while preserving other settings
         self.flat_config = ScenarioConfig(
@@ -114,14 +120,25 @@ class SlabTouchdown:  # pylint: disable=too-many-instance-attributes,too-few-pub
             self.l_BC = self._calc_l_BC()
         except ValueError:
             self.l_BC = self.scenario.L
-        # Assign stage
-        touchdown_mode = "A_free_hanging"
-        if self.scenario.cut_length <= self.l_AB:
+
+        # Assign stage - use forced mode if provided, otherwise calculate from thresholds
+        if self._forced_mode is not None:
+            touchdown_mode = self._forced_mode
+            logger.debug(
+                "Using forced touchdown mode: %s (l_AB=%.2f, l_BC=%.2f, cut_length=%.2f)",
+                touchdown_mode,
+                self.l_AB,
+                self.l_BC,
+                self.scenario.cut_length,
+            )
+        else:
             touchdown_mode = "A_free_hanging"
-        elif self.l_AB < self.scenario.cut_length <= self.l_BC:
-            touchdown_mode = "B_point_contact"
-        elif self.l_BC < self.scenario.cut_length:
-            touchdown_mode = "C_in_contact"
+            if self.scenario.cut_length <= self.l_AB:
+                touchdown_mode = "A_free_hanging"
+            elif self.l_AB < self.scenario.cut_length <= self.l_BC:
+                touchdown_mode = "B_point_contact"
+            elif self.l_BC < self.scenario.cut_length:
+                touchdown_mode = "C_in_contact"
         self.touchdown_mode = touchdown_mode
 
     def _calc_touchdown_distance(self):

weac/core/system_model.py

@@ -28,6 +28,7 @@ from weac.core.eigensystem import Eigensystem
 from weac.core.field_quantities import FieldQuantities
 from weac.core.scenario import Scenario
 from weac.core.slab import Slab
+from weac.components.scenario_config import TouchdownMode
 from weac.core.slab_touchdown import SlabTouchdown
 from weac.core.unknown_constants_solver import UnknownConstantsSolver
 
@@ -164,7 +165,9 @@ class SystemModel:
         if self.config.touchdown:
             logger.info("Solving for Slab Touchdown")
             slab_touchdown = SlabTouchdown(
-                scenario=self.scenario, eigensystem=self.eigensystem
+                scenario=self.scenario,
+                eigensystem=self.eigensystem,
+                forced_mode=self.config.forced_touchdown_mode,
             )
             logger.info(
                 "Original cut_length: %s, touchdown_distance: %s",
@@ -345,6 +348,26 @@ class SystemModel:
             self._invalidate_slab_touchdown()
             self._invalidate_constants()
 
+    def set_forced_touchdown_mode(self, mode: TouchdownMode | None):
+        """
+        Set the forced touchdown mode.
+
+        When set, the touchdown mode calculation will use this mode instead of
+        calculating it from l_AB/l_BC thresholds. This avoids floating-point
+        precision issues when the mode boundaries are recalculated with different
+        scenario parameters.
+
+        Parameters
+        ----------
+        mode : TouchdownMode | None
+            The mode to force ("A_free_hanging", "B_point_contact", "C_in_contact"),
+            or None to use automatic calculation.
+        """
+        if self.config.forced_touchdown_mode != mode:
+            self.config.forced_touchdown_mode = mode
+            self._invalidate_slab_touchdown()
+            self._invalidate_constants()
+
     def _invalidate_eigensystem(self):
         """Invalidate the eigensystem."""
         self.__dict__.pop("eigensystem", None)

{weac-3.1.2.dist-info → weac-3.1.3.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: weac
-Version: 3.1.2
+Version: 3.1.3
 Summary: Weak layer anticrack nucleation model
 Author-email: 2phi GbR <mail@2phi.de>
 License-Expression: MIT

{weac-3.1.2.dist-info → weac-3.1.3.dist-info}/RECORD

@@ -1,32 +1,32 @@
-weac/__init__.py,sha256=CIiVU46-0zhtGbOm2zgB-Pcnh5wzVzn5DojC3CJQmDA,76
+weac/__init__.py,sha256=LEEvHYM1lupYQMMcn2dSLg39Hwlksbx2UAscsPiUt3M,76
 weac/constants.py,sha256=BX8ifZhFciCuzzako1p-2Wh5CWWzR3cPdvyu501UtUs,1194
 weac/logging_config.py,sha256=GJn4Fs80dDRpbPkSPByCmIOLkEI01dbeuMMSv_IDgoM,1103
 weac/analysis/__init__.py,sha256=FDbmQLC5zmT8OJftqNYgihXc8WYB1b9FasFPGzkKVU8,498
 weac/analysis/analyzer.py,sha256=-yoOqsdQvpNvBc-U2nqB4UEKHGTMBu23lVPJqKxUVcs,29095
-weac/analysis/criteria_evaluator.py,sha256=u4nvlkqTxNB0_7xp8ycRPe9RLuTLHP6_VzNlTG9q-nw,48286
+weac/analysis/criteria_evaluator.py,sha256=T4_BdWFHkw8ZpOOWdR0A4Y3GGSZHuMJdtX8bPC1mRdo,48498
 weac/analysis/plotter.py,sha256=ibkH4rzQ2pqQVXPRfek722e8eU3vP63uysb-Q-EzfC0,80556
-weac/components/__init__.py,sha256=94WIUVjPI3U-grN_7v0oi9bZB3Z8rCBcx9wJPOwJBAQ,439
-weac/components/config.py,sha256=tnOnJ0M-_knZBhdr052nDyyFFAZN0f2hQ68XRuXG6d8,869
+weac/components/__init__.py,sha256=zEQCrc9q-yd0B-n1-vLhpZLQD3ZPYhAyivHySnI1zhs,475
+weac/components/config.py,sha256=PTzLLmtWqazcRvkxQt4mCV1gAGAng3oNuQpBbEYc22g,1372
 weac/components/criteria_config.py,sha256=f2agU7nXWURFAw8_68igiSk5aICUxwaB9u3Qasit-Q0,2909
 weac/components/layer.py,sha256=Ut0B9jKW4a6GVPxV43A4VCtAQJYk6ZA3MX8gt2aJBcQ,12491
 weac/components/model_input.py,sha256=oyX4p7PgaAtUXXtpvV8XO_KxwOXtG1TB0x7rXxNfCFA,3330
-weac/components/scenario_config.py,sha256=Tam-m9DQtdjmTm-lKQ8Dcjqje04ttJS2X3v_Nn7AUHQ,2563
+weac/components/scenario_config.py,sha256=kR_YkI1fEKCs8VDQyvc118QCJnTsEFltYuDb-cpSJsc,2641
 weac/components/segment.py,sha256=F279KcAAkRuJKWav_BZ4BanO96WZm4KXtKHinFZki7s,941
 weac/core/__init__.py,sha256=pRyCKD8XD3qXVUWtFG7N3cS91P5x5d8Jpr1hMEgxQ2U,233
 weac/core/eigensystem.py,sha256=b7KXi2weCY9IVlH_7lCTXzKSx_pLdWY-x8BPjHL5nKo,13736
 weac/core/field_quantities.py,sha256=rr0KvM5Eu9AR43CQrrODJJkXkM2b2N3BWk2dE38f64c,9175
 weac/core/scenario.py,sha256=vHHe-JDWc-8ECVfy6lFGVtJmIIKRc0J9iUEGR8zckdw,6050
 weac/core/slab.py,sha256=YwbAf7ogYDn3p0nXfj3IjVqwA6ro3jBSg2wXb_io_gQ,5125
-weac/core/slab_touchdown.py,sha256=fTKrTblYC4i7CDQe_WIReUNdJXES7K5POBb5Fm2kHg4,13980
-weac/core/system_model.py,sha256=rTsQqqKMYcyJCQw3ogsUCKgXmnJlnL9l4tjV5Tfpg1Y,15475
+weac/core/slab_touchdown.py,sha256=gO-nB1RntBdDuxckqnTIoFKnATheX6JeTrzxjo-6wsc,14707
+weac/core/system_model.py,sha256=U7iAcVRZrUbRR-vNhRU1lKWZQwtucMLdrY-K0ZcFcbo,16418
 weac/core/unknown_constants_solver.py,sha256=kzB8351UC03rVcJWez0FAwdawI2tERbpcgqMRKgYbi8,17753
 weac/utils/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 weac/utils/geldsetzer.py,sha256=Rmn_PUyeomPLfxZlj9LDbdBBRDWovavrQenka318TIM,3681
 weac/utils/misc.py,sha256=lGz0IDDJ_3nvYjSkivPJ5Xscl1D_AmvQLSjaL7SUbKs,3674
 weac/utils/snow_types.py,sha256=OKSY8LzfKycSjg7uSUZkSXUzwhcw1D1fZs5S3hVQtUM,1573
 weac/utils/snowpilot_parser.py,sha256=drOLnUg7uzMJMO9qU9x79xZkqQfFPCCQonVpU9S0x8U,13682
-weac-3.1.2.dist-info/licenses/LICENSE,sha256=ojZPWKFHbFGDrlNOvuAKGH9WcKhpLHWZPcQ4SzhK91M,1082
-weac-3.1.2.dist-info/METADATA,sha256=xkJGkSjRAFW5BmqW99JDo9R1G71zq9rlh2gEOBl_P74,25593
-weac-3.1.2.dist-info/WHEEL,sha256=wUyA8OaulRlbfwMtmQsvNngGrxQHAvkKcvRmdizlJi0,92
-weac-3.1.2.dist-info/top_level.txt,sha256=8tyXUHPFU4Ba_5kPtpwvXo5l6GjJmOnODVBJFygpdeE,5
-weac-3.1.2.dist-info/RECORD,,
+weac-3.1.3.dist-info/licenses/LICENSE,sha256=ojZPWKFHbFGDrlNOvuAKGH9WcKhpLHWZPcQ4SzhK91M,1082
+weac-3.1.3.dist-info/METADATA,sha256=D1EzMmxGIS1qFrcq9QnI8P4vKy7Y77GAiTL_LenYvms,25593
+weac-3.1.3.dist-info/WHEEL,sha256=wUyA8OaulRlbfwMtmQsvNngGrxQHAvkKcvRmdizlJi0,92
+weac-3.1.3.dist-info/top_level.txt,sha256=8tyXUHPFU4Ba_5kPtpwvXo5l6GjJmOnODVBJFygpdeE,5
+weac-3.1.3.dist-info/RECORD,,