coordinate-system 2.4.2__cp313-cp313-win_amd64.whl → 2.5.1__cp313-cp313-win_amd64.whl

coordinate_system/__init__.py

@@ -4,65 +4,89 @@ from .coordinate_system import vec3, vec2
 from .coordinate_system import quat
 from .coordinate_system import coord3
 
-# Differential geometry module (v2.2.0+)
+# Differential geometry module (v3.0.0+)
 from .differential_geometry import (
     # Classes
     Surface,
     Sphere,
     Torus,
     MetricTensor,
-    ConnectionOperator,
+    IntrinsicGradientOperator,
+    IntrinsicGradientCurvatureCalculator,
 
     # Main functions
-    compute_metric,
-    compute_connection,
-    compute_curvature_tensor,
-    compute_gaussian_curvature,
-
-    # Aliases
-    compute_frame_derivative,
     compute_intrinsic_gradient,
-    compute_geometric_gradient,
-    compute_metric_tensor,
+    compute_gaussian_curvature,
+    compute_curvature_tensor,
+    compute_all_curvatures,
 )
 
 # High-precision curvature computation module (v2.3.0+)
 from .curvature import (
-    # Main class
+    # Classical method
     CurvatureCalculator,
-
-    # Simplified functions
     gaussian_curvature,
     mean_curvature,
     principal_curvatures,
     all_curvatures,
 
+    # Lie group method
+    LieGroupCurvatureCalculator,
+    gaussian_curvature_lie,
+    compare_methods,
+
     # Utility functions
     derivative_5pt,
     derivative_2nd_5pt,
     richardson_extrapolation,
 )
 
+# Intrinsic Gradient Operator method (v2.4.0+)
+from .intrinsic_gradient_curvature import (
+    IntrinsicGradientOperator,
+    IntrinsicGradientCurvatureCalculator,
+    intrinsic_gradient_gaussian_curvature,
+    intrinsic_gradient_mean_curvature,
+    intrinsic_gradient_principal_curvatures,
+    intrinsic_gradient_all_curvatures,
+)
+
 __all__ = [
+    # Constants
+    'ZERO3','UNITX','UNITY','UNITZ','ONE3','ONE4','ONEC',
+    
     # Core types
     'vec3', 'vec2', 'quat', 'coord3', 'lerp',
 
-    # Differential geometry classes (v2.2.0+)
+    # Differential geometry classes (v3.0.0+)
     'Surface', 'Sphere', 'Torus',
-    'MetricTensor', 'ConnectionOperator',
+    'MetricTensor', 'IntrinsicGradientOperator', 'IntrinsicGradientCurvatureCalculator',
 
-    # Differential geometry functions (v2.2.0+)
-    'compute_metric', 'compute_connection',
-    'compute_curvature_tensor', 'compute_gaussian_curvature',
+    # Differential geometry functions (v3.0.0+)
+    'compute_intrinsic_gradient',
+    'compute_gaussian_curvature',
+    'compute_curvature_tensor',
+    'compute_all_curvatures',
 
-    # Aliases
-    'compute_frame_derivative', 'compute_intrinsic_gradient',
-    'compute_geometric_gradient', 'compute_metric_tensor',
-
-    # High-precision curvature module (v2.3.0+)
+    # High-precision curvature module - Classical method (v2.3.0+)
     'CurvatureCalculator',
     'gaussian_curvature', 'mean_curvature',
     'principal_curvatures', 'all_curvatures',
+
+    # Lie group method (v2.3.0+)
+    'LieGroupCurvatureCalculator',
+    'gaussian_curvature_lie',
+    'compare_methods',
+
+    # Intrinsic Gradient Operator method (v2.4.0+)
+    'IntrinsicGradientOperator',
+    'IntrinsicGradientCurvatureCalculator',
+    'intrinsic_gradient_gaussian_curvature',
+    'intrinsic_gradient_mean_curvature',
+    'intrinsic_gradient_principal_curvatures',
+    'intrinsic_gradient_all_curvatures',
+
+    # Utility functions
     'derivative_5pt', 'derivative_2nd_5pt', 'richardson_extrapolation',
 ]
 

