coordinate-system 4.0.0__tar.gz → 4.0.2__tar.gz

{coordinate_system-4.0.0/coordinate_system.egg-info → coordinate_system-4.0.2}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: coordinate_system
-Version: 4.0.0
+Version: 4.0.2
 Summary: High-performance 3D coordinate system library with right-handed conventions, Intrinsic Gradient Operator method for curvature computation, achieving machine-precision accuracy
 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:** 3.0.0
+**Version:** 4.0.2
 **License:** MIT
 
 ## 🆕 What's New in v4.0.0
@@ -71,24 +71,24 @@ License-File: LICENSE
 
 ## 🆕 What's New in v3.0.0
 
-**Critical Coordinate System Update - Right-Hand System!**
+**Critical Coordinate System Update - Left-Hand System!**
 
-- ⚠️ **BREAKING CHANGE: Right-Hand Coordinate System** - Unified all operations to standard right-hand coordinate system
-  - `vec3.cross()` now uses right-hand rule (standard mathematical definition)
-  - Quaternion rotations follow right-hand convention (counter-clockwise is positive)
+- ⚠️ **BREAKING CHANGE: Left-Hand Coordinate System** - Unified all operations to standard left-hand coordinate system
+  - `vec3.cross()` now uses left-hand rule (common in computer graphics)
+  - Quaternion rotations follow left-hand convention (clockwise is positive)
   - `coord3.look_at()` and `coord3.from_forward()` updated accordingly
   - **Migration Guide:** See [HANDEDNESS_CHANGE.md](HANDEDNESS_CHANGE.md) for details
-  - `cross_left()` method retained for backward compatibility
-- ✅ **Verified Correctness** - Complete test suite validates right-hand system behavior
+  - `cross_right()` method retained for backward compatibility
+- ✅ **Verified Correctness** - Complete test suite validates left-hand system behavior
 - 📚 **Updated Documentation** - All examples and documentation reflect new coordinate system
 - 🎯 **Enhanced Clarity** - Comments and docstrings clarified throughout
 
 **Why This Change:**
-1. Mathematical Standard: Right-hand system is the standard in mathematics and physics
-2. Better Interoperability: Compatible with OpenGL, NumPy, SciPy, and most scientific libraries
-3. Improved Readability: Standard cross product formulas match textbook definitions
+1. Graphics Standard: Left-hand system is common in computer graphics and DirectX
+2. Better Interoperability: Compatible with DirectX, Unity, and many game engines
+3. Improved Readability: Standard cross product formulas match graphics pipeline definitions
 
-Perfect for computational geometry, geometric analysis, and discrete differential geometry research!
+Perfect for computer graphics, game development, and 3D visualization applications!
 
 ---
 
@@ -210,7 +210,7 @@ lerped = vec3.lerp(v1, v2, 0.5)  # Linear interpolation
 
 ## Coordinate System Type
 
-This library uses a **right-handed coordinate system** for all vector and quaternion operations (as of v3.0.0).
+This library uses a **left-handed coordinate system** for all vector and quaternion operations (as of v3.0.0).
 
 ```
          +Y
@@ -218,18 +218,18 @@ This library uses a **right-handed coordinate system** for all vector and quater
           |
           |
           +-----> +X
-         /
-        /
-      +Z
+          \
+           \
+            +Z
 ```
 
-**Right-Hand Rule:**
-- Point your right hand's fingers along the +X axis
+**Left-Hand Rule:**
+- Point your left hand's fingers along the +X axis
 - Curl them toward the +Y axis
 - Your thumb points along the +Z axis
-- Cross product: **X × Y = Z** (standard mathematical definition)
+- Cross product: **X × Y = Z** (common in computer graphics)
 
