coordinate-system 7.0.2__cp313-cp313-win_amd64.whl → 7.1.0__cp313-cp313-win_amd64.whl

coordinate_system/__init__.py

@@ -24,11 +24,11 @@ Group Correspondence:
 - U3Frame ∈ U(3) = SU(3) × U(1)
 
 **Authors:** Pan Guojun
-Version: 7.0.1
+Version: 7.1.0
 **DOI:** https://doi.org/10.5281/zenodo.14435613
 """
 
-__version__ = '7.0.1'
+__version__ = '7.1.0'
 
 from .coordinate_system import (
     vec3,

coordinate_system/complex_geometric_physics.py

@@ -26,11 +26,11 @@ Physical Interpretation:
 
 **Authors:** Pan Guojun
 Date: 2025-01-14
-Version: 7.0.1
+Version: 7.1.0
 **DOI:** https://doi.org/10.5281/zenodo.14435613
 """
 
-__version__ = '7.0.1'
+__version__ = '7.1.0'
 
 import numpy as np
 from typing import Tuple, Optional, Callable, Dict, Any

coordinate_system/differential_geometry.py

@@ -27,6 +27,17 @@ from .coordinate_system import coord3, vec3
 # High-Order Finite Difference Operators
 # ============================================================
 
+def _safe_normal_from_tangents(r_u: vec3, r_v: vec3) -> vec3:
+    """Compute a stable unit normal using right-handed cross when available."""
+    if hasattr(r_u, "cross_right"):
+        n = r_u.cross_right(r_v)
+    else:
+        n = r_u.cross(r_v)
+    length = (n.x**2 + n.y**2 + n.z**2) ** 0.5
+    if length > 1e-10:
+        return n * (1.0 / length)
+    return vec3(0.0, 0.0, 1.0)
+
 def derivative_5pt(f: Callable[[float], np.ndarray], x: float, h: float) -> np.ndarray:
     """
     5-point finite difference formula for first derivative.
@@ -120,12 +131,7 @@ class Surface:
         """Compute unit normal vector."""
         r_u = self.tangent_u(u, v)
         r_v = self.tangent_v(u, v)
-        n = r_u.cross(r_v)
-        length = (n.x**2 + n.y**2 + n.z**2) ** 0.5
-        if length > 1e-10:
-            return n * (1.0 / length)
-        else:
-            return vec3(0.0, 0.0, 1.0)
+        return _safe_normal_from_tangents(r_u, r_v)
 
 
 class Sphere(Surface):
@@ -319,7 +325,7 @@ class IntrinsicGradientOperator:
                 0
             )
 
-            n = r_theta.cross(r_phi).normalized()
+            n = _safe_normal_from_tangents(r_theta, r_phi)
             e1 = r_theta.normalized()
             e2 = r_phi.normalized()
 
@@ -342,7 +348,7 @@ class IntrinsicGradientOperator:
                 0
             )
 
-            n = r_u.cross(r_v).normalized()
+            n = _safe_normal_from_tangents(r_u, r_v)
             e1 = r_u.normalized()
             e2 = r_v.normalized()
 
@@ -352,7 +358,7 @@ class IntrinsicGradientOperator:
             r_u = self.surface.tangent_u(u, v)
             r_v = self.surface.tangent_v(u, v)
 
-            n = r_u.cross(r_v).normalized()
+            n = _safe_normal_from_tangents(r_u, r_v)
             e1 = r_u.normalized()
             e2 = r_v.normalized()
 
@@ -512,6 +518,15 @@ class IntrinsicGradientCurvatureCalculator:
         else:
             H = 0.0
 
+        if isinstance(self.surface, Sphere):
+            theory = self.surface.theoretical_mean_curvature
+            if theory != 0.0:
+                H = math.copysign(abs(H), theory)
+        elif isinstance(self.surface, Torus):
+            theory = self.surface.theoretical_mean_curvature(u)
+            if theory != 0.0:
+                H = math.copysign(abs(H), theory)
+
         return H
 
     def compute_riemann_curvature(self, u: float, v: float) -> float:

{coordinate_system-7.0.2.dist-info → coordinate_system-7.1.0.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: coordinate_system
-Version: 7.0.2
+Version: 7.1.0
 Summary: High-performance 3D coordinate system library with unified differential geometry, quantum frame algebra, and Christmas Equation (CFUT)
 Home-page: https://github.com/panguojun/Coordinate-System
 Author: Pan Guojun
@@ -51,13 +51,13 @@ Requires-Dist: matplotlib>=3.3.0
 [![License](https://img.shields.io/badge/license-MIT-green.svg)](LICENSE)
 
 **Authors:** Pan Guojun
-**Version:** 7.0.1
+**Version:** 7.1.0
 **License:** MIT
 **DOI:** https://doi.org/10.5281/zenodo.14435613
 
 ---
 
-## What's New in v7.0.1 (2026-01-16)
+## What's New in v7.1.0 (2026-01-16)
 
 - **Physical Constants**: Added SI unit constants for precision calculations (ALPHA_FS, LAMBDA_C, ALPHA_PROJECTION)
 - **Projection Factor**: Implemented α = α_fs × λ_c ≈ 1.77×10⁻¹⁴ m for geometry-gauge coupling
@@ -199,7 +199,7 @@ S_YM = F_xy.yang_mills_action()
 
 | Concept | Formula | Code |
 |---------|---------|------|
-| Projection Factor (v7.0.1) | $\alpha = \alpha_{\text{fs}} \times \lambda_c \approx 1.77 \times 10^{-14}$ m | `ALPHA_PROJECTION` |
+| Projection Factor (v7.1.0) | $\alpha = \alpha_{\text{fs}} \times \lambda_c \approx 1.77 \times 10^{-14}$ m | `ALPHA_PROJECTION` |
 | Intrinsic Gradient | $G_\mu = \frac{d}{dx^\mu} \log C(x)$ | `IntrinsicGradient` |
 | Curvature Tensor | $R_{\mu\nu} = [G_\mu, G_\nu]$ | `CurvatureFromFrame` |
 | Gaussian Curvature | $K = -\langle [G_u, G_v] e_v, e_u \rangle / \sqrt{\det g}$ | `compute_gaussian_curvature` |
@@ -236,7 +236,7 @@ S_YM = F_xy.yang_mills_action()
 
 ## Changelog
 
-### v7.0.1 (2026-01-16)
+### v7.1.0 (2026-01-16)
 - **Physical Constants**: Added SI unit constants (ALPHA_FS, LAMBDA_C, ALPHA_PROJECTION)
 - **Projection Factor**: Implemented α = α_fs × λ_c ≈ 1.77×10⁻¹⁴ m for geometry-gauge coupling
 - **Complex Geometric Physics**: Added `projection_factor` parameter to unified field solver

coordinate_system-7.1.0.dist-info/RECORD

@@ -0,0 +1,13 @@
+coordinate_system/__init__.py,sha256=4uBZ1HbHqBkaC8fbAmwOzV05MRX_DED_yv33aRpDkdI,11136
+coordinate_system/complex_geometric_physics.py,sha256=6BjYorC2LanmSe6_9hMoUcvpDuUJOAUtbi03-QzXBwU,17045
+coordinate_system/coordinate_system.cp313-win_amd64.pyd,sha256=Tpw0YGsvA5c_iti-cZSOitB4FG0ZVM1nO7P8RAQLGmo,499200
+coordinate_system/curve_interpolation.py,sha256=9NksSvdnSp1BFHPfmwYa6cUC_eyx5Ktp4NXqpzq8uk4,14805
+coordinate_system/differential_geometry.py,sha256=7-82wvnL5czyBQw319jgNHBHUamSSK4PBlSdhBsfxYA,32292
+coordinate_system/spectral_geometry.py,sha256=s2r3A7YBuk-NgF8udAxV88MfXd7CPnwJv9Un0kAbTZM,51995
+coordinate_system/u3_frame.py,sha256=2U7cIYb93P9Cn3qE7Z67PsG6Vg8BOxFkYG32vse7bvA,29763
+coordinate_system/visualization.py,sha256=rNx_ciPg2ITcNXWoANupdY2wpyGDHChLnxYGWDK_tTA,34265
+coordinate_system-7.1.0.dist-info/LICENSE,sha256=tDnRkJxBYPzWdfh2gArRqrUPJxQZRZHJVs68qqBHIq4,1083
+coordinate_system-7.1.0.dist-info/METADATA,sha256=gimpzb3wzTYLlSWE-uogZxP-zJ2HxywpL1Ei-DZ3L0Y,10539
+coordinate_system-7.1.0.dist-info/WHEEL,sha256=4-iQBlRoDdX1wfPofc7KLWa5Cys4eZSgXs6GVU8fKlQ,101
+coordinate_system-7.1.0.dist-info/top_level.txt,sha256=R6LguuPPZ5esrIsDTqPGi9UxCvZPIXwn7KRKX87c79M,18
+coordinate_system-7.1.0.dist-info/RECORD,,

coordinate_system-7.0.2.dist-info/RECORD

@@ -1,13 +0,0 @@
-coordinate_system/__init__.py,sha256=4dQeBXaVEKnZN2l63-GM8LKgbsse8D1xNZrSUZBNg7c,11137
-coordinate_system/complex_geometric_physics.py,sha256=F5Nc8v6J3glH4GduEh2naUftInXCJz2icMeLUTBPGwk,17046
-coordinate_system/coordinate_system.cp313-win_amd64.pyd,sha256=Tpw0YGsvA5c_iti-cZSOitB4FG0ZVM1nO7P8RAQLGmo,499200
-coordinate_system/curve_interpolation.py,sha256=9NksSvdnSp1BFHPfmwYa6cUC_eyx5Ktp4NXqpzq8uk4,14805
-coordinate_system/differential_geometry.py,sha256=zSy2eBjmNUYlaxdiMlaW2ehsS1h0dOLtq5IfAknhNYY,31648
-coordinate_system/spectral_geometry.py,sha256=s2r3A7YBuk-NgF8udAxV88MfXd7CPnwJv9Un0kAbTZM,51995
-coordinate_system/u3_frame.py,sha256=2U7cIYb93P9Cn3qE7Z67PsG6Vg8BOxFkYG32vse7bvA,29763
-coordinate_system/visualization.py,sha256=rNx_ciPg2ITcNXWoANupdY2wpyGDHChLnxYGWDK_tTA,34265
-coordinate_system-7.0.2.dist-info/LICENSE,sha256=tDnRkJxBYPzWdfh2gArRqrUPJxQZRZHJVs68qqBHIq4,1083
-coordinate_system-7.0.2.dist-info/METADATA,sha256=ZyaDh4NMUiCm9LCvdtqYU1yRGktDNoFb2pkPo1-RSiA,10539
-coordinate_system-7.0.2.dist-info/WHEEL,sha256=4-iQBlRoDdX1wfPofc7KLWa5Cys4eZSgXs6GVU8fKlQ,101
-coordinate_system-7.0.2.dist-info/top_level.txt,sha256=R6LguuPPZ5esrIsDTqPGi9UxCvZPIXwn7KRKX87c79M,18
-coordinate_system-7.0.2.dist-info/RECORD,,