coordinate_system/curvature.py

@@ -12,13 +12,19 @@ Features:
 - Principal directions
 - Convergence analysis
 - Richardson extrapolation
+- Intrinsic Gradient Operator method (NEW)
 
-Method: Classical First/Second Fundamental Forms + 5-point finite differences
-Precision: O(h⁴) convergence, typical error < 0.001%
+Methods Available:
+1. Classical Method: First/Second Fundamental Forms + 5-point finite differences
+   - Precision: O(h⁴) convergence, typical error < 0.001%
 
-Author: PanGuoJun & Claude AI
-Version: 2.3.0
-Date: 2025-10-27
+2. Intrinsic Gradient Operator Method: G_μ = (c(u+h) - c(u))/h
+   - Based on coordinate system changes analysis
+   - Numerically stable and geometrically intuitive
+
+Author: PanGuoJun
+Version: 2.4.0
+Date: 2025-10-30
 """
 
 import numpy as np
@@ -27,9 +33,10 @@ from .differential_geometry import Surface
 from .coordinate_system import coord3, vec3
 from .differential_geometry import (
     MetricTensor,
-    ConnectionOperator,
+    IntrinsicGradientOperator,
+    IntrinsicGradientCurvatureCalculator,
     compute_curvature_tensor as compute_curvature_tensor_lie,
-    compute_metric
+    compute_all_curvatures as compute_all_curvatures_intrinsic
 )
 
 
@@ -588,7 +595,7 @@ class LieGroupCurvatureCalculator:
             )
 
             # Compute metric tensor
-            g = compute_metric(self.surface, u, v)
+            g = MetricTensor.from_surface(self.surface, u, v)
 
             # Extract Gaussian curvature using antisymmetric part
             # K = R_{12}^{antisym} / det(g) = (R_01 - R_10)/2 / det(g)
@@ -644,14 +651,12 @@ class LieGroupCurvatureCalculator:
         self.surface.h = self.h
 
         try:
-            # Compute connection operators
-            conn = ConnectionOperator(
+            # Compute intrinsic gradient operators
+            intrinsic_op = IntrinsicGradientOperator(
                 self.surface,
-                self.scale_factor,
-                self.use_metric_correction,
                 self.h
             )
-            G_u, G_v = conn.compute_both(u, v)
+            G_u, G_v = intrinsic_op.compute_both(u, v)
 
             # Compute curvature tensor
             R_uv = compute_curvature_tensor_lie(
@@ -662,7 +667,7 @@ class LieGroupCurvatureCalculator:
             )
 
             # Compute metric
-            g = compute_metric(self.surface, u, v)
+            g = MetricTensor.from_surface(self.surface, u, v)
 
             # Extract Gaussian curvature using antisymmetric part
             # K = (R_01 - R_10)/2 / det(g)
@@ -778,6 +783,17 @@ def gaussian_curvature_lie(
     return calc.compute_gaussian_curvature(u, v)
 
 
+# ========== Intrinsic Gradient Operator Method ==========
+# Import from the new module
+from .intrinsic_gradient_curvature import (
+    IntrinsicGradientOperator,
+    IntrinsicGradientCurvatureCalculator,
+    intrinsic_gradient_gaussian_curvature,
+    intrinsic_gradient_mean_curvature,
+    intrinsic_gradient_principal_curvatures,
+    intrinsic_gradient_all_curvatures,
+)
+
 # ========== Export ==========
 
 __all__ = [
@@ -793,6 +809,14 @@ __all__ = [
     'gaussian_curvature_lie',
     'compare_methods',
 
+    # Intrinsic Gradient Operator method
+    'IntrinsicGradientOperator',
+    'IntrinsicGradientCurvatureCalculator',
+    'intrinsic_gradient_gaussian_curvature',
+    'intrinsic_gradient_mean_curvature',
+    'intrinsic_gradient_principal_curvatures',
+    'intrinsic_gradient_all_curvatures',
+
     # Utility functions
     'derivative_5pt',
     'derivative_2nd_5pt',