xobjects 0.5.0__tar.gz → 0.5.1__tar.gz

{xobjects-0.5.0/xobjects.egg-info → xobjects-0.5.1}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: xobjects
-Version: 0.5.0
+Version: 0.5.1
 Summary: In-memory serialization and code generator for CPU and GPU
 Home-page: https://xsuite.readthedocs.io/
 Download-URL: https://pypi.python.org/pypi/xobjects

xobjects-0.5.1/xobjects/_version.py

@@ -0,0 +1 @@
+__version__ = "0.5.1"

{xobjects-0.5.0 → xobjects-0.5.1}/xobjects/context_cpu.py

@@ -422,12 +422,6 @@ class ContextCpu(XContext):
             xtr_link_args.append("-fopenmp")
             xtr_compile_args.append("-DXO_CONTEXT_CPU_OPENMP")
             xtr_link_args.append("-DXO_CONTEXT_CPU_OPENMP")
-
-            # https://mac.r-project.org/openmp/
-            # on macos comment the above and uncomment the below flags to compile OpenMP with Xcode clang:
-            # xtr_compile_args.append("-Xclang")
-            # xtr_compile_args.append("-fopenmp")
-            # xtr_link_args.append("-lomp")
         else:
             xtr_compile_args.append("-DXO_CONTEXT_CPU_SERIAL")
             xtr_link_args.append("-DXO_CONTEXT_CPU_SERIAL")

xobjects-0.5.1/xobjects/general.py

@@ -0,0 +1,52 @@
+# copyright ################################# #
+# This file is part of the Xobjects Package.  #
+# Copyright (c) CERN, 2024.                   #
+# ########################################### #
+from numpy.testing import assert_allclose as np_assert_allclose
+import numpy as np
+
+
+class Print:
+    suppress = False
+
+    def __call__(self, *args, **kwargs):
+        if not self.suppress:
+            print(*args, **kwargs)
+
+
+_print = Print()
+
+
+def assert_allclose(a, b, rtol=0, atol=0, max_outliers=0):
+    if hasattr(a, "get"):
+        a = a.get()
+    if hasattr(b, "get"):
+        b = b.get()
+    try:
+        a = np.squeeze(a)
+    except:
+        pass
+    try:
+        b = np.squeeze(b)
+    except:
+        pass
+
+    try:
+        np_assert_allclose(a, b, rtol=rtol, atol=atol)
+    except AssertionError as e:
+        if max_outliers == 0:
+            raise e
+        if not allclose_with_outliers(a, b, rtol, atol, max_outliers):
+            raise AssertionError(
+                "Arrays are not close enough, even with outliers allowed."
+            ) from e
+
+
+def allclose_with_outliers(a, b, rtol=1e-7, atol=0, max_outliers=0):
+    a = np.asanyarray(a)
+    b = np.asanyarray(b)
+    diff = np.abs(a - b)
+    allowed = atol + rtol * np.abs(b)
+    mask = diff > allowed
+    num_outliers = np.count_nonzero(mask)
+    return num_outliers <= max_outliers

{xobjects-0.5.0 → xobjects-0.5.1/xobjects.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: xobjects
-Version: 0.5.0
+Version: 0.5.1
 Summary: In-memory serialization and code generator for CPU and GPU
 Home-page: https://xsuite.readthedocs.io/
 Download-URL: https://pypi.python.org/pypi/xobjects

xobjects-0.5.0/xobjects/_version.py

@@ -1 +0,0 @@
-__version__ = "0.5.0"

xobjects-0.5.0/xobjects/general.py

@@ -1,33 +0,0 @@
-# copyright ################################# #
-# This file is part of the Xobjects Package.  #
-# Copyright (c) CERN, 2024.                   #
-# ########################################### #
-from numpy.testing import assert_allclose as np_assert_allclose
-import numpy as np
-
-
-class Print:
-    suppress = False
-
-    def __call__(self, *args, **kwargs):
-        if not self.suppress:
-            print(*args, **kwargs)
-
-
-_print = Print()
-
-
-def assert_allclose(a, b, rtol=1e-7, atol=1e-7):
-    if hasattr(a, "get"):
-        a = a.get()
-    if hasattr(b, "get"):
-        b = b.get()
-    try:
-        a = np.squeeze(a)
-    except:
-        pass
-    try:
-        b = np.squeeze(b)
-    except:
-        pass
-    np_assert_allclose(a, b, rtol=rtol, atol=atol)