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

coordinate_system/__init__.py

@@ -23,13 +23,23 @@ Group Correspondence:
 - FourierFrame ∈ GL(1,C) = U(1) × R⁺
 - U3Frame ∈ U(3) = SU(3) × U(1)
 
-Version: 7.0.1
-DOI: https://doi.org/10.5281/zenodo.18217542
+**Authors:** Pan Guojun
+Version: 7.0.1
+**DOI:** https://doi.org/10.5281/zenodo.14435613
 """
 
 __version__ = '7.0.1'
 
-from .coordinate_system import vec3, vec2
+from .coordinate_system import (
+    vec3,
+    vec2,
+    cross,
+    cross_right,
+    set_handedness,
+    get_handedness,
+    is_left_handed,
+    is_right_handed,
+)
 from .coordinate_system import quat
 from .coordinate_system import coord3
 
@@ -41,24 +51,32 @@ from .differential_geometry import (
     Torus,
 
     # Core classes
-    MetricTensor,
-    GradientResult,
-    IntrinsicGradientOperator,
-    IntrinsicGradientCurvatureCalculator,
-    CurvatureCalculator,
+    MetricTensor,
+    GradientResult,
+    IntrinsicGradientOperator,
+    IntrinsicGradientCurvatureCalculator,
+    LieGroupCurvatureCalculator,
+    CurvatureCalculator,
 
     # Intrinsic gradient method functions (default)
-    compute_gaussian_curvature,
-    compute_mean_curvature,
-    compute_riemann_curvature,
-    compute_all_curvatures,
-    compute_intrinsic_gradient,
-
-    # Classical method functions
-    gaussian_curvature_classical,
-    mean_curvature_classical,
-    principal_curvatures_classical,
-    all_curvatures_classical,
+    compute_gaussian_curvature,
+    compute_mean_curvature,
+    compute_riemann_curvature,
+    compute_curvature_tensor,
+    compute_all_curvatures,
+    compute_intrinsic_gradient,
+
+    # Intrinsic gradient legacy aliases
+    intrinsic_gradient_gaussian_curvature,
+    intrinsic_gradient_mean_curvature,
+    intrinsic_gradient_principal_curvatures,
+    intrinsic_gradient_all_curvatures,
+
+    # Classical method functions
+    gaussian_curvature_classical,
+    mean_curvature_classical,
+    principal_curvatures_classical,
+    all_curvatures_classical,
 
     # Backward compatibility aliases
     gaussian_curvature,
@@ -66,8 +84,9 @@ from .differential_geometry import (
     principal_curvatures,
     all_curvatures,
 
-    # Method comparison
-    compare_methods,
+    # Method comparison
+    compare_methods,
+    gaussian_curvature_lie,
 
     # Utility functions
     derivative_5pt,
@@ -115,10 +134,11 @@ from .u3_frame import (
 )
 
 # Complex Geometric Physics module (Christmas Equation, CFUT)
-from .complex_geometric_physics import (
-    # Core classes (U3Frame imported from u3_frame module above)
-    EnergyMomentumTensor,
-    ChristmasEquation,
+from .complex_geometric_physics import (
+    # Core classes (U3Frame imported from u3_frame module above)
+    ComplexFrame,
+    EnergyMomentumTensor,
+    ChristmasEquation,
 
     # Utility functions
     create_flat_spacetime_frame,
@@ -158,22 +178,28 @@ __all__ = [
     # Constants
     'ZERO3', 'UNITX', 'UNITY', 'UNITZ', 'ONE3', 'ONE4', 'ONEC',
 
-    # Core types
-    'vec3', 'vec2', 'quat', 'coord3', 'lerp',
+    # Core types
+    'vec3', 'vec2', 'quat', 'coord3', 'lerp', 'cross', 'cross_right',
+    'set_handedness', 'get_handedness', 'is_left_handed', 'is_right_handed',
 
     # Differential geometry - Surface classes
     'Surface', 'Sphere', 'Torus',
 
     # Differential geometry - Core classes
     'MetricTensor', 'GradientResult',
-    'IntrinsicGradientOperator', 'IntrinsicGradientCurvatureCalculator',
-    'CurvatureCalculator',
+    'IntrinsicGradientOperator', 'IntrinsicGradientCurvatureCalculator',
+    'LieGroupCurvatureCalculator',
+    'CurvatureCalculator',
 
     # Differential geometry - Intrinsic gradient method (default)
-    'compute_gaussian_curvature', 'compute_mean_curvature',
-    'compute_riemann_curvature', 'compute_all_curvatures',
-    'compute_intrinsic_gradient',
-
+    'compute_gaussian_curvature', 'compute_mean_curvature',
+    'compute_riemann_curvature', 'compute_curvature_tensor', 'compute_all_curvatures',
+    'compute_intrinsic_gradient',
+
+    # Differential geometry - Intrinsic gradient legacy aliases
+    'intrinsic_gradient_gaussian_curvature', 'intrinsic_gradient_mean_curvature',
+    'intrinsic_gradient_principal_curvatures', 'intrinsic_gradient_all_curvatures',
+
     # Differential geometry - Classical method
     'gaussian_curvature_classical', 'mean_curvature_classical',
     'principal_curvatures_classical', 'all_curvatures_classical',
@@ -183,8 +209,8 @@ __all__ = [
     'principal_curvatures', 'all_curvatures',
 
     # Differential geometry - Comparison and utilities
-    'compare_methods',
-    'derivative_5pt', 'derivative_2nd_5pt', 'richardson_extrapolation',
+    'compare_methods', 'gaussian_curvature_lie',
+    'derivative_5pt', 'derivative_2nd_5pt', 'richardson_extrapolation',
 
     # Spectral geometry module (FourierFrame, GL(1,C))
     'FourierFrame', 'FourierFrameSpectrum',
@@ -193,11 +219,17 @@ __all__ = [
     'spectral_transform', 'inverse_spectral_transform',
     'HBAR', 'GPU_AVAILABLE',
 
-    # U(3) Frame module (Gauge theory)
-    'U3Frame', 'SU3Component',
-    'GaugeConnection', 'FieldStrength',
-    'SymmetryBreakingPotential',
-
+    # U(3) Frame module (Gauge theory)
+    'U3Frame', 'SU3Component',
+    'GaugeConnection', 'FieldStrength',
+    'SymmetryBreakingPotential',
+
+    # Complex Geometric Physics
+    'ComplexFrame', 'EnergyMomentumTensor', 'ChristmasEquation',
+    'create_flat_spacetime_frame', 'create_curved_spacetime_frame',
+    'create_gauge_field_frame', 'M_PLANCK', 'LAMBDA_TOPO',
+    'ALPHA_FS', 'LAMBDA_C', 'ALPHA_PROJECTION',
+
     # Visualization
     'CoordinateSystemVisualizer', 'CurveVisualizer', 'ParametricCurve',
     'visualize_coord_system', 'visualize_curve',
@@ -284,16 +316,26 @@ class CoordTuple(tuple):
         return (result.x, result.y, result.z)
 
 
-# Store original coord3 operators
-_original_coord3_mul = coord3.__mul__
-_original_coord3_rmul = coord3.__rmul__
-_original_coord3_truediv = getattr(coord3, '__truediv__', None)
-
-
-def _new_coord3_mul(self, other):
-    """Enhanced multiplication operator for coord3."""
-    if isinstance(other, tuple):
-        other = CoordTuple(other)
+# Store original coord3 operators
+_original_coord3_init = coord3.__init__
+_original_coord3_mul = coord3.__mul__
+_original_coord3_rmul = coord3.__rmul__
+_original_coord3_truediv = getattr(coord3, '__truediv__', None)
+
+
+def _new_coord3_init(self, *args, **kwargs):
+    """Extended coord3 initializer with two-axis convenience."""
+    if len(args) == 2 and all(isinstance(arg, vec3) for arg in args):
+        ux, uy = args
+        uz = ux.cross(uy)
+        return _original_coord3_init(self, ux, uy, uz)
+    return _original_coord3_init(self, *args, **kwargs)
+
+
+def _new_coord3_mul(self, other):
+    """Enhanced multiplication operator for coord3."""
+    if isinstance(other, tuple):
+        other = CoordTuple(other)
         return other * self
     return _original_coord3_mul(self, other)
 
@@ -316,7 +358,8 @@ def _new_coord3_truediv(self, other):
     raise TypeError(f"unsupported operand type(s) for /: 'coord3' and {type(other).__name__}")
 
 
-# Apply enhancements to coord3 operators
-coord3.__mul__ = _new_coord3_mul
-coord3.__rmul__ = _new_coord3_rmul
-coord3.__truediv__ = _new_coord3_truediv
+# Apply enhancements to coord3 operators
+coord3.__init__ = _new_coord3_init
+coord3.__mul__ = _new_coord3_mul
+coord3.__rmul__ = _new_coord3_rmul
+coord3.__truediv__ = _new_coord3_truediv

coordinate_system/complex_geometric_physics.py

@@ -1,33 +1,34 @@
-"""
-Complex Geometric Physics - Unified Framework
-================================================================================
-
-Implementation of the "Christmas Equation" and Complex Frame Unified Theory (CFUT)
-
-The Christmas Equation (Complete complex form):
-    M_P²/2 Ĝ_μν[U] + λ/(32π²) ∇̂_(μ K̄_ν)[U] = T̂_μν^(top)[U] + T̂_μν^(mat)
-
-Where:
-    - Ĝ_μν: Einstein tensor from complex frame U(x)
-    - K̄_μ: Chern-Simons current (topological)
-    - T̂_μν^(top): Topological energy-momentum tensor
-    - T̂_μν^(mat): Matter energy-momentum tensor
-
-Core Theory:
-    - Universe as U(3) complex frame field: U(x) ∈ U(3)
-    - Real-imaginary decomposition: U = U^(R) + iU^(I)
-    - Real part: geometric properties (metric, curvature, spacetime)
-    - Imaginary part: topological properties (phase winding, gauge symmetry)
-
-Physical Interpretation:
-    - Geometry + Topology = Complex Matter + Topological Force
-    - Unifies gravity, gauge fields, dark matter, and topology
-
-Author: Enhanced by AI following theoretical framework
-Date: 2025-01-14
-Version: 7.0.1
-DOI: https://doi.org/10.5281/zenodo.18217542
-"""
+"""
+Complex Geometric Physics - Unified Framework
+================================================================================
+
+Implementation of the "Christmas Equation" and Complex Frame Unified Theory (CFUT).
+
+The Christmas Equation (complete complex form):
+    (M_P^2 / 2) Ĝ_{μν}[U] + (λ / (32π^2)) ∇̂_{(μ} K̄_{ν)}[U]
+        = T̂_{μν}^{(top)}[U] + T̂_{μν}^{(mat)}
+
+Where:
+    - Ĝ_{μν}: Einstein tensor from complex frame U(x)
+    - K̄_μ: Chern-Simons current (topological)
+    - T̂_{μν}^{(top)}: Topological energy-momentum tensor
+    - T̂_{μν}^{(mat)}: Matter energy-momentum tensor
+
+Core Theory:
+    - Universe as U(3) complex frame field: U(x) ∈ U(3)
+    - Real-imaginary decomposition: U = U^(R) + i U^(I)
+    - Real part: geometric properties (metric, curvature, spacetime)
+    - Imaginary part: topological properties (phase winding, gauge symmetry)
+
+Physical Interpretation:
+    - Geometry + Topology = Complex Matter + Topological Force
+    - Unifies gravity, gauge fields, dark matter, and topology
+
+**Authors:** Pan Guojun
+Date: 2025-01-14
+Version: 7.0.1
+**DOI:** https://doi.org/10.5281/zenodo.14435613
+"""
 
 __version__ = '7.0.1'
 
@@ -36,20 +37,20 @@ from typing import Tuple, Optional, Callable, Dict, Any
 from dataclasses import dataclass
 import warnings
 
-# Physical constants (SI units for precision calculations)
-HBAR_SI = 1.054571817e-34  # Reduced Planck constant [J·s]
-H_PLANCK_SI = 6.62607015e-34  # Planck constant [J·s]
+# Physical constants (SI units for precision calculations)
+HBAR_SI = 1.054571817e-34  # Reduced Planck constant [J·s]
+H_PLANCK_SI = 6.62607015e-34  # Planck constant [J·s]
 C_LIGHT_SI = 2.99792458e8  # Speed of light [m/s]
 E_CHARGE_SI = 1.602176634e-19  # Elementary charge [C]
 EPSILON_0_SI = 8.8541878128e-12  # Vacuum permittivity [F/m]
 M_ELECTRON_SI = 9.1093837015e-31  # Electron mass [kg]
 
 # Derived constants
-ALPHA_FS = E_CHARGE_SI**2 / (4 * np.pi * EPSILON_0_SI * HBAR_SI * C_LIGHT_SI)  # Fine structure constant ≈ 1/137
-LAMBDA_C = H_PLANCK_SI / (M_ELECTRON_SI * C_LIGHT_SI)  # Compton wavelength [m]
-ALPHA_PROJECTION = ALPHA_FS * LAMBDA_C  # Projection factor α ≈ 1.77×10⁻¹⁴ m
-
-# Physical constants (natural units: ℏ = c = 1)
+ALPHA_FS = E_CHARGE_SI**2 / (4 * np.pi * EPSILON_0_SI * HBAR_SI * C_LIGHT_SI)  # Fine structure constant ≈ 1/137
+LAMBDA_C = H_PLANCK_SI / (M_ELECTRON_SI * C_LIGHT_SI)  # Compton wavelength [m]
+ALPHA_PROJECTION = ALPHA_FS * LAMBDA_C  # Projection factor α ≈ 1.77×10⁻¹⁴ m
+
+# Physical constants (natural units: ℏ = c = 1)
 M_PLANCK = 2.435e18  # Planck mass [GeV]
 LAMBDA_TOPO = 0.1008  # Topological coupling constant (from theory)
 HBAR = 1.0  # Reduced Planck constant (natural units)
@@ -67,30 +68,33 @@ except ImportError:
         FieldStrength = None
         warnings.warn("U3Frame not available. Some features will be limited.")
 
+# Compatibility alias for older APIs.
+ComplexFrame = U3Frame
+
 
 # ============================================================
 # Core Data Structures
 # ============================================================
 
 @dataclass
-class EnergyMomentumTensor:
-    """
-    Energy-momentum tensor T_μν
-
-    Decomposition:
-        T_μν = T_μν^(R) + iT_μν^(I)
-        - Real part: mass-energy density
-        - Imaginary part: charge current
-    """
-    real_part: np.ndarray  # 4×4 real symmetric tensor
-    imag_part: np.ndarray  # 4×4 real tensor
+class EnergyMomentumTensor:
+    """
+    Energy-momentum tensor T_{μν}
+
+    Decomposition:
+        T_{μν} = T_{μν}^{(R)} + i T_{μν}^{(I)}
+        - Real part: mass-energy density
+        - Imaginary part: charge current
+    """
+    real_part: np.ndarray  # 4×4 real symmetric tensor
+    imag_part: np.ndarray  # 4×4 real tensor
 
     def __post_init__(self):
-        """Validate tensor structure"""
-        if self.real_part.shape != (4, 4):
-            raise ValueError(f"Real part must be 4×4, got {self.real_part.shape}")
-        if self.imag_part.shape != (4, 4):
-            raise ValueError(f"Imaginary part must be 4×4, got {self.imag_part.shape}")
+        """Validate tensor structure"""
+        if self.real_part.shape != (4, 4):
+            raise ValueError(f"Real part must be 4×4, got {self.real_part.shape}")
+        if self.imag_part.shape != (4, 4):
+            raise ValueError(f"Imaginary part must be 4×4, got {self.imag_part.shape}")
 
     @property
     def complex_tensor(self) -> np.ndarray:
@@ -110,46 +114,47 @@ class EnergyMomentumTensor:
 # Christmas Equation Implementation
 # ============================================================
 
-class ChristmasEquation:
-    """
-    The Christmas Equation - Unified field equation
-
-    Complete form:
-        M_P²/2 Ĝ_μν[U] + λ/(32π²) ∇̂_(μ K̄_ν)[U] = T̂_μν^(top)[U] + T̂_μν^(mat)
-
-    Components:
-        - Left side: Geometry + Topology
-        - Right side: Topological source + Matter source
-    """
-
-    def __init__(self,
-                 planck_mass: float = M_PLANCK,
-                 topo_coupling: float = LAMBDA_TOPO,
-                 projection_factor: float = ALPHA_PROJECTION):
-        """
-        Initialize Christmas Equation solver
-
-        Args:
-            planck_mass: Planck mass M_P [GeV]
-            topo_coupling: Topological coupling constant λ
-            projection_factor: Projection factor α = α_fs × λ_c [m]
-        """
+class ChristmasEquation:
+    """
+    The Christmas Equation - unified field equation
+
+    Complete form:
+        (M_P^2 / 2) Ĝ_{μν}[U] + (λ / (32π^2)) ∇̂_{(μ} K̄_{ν)}[U]
+            = T̂_{μν}^{(top)}[U] + T̂_{μν}^{(mat)}
+
+    Components:
+        - Left side: Geometry + Topology
+        - Right side: Topological source + Matter source
+    """
+
+    def __init__(self,
+                 planck_mass: float = M_PLANCK,
+                 topo_coupling: float = LAMBDA_TOPO,
+                 projection_factor: float = ALPHA_PROJECTION):
+        """
+        Initialize Christmas Equation solver.
+
+        Args:
+            planck_mass: Planck mass M_P [GeV]
+            topo_coupling: Topological coupling constant λ
+            projection_factor: Projection factor α = α_fs × λ_c [m]
+        """
         self.M_P = planck_mass
         self.lambda_topo = topo_coupling
         self.alpha_proj = projection_factor
 
-    def einstein_tensor(self, frame: U3Frame) -> np.ndarray:
-        """
-        Compute Einstein tensor Ĝ_μν from complex frame
-
-        Ĝ_μν = R_μν - (1/2)g_μν R
-
-        Args:
-            frame: Complex frame field U(x)
-
-        Returns:
-            3×3 Einstein tensor (spatial part)
-        """
+    def einstein_tensor(self, frame: U3Frame) -> np.ndarray:
+        """
+        Compute Einstein tensor Ĝ_{μν} from complex frame.
+
+        Ĝ_{μν} = R_{μν} - (1/2) g_{μν} R
+
+        Args:
+            frame: Complex frame field U(x)
+
+        Returns:
+            3×3 Einstein tensor (spatial part)
+        """
         # Extract metric from real part
         g = frame.metric_tensor
 
@@ -166,14 +171,14 @@ class ChristmasEquation:
 
         return G_tensor
 
-    def chern_simons_current(self, frame: U3Frame) -> np.ndarray:
-        """
-        Compute Chern-Simons current K̄_μ from imaginary part
-
-        K̄_μ = ε_μνρσ Tr(A^ν F^ρσ - (2/3)A^ν A^ρ A^σ)
-
-        Args:
-            frame: Complex frame field U(x)
+    def chern_simons_current(self, frame: U3Frame) -> np.ndarray:
+        """
+        Compute Chern-Simons current K̄_μ from imaginary part.
+
+        K̄_μ = ε_{μνρσ} Tr(A^ν F^{ρσ} - (2/3) A^ν A^ρ A^σ)
+
+        Args:
+            frame: Complex frame field U(x)
 
         Returns:
             4-vector Chern-Simons current
@@ -193,9 +198,9 @@ class ChristmasEquation:
 
         return K_current
 
-    def topological_energy_momentum(self, frame: U3Frame) -> EnergyMomentumTensor:
-        """
-        Compute topological energy-momentum tensor T̂_μν^(top)
+    def topological_energy_momentum(self, frame: U3Frame) -> EnergyMomentumTensor:
+        """
+        Compute topological energy-momentum tensor T̂_{μν}^{(top)}.
 
         From topological defects (instantons, vortices)
 
@@ -222,13 +227,14 @@ class ChristmasEquation:
 
         return EnergyMomentumTensor(T_real, T_imag)
 
-    def solve_christmas_equation(self,
-                                 frame: U3Frame,
-                                 matter_tensor: EnergyMomentumTensor) -> Dict[str, Any]:
-        """
-        Solve the Christmas Equation
-
-        M_P²/2 Ĝ_μν + λ/(32π²) ∇̂_(μ K̄_ν) = T̂_μν^(top) + T̂_μν^(mat)
+    def solve_christmas_equation(self,
+                                 frame: U3Frame,
+                                 matter_tensor: EnergyMomentumTensor) -> Dict[str, Any]:
+        """
+        Solve the Christmas Equation.
+
+        (M_P^2 / 2) Ĝ_{μν} + (λ / (32π^2)) ∇̂_{(μ} K̄_{ν)}
+            = T̂_{μν}^{(top)} + T̂_{μν}^{(mat)}
 
         Args:
             frame: Complex frame field U(x)
@@ -275,8 +281,8 @@ class ChristmasEquation:
         """
         Compute Ricci tensor from frame (simplified)
 
