numpy-algebras 0.1.0__tar.gz

numpy_algebras-0.1.0/LICENSE

@@ -0,0 +1,21 @@
+MIT License
+
+Copyright (c) 2026 Daniel Filonik
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in all
+copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+SOFTWARE.

numpy_algebras-0.1.0/PKG-INFO

@@ -0,0 +1,8 @@
+Metadata-Version: 2.4
+Name: numpy-algebras
+Version: 0.1.0
+Summary: Custom number types for numpy.
+Requires-Python: >=3.10
+License-File: LICENSE
+Requires-Dist: numpy
+Dynamic: license-file

numpy_algebras-0.1.0/README.md

@@ -0,0 +1,65 @@
+# algpy
+
+## Build
+
+    pip install .
+
+For development use `pip install -e .`.
+
+## Algebras
+
+The algebra constructors create scalar types that can also be used as NumPy
+dtypes. Parameters specify the squares of the algebra's generators.
+
+### Multicomplex
+
+The multicomplex algebra with parameters `k0, ..., k(n-1)` is generated over
+the reals by commuting elements `e_i` satisfying `e_i² = k_i`. Its basis
+consists of products of distinct generators, giving `2ⁿ` real components.
+`multicomplex` accepts one to four parameters; the numbered
+`multicomplex1`–`multicomplex4` constructors are also available.
+
+```python
+from algpy.algebras import multicomplex
+
+mc = multicomplex(-1, 1)
+```
+
+The two-component specializations are `complex1`, given by
+`multicomplex1(-1)`, `dual1`, given by `multicomplex1(0)`, and `splitcomplex1`,
+given by `multicomplex1(+1)`.
+
+### Cayley-Dickson
+
+The Cayley-Dickson construction builds successively doubled real algebras.
+At each step it adjoins a new generator whose square is set by the next
+parameter; the multiplication uses the Cayley-Dickson conjugation rule and is
+generally noncommutative (and, at higher dimensions, may be nonassociative).
+`cayleydickson` accepts one to four parameters, with numbered
+`cayleydickson1`–`cayleydickson4` constructors for the corresponding
+dimensions.
+
+The `quaternion` type is the Hamilton quaternion specialization
+`cayleydickson2(-1, -1)`. Its components are ordered `(w, x, y, z)`. It
+supports arithmetic, `exp`, `log`, `norm2`, and equality for scalars and NumPy
+arrays.
+
+```python
+from algpy.algebras import cayleydickson, quaternion
+
+cd = cayleydickson(-1, -1)
+q = quaternion(1, 2, 3, 4)
+```
+
+### Clifford
+
+The Clifford algebra with parameters `k0, ..., k(n-1)` is generated by elements
+`e_i` with `e_i² = k_i` and `e_i e_j = -e_j e_i` for distinct generators.
+Products of distinct generators form a `2ⁿ`-element basis. `clifford` accepts
+one or two parameters; `clifford1` and `clifford2` are also available.
+
+```python
+from algpy.algebras import clifford
+
+cl = clifford(-1, 1)
+```

numpy_algebras-0.1.0/pyproject.toml

@@ -0,0 +1,13 @@
+[build-system]
+requires = ["setuptools>=61", "numpy"]
+build-backend = "setuptools.build_meta"
+
+[project]
+name = "numpy-algebras"
+version = "0.1.0"
+description = "Custom number types for numpy."
+requires-python = ">=3.10"
+dependencies = ["numpy"]
+
+[tool.setuptools.packages.find]
+where = ["src"]

numpy_algebras-0.1.0/setup.cfg

@@ -0,0 +1,4 @@
+[egg_info]
+tag_build = 
+tag_date = 0
+

numpy_algebras-0.1.0/setup.py

