coordinate-system 2.3.1__cp313-cp313-win_amd64.whl → 2.3.2__cp313-cp313-win_amd64.whl

coordinate_system/curvature.py

@@ -532,7 +532,7 @@ class LieGroupCurvatureCalculator:
     def __init__(
         self,
         surface: Surface,
-        step_size: float = 1e-5,
+        step_size: float = 1e-3,
         scale_factor: float = 2.0,
         use_metric_correction: bool = True,
         use_lie_derivative: bool = True
@@ -542,9 +542,8 @@ class LieGroupCurvatureCalculator:
 
         Args:
             surface: Surface object
-            step_size: Finite difference step size (default: 1e-5)
-                      Note: Smaller than classical method due to frame normalization
-            scale_factor: Scaling factor for connection correction (default: 2.0)
+            step_size: Finite difference step size (default: 1e-3, optimal for O(h²) convergence)
+            scale_factor: Scaling factor for connection correction (default: 2.0, precomputed correction for self-referencing terms)
             use_metric_correction: Whether to apply metric correction (default: True)
             use_lie_derivative: Whether to include Lie derivative term (default: True)
                               Highly recommended for accuracy!

coordinate_system/differential_geometry.py

@@ -344,7 +344,7 @@ class ConnectionOperator:
 
         Args:
             surface: Surface object
-            scale_factor: Scaling factor for correction (default: 2.0, optimal for constant curvature)
+            scale_factor: Scaling factor for correction (default: 2.0, precomputed numerical correction from curvature estimation)
             use_metric_correction: Whether to apply metric correction (default: True)
             step_size: Numerical differentiation step size (default: surface.h)
         """
@@ -382,7 +382,8 @@ class ConnectionOperator:
             t_v = UNITY
 
         # Create coord3 from basis vectors
-        return coord3(ZERO3, t_u, t_v, n)
+        # Use 3-parameter constructor "From three axes" (normalized vectors only, no origin)
+        return coord3(t_u, t_v, n)
 
     def _compute_embedding_frame(self, u: float, v: float) -> coord3:
         """
@@ -589,7 +590,7 @@ def compute_connection(
         surface: Surface object
         u, v: Parameter coordinates
         direction: 'u' or 'v'
-        scale_factor: Scaling factor (default: 2.0)
+        scale_factor: Scaling factor (default: 2.0, precomputed numerical correction)
         use_metric_correction: Whether to apply metric correction (default: True)
 
     Returns:
@@ -629,7 +630,7 @@ def compute_curvature_tensor(
     Args:
         surface: Surface object
         u, v: Parameter coordinates
-        scale_factor: Scaling factor for connection operators
+        scale_factor: Scaling factor for connection operators (default: 2.0, precomputed numerical correction)
         use_metric_correction: Whether to apply metric correction
         use_lie_derivative: Whether to include Lie derivative term (default: True)
                            For spheres, this improves accuracy by 24×!
@@ -701,14 +702,14 @@ def compute_gaussian_curvature(
     """
     Compute Gaussian curvature K at a point
 
-    K = R_12 / det(g)
+    K = R_12^antisym / det(g) = (R_01 - R_10) / (2 × det(g))
 
-    where R_12 is extracted from the curvature tensor R_uv
+    where R_12^antisym is the antisymmetric part of the curvature tensor R_uv
 
     Args:
         surface: Surface object
         u, v: Parameter coordinates
-        scale_factor: Scaling factor
+        scale_factor: Scaling factor (default: 2.0, precomputed correction for numerical method's self-referencing terms)
         use_metric_correction: Whether to apply metric correction
         use_lie_derivative: Whether to include Lie derivative term
 
@@ -732,14 +733,17 @@ def compute_gaussian_curvature(
     # Compute metric
     g = compute_metric(surface, u, v)
 
-    # Extract R_12 (antisymmetric part of tangent-tangent block)
-    # R_12 is the (1,2) element of R_uv matrix
-    # In coord3: uy.x represents the second row, first column
-    R_12 = R_uv.uy.x
+    # Extract antisymmetric part of curvature tensor
+    # R_01 is the (0,1) element (first row, second column)
+    # R_10 is the (1,0) element (second row, first column)
+    # In coord3: ux.y represents R_01, uy.x represents R_10
+    R_01 = R_uv.ux.y
+    R_10 = R_uv.uy.x
+    R_12_antisym = (R_01 - R_10) / 2.0
 
     # Gaussian curvature
     if abs(g.det) > 1e-10:
-        K = R_12 / g.det
+        K = R_12_antisym / g.det
     else:
         K = 0.0
 

{coordinate_system-2.3.1.dist-info → coordinate_system-2.3.2.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: coordinate_system
-Version: 2.3.1
+Version: 2.3.2
 Summary: High-performance 3D coordinate system library with C++ optimized differential geometry and high-precision discrete curvature computation (Gaussian/Mean/Principal curvatures)
 Home-page: https://github.com/panguojun/Coordinate-System
 Author: PanGuoJun

coordinate_system-2.3.2.dist-info/RECORD

@@ -0,0 +1,9 @@
+coordinate_system/__init__.py,sha256=S-h_uoibcICCh5Tk_-Q6_aMwiCQ2vvKUgE4BzqCHpJ0,7086
+coordinate_system/coordinate_system.cp313-win_amd64.pyd,sha256=P0KTnkBHLBixO2-g7WRc608Iw-mvWZS8p1EnvcSYZ7A,490496
+coordinate_system/curvature.py,sha256=8KLC3xupF69W5-DHRy6lfDUgSErc7SbGA3gPiim0CUo,25232
+coordinate_system/differential_geometry.py,sha256=0dPYkiFgAg1UH-O7PQmAJ2Mmm0VQLQSU7bZThy4vZb8,24111
+coordinate_system-2.3.2.dist-info/LICENSE,sha256=tDnRkJxBYPzWdfh2gArRqrUPJxQZRZHJVs68qqBHIq4,1083
+coordinate_system-2.3.2.dist-info/METADATA,sha256=BjFU1EdnAkoZEelXOf4bpFhA3m220Lia4fbcj0yR5yA,22632
+coordinate_system-2.3.2.dist-info/WHEEL,sha256=4-iQBlRoDdX1wfPofc7KLWa5Cys4eZSgXs6GVU8fKlQ,101
+coordinate_system-2.3.2.dist-info/top_level.txt,sha256=R6LguuPPZ5esrIsDTqPGi9UxCvZPIXwn7KRKX87c79M,18
+coordinate_system-2.3.2.dist-info/RECORD,,

coordinate_system-2.3.1.dist-info/RECORD

@@ -1,9 +0,0 @@
-coordinate_system/__init__.py,sha256=S-h_uoibcICCh5Tk_-Q6_aMwiCQ2vvKUgE4BzqCHpJ0,7086
-coordinate_system/coordinate_system.cp313-win_amd64.pyd,sha256=t1odSp5QdNMLQ_W5WW-RRR1ItKzwRfB53TUSCE1AKTI,490496
-coordinate_system/curvature.py,sha256=-KVkyqeuySr2wJaCIbxSPdZ351VuNRMaeyoY32-czd8,25235
-coordinate_system/differential_geometry.py,sha256=ZDYJzg7uLYpak4JWHMn7oufnhpF-Lh2dh8x7UpRQl2w,23645
-coordinate_system-2.3.1.dist-info/LICENSE,sha256=tDnRkJxBYPzWdfh2gArRqrUPJxQZRZHJVs68qqBHIq4,1083
-coordinate_system-2.3.1.dist-info/METADATA,sha256=2Rh4Is9rUKm8hhyrqZDdQ-bGVRu_cdue2LZQIzxsoD0,22632
-coordinate_system-2.3.1.dist-info/WHEEL,sha256=4-iQBlRoDdX1wfPofc7KLWa5Cys4eZSgXs6GVU8fKlQ,101
-coordinate_system-2.3.1.dist-info/top_level.txt,sha256=R6LguuPPZ5esrIsDTqPGi9UxCvZPIXwn7KRKX87c79M,18
-coordinate_system-2.3.1.dist-info/RECORD,,