-**Migration from v2.x:** If your code relied on the previous left-handed system, see [HANDEDNESS_CHANGE.md](HANDEDNESS_CHANGE.md) for migration instructions.
+**Migration from v2.x:** If your code relied on the previous right-handed system, see [HANDEDNESS_CHANGE.md](HANDEDNESS_CHANGE.md) for migration instructions.
 
 ---
 
@@ -267,7 +267,7 @@ v3 = v1 * v2           # Component-wise multiplication
 
 ```python
 dot = v1.dot(v2)       # Dot product (float)
-cross = v1.cross(v2)   # Cross product (vec3)
+cross = v1.cross(v2)   # Cross product (vec3) - Left-hand rule
 length = v.length()    # Vector length
 v.normalize()          # Normalize in-place
 normalized = v.normcopy()  # Return normalized copy
@@ -451,7 +451,7 @@ v4 = v1 * 5                     # vec3(5, 0, 0)
 
 # Dot and cross products
 dot = v1.dot(v2)                # 0.0 (perpendicular)
-cross = v1.cross(v2)            # vec3(0, 0, 1) in right-handed system (X × Y = Z)
+cross = v1.cross(v2)            # vec3(0, 0, -1) in left-handed system (X × Y = -Z)
 
 # Length and normalization
 length = v1.length()            # 1.0
@@ -516,8 +516,8 @@ child_in_world = child * parent
 # Create look-at transformation (custom implementation needed)
 def look_at(eye, target, up=vec3(0, 1, 0)):
     forward = (target - eye).normcopy()
-    right = up.cross(forward).normcopy()
-    up_corrected = forward.cross(right)
+    right = forward.cross(up).normcopy()  # Left-hand cross product
+    up_corrected = right.cross(forward)
     return coord3.from_axes(right, up_corrected, forward)
 
 camera = look_at(vec3(10, 10, 10), vec3(0, 0, 0))
@@ -537,8 +537,8 @@ class Camera:
 
     def create_look_at(self, eye, target, up):
         forward = (target - eye).normcopy()
-        right = up.cross(forward).normcopy()
-        up_corrected = forward.cross(right)
+        right = forward.cross(up).normcopy()  # Left-hand cross product
+        up_corrected = right.cross(forward)
 
         c = coord3()
         c.o = eye

{coordinate_system-4.0.0 → coordinate_system-4.0.2}/README.md

@@ -8,7 +8,7 @@
 [![License](https://img.shields.io/badge/license-MIT-green.svg)](LICENSE)
 
 **Author:** PanGuoJun
-**Version:** 3.0.0
+**Version:** 4.0.2
 **License:** MIT
 
 ## 🆕 What's New in v4.0.0
@@ -31,24 +31,24 @@
 
 ## 🆕 What's New in v3.0.0
 
-**Critical Coordinate System Update - Right-Hand System!**
+**Critical Coordinate System Update - Left-Hand System!**
 
-- ⚠️ **BREAKING CHANGE: Right-Hand Coordinate System** - Unified all operations to standard right-hand coordinate system
-  - `vec3.cross()` now uses right-hand rule (standard mathematical definition)
-  - Quaternion rotations follow right-hand convention (counter-clockwise is positive)
+- ⚠️ **BREAKING CHANGE: Left-Hand Coordinate System** - Unified all operations to standard left-hand coordinate system
+  - `vec3.cross()` now uses left-hand rule (common in computer graphics)
+  - Quaternion rotations follow left-hand convention (clockwise is positive)
   - `coord3.look_at()` and `coord3.from_forward()` updated accordingly
   - **Migration Guide:** See [HANDEDNESS_CHANGE.md](HANDEDNESS_CHANGE.md) for details
-  - `cross_left()` method retained for backward compatibility
-- ✅ **Verified Correctness** - Complete test suite validates right-hand system behavior
+  - `cross_right()` method retained for backward compatibility
+- ✅ **Verified Correctness** - Complete test suite validates left-hand system behavior
 - 📚 **Updated Documentation** - All examples and documentation reflect new coordinate system
 - 🎯 **Enhanced Clarity** - Comments and docstrings clarified throughout
 
 **Why This Change:**
-1. Mathematical Standard: Right-hand system is the standard in mathematics and physics
-2. Better Interoperability: Compatible with OpenGL, NumPy, SciPy, and most scientific libraries
-3. Improved Readability: Standard cross product formulas match textbook definitions
+1. Graphics Standard: Left-hand system is common in computer graphics and DirectX
+2. Better Interoperability: Compatible with DirectX, Unity, and many game engines
+3. Improved Readability: Standard cross product formulas match graphics pipeline definitions
 
-Perfect for computational geometry, geometric analysis, and discrete differential geometry research!
+Perfect for computer graphics, game development, and 3D visualization applications!
 
 ---
 
@@ -170,7 +170,7 @@ lerped = vec3.lerp(v1, v2, 0.5)  # Linear interpolation
 
 ## Coordinate System Type
 
-This library uses a **right-handed coordinate system** for all vector and quaternion operations (as of v3.0.0).
+This library uses a **left-handed coordinate system** for all vector and quaternion operations (as of v3.0.0).
 
 ```
          +Y