@@ -0,0 +1,194 @@
+import numpy
+from setuptools import Extension, setup
+
+extensions = [
+    Extension(
+        "algpy.algebras.quaternion._quaternion",
+        ["src/algpy/algebras/quaternion/_quaternion.c"],
+        depends=[
+            "src/algpy/algebras/quaternion/quaternion.h",
+            "src/algpy/vectorspaces/vec4.h",
+            "src/algpy/vectorspaces/operations.h",
+        ],
+        include_dirs=[numpy.get_include()],
+    ),
+    Extension(
+        "algpy.algebras.complex._complex1",
+        ["src/algpy/algebras/complex/_complex1.c"],
+        depends=[
+            "src/algpy/algebras/complex/complex1.h",
+            "src/algpy/vectorspaces/vec2.h",
+            "src/algpy/vectorspaces/operations.h",
+        ],
+        include_dirs=[numpy.get_include()],
+    ),
+    Extension(
+        "algpy.algebras.dual._dual1",
+        ["src/algpy/algebras/dual/_dual1.c"],
+        depends=[
+            "src/algpy/algebras/dual/dual1.h",
+            "src/algpy/vectorspaces/vec2.h",
+            "src/algpy/vectorspaces/operations.h",
+        ],
+        include_dirs=[numpy.get_include()],
+    ),
+    Extension(
+        "algpy.algebras.splitcomplex._splitcomplex1",
+        ["src/algpy/algebras/splitcomplex/_splitcomplex1.c"],
+        depends=[
+            "src/algpy/algebras/splitcomplex/splitcomplex1.h",
+            "src/algpy/vectorspaces/vec2.h",
+            "src/algpy/vectorspaces/operations.h",
+        ],
+        include_dirs=[numpy.get_include()],
+    ),
+    Extension(
+        "algpy.algebras.general.multicomplex._multicomplex1",
+        ["src/algpy/algebras/general/multicomplex/_multicomplex1.c"],
+        depends=[
+            "src/algpy/algebras/complex/complex1.h",
+            "src/algpy/algebras/general/multicomplex/multicomplex1.h",
+            "src/algpy/vectorspaces/vec1.h",
+            "src/algpy/vectorspaces/vec2.h",
+            "src/algpy/vectorspaces/operations.h",
+        ],
+        include_dirs=[numpy.get_include()],
+    ),
+    Extension(
+        "algpy.algebras.general.clifford._clifford1",
+        ["src/algpy/algebras/general/clifford/_clifford1.c"],
+        depends=[
+            "src/algpy/algebras/general/clifford/clifford1.h",
+            "src/algpy/algebras/general/multicomplex/multicomplex1.h",
+            "src/algpy/algebras/general/multicomplex/_multicomplex1.c",
+            "src/algpy/vectorspaces/vec1.h",
+            "src/algpy/vectorspaces/vec2.h",
+            "src/algpy/vectorspaces/operations.h",
+            "src/algpy/algebras/complex/complex1.h",
+        ],
+        include_dirs=[numpy.get_include()],
+    ),
+    Extension(
+        "algpy.algebras.general.clifford._clifford2",
+        ["src/algpy/algebras/general/clifford/_clifford2.c"],
+        depends=[
+            "src/algpy/algebras/general/clifford/clifford2.h",
+            "src/algpy/algebras/general/cayleydickson/cayleydickson2.h",
+            "src/algpy/algebras/complex/complex1.h",
+            "src/algpy/vectorspaces/vec2.h",
+            "src/algpy/vectorspaces/vec4.h",
+            "src/algpy/vectorspaces/operations.h",
+        ],
+        include_dirs=[numpy.get_include()],
+    ),
+    Extension(
+        "algpy.algebras.general.cayleydickson._cayleydickson1",
+        ["src/algpy/algebras/general/cayleydickson/_cayleydickson1.c"],
+        depends=[
+            "src/algpy/algebras/general/cayleydickson/cayleydickson1.h",
+            "src/algpy/algebras/general/multicomplex/multicomplex1.h",
+            "src/algpy/algebras/general/multicomplex/_multicomplex1.c",
+            "src/algpy/vectorspaces/vec1.h",
+            "src/algpy/vectorspaces/vec2.h",
+            "src/algpy/vectorspaces/operations.h",
+            "src/algpy/algebras/complex/complex1.h",
+        ],
+        include_dirs=[numpy.get_include()],
+    ),
+    Extension(
+        "algpy.algebras.general.cayleydickson._cayleydickson2",
+        ["src/algpy/algebras/general/cayleydickson/_cayleydickson2.c"],
+        depends=[
+            "src/algpy/algebras/general/cayleydickson/cayleydickson2.h",
+            "src/algpy/algebras/complex/complex1.h",
+            "src/algpy/vectorspaces/vec2.h",
+            "src/algpy/vectorspaces/vec4.h",
+            "src/algpy/vectorspaces/operations.h",
+        ],
+        include_dirs=[numpy.get_include()],
+    ),
+    Extension(
+        "algpy.algebras.general.cayleydickson._cayleydickson3",
+        ["src/algpy/algebras/general/cayleydickson/_cayleydickson3.c"],
+        depends=[
+            "src/algpy/algebras/general/cayleydickson/cayleydickson3.h",
+            "src/algpy/algebras/general/cayleydickson/cayleydickson2.h",
+            "src/algpy/algebras/complex/complex1.h",
+            "src/algpy/vectorspaces/vec2.h",
+            "src/algpy/vectorspaces/vec3.h",
+            "src/algpy/vectorspaces/vec4.h",
+            "src/algpy/vectorspaces/vec8.h",
+            "src/algpy/vectorspaces/operations.h",
+        ],
+        include_dirs=[numpy.get_include()],
+    ),
+    Extension(
+        "algpy.algebras.general.cayleydickson._cayleydickson4",
+        ["src/algpy/algebras/general/cayleydickson/_cayleydickson4.c"],
+        depends=[
+            "src/algpy/algebras/general/cayleydickson/cayleydickson4.h",
+            "src/algpy/algebras/general/cayleydickson/cayleydickson3.h",
+            "src/algpy/algebras/general/cayleydickson/cayleydickson2.h",
+            "src/algpy/algebras/complex/complex1.h",
+            "src/algpy/vectorspaces/vec2.h",
+            "src/algpy/vectorspaces/vec3.h",
+            "src/algpy/vectorspaces/vec4.h",
+            "src/algpy/vectorspaces/vec8.h",
+            "src/algpy/vectorspaces/vec16.h",
+            "src/algpy/vectorspaces/operations.h",
+        ],
+        include_dirs=[numpy.get_include()],
+    ),
+    Extension(
+        "algpy.algebras.general.multicomplex._multicomplex2",
+        ["src/algpy/algebras/general/multicomplex/_multicomplex2.c"],
+        depends=[
+            "src/algpy/algebras/complex/complex1.h",
+            "src/algpy/algebras/general/multicomplex/multicomplex1.h",
+            "src/algpy/algebras/general/multicomplex/multicomplex2.h",
+            "src/algpy/vectorspaces/vec1.h",
+            "src/algpy/vectorspaces/vec2.h",
+            "src/algpy/vectorspaces/vec4.h",
+            "src/algpy/vectorspaces/operations.h",
+        ],
+        include_dirs=[numpy.get_include()],
+    ),
+    Extension(
+        "algpy.algebras.general.multicomplex._multicomplex3",
+        ["src/algpy/algebras/general/multicomplex/_multicomplex3.c"],
+        depends=[
+            "src/algpy/algebras/complex/complex1.h",
+            "src/algpy/algebras/general/multicomplex/multicomplex1.h",
+            "src/algpy/algebras/general/multicomplex/multicomplex2.h",
+            "src/algpy/algebras/general/multicomplex/multicomplex3.h",
+            "src/algpy/vectorspaces/vec1.h",
+            "src/algpy/vectorspaces/vec2.h",
+            "src/algpy/vectorspaces/vec3.h",
+            "src/algpy/vectorspaces/vec4.h",
+            "src/algpy/vectorspaces/vec8.h",
+            "src/algpy/vectorspaces/operations.h",
+        ],
+        include_dirs=[numpy.get_include()],
+    ),
+    Extension(
+        "algpy.algebras.general.multicomplex._multicomplex4",
+        ["src/algpy/algebras/general/multicomplex/_multicomplex4.c"],
+        depends=[
+            "src/algpy/algebras/complex/complex1.h",
+            "src/algpy/algebras/general/multicomplex/multicomplex1.h",
+            "src/algpy/algebras/general/multicomplex/multicomplex2.h",
+            "src/algpy/algebras/general/multicomplex/multicomplex3.h",
+            "src/algpy/algebras/general/multicomplex/multicomplex4.h",
+            "src/algpy/vectorspaces/vec1.h",
+            "src/algpy/vectorspaces/vec2.h",
+            "src/algpy/vectorspaces/vec3.h",
+            "src/algpy/vectorspaces/vec4.h",
+            "src/algpy/vectorspaces/vec8.h",
+            "src/algpy/vectorspaces/vec16.h",
+            "src/algpy/vectorspaces/operations.h",
+        ],
+        include_dirs=[numpy.get_include()],
+    ),
+]
+
+setup(ext_modules=extensions)

