coordinate-system 7.0.0__tar.gz → 7.0.2__tar.gz

{coordinate_system-7.0.0/coordinate_system.egg-info → coordinate_system-7.0.2}/PKG-INFO

@@ -1,16 +1,16 @@
 Metadata-Version: 2.1
 Name: coordinate_system
-Version: 7.0.0
+Version: 7.0.2
 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: PanGuoJun
-Author-email: PanGuoJun <18858146@qq.com>
+Author: Pan Guojun
+Author-email: Pan Guojun <18858146@qq.com>
 License: MIT
 Project-URL: Homepage, https://github.com/panguojun/Coordinate-System
 Project-URL: Documentation, https://github.com/panguojun/Coordinate-System/blob/main/README.md
 Project-URL: Repository, https://github.com/panguojun/Coordinate-System
 Project-URL: Bug Reports, https://github.com/panguojun/Coordinate-System/issues
-Project-URL: DOI, https://zenodo.org/records/14435613
+Project-URL: DOI, https://doi.org/10.5281/zenodo.14435613
 Keywords: 3d,math,vector,quaternion,coordinate-system,geometry,differential-geometry,curvature,spectral-analysis,christmas-equation,cfut,unified-field-theory,complex-frame
 Platform: Windows
 Platform: Linux
@@ -50,17 +50,26 @@ Requires-Dist: matplotlib>=3.3.0
 [![Python](https://img.shields.io/pypi/pyversions/coordinate-system.svg)](https://pypi.org/project/coordinate-system/)
 [![License](https://img.shields.io/badge/license-MIT-green.svg)](LICENSE)
 
-**Author:** PanGuoJun
-**Version:** 7.0.0-alpha
+**Authors:** Pan Guojun
+**Version:** 7.0.1
 **License:** MIT
-**DOI:** [10.5281/zenodo.14435613](https://zenodo.org/records/14435613)
+**DOI:** https://doi.org/10.5281/zenodo.14435613
 
 ---
 
+## What's New in v7.0.1 (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
+- **Numerical Verification**: New `verify_cfut_theory.py` script with complete experimental validation
+- **Equation Decomposition**: Real/imaginary parts analysis with classical limit verification
+- **Experimental Comparison**: All calculations validated against CODATA 2018 data (error < 10⁻⁹)
+- **Documentation**: Added `VERIFICATION_REPORT.md` with detailed numerical results
+
 ## What's New in v7.0.0-alpha (2026-01-14)
 
-- **🎄 Christmas Equation**: New `complex_geometric_physics` module implementing unified field theory
-- **Complex Frame Unified Theory (CFUT)**: Geometry + Topology = Complex Matter + Topological Force
+- **Complex Geometric Physics Module**: New `complex_geometric_physics` module for field theory calculations
+- **U(3) Frame Implementation**: U(3) complex field with real-imaginary decomposition
 - **Einstein Tensor**: Compute Ĝ_μν from complex frame field U(x)
 - **Chern-Simons Current**: Topological current K̄_μ for gauge field analysis
 - **Energy-Momentum Tensor**: Real-imaginary decomposition for matter and topology
@@ -73,7 +82,7 @@ Requires-Dist: matplotlib>=3.3.0
 ```
 coordinate_system/
 ├── coordinate_system.pyd/.so      # C++ core (vec3, quat, coord3)
-├── complex_geometric_physics.py   # 🎄 Christmas Equation, CFUT unified theory
+├── complex_geometric_physics.py   # Complex geometric physics module
 ├── spectral_geometry.py           # FourierFrame [GL(1,C)], spectral analysis
 ├── u3_frame.py                    # U3Frame [U(3)], gauge field theory
 ├── differential_geometry.py       # Surface curvature calculation
@@ -86,9 +95,8 @@ coordinate_system/
 | Class | Group | DOF | Use Case |
 |-------|-------|-----|----------|
 | `coord3` | Sim(3) = R³ ⋊ (SO(3) × R⁺) | 10 | 3D coordinate transform |
-| `ComplexFrame` | U(3) complex field | 18 | 🎄 Unified field theory (CFUT) |
 | `FourierFrame` | GL(1,C) = U(1) × R⁺ | 2 | Spectral geometry, heat kernel |
-| `U3Frame` | U(3) = SU(3) × U(1) | 9 | Gauge field theory |
+| `U3Frame` | U(3) = SU(3) × U(1) | 9 | Gauge field theory, gauge transforms |
 
 ---
 
@@ -184,53 +192,6 @@ F_xy = conn_x.field_strength(conn_y)
 S_YM = F_xy.yang_mills_action()
 ```
 
-### 🎄 Complex Geometric Physics (Christmas Equation)
-
-```python
-from coordinate_system import (
-    ComplexFrame,
-    EnergyMomentumTensor,
-    ChristmasEquation,
-    create_flat_spacetime_frame,
-    create_curved_spacetime_frame,
-    create_gauge_field_frame,
-    M_PLANCK,
-    LAMBDA_TOPO
-)
-import numpy as np
-
-# Create complex frames
-flat_frame = create_flat_spacetime_frame()
-curved_frame = create_curved_spacetime_frame(curvature=0.1)
-gauge_frame = create_gauge_field_frame(field_strength=0.1)
-
-# Initialize Christmas Equation solver
-solver = ChristmasEquation()
-print(f"Planck mass: {solver.M_P:.3e} GeV")
-print(f"Topological coupling: {solver.lambda_topo:.4f}")
-
-# Compute geometric quantities
-G_tensor = solver.einstein_tensor(curved_frame)
-K_current = solver.chern_simons_current(gauge_frame)
-
-# Create matter energy-momentum tensor
-matter_real = np.diag([1.0, 0.1, 0.1, 0.1])
-matter_imag = np.zeros((4, 4))
-T_matter = EnergyMomentumTensor(matter_real, matter_imag)
-
-# Solve the Christmas Equation
-# M_P²/2 Ĝ_μν[U] + λ/(32π²) ∇̂_(μ K̄_ν)[U] = T̂_μν^(top)[U] + T̂_μν^(mat)
-solution = solver.solve_christmas_equation(gauge_frame, T_matter)
-print(f"Geometric term norm: {np.linalg.norm(solution['geometric_term']):.6e}")
-print(f"Topological term norm: {np.linalg.norm(solution['topological_term']):.6e}")
-print(f"Equation balanced: {solution['balanced']}")
-```
-
-**The Christmas Equation** unifies geometry and topology:
-- **Left side**: Geometry (Einstein tensor) + Topology (Chern-Simons current)
-- **Right side**: Topological source + Matter source
-- **Real part U^(R)**: Geometric properties (metric, curvature, spacetime)
-- **Imaginary part U^(I)**: Topological properties (phase winding, gauge symmetry)
 
 ---
 
@@ -238,9 +199,7 @@ print(f"Equation balanced: {solution['balanced']}")
 
 | Concept | Formula | Code |
 |---------|---------|------|
-| 🎄 **Christmas Equation** | $\frac{M_P^2}{2} \hat{G}_{\mu\nu}[U] + \frac{\lambda}{32\pi^2} \hat{\nabla}_{(\mu} \bar{K}_{\nu)}[U] = \hat{T}_{\mu\nu}^{(\text{top})}[U] + \hat{T}_{\mu\nu}^{(\text{mat})}$ | `ChristmasEquation.solve_christmas_equation()` |
-| Einstein Tensor | $\hat{G}_{\mu\nu} = R_{\mu\nu} - \frac{1}{2}g_{\mu\nu} R$ | `ChristmasEquation.einstein_tensor()` |
-| Chern-Simons Current | $\bar{K}_\mu = \varepsilon_{\mu\nu\rho\sigma} \text{Tr}(A^\nu F^{\rho\sigma})$ | `ChristmasEquation.chern_simons_current()` |
+| Projection Factor (v7.0.1) | $\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` |
@@ -277,14 +236,24 @@ print(f"Equation balanced: {solution['balanced']}")
 
 ## Changelog
 
+### v7.0.1 (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
+- **Numerical Verification**: New `verify_cfut_theory.py` script with complete validation
+- **Equation Decomposition**: Real/imaginary parts analysis with classical limit verification
+- **Experimental Validation**: All calculations compared with CODATA 2018 data (error < 10⁻⁹)
+- **Documentation**: Added `VERIFICATION_REPORT.md` with detailed numerical results
+- **Bug Fix**: Corrected Compton wavelength calculation (h instead of ℏ)
+
 ### v7.0.0-alpha (2026-01-14)
-- 🎄 **Christmas Equation**: New `complex_geometric_physics.py` module
-- **Complex Frame Unified Theory (CFUT)**: Unified field equation implementation
-- **ComplexFrame**: U(x) = U^(R)(x) + iU^(I)(x) decomposition
+- **Complex Geometric Physics Module**: New `complex_geometric_physics.py` module
+- **U(3) Frame Implementation**: U(3) complex field with real-imaginary decomposition
+- **U3Frame**: U(x) = U^(R)(x) + iU^(I)(x) decomposition
 - **EnergyMomentumTensor**: Real-imaginary tensor decomposition
-- **ChristmasEquation**: Einstein tensor, Chern-Simons current, topological energy-momentum
+- **Unified Field Solver**: Einstein tensor, Chern-Simons current, topological energy-momentum
 - **Complete English Translation**: All documentation and code comments in English
-- **DOI**: Added Zenodo DOI 10.5281/zenodo.14435613
+- **DOI:** https://doi.org/10.5281/zenodo.14435613
 
 ### v6.0.4 (2025-12-08)
 - `frames.py` → `spectral_geometry.py`

{coordinate_system-7.0.0 → coordinate_system-7.0.2}/README.md

@@ -6,17 +6,26 @@
 [![Python](https://img.shields.io/pypi/pyversions/coordinate-system.svg)](https://pypi.org/project/coordinate-system/)
 [![License](https://img.shields.io/badge/license-MIT-green.svg)](LICENSE)
 
-**Author:** PanGuoJun
-**Version:** 7.0.0-alpha
+**Authors:** Pan Guojun
+**Version:** 7.0.1
 **License:** MIT
-**DOI:** [10.5281/zenodo.14435613](https://zenodo.org/records/14435613)
+**DOI:** https://doi.org/10.5281/zenodo.14435613
 
 ---
 
+## What's New in v7.0.1 (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
+- **Numerical Verification**: New `verify_cfut_theory.py` script with complete experimental validation
+- **Equation Decomposition**: Real/imaginary parts analysis with classical limit verification
+- **Experimental Comparison**: All calculations validated against CODATA 2018 data (error < 10⁻⁹)
+- **Documentation**: Added `VERIFICATION_REPORT.md` with detailed numerical results
+
 ## What's New in v7.0.0-alpha (2026-01-14)
 
-- **🎄 Christmas Equation**: New `complex_geometric_physics` module implementing unified field theory
-- **Complex Frame Unified Theory (CFUT)**: Geometry + Topology = Complex Matter + Topological Force
+- **Complex Geometric Physics Module**: New `complex_geometric_physics` module for field theory calculations
+- **U(3) Frame Implementation**: U(3) complex field with real-imaginary decomposition
 - **Einstein Tensor**: Compute Ĝ_μν from complex frame field U(x)
 - **Chern-Simons Current**: Topological current K̄_μ for gauge field analysis
 - **Energy-Momentum Tensor**: Real-imaginary decomposition for matter and topology
@@ -29,7 +38,7 @@
 ```
 coordinate_system/
 ├── coordinate_system.pyd/.so      # C++ core (vec3, quat, coord3)
-├── complex_geometric_physics.py   # 🎄 Christmas Equation, CFUT unified theory
+├── complex_geometric_physics.py   # Complex geometric physics module
 ├── spectral_geometry.py           # FourierFrame [GL(1,C)], spectral analysis
 ├── u3_frame.py                    # U3Frame [U(3)], gauge field theory
 ├── differential_geometry.py       # Surface curvature calculation
@@ -42,9 +51,8 @@ coordinate_system/
 | Class | Group | DOF | Use Case |
 |-------|-------|-----|----------|
 | `coord3` | Sim(3) = R³ ⋊ (SO(3) × R⁺) | 10 | 3D coordinate transform |
-| `ComplexFrame` | U(3) complex field | 18 | 🎄 Unified field theory (CFUT) |
 | `FourierFrame` | GL(1,C) = U(1) × R⁺ | 2 | Spectral geometry, heat kernel |
-| `U3Frame` | U(3) = SU(3) × U(1) | 9 | Gauge field theory |
+| `U3Frame` | U(3) = SU(3) × U(1) | 9 | Gauge field theory, gauge transforms |
 
 ---
 
@@ -140,53 +148,6 @@ F_xy = conn_x.field_strength(conn_y)
 S_YM = F_xy.yang_mills_action()
 ```
 
-### 🎄 Complex Geometric Physics (Christmas Equation)
-
-```python
-from coordinate_system import (
-    ComplexFrame,
-    EnergyMomentumTensor,
-    ChristmasEquation,
-    create_flat_spacetime_frame,
-    create_curved_spacetime_frame,
-    create_gauge_field_frame,
-    M_PLANCK,
-    LAMBDA_TOPO
-)
-import numpy as np
-
-# Create complex frames
-flat_frame = create_flat_spacetime_frame()
-curved_frame = create_curved_spacetime_frame(curvature=0.1)
-gauge_frame = create_gauge_field_frame(field_strength=0.1)
-
-# Initialize Christmas Equation solver
-solver = ChristmasEquation()
-print(f"Planck mass: {solver.M_P:.3e} GeV")
-print(f"Topological coupling: {solver.lambda_topo:.4f}")
-
-# Compute geometric quantities
-G_tensor = solver.einstein_tensor(curved_frame)
-K_current = solver.chern_simons_current(gauge_frame)
-
-# Create matter energy-momentum tensor
-matter_real = np.diag([1.0, 0.1, 0.1, 0.1])
-matter_imag = np.zeros((4, 4))
-T_matter = EnergyMomentumTensor(matter_real, matter_imag)
-
-# Solve the Christmas Equation
-# M_P²/2 Ĝ_μν[U] + λ/(32π²) ∇̂_(μ K̄_ν)[U] = T̂_μν^(top)[U] + T̂_μν^(mat)
-solution = solver.solve_christmas_equation(gauge_frame, T_matter)
-print(f"Geometric term norm: {np.linalg.norm(solution['geometric_term']):.6e}")
-print(f"Topological term norm: {np.linalg.norm(solution['topological_term']):.6e}")
-print(f"Equation balanced: {solution['balanced']}")
-```
-
-**The Christmas Equation** unifies geometry and topology:
-- **Left side**: Geometry (Einstein tensor) + Topology (Chern-Simons current)
-- **Right side**: Topological source + Matter source
-- **Real part U^(R)**: Geometric properties (metric, curvature, spacetime)
-- **Imaginary part U^(I)**: Topological properties (phase winding, gauge symmetry)
 
 ---
 
@@ -194,9 +155,7 @@ print(f"Equation balanced: {solution['balanced']}")
 
 | Concept | Formula | Code |
 |---------|---------|------|
-| 🎄 **Christmas Equation** | $\frac{M_P^2}{2} \hat{G}_{\mu\nu}[U] + \frac{\lambda}{32\pi^2} \hat{\nabla}_{(\mu} \bar{K}_{\nu)}[U] = \hat{T}_{\mu\nu}^{(\text{top})}[U] + \hat{T}_{\mu\nu}^{(\text{mat})}$ | `ChristmasEquation.solve_christmas_equation()` |
-| Einstein Tensor | $\hat{G}_{\mu\nu} = R_{\mu\nu} - \frac{1}{2}g_{\mu\nu} R$ | `ChristmasEquation.einstein_tensor()` |
-| Chern-Simons Current | $\bar{K}_\mu = \varepsilon_{\mu\nu\rho\sigma} \text{Tr}(A^\nu F^{\rho\sigma})$ | `ChristmasEquation.chern_simons_current()` |
+| Projection Factor (v7.0.1) | $\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` |
@@ -233,14 +192,24 @@ print(f"Equation balanced: {solution['balanced']}")
 
 ## Changelog
 
+### v7.0.1 (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
+- **Numerical Verification**: New `verify_cfut_theory.py` script with complete validation
+- **Equation Decomposition**: Real/imaginary parts analysis with classical limit verification
+- **Experimental Validation**: All calculations compared with CODATA 2018 data (error < 10⁻⁹)
+- **Documentation**: Added `VERIFICATION_REPORT.md` with detailed numerical results
+- **Bug Fix**: Corrected Compton wavelength calculation (h instead of ℏ)
+
 ### v7.0.0-alpha (2026-01-14)
-- 🎄 **Christmas Equation**: New `complex_geometric_physics.py` module
-- **Complex Frame Unified Theory (CFUT)**: Unified field equation implementation
-- **ComplexFrame**: U(x) = U^(R)(x) + iU^(I)(x) decomposition
+- **Complex Geometric Physics Module**: New `complex_geometric_physics.py` module
+- **U(3) Frame Implementation**: U(3) complex field with real-imaginary decomposition
+- **U3Frame**: U(x) = U^(R)(x) + iU^(I)(x) decomposition
 - **EnergyMomentumTensor**: Real-imaginary tensor decomposition
-- **ChristmasEquation**: Einstein tensor, Chern-Simons current, topological energy-momentum
+- **Unified Field Solver**: Einstein tensor, Chern-Simons current, topological energy-momentum
 - **Complete English Translation**: All documentation and code comments in English
-- **DOI**: Added Zenodo DOI 10.5281/zenodo.14435613
+- **DOI:** https://doi.org/10.5281/zenodo.14435613
 
 ### v6.0.4 (2025-12-08)
 - `frames.py` → `spectral_geometry.py`

{coordinate_system-7.0.0 → coordinate_system-7.0.2}/coordinate_system/__init__.py

@@ -23,12 +23,23 @@ Group Correspondence:
 - FourierFrame ∈ GL(1,C) = U(1) × R⁺
 - U3Frame ∈ U(3) = SU(3) × U(1)
 
-Version: 7.0.0-alpha
+**Authors:** Pan Guojun
+Version: 7.0.1
+**DOI:** https://doi.org/10.5281/zenodo.14435613
 """
 
-__version__ = '7.0.0-alpha'
-
-from .coordinate_system import vec3, vec2
+__version__ = '7.0.1'
+
+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
 
@@ -40,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,
@@ -65,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,
@@ -114,11 +134,11 @@ from .u3_frame import (
 )
 
 # Complex Geometric Physics module (Christmas Equation, CFUT)
-from .complex_geometric_physics import (
-    # Core classes
-    ComplexFrame,
-    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,
@@ -128,6 +148,9 @@ from .complex_geometric_physics import (
     # Constants
     M_PLANCK,
     LAMBDA_TOPO,
+    ALPHA_FS,
+    LAMBDA_C,
+    ALPHA_PROJECTION,
 )
 
 # Visualization module
@@ -155,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',
@@ -180,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',
@@ -190,16 +219,17 @@ __all__ = [
     'spectral_transform', 'inverse_spectral_transform',
     'HBAR', 'GPU_AVAILABLE',
 
-    # U(3) Frame module (Gauge theory)
-    'U3Frame', 'SU3Component',
-    'GaugeConnection', 'FieldStrength',
-    'SymmetryBreakingPotential',
-
-    # Complex Geometric Physics (Christmas Equation, CFUT)
-    'ComplexFrame', 'EnergyMomentumTensor', 'ChristmasEquation',
-    'create_flat_spacetime_frame', 'create_curved_spacetime_frame', 'create_gauge_field_frame',
-    'M_PLANCK', 'LAMBDA_TOPO',
-
+    # 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',
@@ -286,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)
 
@@ -318,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