coordinate-system 2.4.0__tar.gz → 2.4.2__tar.gz

{coordinate_system-2.4.0/coordinate_system.egg-info → coordinate_system-2.4.2}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: coordinate_system
-Version: 2.4.0
+Version: 2.4.2
 Summary: High-performance 3D coordinate system library with machine-precision Gaussian curvature computation (0.000% error) using corrected antisymmetric extraction formula
 Home-page: https://github.com/panguojun/Coordinate-System
 Author: PanGuoJun
@@ -48,7 +48,7 @@ License-File: LICENSE
 [![License](https://img.shields.io/badge/license-MIT-green.svg)](LICENSE)
 
 **Author:** PanGuoJun
-**Version:** 2.4.0
+**Version:** 2.4.1
 **License:** MIT
 
 ## 🆕 What's New in v2.4.0
@@ -99,154 +99,6 @@ Perfect for computational geometry, geometric analysis, and discrete differentia
 
 ---
 
-## 🆕 Differential Geometry (NEW in v2.2.0)
-
-### Surface Representation
-
-```python
-from coordinate_system import Sphere, Torus, Surface
-
-# Create a sphere
-sphere = Sphere(radius=2.0)
-
-# Create a torus
-torus = Torus(major_radius=3.0, minor_radius=1.0)
-
-# Or define your own surface by subclassing Surface
-class MySurface(Surface):
-    def position(self, u, v):
-        # Return vec3(x, y, z) for parameters (u, v)
-        pass
-```
-
-### Metric Tensor (First Fundamental Form)
-
-```python
-from coordinate_system import compute_metric
-import math
-
-# Compute metric tensor at a point
-g = compute_metric(sphere, u=math.pi/4, v=math.pi/3)
-
-print(f"E = {g.E:.6f}")   # u-direction metric
-print(f"F = {g.F:.6f}")   # cross term (0 if orthogonal)
-print(f"G = {g.G:.6f}")   # v-direction metric
-print(f"det(g) = {g.det:.6f}")
-print(f"Metric correction = {g.correction_factor():.6f}")
-```
-
-### Connection Operator (Intrinsic Gradient)
-
-```python
-from coordinate_system import compute_connection
-
-# Compute connection operator (frame derivative)
-G_u = compute_connection(sphere, u=math.pi/4, v=math.pi/3, direction='u')
-G_v = compute_connection(sphere, u=math.pi/4, v=math.pi/3, direction='v')
-
-print(f"G_u norm: {G_u.norm():.8f}")
-print(f"G_v norm: {G_v.norm():.8f}")
-```
-
-**Aliases available:**
-- `compute_frame_derivative` (standard mathematical term)
-- `compute_intrinsic_gradient` (intuitive name)
-- `compute_geometric_gradient` (alternative)
-
-### Curvature Tensor
-
-```python
-from coordinate_system import compute_curvature_tensor
-
-# Compute complete curvature tensor R_uv
-R_uv = compute_curvature_tensor(
-    sphere,
-    u=math.pi/4,
-    v=math.pi/3,
-    use_lie_derivative=True  # CRITICAL for accuracy!
-)
-
-print(f"Curvature tensor norm: {R_uv.norm():.8f}")
-
-# R_uv is a 3×3 coord3 object with structure:
-# [[R_11, R_12, R_13],   ← Tangent-tangent (intrinsic)
-#  [R_21, R_22, R_23],   ← Tangent-normal (second fundamental form)
-#  [R_31, R_32, R_33]]   ← Normal-normal (extrinsic)
-```
-
-### Gaussian Curvature (Quickest Method)
-
-```python
-from coordinate_system import Sphere, compute_gaussian_curvature
-import math
-
-# One-line curvature computation!
-sphere = Sphere(radius=2.0)
-K = compute_gaussian_curvature(sphere, u=math.pi/4, v=math.pi/3)
-
-print(f"Gaussian curvature K = {K:.6f}")
-# Expected: K = 1/R² = 0.25 for a sphere
-
-# Theoretical value
-K_theory = 1.0 / (2.0 ** 2)
-error = abs(K - K_theory) / K_theory * 100
-print(f"Relative error: {error:.2f}%")  # Typically < 3% !
-```
-
-### Accuracy Comparison
-
-**On a sphere with radius R=2 (24 test points across different angles):**
-
-| Method | Antisymmetric Extraction | Accuracy | Notes |
-|--------|-------------------------|----------|-------|
-| Old (R_10 only) | ❌ | ~2.5% error | Previous versions |
-| **v2.4.0 (Corrected)** | ✓ | **0.000% error** | **Machine precision!** 🎉 |
-
-**Key Formula:** `K = (R_01 - R_10) / (2 · det(g))`
-
-The antisymmetric part `(R_01 - R_10)/2` correctly eliminates metric influence and extracts pure curvature. This formula is verified against traditional Christoffel symbol methods and achieves identical results at machine precision (< 10⁻¹⁵ error).
-
-### Key Parameters
-
-```python
-# For constant curvature surfaces (spheres):
-K = compute_gaussian_curvature(
-    surface,
-    u, v,
-    scale_factor=2.0,           # Optimal for spheres
-    use_metric_correction=True,  # Essential!
-    use_lie_derivative=True      # Critical for accuracy!
-)
-
-# These are the default values, optimized for best accuracy
-```
-
-### Complete Example
-
-See `examples/curvature_computation.py` for comprehensive demonstrations, or `examples/quickstart.py` for the simplest possible usage.
-
-```python
-# examples/quickstart.py
-import math
-from coordinate_system import Sphere, compute_gaussian_curvature
-
-# Step 1: Create a sphere
-sphere = Sphere(radius=2.0)
-
-# Step 2: Choose a point
-u, v = math.pi/4, math.pi/3
-
-# Step 3: Compute curvature (one line!)
-K = compute_gaussian_curvature(sphere, u, v)
-
-# Step 4: Compare with theory
-K_theory = 1.0 / (2.0 ** 2)  # K = 1/R²
-print(f"Computed: K = {K:.6f}")
-print(f"Theory:   K = {K_theory:.6f}")
-print(f"Error: {abs(K-K_theory)/K_theory*100:.2f}%")
-```
-
----
 
 ## 📚 Documentation
 

{coordinate_system-2.4.0 → coordinate_system-2.4.2}/README.md

@@ -8,7 +8,7 @@
 [![License](https://img.shields.io/badge/license-MIT-green.svg)](LICENSE)
 
 **Author:** PanGuoJun
-**Version:** 2.4.0
+**Version:** 2.4.1
 **License:** MIT
 
 ## 🆕 What's New in v2.4.0
@@ -59,154 +59,6 @@ Perfect for computational geometry, geometric analysis, and discrete differentia
 
 ---
 
-## 🆕 Differential Geometry (NEW in v2.2.0)
-
-### Surface Representation
-
-```python
-from coordinate_system import Sphere, Torus, Surface
-
-# Create a sphere
-sphere = Sphere(radius=2.0)
-
-# Create a torus
-torus = Torus(major_radius=3.0, minor_radius=1.0)
-
-# Or define your own surface by subclassing Surface
-class MySurface(Surface):
-    def position(self, u, v):
-        # Return vec3(x, y, z) for parameters (u, v)
-        pass
-```
-
-### Metric Tensor (First Fundamental Form)
-
-```python
-from coordinate_system import compute_metric
-import math
-
-# Compute metric tensor at a point
-g = compute_metric(sphere, u=math.pi/4, v=math.pi/3)
-
-print(f"E = {g.E:.6f}")   # u-direction metric
-print(f"F = {g.F:.6f}")   # cross term (0 if orthogonal)
-print(f"G = {g.G:.6f}")   # v-direction metric
-print(f"det(g) = {g.det:.6f}")
-print(f"Metric correction = {g.correction_factor():.6f}")
-```
-
-### Connection Operator (Intrinsic Gradient)
-
-```python
-from coordinate_system import compute_connection
-
-# Compute connection operator (frame derivative)
-G_u = compute_connection(sphere, u=math.pi/4, v=math.pi/3, direction='u')
-G_v = compute_connection(sphere, u=math.pi/4, v=math.pi/3, direction='v')
-
-print(f"G_u norm: {G_u.norm():.8f}")
-print(f"G_v norm: {G_v.norm():.8f}")
-```
-
-**Aliases available:**
-- `compute_frame_derivative` (standard mathematical term)
-- `compute_intrinsic_gradient` (intuitive name)
-- `compute_geometric_gradient` (alternative)
-
-### Curvature Tensor
-
-```python
-from coordinate_system import compute_curvature_tensor
-
-# Compute complete curvature tensor R_uv
-R_uv = compute_curvature_tensor(
-    sphere,
-    u=math.pi/4,
-    v=math.pi/3,
-    use_lie_derivative=True  # CRITICAL for accuracy!
-)
-
-print(f"Curvature tensor norm: {R_uv.norm():.8f}")
-
-# R_uv is a 3×3 coord3 object with structure:
-# [[R_11, R_12, R_13],   ← Tangent-tangent (intrinsic)
-#  [R_21, R_22, R_23],   ← Tangent-normal (second fundamental form)
-#  [R_31, R_32, R_33]]   ← Normal-normal (extrinsic)
-```
-
-### Gaussian Curvature (Quickest Method)
-
-```python
-from coordinate_system import Sphere, compute_gaussian_curvature
-import math
-
-# One-line curvature computation!
-sphere = Sphere(radius=2.0)
-K = compute_gaussian_curvature(sphere, u=math.pi/4, v=math.pi/3)
-
-print(f"Gaussian curvature K = {K:.6f}")
-# Expected: K = 1/R² = 0.25 for a sphere
-
-# Theoretical value
-K_theory = 1.0 / (2.0 ** 2)
-error = abs(K - K_theory) / K_theory * 100
-print(f"Relative error: {error:.2f}%")  # Typically < 3% !
-```
-
-### Accuracy Comparison
-
-**On a sphere with radius R=2 (24 test points across different angles):**
-
-| Method | Antisymmetric Extraction | Accuracy | Notes |
-|--------|-------------------------|----------|-------|
-| Old (R_10 only) | ❌ | ~2.5% error | Previous versions |
-| **v2.4.0 (Corrected)** | ✓ | **0.000% error** | **Machine precision!** 🎉 |
-
-**Key Formula:** `K = (R_01 - R_10) / (2 · det(g))`
-
-The antisymmetric part `(R_01 - R_10)/2` correctly eliminates metric influence and extracts pure curvature. This formula is verified against traditional Christoffel symbol methods and achieves identical results at machine precision (< 10⁻¹⁵ error).
-
-### Key Parameters
-
-```python
-# For constant curvature surfaces (spheres):
-K = compute_gaussian_curvature(
-    surface,
-    u, v,
-    scale_factor=2.0,           # Optimal for spheres
-    use_metric_correction=True,  # Essential!
-    use_lie_derivative=True      # Critical for accuracy!
-)
-
-# These are the default values, optimized for best accuracy
-```
-
-### Complete Example
-
-See `examples/curvature_computation.py` for comprehensive demonstrations, or `examples/quickstart.py` for the simplest possible usage.
-
-```python
-# examples/quickstart.py
-import math
-from coordinate_system import Sphere, compute_gaussian_curvature
-
-# Step 1: Create a sphere
-sphere = Sphere(radius=2.0)
-
-# Step 2: Choose a point
-u, v = math.pi/4, math.pi/3
-
-# Step 3: Compute curvature (one line!)
-K = compute_gaussian_curvature(sphere, u, v)
-
-# Step 4: Compare with theory
-K_theory = 1.0 / (2.0 ** 2)  # K = 1/R²
-print(f"Computed: K = {K:.6f}")
-print(f"Theory:   K = {K_theory:.6f}")
-print(f"Error: {abs(K-K_theory)/K_theory*100:.2f}%")
-```
-
----
 
 ## 📚 Documentation
 

{coordinate_system-2.4.0 → coordinate_system-2.4.2/coordinate_system.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: coordinate_system
-Version: 2.4.0
+Version: 2.4.2
 Summary: High-performance 3D coordinate system library with machine-precision Gaussian curvature computation (0.000% error) using corrected antisymmetric extraction formula
 Home-page: https://github.com/panguojun/Coordinate-System
 Author: PanGuoJun
@@ -48,7 +48,7 @@ License-File: LICENSE
 [![License](https://img.shields.io/badge/license-MIT-green.svg)](LICENSE)
 
 **Author:** PanGuoJun
-**Version:** 2.4.0
+**Version:** 2.4.1
 **License:** MIT
 
 ## 🆕 What's New in v2.4.0
@@ -99,154 +99,6 @@ Perfect for computational geometry, geometric analysis, and discrete differentia
 
 ---
 
-## 🆕 Differential Geometry (NEW in v2.2.0)
-
-### Surface Representation
-
-```python
-from coordinate_system import Sphere, Torus, Surface
-
-# Create a sphere
-sphere = Sphere(radius=2.0)
-
-# Create a torus
-torus = Torus(major_radius=3.0, minor_radius=1.0)
-
-# Or define your own surface by subclassing Surface
-class MySurface(Surface):
-    def position(self, u, v):
-        # Return vec3(x, y, z) for parameters (u, v)
-        pass
-```
-
-### Metric Tensor (First Fundamental Form)
-
-```python
-from coordinate_system import compute_metric
-import math
-
-# Compute metric tensor at a point
-g = compute_metric(sphere, u=math.pi/4, v=math.pi/3)
-
-print(f"E = {g.E:.6f}")   # u-direction metric
-print(f"F = {g.F:.6f}")   # cross term (0 if orthogonal)
-print(f"G = {g.G:.6f}")   # v-direction metric
-print(f"det(g) = {g.det:.6f}")
-print(f"Metric correction = {g.correction_factor():.6f}")
-```
-
-### Connection Operator (Intrinsic Gradient)
-
-```python
-from coordinate_system import compute_connection
-
-# Compute connection operator (frame derivative)
-G_u = compute_connection(sphere, u=math.pi/4, v=math.pi/3, direction='u')
-G_v = compute_connection(sphere, u=math.pi/4, v=math.pi/3, direction='v')
-
-print(f"G_u norm: {G_u.norm():.8f}")
-print(f"G_v norm: {G_v.norm():.8f}")
-```
-
-**Aliases available:**
-- `compute_frame_derivative` (standard mathematical term)
-- `compute_intrinsic_gradient` (intuitive name)
-- `compute_geometric_gradient` (alternative)
-
-### Curvature Tensor
-
-```python
-from coordinate_system import compute_curvature_tensor
-
-# Compute complete curvature tensor R_uv
-R_uv = compute_curvature_tensor(
-    sphere,
-    u=math.pi/4,
-    v=math.pi/3,
-    use_lie_derivative=True  # CRITICAL for accuracy!
-)
-
-print(f"Curvature tensor norm: {R_uv.norm():.8f}")
-
-# R_uv is a 3×3 coord3 object with structure:
-# [[R_11, R_12, R_13],   ← Tangent-tangent (intrinsic)
-#  [R_21, R_22, R_23],   ← Tangent-normal (second fundamental form)
-#  [R_31, R_32, R_33]]   ← Normal-normal (extrinsic)
-```
-
-### Gaussian Curvature (Quickest Method)
-
-```python
-from coordinate_system import Sphere, compute_gaussian_curvature
-import math
-
-# One-line curvature computation!
-sphere = Sphere(radius=2.0)
-K = compute_gaussian_curvature(sphere, u=math.pi/4, v=math.pi/3)
-
-print(f"Gaussian curvature K = {K:.6f}")
-# Expected: K = 1/R² = 0.25 for a sphere
-
-# Theoretical value
-K_theory = 1.0 / (2.0 ** 2)
-error = abs(K - K_theory) / K_theory * 100
-print(f"Relative error: {error:.2f}%")  # Typically < 3% !
-```
-
-### Accuracy Comparison
-
-**On a sphere with radius R=2 (24 test points across different angles):**
-
-| Method | Antisymmetric Extraction | Accuracy | Notes |
-|--------|-------------------------|----------|-------|
-| Old (R_10 only) | ❌ | ~2.5% error | Previous versions |
-| **v2.4.0 (Corrected)** | ✓ | **0.000% error** | **Machine precision!** 🎉 |
-
-**Key Formula:** `K = (R_01 - R_10) / (2 · det(g))`
-
-The antisymmetric part `(R_01 - R_10)/2` correctly eliminates metric influence and extracts pure curvature. This formula is verified against traditional Christoffel symbol methods and achieves identical results at machine precision (< 10⁻¹⁵ error).
-
-### Key Parameters
-
-```python
-# For constant curvature surfaces (spheres):
-K = compute_gaussian_curvature(
-    surface,
-    u, v,
-    scale_factor=2.0,           # Optimal for spheres
-    use_metric_correction=True,  # Essential!
-    use_lie_derivative=True      # Critical for accuracy!
-)
-
-# These are the default values, optimized for best accuracy
-```
-
-### Complete Example
-
-See `examples/curvature_computation.py` for comprehensive demonstrations, or `examples/quickstart.py` for the simplest possible usage.
-
-```python
-# examples/quickstart.py
-import math
-from coordinate_system import Sphere, compute_gaussian_curvature
-
-# Step 1: Create a sphere
-sphere = Sphere(radius=2.0)
-
-# Step 2: Choose a point
-u, v = math.pi/4, math.pi/3
-
-# Step 3: Compute curvature (one line!)
-K = compute_gaussian_curvature(sphere, u, v)
-
-# Step 4: Compare with theory
-K_theory = 1.0 / (2.0 ** 2)  # K = 1/R²
-print(f"Computed: K = {K:.6f}")
-print(f"Theory:   K = {K_theory:.6f}")
-print(f"Error: {abs(K-K_theory)/K_theory*100:.2f}%")
-```
-
----
 
 ## 📚 Documentation
 

{coordinate_system-2.4.0 → coordinate_system-2.4.2}/coordinate_system_binding.cpp

@@ -397,51 +397,27 @@ PYBIND11_MODULE(coordinate_system, m) {
             return coord3(vec3::ZERO, q);
         }, "Create coordinate system at origin with rotation", py::arg("rotation"))
 
-        // Differential Geometry - Connection Operator (Intrinsic Gradient)
-        .def("compute_metric_det", &coord3::compute_metric_det,
+        // Differential Geometry
+		.def("metric", &coord3::metric,
             "Compute metric tensor determinant: det(g) = E*G - F²\n\n"
             "For a parametric surface, the first fundamental form is g = [[E, F], [F, G]].\n"
-            "Returns det(g) = E*G - F² which appears in Gaussian curvature formula.")
-        .def_static("connection_operator_base", &coord3::connection_operator_base,
-            "Compute base connection operator: G' = (c2*c1^-1)/C2 - I/C1\n\n"
-            "This is the fundamental discrete connection operator formula.\n"
-            "G' represents the discrete covariant derivative of the frame field.\n\n"
-            "Args:\n"
-            "    c1: coordinate frame at point (u, v)\n"
-            "    c2: coordinate frame at point (u+h, v) or (u, v+h)\n"
-            "    C1: scale factor at c1\n"
-            "    C2: scale factor at c2\n\n"
-            "Returns:\n"
-            "    G' (base connection operator as coord3)",
-            py::arg("c1"), py::arg("c2"), py::arg("C1"), py::arg("C2"))
-        .def_static("connection_operator", &coord3::connection_operator,
-            "Compute connection operator with corrections: G = f(G')\n\n"
-            "Applies scale correction and optional metric correction to G'.\n\n"
-            "Args:\n"
-            "    c1: coordinate frame at point (u, v)\n"
-            "    c2: coordinate frame at point (u+h, v) or (u, v+h)\n"
-            "    C1: scale factor at c1\n"
-            "    C2: scale factor at c2\n"
-            "    scale_factor: overall scale multiplier (optimal: 2.0 for spheres)\n"
-            "    use_metric_correction: if true, applies 1/√det(g) correction\n\n"
-            "Returns:\n"
-            "    G (corrected connection operator)\n\n"
-            "For constant curvature surfaces (spheres), optimal parameters are:\n"
-            "  - scale_factor = 2.0\n"
-            "  - use_metric_correction = true\n"
-            "  - This achieves ~2.22% accuracy in Gaussian curvature computation",
-            py::arg("c1"), py::arg("c2"), py::arg("C1"), py::arg("C2"),
-            py::arg("scale_factor") = 2.0, py::arg("use_metric_correction") = true)
+            "Returns det(g) = E*G - F² which appears in Gaussian curvature formula.")  
+			
+        .def("metric_det", &coord3::metric_det,
+            "Compute metric tensor determinant: det(g) = E*G - F²\n\n"
+            "For a parametric surface, the first fundamental form is g = [[E, F], [F, G]].\n"
+            "Returns det(g) = E*G - F² which appears in Gaussian curvature formula.")        
+        
         .def_static("lie_bracket", &coord3::lie_bracket,
-            "Compute Lie bracket [G_u, G_v] = G_u * G_v - G_v * G_u\n\n"
+            "Compute Lie bracket [A, B] = A * B - B * A\n\n"
             "The Lie bracket measures the non-commutativity of the connection operators.\n"
             "It appears in the curvature tensor formula as the Lie derivative term.\n\n"
             "Args:\n"
-            "    G_u: connection operator in u-direction\n"
-            "    G_v: connection operator in v-direction\n\n"
+            "    A: connection operator in u-direction\n"
+            "    B: connection operator in v-direction\n\n"
             "Returns:\n"
-            "    [G_u, G_v] (Lie bracket as coord3)",
-            py::arg("G_u"), py::arg("G_v"));
+            "    [A, B] (Lie bracket as coord3)",
+            py::arg("A"), py::arg("B"));
 
     // ===========================================================================
     // Utility Functions

{coordinate_system-2.4.0 → coordinate_system-2.4.2}/pmsys_minimal.hpp

@@ -3762,15 +3762,9 @@ struct  coord3 : vcoord3
 	}
 
 	// ============================================================================
-	// 几何运算：梯度、插值、距离等
+	// 几何运算：插值、距离等
 	// ============================================================================
-	// 梯度坐标系：用于计算联络
-	// V2 - V1 = V1 * G, 其中 G = C2 / C1 - I
-	static coord3 grad(const coord3& c1, const coord3& c2)
-	{
-		return c2 / c1 - ONE;
-	}
-
+	
 	// 线性插值：位置和旋转的独立插值
 	static coord3 lerp(const coord3& c1, const coord3& c2, real t)
 	{
@@ -3802,6 +3796,48 @@ struct  coord3 : vcoord3
 	{
 		return Q().angle_to(other.Q());
 	}
+
+	// ============================================================================
+	// 曲率计算
+	// ============================================================================
+	// 梯度坐标系
+	// V2 - V1 = V1 * G, 其中 G = C2 / C1 - I
+	static coord3 grad(const coord3& c1, const coord3& c2)
+	{
+		return c2 / c1 - ONE;
+	}
+	
+	vec3 metric() const
+	{
+		// 对于曲面，我们通常只关心 u,v 方向的切向量
+		vec3 vx = VX();  // ∂r/∂u 方向
+		vec3 vy = VY();  // ∂r/∂v 方向
+
+		real E = vx.dot(vx);  // g_uu = ∂r/∂u · ∂r/∂u
+		real F = vx.dot(vy);  // g_uv = ∂r/∂u · ∂r/∂v  
+		real G = vy.dot(vy);  // g_vv = ∂r/∂v · ∂r/∂v
+
+		return vec3(E, F, G);
+	}
+	// 计算第一基本形式的行列式
+	real metric_det() const 
+	{
+		vec3 vx = VX();
+		vec3 vy = VY();
+		vec3 vz = VZ();  // 这个vz不应该参与2D曲面计算！
+
+		real E = vx.dot(vx);
+		real F = vx.dot(vy);
+		real G = vy.dot(vy);
+
+		return E * G - F * F;
+	}
+
+	static coord3 lie_bracket(const coord3& A, const coord3& B) {
+		// [A, B] = A * B - B * A
+		return A * B - B * A;  // 需要正确定义这个运算
+	}
+	
 	std::string serialise() const
 	{
 		vec3 eu = coord2eulers();
@@ -3847,134 +3883,7 @@ struct  coord3 : vcoord3
 		}
 		return hash;
 	}
-
-	// ============================================================================
-	// Differential Geometry: Connection Operator (Intrinsic Gradient)
-	// ============================================================================
-
-	/**
-	 * Compute metric tensor determinant: det(g) = E*G - F²
-	 *
-	 * For a parametric surface, the first fundamental form is:
-	 *   g = [[E, F], [F, G]]
-	 * where:
-	 *   E = VX · VX (squared length of X-axis basis vector)
-	 *   F = VX · VY (dot product of X and Y basis vectors)
-	 *   G = VY · VY (squared length of Y-axis basis vector)
-	 *
-	 * Returns: det(g) = E*G - F²
-	 */
-	real compute_metric_det() const
-	{
-		vec3 vx = VX();
-		vec3 vy = VY();
-		real E = vx.dot(vx);
-		real F = vx.dot(vy);
-		real G = vy.dot(vy);
-		return E * G - F * F;
-	}
-
-	/**
-	 * Compute base connection operator: G' = (c2*c1^-1)/C2 - I/C1
-	 *
-	 * This is the fundamental discrete connection operator formula.
-	 *
-	 * Parameters:
-	 *   c1 - coordinate frame at point (u, v)
-	 *   c2 - coordinate frame at point (u+h, v) or (u, v+h)
-	 *   C1 - scale factor at c1 (typically sqrt(det(g1)) or c1.s.mean())
-	 *   C2 - scale factor at c2 (typically sqrt(det(g2)) or c2.s.mean())
-	 *
-	 * Returns: G' (base connection operator as coord3)
-	 *
-	 * Mathematical interpretation:
-	 *   G' represents the discrete covariant derivative of the frame field.
-	 *   It measures how the coordinate frame changes as we move along the surface.
-	 */
-	static coord3 connection_operator_base(
-		const coord3& c1,
-		const coord3& c2,
-		real C1,
-		real C2
-	)
-	{
-		// G' = (c2*c1^-1)/C2 - I/C1
-		// where c1^-1 is computed as c2/c1
-		coord3 identity; // Identity coordinate frame
-		coord3 term1 = (c2 / c1) / C2;
-		coord3 term2 = identity / C1;
-		return term1 - term2;
-	}
-
-	/**
-	 * Compute connection operator with corrections: G = f(G')
-	 *
-	 * Applies scale correction and optional metric correction to G'.
-	 *
-	 * Parameters:
-	 *   c1, c2 - coordinate frames (see connection_operator_base)
-	 *   C1, C2 - scale factors
-	 *   scale_factor - overall scale multiplier (optimal: 2.0 for spheres)
-	 *   use_metric_correction - if true, applies 1/√det(g) correction
-	 *
-	 * Returns: G (corrected connection operator)
-	 *
-	 * Formula:
-	 *   If use_metric_correction:
-	 *     G = scale_factor * G' / √det(g)
-	 *   Else:
-	 *     G = scale_factor * G'
-	 *
-	 * For constant curvature surfaces (spheres), optimal parameters are:
-	 *   - scale_factor = 2.0
-	 *   - use_metric_correction = true
-	 *   - This achieves ~2.22% accuracy in Gaussian curvature computation
-	 */
-	static coord3 connection_operator(
-		const coord3& c1,
-		const coord3& c2,
-		real C1,
-		real C2,
-		real scale_factor = 2.0,
-		bool use_metric_correction = true
-	)
-	{
-		// Compute base operator
-		coord3 G_prime = connection_operator_base(c1, c2, C1, C2);
-
-		// Apply scale correction
-		coord3 G = G_prime * scale_factor;
-
-		// Apply metric correction if requested
-		if (use_metric_correction)
-		{
-			real det_g = c1.compute_metric_det();
-			if (det_g > 1e-10)
-			{
-				real metric_factor = 1.0 / sqrt(det_g);
-				G = G * metric_factor;
-			}
-		}
-
-		return G;
-	}
-
-	/**
-	 * Compute Lie bracket [G_u, G_v] = G_u * G_v - G_v * G_u
-	 *
-	 * The Lie bracket measures the non-commutativity of the connection operators.
-	 * It appears in the curvature tensor formula as the Lie derivative term.
-	 *
-	 * Parameters:
-	 *   G_u - connection operator in u-direction
-	 *   G_v - connection operator in v-direction
-	 *
-	 * Returns: [G_u, G_v] = G_u * G_v - G_v * G_u
-	 */
-	static coord3 lie_bracket(const coord3& G_u, const coord3& G_v)
-	{
-		return G_u * G_v - G_v * G_u;
-	}
+	
 };
 
 inline const coord3 coord3::ZERO = {ucoord3::ONE, vec3::ZERO, vec3::ZERO };

{coordinate_system-2.4.0 → coordinate_system-2.4.2}/setup.py

@@ -2,7 +2,7 @@
 setup.py - Cross-platform setup for coordinate_system package
 
 Author: PanGuoJun
-Version: 2.4.0
+Version: 2.4.2
 License: MIT
 """
 
@@ -69,7 +69,7 @@ ext_modules = [
 
 setup(
     name='coordinate_system',
-    version='2.4.0',
+    version='2.4.2',
     packages=find_packages(),
     ext_modules=ext_modules,  # Add extension modules