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

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

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: coordinate_system
-Version: 4.0.0
+Version: 4.0.1
 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

{coordinate_system-4.0.0 → coordinate_system-4.0.1}/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.1/coordinate_system.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: coordinate_system
-Version: 4.0.0
+Version: 4.0.1
 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

{coordinate_system-4.0.0 → coordinate_system-4.0.1}/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.1}/pmsys_minimal.hpp

@@ -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:
+*			 	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:
-*				 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
+*      			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.1}/setup.py

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