-        In full implementation, use intrinsic gradient method:
-        R_μν = [G_μ, G_ν] where G_μ = ∂_μ log U
+        In full implementation, use intrinsic gradient method:
+        R_{μν} = [G_μ, G_ν] where G_μ = ∂_μ log U
 
         Args:
             frame: Complex frame field
@@ -295,7 +301,7 @@ class ChristmasEquation:
 
     def _compute_field_strength(self, frame: U3Frame) -> np.ndarray:
         """
-        Compute field strength tensor F_μν = ∂_μ A_ν - ∂_ν A_μ + [A_μ, A_ν]
+        Compute field strength tensor F_{μν} = ∂_μ A_ν - ∂_ν A_μ + [A_μ, A_ν]
 
         Args:
             frame: Complex frame field
@@ -305,14 +311,14 @@ class ChristmasEquation:
         """
         A = frame.gauge_potential
 
-        # Simplified: F ≈ [A, A]
+        # Simplified: F ~= [A, A]
         F = A @ A - A.T @ A.T
 
         return F
 
     def _compute_topological_charge(self, frame: U3Frame) -> float:
         """
-        Compute topological charge Q = (1/32π²) ∫ Tr(F ∧ F)
+        Compute topological charge Q = (1/(32π^2)) ∫ Tr(F ∧ F)
 
         Args:
             frame: Complex frame field
@@ -412,6 +418,7 @@ def create_gauge_field_frame(field_strength: float = 0.1,
 # ============================================================
 
 __all__ = [
+    'ComplexFrame',
     'U3Frame',
     'EnergyMomentumTensor',
     'ChristmasEquation',
@@ -458,7 +465,7 @@ def demonstrate():
 
     solver = ChristmasEquation()
     print(f"   Planck mass M_P: {solver.M_P:.3e} GeV")
-    print(f"   Topological coupling λ: {solver.lambda_topo:.4f}")
+    print(f"   Topological coupling λ: {solver.lambda_topo:.4f}")
 
     # 3. Compute geometric quantities
     print("\n3. Geometric Quantities")
@@ -497,12 +504,13 @@ def demonstrate():
 
     # 6. Theory summary
     print("\n" + "=" * 80)
-    print("Core Theory Summary:")
-    print("  • Christmas Equation: M_P²/2 Ĝ_μν + λ/(32π²) ∇̂_(μ K̄_ν) = T̂_μν^(top) + T̂_μν^(mat)")
-    print("  • U(x) = U^(R)(x) + iU^(I)(x)  [Complex frame decomposition]")
-    print("  • Real part: Geometry (metric, curvature, spacetime)")
-    print("  • Imaginary part: Topology (phase, gauge field, winding)")
-    print("  • Unifies: Gravity + Gauge fields + Dark matter + Topology")
+    print("Core Theory Summary:")
+    print("  - Christmas Equation: (M_P^2 / 2) Ĝ_{μν} + (λ / (32π^2)) ∇̂_{(μ} K̄_{ν)}"
+          " = T̂_{μν}^{(top)} + T̂_{μν}^{(mat)}")
+    print("  - U(x) = U^(R)(x) + iU^(I)(x)  [Complex frame decomposition]")
+    print("  - Real part: Geometry (metric, curvature, spacetime)")
+    print("  - Imaginary part: Topology (phase, gauge field, winding)")
+    print("  - Unifies: Gravity + Gauge fields + Dark matter + Topology")
     print("=" * 80)
 
 

coordinate_system/curve_interpolation.py

@@ -12,9 +12,10 @@ Main Features:
 - B-spline and frame field hybrid interpolation
 - Curvature distribution analysis
 
-Author: PanGuoJun
-Date: 2025-12-01
-"""
+**Authors:** Pan Guojun
+Date: 2025-12-01
+**DOI:** https://doi.org/10.5281/zenodo.14435613
+"""
 
 import math
 import numpy as np