crc-framework 0.2.0__tar.gz → 0.2.4__tar.gz

{crc_framework-0.2.0 → crc_framework-0.2.4}/.github/workflows/release.yaml

@@ -69,7 +69,6 @@ jobs:
             3.11
             3.12
             3.13
-            3.13t
             3.14
             3.14t
       - name: Build wheels
@@ -79,12 +78,11 @@ jobs:
         with:
           target: ${{ matrix.target }}
           args: >-
-            --release 
-            --out dist 
+            --release
+            --out dist
             --find-interpreter
           sccache: "true"
           container: "off"
-          rust-toolchain: nightly-2026-01-27
       - name: Upload wheels
         uses: actions/upload-artifact@v4
         with:
@@ -102,7 +100,6 @@ jobs:
         with:
           command: sdist
           args: --out dist
-          rust-toolchain: nightly-2026-01-27
       - name: Upload sdist
         uses: actions/upload-artifact@v4
         with:

{crc_framework-0.2.0 → crc_framework-0.2.4}/Cargo.lock

@@ -15,12 +15,6 @@ dependencies = [
  "zerocopy",
 ]
 
-[[package]]
-name = "autocfg"
-version = "1.5.1"
-source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "f2032f911046de80f0a198e0901378627c33f59ea0ac00e363d481118bd70a53"
-
 [[package]]
 name = "cfg-if"
 version = "1.0.4"