numpy_algebras-0.1.0/src/algpy/algebras/__init__.py

@@ -0,0 +1,34 @@
+"""Custom number types for numpy, implemented in Cython."""
+from .general.multicomplex import multicomplex, multicomplex1, multicomplex2, multicomplex3, multicomplex4
+from .general.clifford import clifford, clifford1, clifford2
+from .general.cayleydickson import (
+    cayleydickson,
+    cayleydickson1,
+    cayleydickson2,
+    cayleydickson3,
+    cayleydickson4,
+)
+from .quaternion import quaternion
+from .complex import complex1
+from .dual import dual1
+from .splitcomplex import splitcomplex1
+
+__all__ = [
+    "multicomplex",
+    "multicomplex1",
+    "multicomplex2",
+    "multicomplex3",
+    "multicomplex4",
+    "clifford",
+    "clifford1",
+    "clifford2",
+    "cayleydickson",
+    "cayleydickson1",
+    "cayleydickson2",
+    "cayleydickson3",
+    "cayleydickson4",
+    "quaternion",
+    "complex1",
+    "dual1",
+    "splitcomplex1",
+]

numpy_algebras-0.1.0/src/algpy/algebras/_conversions.py

@@ -0,0 +1,19 @@
+import numpy as _np
+
+
+def attach(scalar_type, n):
+    def _to_float(a, dtype=float):
+        a = _np.asarray(a, dtype=scalar_type)
+        flat = _np.ascontiguousarray(a.reshape(-1))
+        return flat.view(float).reshape(a.shape + (n,)).astype(dtype, copy=False)
+
+    def _from_float(a):
+        a = _np.asarray(a, dtype=float)
+        if a.ndim == 0 or a.shape[-1] != n:
+            raise ValueError(f"last axis must have length {n}")
+        flat = _np.ascontiguousarray(a.reshape(-1))
+        return flat.view(scalar_type).reshape(a.shape[:-1])
+
+    scalar_type.N = n
+    scalar_type.to_float = staticmethod(_to_float)
+    scalar_type.from_float = staticmethod(_from_float)

numpy_algebras-0.1.0/src/algpy/algebras/complex/__init__.py

@@ -0,0 +1,7 @@
+"""Complex1 number type registered as a NumPy user-defined dtype."""
+from ._complex1 import complex1
+from .._conversions import attach as _attach_conversions
+
+_attach_conversions(complex1, 2)
+
+__all__ = ["complex1"]