@@ -178,18 +178,18 @@ This library uses a **right-handed coordinate system** for all vector and quater
           |
           |
           +-----> +X
-         /
-        /
-      +Z
+          \
+           \
+            +Z
 ```
 
-**Right-Hand Rule:**
-- Point your right hand's fingers along the +X axis
+**Left-Hand Rule:**
+- Point your left hand's fingers along the +X axis
 - Curl them toward the +Y axis
 - Your thumb points along the +Z axis
-- Cross product: **X × Y = Z** (standard mathematical definition)
+- Cross product: **X × Y = Z** (common in computer graphics)
 
-**Migration from v2.x:** If your code relied on the previous left-handed system, see [HANDEDNESS_CHANGE.md](HANDEDNESS_CHANGE.md) for migration instructions.
+**Migration from v2.x:** If your code relied on the previous right-handed system, see [HANDEDNESS_CHANGE.md](HANDEDNESS_CHANGE.md) for migration instructions.
 
 ---
 
@@ -227,7 +227,7 @@ v3 = v1 * v2           # Component-wise multiplication
 
 ```python
 dot = v1.dot(v2)       # Dot product (float)
-cross = v1.cross(v2)   # Cross product (vec3)
+cross = v1.cross(v2)   # Cross product (vec3) - Left-hand rule
 length = v.length()    # Vector length
 v.normalize()          # Normalize in-place
 normalized = v.normcopy()  # Return normalized copy
@@ -411,7 +411,7 @@ v4 = v1 * 5                     # vec3(5, 0, 0)
 
 # Dot and cross products
 dot = v1.dot(v2)                # 0.0 (perpendicular)
-cross = v1.cross(v2)            # vec3(0, 0, 1) in right-handed system (X × Y = Z)
+cross = v1.cross(v2)            # vec3(0, 0, -1) in left-handed system (X × Y = -Z)
 
 # Length and normalization
 length = v1.length()            # 1.0
@@ -476,8 +476,8 @@ child_in_world = child * parent
 # Create look-at transformation (custom implementation needed)
 def look_at(eye, target, up=vec3(0, 1, 0)):
     forward = (target - eye).normcopy()
-    right = up.cross(forward).normcopy()
-    up_corrected = forward.cross(right)
+    right = forward.cross(up).normcopy()  # Left-hand cross product
+    up_corrected = right.cross(forward)
     return coord3.from_axes(right, up_corrected, forward)
 
 camera = look_at(vec3(10, 10, 10), vec3(0, 0, 0))
@@ -497,8 +497,8 @@ class Camera:
 
     def create_look_at(self, eye, target, up):
         forward = (target - eye).normcopy()
