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

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