o reverting to max(-1.0, min(1.0, d_dot_norm)) as approach is the most reliable solution since it's simple, fast, and fully supported by Numba. Numba doesn't fully support np.clip in nopython mode.

Author: rajeeja

uxarray/grid/dual.py

@@ -238,7 +238,7 @@ def _order_nodes(
             d_dot_norm = np.dot(node_zero, node_diff) / (node_zero_mag * node_diff_mag)
 
             # Clamp to valid range for arccos to avoid numerical errors
-            d_dot_norm = np.clip(d_dot_norm, -1.0, 1.0)
+            d_dot_norm = max(-1.0, min(1.0, d_dot_norm))
 
             d_angles[j] = np.arccos(d_dot_norm)