@@ -49,7 +43,7 @@ dependencies = [
 
 [[package]]
 name = "crc-framework-core"
-version = "0.2.0"
+version = "0.2.4"
 dependencies = [
  "h3o",
  "phf",
@@ -58,7 +52,7 @@ dependencies = [
 
 [[package]]
 name = "crc-framework-rs"
-version = "0.2.0"
+version = "0.2.4"
 dependencies = [
  "crc-framework-core",
  "h3o",
@@ -137,15 +131,6 @@ version = "0.5.0"
 source = "registry+https://github.com/rust-lang/crates.io-index"
 checksum = "2304e00983f87ffb38b55b444b5e3b60a884b5d30c0fca7d82fe33449bbe55ea"
 
-[[package]]
-name = "indoc"
-version = "2.0.7"
-source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "79cf5c93f93228cf8efb3ba362535fb11199ac548a09ce117c9b1adc3030d706"
-dependencies = [
- "rustversion",
-]
-
 [[package]]
 name = "libc"
 version = "0.2.188"
@@ -158,15 +143,6 @@ version = "0.2.16"
 source = "registry+https://github.com/rust-lang/crates.io-index"
 checksum = "b6d2cec3eae94f9f509c767b45932f1ada8350c4bdb85af2fcab4a3c14807981"
 
-[[package]]
-name = "memoffset"
-version = "0.9.1"
-source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "488016bfae457b036d996092f6cb448677611ce4449e970ceaf42695203f218a"
-dependencies = [
- "autocfg",
-]
-
 [[package]]
 name = "once_cell"
 version = "1.21.4"
@@ -242,35 +218,32 @@ dependencies = [
 
 [[package]]
 name = "pyo3"
-version = "0.27.2"
+version = "0.29.0"
 source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "ab53c047fcd1a1d2a8820fe84f05d6be69e9526be40cb03b73f86b6b03e6d87d"
+checksum = "cd274650b21d4bfc26a0a47587962c1edb425f69287324355cd040c3ea66071c"
 dependencies = [
- "indoc",
  "libc",
- "memoffset",
  "once_cell",
  "portable-atomic",
  "pyo3-build-config",
  "pyo3-ffi",
  "pyo3-macros",
- "unindent",
 ]
 
 [[package]]
 name = "pyo3-build-config"
-version = "0.27.2"
+version = "0.29.0"
 source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "b455933107de8642b4487ed26d912c2d899dec6114884214a0b3bb3be9261ea6"
+checksum = "c5e2a7d2f0d013342f295c048ad19237add5154a55b1c5a254c0ec93d4109078"
 dependencies = [
  "target-lexicon",
 ]
 
 [[package]]
 name = "pyo3-ffi"
-version = "0.27.2"
+version = "0.29.0"
 source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "1c85c9cbfaddf651b1221594209aed57e9e5cff63c4d11d1feead529b872a089"
+checksum = "ca85c467da1bbc8d866eea5deff9cf29ea5f7785054a17da36e65bda9c05845b"
 dependencies = [
  "libc",
  "pyo3-build-config",
@@ -278,9 +251,9 @@ dependencies = [
 
 [[package]]
 name = "pyo3-macros"
-version = "0.27.2"
+version = "0.29.0"
 source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "0a5b10c9bf9888125d917fb4d2ca2d25c8df94c7ab5a52e13313a07e050a3b02"
+checksum = "9ac53762fd065daa3194dd09337a38bd793a188100fd1a9304c4ab312d901771"
 dependencies = [
  "proc-macro2",
  "pyo3-macros-backend",
@@ -290,13 +263,12 @@ dependencies = [
 
 [[package]]
 name = "pyo3-macros-backend"
-version = "0.27.2"
+version = "0.29.0"
 source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "03b51720d314836e53327f5871d4c0cfb4fb37cc2c4a11cc71907a86342c40f9"
+checksum = "4ca3a1557399783172dc5bf39cfca835157732532cba56b71d2292161e53b362"
 dependencies = [
  "heck",
  "proc-macro2",
- "pyo3-build-config",
  "quote",
  "syn 2.0.119",
 ]
@@ -345,12 +317,6 @@ dependencies = [
  "getrandom 0.3.4",
 ]
 
-[[package]]
-name = "rustversion"
-version = "1.0.23"
-source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "cf54715a573b99ac80df0bc206da022bcd442c974952c7b9720069370852e21f"
-
 [[package]]
 name = "serde"
 version = "1.0.229"
@@ -429,12 +395,6 @@ version = "1.0.24"
 source = "registry+https://github.com/rust-lang/crates.io-index"
 checksum = "e6e4313cd5fcd3dad5cafa179702e2b244f760991f45397d14d4ebf38247da75"
 
-[[package]]
-name = "unindent"
-version = "0.2.4"
-source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "7264e107f553ccae879d21fbea1d6724ac785e8c3bfc762137959b5802826ef3"
-
 [[package]]
 name = "version_check"
 version = "0.9.5"

{crc_framework-0.2.0 → crc_framework-0.2.4}/Cargo.toml

@@ -1,6 +1,6 @@
 [package]
 name = "crc-framework-rs"
-version = "0.2.0"
+version = "0.2.4"
 edition = "2024"
 rust-version = "1.85"
 description = "Rust-powered Python bindings for the Climate Risk Commons Framework"
@@ -18,4 +18,4 @@ crate-type = ["cdylib"]
 [dependencies]
 crc-framework-core = { path = "crates/core" }
 h3o = "=0.9.3"
-pyo3 = { version = "0.27", features = ["extension-module"] }
+pyo3 = { version = "0.29", features = ["extension-module"] }

{crc_framework-0.2.0 → crc_framework-0.2.4}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: crc-framework
-Version: 0.2.0
+Version: 0.2.4
 Classifier: Development Status :: 3 - Alpha
 Classifier: License :: OSI Approved :: GNU Affero General Public License v3 or later (AGPLv3+)
 Classifier: Operating System :: OS Independent
@@ -20,6 +20,9 @@ Author: Riskthinking.AI
 License-Expression: AGPL-3.0-or-later
 Requires-Python: >=3.9
 Description-Content-Type: text/markdown; charset=UTF-8; variant=GFM
+Project-URL: Homepage, https://github.com/RiskThinking/crc-framework-rs
+Project-URL: Issues, https://github.com/RiskThinking/crc-framework-rs/issues
+Project-URL: Repository, https://github.com/RiskThinking/crc-framework-rs
 
 # Climate Risk Commons Framework
 
@@ -92,6 +95,7 @@ Choose a distribution according to the form of the available data:
 - `FittedDistribution` for a known parametric family and its parameters.
 - `HurdleDistribution` for a point mass followed by a truncated parametric
   tail.
+- `PointMassDistribution` for a value that is constant at every probability.
 
 `TabulatedDistribution` does not extrapolate by default. Querying outside its
 declared probability range raises an error unless `extrapolate=True` is set.
@@ -307,8 +311,9 @@ continent and intersecting ISO3 country codes.
 The package exports the following primary interfaces:
 
 - Distributions: `EmpiricalDistribution`, `TabulatedDistribution`,
-  `FittedDistribution`, `HurdleDistribution`, `fit_distribution`,
-  `fit_quantiles`, `fit_hurdle_quantiles`, `fit_all`, and `quality_metrics`.
+  `FittedDistribution`, `HurdleDistribution`, `PointMassDistribution`,
+  `fit_distribution`, `fit_quantiles`, `fit_hurdle_quantiles`, `fit_all`, and
+  `quality_metrics`.
 - Transforms: `impacts`, `ImpactRegistry`, `ImpactFunction`, `ClimateImpact`,
   `LinearImpact`, `SigmoidImpact`, `PiecewiseLinearImpact`, `CallableImpact`,
   and `CallableTransform`.

{crc_framework-0.2.0 → crc_framework-0.2.4}/README.md

@@ -69,6 +69,7 @@ Choose a distribution according to the form of the available data:
 - `FittedDistribution` for a known parametric family and its parameters.
 - `HurdleDistribution` for a point mass followed by a truncated parametric
   tail.
+- `PointMassDistribution` for a value that is constant at every probability.
 
 `TabulatedDistribution` does not extrapolate by default. Querying outside its
 declared probability range raises an error unless `extrapolate=True` is set.
@@ -284,8 +285,9 @@ continent and intersecting ISO3 country codes.
 The package exports the following primary interfaces:
 
 - Distributions: `EmpiricalDistribution`, `TabulatedDistribution`,
-  `FittedDistribution`, `HurdleDistribution`, `fit_distribution`,
-  `fit_quantiles`, `fit_hurdle_quantiles`, `fit_all`, and `quality_metrics`.
+  `FittedDistribution`, `HurdleDistribution`, `PointMassDistribution`,
+  `fit_distribution`, `fit_quantiles`, `fit_hurdle_quantiles`, `fit_all`, and
+  `quality_metrics`.
 - Transforms: `impacts`, `ImpactRegistry`, `ImpactFunction`, `ClimateImpact`,
   `LinearImpact`, `SigmoidImpact`, `PiecewiseLinearImpact`, `CallableImpact`,
   and `CallableTransform`.

{crc_framework-0.2.0 → crc_framework-0.2.4}/crates/core/Cargo.toml

@@ -1,6 +1,6 @@
 [package]
 name = "crc-framework-core"
-version = "0.2.0"
+version = "0.2.4"
 edition = "2024"
 rust-version = "1.85"
 description = "Core distribution and climate-risk algorithms for CRC Framework"

{crc_framework-0.2.0 → crc_framework-0.2.4}/crates/core/src/distribution.rs

@@ -256,6 +256,47 @@ impl Distribution for EmpiricalDistribution {
     }
 }
 
+/// A degenerate distribution whose entire probability mass is at one value.
+#[derive(Debug, Clone, Copy)]
+pub struct PointMassDistribution {
+    location: f64,
+}
+
+impl PointMassDistribution {
+    pub fn new(location: f64) -> Result<Self> {
+        if !location.is_finite() {
+            return Err(CrcError::InvalidInput(
+                "point-mass location must be finite".into(),
+            ));
+        }
+        Ok(Self { location })
+    }
+
+    pub fn location(&self) -> f64 {
+        self.location
+    }
+}
+
+impl Distribution for PointMassDistribution {
+    fn pdf(&self, _x: f64) -> f64 {
+        // A Dirac mass has no finite density with respect to Lebesgue measure.
+        0.0
+    }
+
+    fn cdf(&self, x: f64) -> f64 {
+        if x < self.location { 0.0 } else { 1.0 }
+    }
+
+    fn ppf(&self, probability: f64) -> Result<f64> {
+        validate_probability(probability)?;
+        Ok(self.location)
+    }
+
+    fn sample(&self, size: usize, _seed: Option<u64>) -> Result<Vec<f64>> {
+        Ok(vec![self.location; size])
+    }
+}
+
 #[derive(Debug, Clone, Copy, PartialEq, Eq)]
 pub enum DistributionFamily {
     GenExtreme,
@@ -1129,6 +1170,18 @@ fn erf(x: f64) -> f64 {
 mod tests {
     use super::*;
 
+    #[test]
+    fn point_mass_is_exact_at_every_probability() {
+        let distribution = PointMassDistribution::new(2.5).unwrap();
+        assert_eq!(
+            distribution.quantiles(&[0.0, 0.5, 1.0]).unwrap(),
+            vec![2.5; 3]
+        );
+        assert_eq!(distribution.cdf(2.49), 0.0);
+        assert_eq!(distribution.cdf(2.5), 1.0);
+        assert_eq!(distribution.sample(3, Some(42)).unwrap(), vec![2.5; 3]);
+    }
+
     #[test]
     fn return_periods_convert_to_upper_tail_probabilities() {
         let distribution = TabulatedDistribution::from_return_periods(

{crc_framework-0.2.0 → crc_framework-0.2.4}/crates/core/src/lib.rs

@@ -10,8 +10,8 @@ pub mod transform;
 pub use distribution::{
     DiagnosticMetrics, Distribution, DistributionFamily, EmpiricalDistribution, FitResult,
     FittedDistribution, HurdleDistribution, HurdleQuantileFitDiagnostics, HurdleQuantileFitResult,
-    Interpolation, QuantileFitDiagnostics, QuantileFitResult, TabulatedDistribution, Tail,
-    fit_hurdle_quantiles, fit_quantiles,
+    Interpolation, PointMassDistribution, QuantileFitDiagnostics, QuantileFitResult,
+    TabulatedDistribution, Tail, fit_hurdle_quantiles, fit_quantiles,
 };
 pub use error::{CrcError, Result};
 pub use metrics::{

{crc_framework-0.2.0 → crc_framework-0.2.4}/crates/core/src/reference/curve_fitting.rs

@@ -469,6 +469,17 @@ where
     let mut values: Vec<f64> = simplex.iter().map(|point| objective(point)).collect();
     let mut evaluations = values.len();
 
+    // Multi-start callers (e.g. GenExtreme's 5 shape starts) run this to
+    // completion for every start even when a start is clearly not
+    // converging. Nelder-Mead's improvement curve is front-loaded, so a
+    // start that hasn't beaten its best objective in STALL_LIMIT iterations
+    // is treated the same as one that exhausted max_iterations: it returns
+    // early as not-converged, which the caller already ranks below any
+    // converged candidate.
+    const STALL_LIMIT: usize = 300;
+    let mut best_seen = f64::INFINITY;
+    let mut stall_count = 0usize;
+
     for iteration in 0..max_iterations {
         let mut order: Vec<usize> = (0..simplex.len()).collect();
         order.sort_by(|&left, &right| values[left].total_cmp(&values[right]));
@@ -498,6 +509,22 @@ where
             };
         }
 
+        if values[best] < best_seen - objective_tolerance {
+            best_seen = values[best];
+            stall_count = 0;
+        } else {
+            stall_count += 1;
+            if stall_count >= STALL_LIMIT {
+                return SimplexResult {
+                    parameters: simplex[best].clone(),
+                    objective: values[best],
+                    converged: false,
+                    iterations: iteration,
+                    evaluations,
+                };
+            }
+        }
+
         let mut centroid = vec![0.0; dimensions];
         for &index in &order[..dimensions] {
             for (dimension, value) in centroid.iter_mut().enumerate() {

{crc_framework-0.2.0 → crc_framework-0.2.4}/pyproject.toml

@@ -28,6 +28,11 @@ classifiers = [
 ]
 dynamic = ["version"]
 
+[project.urls]
+Homepage = "https://github.com/RiskThinking/crc-framework-rs"
+Repository = "https://github.com/RiskThinking/crc-framework-rs"
+Issues = "https://github.com/RiskThinking/crc-framework-rs/issues"
+
 [project.optional-dependencies]
 test = ["mypy", "pytest"]
 

{crc_framework-0.2.0 → crc_framework-0.2.4}/python/crc_framework/__init__.py

@@ -11,6 +11,7 @@ from .distributions import (
     HurdleDistribution,
     HurdleQuantileFitDiagnostics,
     HurdleQuantileFitResult,
+    PointMassDistribution,
     QuantileFitDiagnostics,
     QuantileFitResult,
     TabulatedDistribution,
@@ -44,34 +45,35 @@ from .transforms import (
 )
 
 __all__ = [
+    "HORIZONS",
+    "PATHWAYS",
+    "RISK_FACTORS",
     "BinaryOutcome",
     "CallableImpact",
     "CallableTransform",
     "ClimateImpact",
     "Distribution",
     "EmpiricalDistribution",
-    "FitDiagnostics",
     "FitConstraints",
+    "FitDiagnostics",
     "FitResult",
     "FittedDistribution",
+    "Geography",
     "HurdleDistribution",
     "HurdleQuantileFitDiagnostics",
     "HurdleQuantileFitResult",
-    "Geography",
-    "HORIZONS",
     "ImpactFunction",
     "ImpactRegistry",
     "LinearImpact",
     "Microscore",
     "MicroscoreSuite",
-    "PATHWAYS",
     "Pathway",
     "PiecewiseLinearImpact",
-    "RISK_FACTORS",
-    "RiskFactor",
-    "RiskResult",
+    "PointMassDistribution",
     "QuantileFitDiagnostics",
     "QuantileFitResult",
+    "RiskFactor",
+    "RiskResult",
     "ScenarioMetadata",
     "SigmoidImpact",
     "TabulatedDistribution",

{crc_framework-0.2.0 → crc_framework-0.2.4}/python/crc_framework/_core.pyi

@@ -100,6 +100,7 @@ class NativeRiskResult:
     branch_count: int
 
 def empirical_distribution(samples: Sequence[float]) -> NativeDistribution: ...
+def point_mass_distribution(location: float) -> NativeDistribution: ...
 def tabulated_distribution(
     probabilities: Sequence[float],
     values: Sequence[float],

{crc_framework-0.2.0 → crc_framework-0.2.4}/python/crc_framework/distributions/__init__.py

@@ -16,24 +16,26 @@ from .fitting import (
     quality_metrics,
 )
 from .hurdle import HurdleDistribution
+from .point_mass import PointMassDistribution
 from .tabulated import Interpolation, TabulatedDistribution, Tail
 
 __all__ = [
     "Distribution",
     "DistributionFamily",
     "EmpiricalDistribution",
-    "FitDiagnostics",
     "FitConstraints",
+    "FitDiagnostics",
     "FitResult",
     "FittedDistribution",
     "HurdleDistribution",
     "HurdleQuantileFitDiagnostics",
     "HurdleQuantileFitResult",
     "Interpolation",
-    "TabulatedDistribution",
-    "Tail",
+    "PointMassDistribution",
     "QuantileFitDiagnostics",
     "QuantileFitResult",
+    "TabulatedDistribution",
+    "Tail",
     "fit_all",
     "fit_distribution",
     "fit_hurdle_quantiles",

crc_framework-0.2.4/python/crc_framework/distributions/point_mass.py

@@ -0,0 +1,55 @@
+from __future__ import annotations
+
+import numpy as np
+import numpy.typing as npt
+
+from crc_framework import _core
+
+from .base import ArrayLike
+from .empirical import _evaluate
+
+
+class PointMassDistribution:
+    """A distribution concentrated exactly at one finite value."""
+
+    def __init__(
+        self,
+        location: float,
+        *,
+        _native: _core.NativeDistribution | None = None,
+    ) -> None:
+        self.location = float(location)
+        self._native = _native or _core.point_mass_distribution(self.location)
+
+    def point_mass(self, x: ArrayLike) -> float | npt.NDArray[np.float64]:
+        values = np.asarray(x, dtype=np.float64)
+        masses = np.where(values == self.location, 1.0, 0.0)
+        if values.ndim == 0:
+            return float(masses)
+        return np.asarray(masses, dtype=np.float64)
+
+    def pdf(self, x: ArrayLike) -> float | npt.NDArray[np.float64]:
+        return _evaluate(self._native, "pdf", x)
+
+    def cdf(self, x: ArrayLike) -> float | npt.NDArray[np.float64]:
+        return _evaluate(self._native, "cdf", x)
+
+    def ppf(self, probability: ArrayLike) -> float | npt.NDArray[np.float64]:
+        return _evaluate(self._native, "ppf", probability)
+
+    def quantiles(self, probabilities: ArrayLike) -> npt.NDArray[np.float64]:
+        values = np.asarray(probabilities, dtype=np.float64)
+        return np.asarray(
+            self._native.quantiles(values.reshape(-1).tolist()), dtype=np.float64
+        ).reshape(values.shape)
+
+    def sample(
+        self,
+        size: int,
+        *,
+        rng: np.random.Generator | None = None,
+        seed: int | None = None,
+    ) -> npt.NDArray[np.float64]:
+        if rng is not None and seed is not None:
+            raise ValueError("provide either rng or seed, not both")
+        return np.full(size, self.location, dtype=np.float64)

{crc_framework-0.2.0 → crc_framework-0.2.4}/src/lib.rs

@@ -2,8 +2,8 @@ use std::{collections::HashMap, sync::Arc};
 
 use crc_framework_core::{
     BinaryOutcome, Distribution, DistributionFamily, EmpiricalDistribution, FittedDistribution,
-    HurdleDistribution, ImpactContext, ImpactRegistry, Interpolation, ScenarioMetadata,
-    TabulatedDistribution, Tail, Transform, compute_risk, compute_spanning_set,
+    HurdleDistribution, ImpactContext, ImpactRegistry, Interpolation, PointMassDistribution,
+    ScenarioMetadata, TabulatedDistribution, Tail, Transform, compute_risk, compute_spanning_set,
     fit_hurdle_quantiles as core_fit_hurdle_quantiles, fit_quantiles as core_fit_quantiles,
     generate_microscores,
 };
@@ -246,6 +246,18 @@ fn empirical_distribution(samples: Vec<f64>) -> PyResult<PyDistribution> {
     })
 }
 
+#[pyfunction]
+fn point_mass_distribution(location: f64) -> PyResult<PyDistribution> {
+    let distribution = PointMassDistribution::new(location).map_err(py_error)?;
+    Ok(PyDistribution {
+        inner: Arc::new(distribution),
+        family: None,
+        shape: None,
+        location: Some(location),
+        scale: None,
+    })
+}
+
 fn interpolation(value: &str) -> PyResult<Interpolation> {
     match value {
         "linear_probability" => Ok(Interpolation::LinearProbability),
@@ -365,6 +377,7 @@ fn fit_candidates(samples: Vec<f64>) -> PyResult<Vec<PyFitResult>> {
 #[pyfunction]
 #[pyo3(signature = (probabilities, values, family, weights=None))]
 fn fit_quantiles(
+    py: Python<'_>,
     probabilities: Vec<f64>,
     values: Vec<f64>,
     family: &str,
@@ -372,7 +385,7 @@ fn fit_quantiles(
 ) -> PyResult<PyQuantileFitResult> {
     let family = DistributionFamily::from_name(family)
         .ok_or_else(|| PyValueError::new_err(format!("unknown distribution family {family}")))?;
-    core_fit_quantiles(&probabilities, &values, weights.as_deref(), family)
+    py.detach(move || core_fit_quantiles(&probabilities, &values, weights.as_deref(), family))
         .map(quantile_fit_result_py)
         .map_err(py_error)
 }
@@ -387,6 +400,7 @@ fn fit_quantiles(
     weights=None
 ))]
 fn fit_hurdle_quantiles(
+    py: Python<'_>,
     probabilities: Vec<f64>,
     values: Vec<f64>,
     family: &str,
@@ -396,14 +410,16 @@ fn fit_hurdle_quantiles(
 ) -> PyResult<PyHurdleQuantileFitResult> {
     let family = DistributionFamily::from_name(family)
         .ok_or_else(|| PyValueError::new_err(format!("unknown distribution family {family}")))?;
-    core_fit_hurdle_quantiles(
-        &probabilities,
-        &values,
-        weights.as_deref(),
-        family,
-        atom_location,
-        atom_probability,
-    )
+    py.detach(move || {
+        core_fit_hurdle_quantiles(
+            &probabilities,
+            &values,
+            weights.as_deref(),
+            family,
+            atom_location,
+            atom_probability,
+        )
+    })
     .map(hurdle_quantile_fit_result_py)
     .map_err(py_error)
 }
@@ -869,6 +885,7 @@ fn _core(m: &Bound<'_, PyModule>) -> PyResult<()> {
     m.add_class::<PyRiskLevel>()?;
     m.add_class::<PyRiskResult>()?;
     m.add_function(wrap_pyfunction!(empirical_distribution, m)?)?;
+    m.add_function(wrap_pyfunction!(point_mass_distribution, m)?)?;
     m.add_function(wrap_pyfunction!(tabulated_distribution, m)?)?;
     m.add_function(wrap_pyfunction!(return_period_distribution, m)?)?;
     m.add_function(wrap_pyfunction!(fitted_distribution, m)?)?;

{crc_framework-0.2.0 → crc_framework-0.2.4}/tests/test_sdk.py

@@ -1,13 +1,13 @@
 import unittest
 
 import numpy as np
-
 from crc_framework import (
     BinaryOutcome,
     EmpiricalDistribution,
     FitConstraints,
     FittedDistribution,
     HurdleDistribution,
+    PointMassDistribution,
     RiskFactor,
     ScenarioMetadata,
     TabulatedDistribution,
@@ -27,6 +27,18 @@ from crc_framework import (
 
 
 class DistributionTests(unittest.TestCase):
+    def test_point_mass_is_exact_at_every_probability(self):
+        distribution = PointMassDistribution(0.0)
+
+        np.testing.assert_array_equal(
+            distribution.quantiles([0.0, 0.5, 0.999, 1.0]),
+            np.zeros(4),
+        )
+        self.assertEqual(distribution.cdf(-0.1), 0.0)
+        self.assertEqual(distribution.cdf(0.0), 1.0)
+        self.assertEqual(distribution.point_mass(0.0), 1.0)
+        np.testing.assert_array_equal(distribution.sample(3), np.zeros(3))
+
     def test_return_period_input_uses_canonical_probability(self):
         distribution = TabulatedDistribution.from_return_periods(
             [10, 100, 1000],
@@ -44,7 +56,9 @@ class DistributionTests(unittest.TestCase):
         self.assertEqual(result.distribution.quantiles([0.1, 0.5, 0.9]).shape, (3,))
         self.assertEqual(result.distribution.sample(4, seed=42).shape, (4,))
         diagnostics = quality_metrics(source, result.distribution)
-        self.assertAlmostEqual(diagnostics.ks_statistic, result.diagnostics.ks_statistic)
+        self.assertAlmostEqual(
+            diagnostics.ks_statistic, result.diagnostics.ks_statistic
+        )
         constrained = fit_distribution(
             source,
             constraints=FitConstraints(probability=0.99, maximum_value=100.0),
@@ -72,16 +86,12 @@ class DistributionTests(unittest.TestCase):
         )
 
     def test_sample_fit_rejects_probability_knots(self):
-        tabulated = TabulatedDistribution(
-            [0.1, 0.3, 0.6, 0.9], [0.0, 1.0, 2.0, 3.0]
-        )
+        tabulated = TabulatedDistribution([0.1, 0.3, 0.6, 0.9], [0.0, 1.0, 2.0, 3.0])
         with self.assertRaisesRegex(ValueError, "fit_quantiles"):
             fit_distribution(tabulated)  # type: ignore[arg-type]
 
     def test_hurdle_constructor_preserves_atom_and_truncated_tail(self):
-        base = FittedDistribution.from_parameters(
-            "gumbel_r", location=0.5, scale=1.2
-        )
+        base = FittedDistribution.from_parameters("gumbel_r", location=0.5, scale=1.2)
         hurdle = HurdleDistribution(base, atom_probability=0.4)
         self.assertEqual(hurdle.ppf(0.4), 0.0)
         self.assertEqual(hurdle.point_mass(0.0), 0.4)
@@ -158,9 +168,7 @@ class PipelineTests(unittest.TestCase):
 class SpatialAndConstantTests(unittest.TestCase):
     def test_reference_spatial_tables_are_exposed_without_silent_defaults(self):
         self.assertEqual(lookup_continent(594475228122316799), "Europe")
-        self.assertEqual(
-            lookup_ipcc_region(594479566039285759), "Greenland/Iceland"
-        )
+        self.assertEqual(lookup_ipcc_region(594479566039285759), "Greenland/Iceland")
         geography = lookup_geography(600550049193132031)
         self.assertIsNotNone(geography)
         self.assertEqual(geography.continent, "Africa")