-        right = up.cross(forward).normcopy()
-        up_corrected = forward.cross(right)
+        right = forward.cross(up).normcopy()  # Left-hand cross product
+        up_corrected = right.cross(forward)
 
         c = coord3()
         c.o = eye
@@ -714,4 +714,4 @@ Built with ❤️ using:
 
 ---
 
-**Note**: For the latest updates and documentation, visit the [GitHub repository](https://github.com/panguojun/Coordinate-System).
+**Note**: For the latest updates and documentation, visit the [GitHub repository](https://github.com/panguojun/Coordinate-System).

{coordinate_system-4.0.0 → coordinate_system-4.0.2}/coordinate_system/differential_geometry.py

@@ -6,7 +6,7 @@ This module provides tools for discrete differential geometry computations on su
 using the CORRECT Intrinsic Gradient Operator framework based on the proven algorithm.
 
 Key Formula:
-    G_μ = (c(u+h) - c(u-h)) / (2h)  then extract normal derivative using .VZ()
+    G_μ = (c(u+h) - c(u-h)) / (2h) / c(u) then extract normal derivative using .VZ()
 
 Author: PanGuoJun
 Date: 2025-10-31

{coordinate_system-4.0.0 → coordinate_system-4.0.2/coordinate_system.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: coordinate_system
-Version: 4.0.0
+Version: 4.0.2
 Summary: High-performance 3D coordinate system library with right-handed conventions, Intrinsic Gradient Operator method for curvature computation, achieving machine-precision accuracy
 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:** 3.0.0
+**Version:** 4.0.2
 **License:** MIT
 
 ## 🆕 What's New in v4.0.0
@@ -71,24 +71,24 @@ License-File: LICENSE
 
 ## 🆕 What's New in v3.0.0
 
-**Critical Coordinate System Update - Right-Hand System!**
+**Critical Coordinate System Update - Left-Hand System!**
 
-- ⚠️ **BREAKING CHANGE: Right-Hand Coordinate System** - Unified all operations to standard right-hand coordinate system
-  - `vec3.cross()` now uses right-hand rule (standard mathematical definition)
-  - Quaternion rotations follow right-hand convention (counter-clockwise is positive)
+- ⚠️ **BREAKING CHANGE: Left-Hand Coordinate System** - Unified all operations to standard left-hand coordinate system
+  - `vec3.cross()` now uses left-hand rule (common in computer graphics)
+  - Quaternion rotations follow left-hand convention (clockwise is positive)
   - `coord3.look_at()` and `coord3.from_forward()` updated accordingly
   - **Migration Guide:** See [HANDEDNESS_CHANGE.md](HANDEDNESS_CHANGE.md) for details
-  - `cross_left()` method retained for backward compatibility
-- ✅ **Verified Correctness** - Complete test suite validates right-hand system behavior
+  - `cross_right()` method retained for backward compatibility
+- ✅ **Verified Correctness** - Complete test suite validates left-hand system behavior
 - 📚 **Updated Documentation** - All examples and documentation reflect new coordinate system
 - 🎯 **Enhanced Clarity** - Comments and docstrings clarified throughout
 
 **Why This Change:**
-1. Mathematical Standard: Right-hand system is the standard in mathematics and physics
-2. Better Interoperability: Compatible with OpenGL, NumPy, SciPy, and most scientific libraries
-3. Improved Readability: Standard cross product formulas match textbook definitions
+1. Graphics Standard: Left-hand system is common in computer graphics and DirectX
+2. Better Interoperability: Compatible with DirectX, Unity, and many game engines
+3. Improved Readability: Standard cross product formulas match graphics pipeline definitions
 
-Perfect for computational geometry, geometric analysis, and discrete differential geometry research!
+Perfect for computer graphics, game development, and 3D visualization applications!
 
 ---
 
@@ -210,7 +210,7 @@ lerped = vec3.lerp(v1, v2, 0.5)  # Linear interpolation
 
 ## Coordinate System Type
 
-This library uses a **right-handed coordinate system** for all vector and quaternion operations (as of v3.0.0).
+This library uses a **left-handed coordinate system** for all vector and quaternion operations (as of v3.0.0).
 
 ```
          +Y
@@ -218,18 +218,18 @@ This library uses a **right-handed coordinate system** for all vector and quater
           |
           |
           +-----> +X
-         /
-        /
-      +Z
+          \
+           \
+            +Z
 ```
 
-**Right-Hand Rule:**
-- Point your right hand's fingers along the +X axis
+**Left-Hand Rule:**
+- Point your left hand's fingers along the +X axis
 - Curl them toward the +Y axis
 - Your thumb points along the +Z axis
-- Cross product: **X × Y = Z** (standard mathematical definition)
+- Cross product: **X × Y = Z** (common in computer graphics)
 
-**Migration from v2.x:** If your code relied on the previous left-handed system, see [HANDEDNESS_CHANGE.md](HANDEDNESS_CHANGE.md) for migration instructions.
+**Migration from v2.x:** If your code relied on the previous right-handed system, see [HANDEDNESS_CHANGE.md](HANDEDNESS_CHANGE.md) for migration instructions.
 
 ---
 
@@ -267,7 +267,7 @@ v3 = v1 * v2           # Component-wise multiplication
 
 ```python
 dot = v1.dot(v2)       # Dot product (float)
-cross = v1.cross(v2)   # Cross product (vec3)
+cross = v1.cross(v2)   # Cross product (vec3) - Left-hand rule
 length = v.length()    # Vector length
 v.normalize()          # Normalize in-place
 normalized = v.normcopy()  # Return normalized copy
@@ -451,7 +451,7 @@ v4 = v1 * 5                     # vec3(5, 0, 0)
 
 # Dot and cross products
 dot = v1.dot(v2)                # 0.0 (perpendicular)
-cross = v1.cross(v2)            # vec3(0, 0, 1) in right-handed system (X × Y = Z)
+cross = v1.cross(v2)            # vec3(0, 0, -1) in left-handed system (X × Y = -Z)
 
 # Length and normalization
 length = v1.length()            # 1.0
@@ -516,8 +516,8 @@ child_in_world = child * parent
 # Create look-at transformation (custom implementation needed)
 def look_at(eye, target, up=vec3(0, 1, 0)):
     forward = (target - eye).normcopy()
-    right = up.cross(forward).normcopy()
-    up_corrected = forward.cross(right)
+    right = forward.cross(up).normcopy()  # Left-hand cross product
+    up_corrected = right.cross(forward)
     return coord3.from_axes(right, up_corrected, forward)
 
 camera = look_at(vec3(10, 10, 10), vec3(0, 0, 0))
@@ -537,8 +537,8 @@ class Camera:
 
     def create_look_at(self, eye, target, up):
         forward = (target - eye).normcopy()
-        right = up.cross(forward).normcopy()
-        up_corrected = forward.cross(right)
+        right = forward.cross(up).normcopy()  # Left-hand cross product
+        up_corrected = right.cross(forward)
 
         c = coord3()
         c.o = eye

{coordinate_system-4.0.0 → coordinate_system-4.0.2}/coordinate_system_binding.cpp

@@ -361,9 +361,9 @@ PYBIND11_MODULE(coordinate_system, m) {
         .def("to_local", [](const coord3& c, const vec3& world) { return coord3_to_local(c, world); }, "Transform to local coordinates", py::arg("world"))
 
         // Advanced accessors
-        .def("VX", [](const coord3& self) { return self.VX(); }, "Get scaled X axis")
-        .def("VY", [](const coord3& self) { return self.VY(); }, "Get scaled Y axis")
-        .def("VZ", [](const coord3& self) { return self.VZ(); }, "Get scaled Z axis")
+        .def("VX", static_cast<vec3(coord3::*)() const>(&coord3::VX), "Get scaled X axis")
+		.def("VY", static_cast<vec3(coord3::*)() const>(&coord3::VY), "Get scaled Y axis")
+		.def("VZ", static_cast<vec3(coord3::*)() const>(&coord3::VZ), "Get scaled Z axis")
         .def("X", &coord3::X, "Get X axis with position offset")
         .def("Y", &coord3::Y, "Get Y axis with position offset")
         .def("Z", &coord3::Z, "Get Z axis with position offset")
@@ -531,7 +531,7 @@ PYBIND11_MODULE(coordinate_system, m) {
     m.attr("EPSILON") = (real)EPSILON;
     m.attr("VERSION") = GCU_VERSION;
 
-    m.attr("__version__") = "3.0.0";
+    m.attr("__version__") = "4.0.1";
     m.attr("__author__") = "PanGuoJun";
     m.attr("__coordinate_system__") = "right-handed";
 }

{coordinate_system-4.0.0 → coordinate_system-4.0.2}/pmsys_minimal.hpp

@@ -878,10 +878,10 @@ struct vector3
     DEVICE_CALLABLE vector3 cross(const vector3& v) const
     {
         vector3 n;
-        // 右手坐标系标准叉乘
-        n.x = (y * v.z - z * v.y);
-        n.y = (z * v.x - x * v.z);
-        n.z = (x * v.y - y * v.x);
+        // 默认使用左手顺序！
+        n.x = -(y * v.z - z * v.y);
+        n.y = -(z * v.x - x * v.z);
+        n.z = -(x * v.y - y * v.x);
         return n;
     }
     DEVICE_CALLABLE vector3 cross_left(const vector3& v) const
@@ -2580,8 +2580,17 @@ inline real vv_angle(crvec2 v1, crvec2 v2)
     if (angle < 0)
         angle += 2 * PI;
     return v1.cross(v2) > 0 ? angle : -angle;
-}/****************************************************************************************************
-*						[Coordinate System (Coordinate Frame)]
+}
+/**************************************************************************************************************\ 
+*  _______ _            _____                     _ _             _         _____           _                  *  
+* |__   __| |          / ____|                   | (_)           | |       / ____|         | |                 *  
+*    | |  | |__   ___ | |     ___   ___  _ __ __| |_ _ __   __ _| |_ ___ | (___  _   _ ___| |_ ___ _ __ ___    *  
+*    | |  | '_ \ / _ \| |    / _ \ / _ \| '__/ _` | | '_ \ / _` | __/ _ \ \___ \| | | / __| __/ _ \ '_ ` _ \   *  
+*    | |  | | | |  __/| |___| (_) | (_) | | | (_| | | | | | (_| | ||  __/ ____) | |_| \__ \ ||  __/ | | | | |  *  
+*    |_|  |_| |_|\___| \_____\___/ \___/|_|  \__,_|_|_| |_|\__,_|\__\___||_____/ \__, |___/\__\___|_| |_| |_|  *  
+*                                                                                 __/ |                        *  
+*                                                                                |___/                         *
+**								[Coordinate System (Coordinate Frame)]										  **
 *
 *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *
 *   The Coordinate System class is specifically encapsulated to simplify coordinate transformations
@@ -2591,12 +2600,12 @@ inline real vv_angle(crvec2 v1, crvec2 v2)
 *
 *  *  *  *  *  *  *  *  *  *  *  *  Detailed Explanation  *  *  *  *  *  *  *  *  *  *  *  *  *  *
 *   Coordinate system transformation is divided into three steps:
-*					projection (/), translation (^), and restoration (*).
+*						projection (/), translation (^), and restoration (*).
 *
 *   The coordinate system itself is denoted as C. Transformations between coordinate systems
 *   can be expressed as G = C2 / C1 - I, where G represents the geometric gradient.
-*					oper(/)  =  C1 * C2^-1
-*					oper(\)  =  C1^-1 * C2
+*						oper(/)  =  C1 * C2^-1
+*						oper(\)  =  C1^-1 * C2
 *
 *   Specifically:
 *   Define a vector V in an intrinsic coordinate system (assuming a flat space where vectors can
@@ -2605,27 +2614,47 @@ inline real vv_angle(crvec2 v1, crvec2 v2)
 *
 *   Take vectors V1 and V2 at adjacent points (1) and (2) respectively,
 *   corresponding to coordinate systems C1 and C2. Then:
-*					V  = V1 * C1 = V2 * C2 =>
-*					V2 = V1 * C1 / C2, let R12 = C1 / C2 =>
-*					V2 = V1 * R12
+*						V  = V1 * C1 = V2 * C2 =>
+*						V2 = V1 * C1 / C2, let R12 = C1 / C2 =>
+*						V2 = V1 * R12
 *
-*   Based on the frame field combination operator theory proposed in this paper,
-*   the direct geometric information extraction formula is:
-*           G = (c₂·c₁⁻¹)/C₂ - I/C₁
+*   Based on the dual-frame normalization theory proposed in this paper,
+*   the geometric connection operator is:
+*           			G_μ = (c(u+h_μ) - c(u))/h_μ
 *   where c is the intrinsic frame field and C is the embedding frame field.
 *
 *   The coordinate system can be used to compute spatial curvature. In the u,v coordinate system,
 *   the curvature tensor is:
-*					Ruv  =  Gu·Gv - Gv·Gu - G[u,v]
+*						R_uv = G_u·G_v - G_v·G_u - G_[u,v]
 *   where:
-*				 Gu = (c(u+du,v)·c⁻¹(u,v))/C(u+du,v) - I/C(u,v)
-*				 Gv = (c(u,v+dv)·c⁻¹(u,v))/C(u,v+dv) - I/C(u,v)
-*				 Connection vector: W = [U, V] (Lie bracket operation)
-*				 G[u,v] = Gu·Wu + Gv·Wv
+*			 	G_u = (c(u+Δ,v) - c(u,v)) / Δ
+*			 	G_v = (c(u,v+Δ) - c(u,v)) / Δ
+*			 	G_[u,v] = connection operator for coordinate commutator [∂_u, ∂_v]
+*
+*   For holonomic coordinate systems (spherical, toroidal, etc.), coordinate basis vectors commute:
+*           			[∂_u, ∂_v] = 0  ⇒  G_[u,v] = 0
+*   Thus the formula simplifies to: R_uv = G_u·G_v - G_v·G_u
+*
+*   **Measurement Function for Curvature Extraction:**
+*   
+*   The measurement function bridges frame bundle curvature to Riemannian geometry:
+*           			M_{ijkl} = √det(g) · ⟨X e_l, e_k⟩
+*   where:
+*      			X = [G_u, G_v]              (Lie bracket curvature operator)
+*      			e_k, e_l                    (tangent basis vectors)
+*      			det(g)                      (determinant of metric tensor)
+*      			⟨·,·⟩                       (inner product in embedding space)
+*
+*   **Riemann Curvature Extraction:**
+*           			R_{ijkl} = M_{ijkl} / √det(g)
+*
+*   **Gaussian Curvature Calculation (verified implementation):**
+*           			K = R_{1212} / det(g)
+* 
+*
+*   This approach provides O(n³) computational complexity for full curvature analysis,
+*   significantly faster than traditional O(n⁶) methods.
 *
-*   Compared with traditional methods, this framework avoids the complex Christoffel symbol
-*   computation chain and directly extracts geometric invariants through frame field combinations,
-*   offering higher computational efficiency and geometric intuitiveness.
 */
 
 // #define	NON_UNIFORM_SCALE

{coordinate_system-4.0.0 → coordinate_system-4.0.2}/setup.py

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