nulapack 0.1.0__cp312-cp312-win_amd64.whl → 0.2.0__cp312-cp312-win_amd64.whl

nulapack/__init__.py

@@ -1,6 +1,7 @@
 from importlib.metadata import PackageNotFoundError, version
 
 from .cholesky import cholesky
+from .crout import crout
 from .doolittle import doolittle
 from .gauss_seidel import gauss_seidel
 from .jacobi import jacobi
@@ -13,4 +14,4 @@ except PackageNotFoundError:
     # package is not installed
     pass
 
-__all__ = ["__version__", "cholesky", "doolittle", "gauss_seidel", "jacobi", "thomas"]
+__all__ = ["__version__", "cholesky", "crout", "doolittle", "gauss_seidel", "jacobi", "thomas"]

nulapack/crout.py

@@ -0,0 +1,44 @@
+import _nulapack
+import numpy as np
+
+
+def crout(a: np.ndarray):
+    """
+    Compute the LU Crout decomposition of a general matrix A.
+
+    Parameters
+    ----------
+    a : ndarray
+        Coefficient matrix (n x n) stored as a full matrix.
+
+    Returns
+    -------
+    L : ndarray
+        Lower triangular matrix from the factorization.
+    U : ndarray
+        Upper triangular matrix from the factorization (with ones on diagonal).
+    info : int
+        0 if success, <0 if a zero diagonal in L was detected.
+    """
+    a = np.ascontiguousarray(a)
+    n = a.shape[0]
+
+    a_flat = a.ravel(order="C")
+    l_flat = np.zeros_like(a_flat)
+    u_flat = np.zeros_like(a_flat)
+    info = np.zeros(1, dtype=np.int32)
+
+    if np.issubdtype(a.dtype, np.floating):
+        if a.dtype == np.float32:
+            _nulapack.sgectrf(n, a_flat, l_flat, u_flat, info)
+        else:  # float64
+            _nulapack.dgectrf(n, a_flat, l_flat, u_flat, info)
+    elif np.issubdtype(a.dtype, np.complexfloating):
+        if a.dtype == np.complex64:
+            _nulapack.cgectrf(n, a_flat, l_flat, u_flat, info)
+        else:  # complex128
+            _nulapack.zgectrf(n, a_flat, l_flat, u_flat, info)
+    else:
+        raise TypeError(f"Unsupported array dtype: {a.dtype}")
+
+    return l_flat.reshape(n, n, order="C"), u_flat.reshape(n, n, order="C"), int(info[0])

{nulapack-0.1.0.dist-info → nulapack-0.2.0.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: nulapack
-Version: 0.1.0
+Version: 0.2.0
 Summary: NUmerical Linear Algebra PACKage
 Keywords: numerical,linear,algebra
 Author-Email: Saud Zahir <m.saud.zahir@gmail.com>

nulapack-0.2.0.dist-info/RECORD

@@ -0,0 +1,13 @@
+nulapack-0.2.0.dist-info/METADATA,sha256=DQjLWgONIz7g2ey4AN6FNQ0NC8K1fJBzueSQwaFAJ9U,45248
+nulapack-0.2.0.dist-info/WHEEL,sha256=50PeAbplA6PkI0hYOYoeacB9US1R6EguyfOnsccH0WU,85
+nulapack-0.2.0.dist-info/LICENSE,sha256=IwGE9guuL-ryRPEKi6wFPI_zOhg7zDZbTYuHbSt_SAk,35823
+_nulapack.cp312-win_amd64.pyd,sha256=NSXU-lFo0gHPrlcgt1NhxxqqWX4mu8Fe7y9-ezaC-Gk,188166
+_nulapack.cp312-win_amd64.dll.a,sha256=4extQMPR95OK25Z1J7lDOs7PavUkMbMWkftHz-5q0Ak,1752
+nulapack/cholesky.py,sha256=q22LfUITIW4neCoND2EVdNEKPcf5zu12XK3mLYXjZS8,1397
+nulapack/crout.py,sha256=xqLzhYx_vImxvHlIoZJWS9QUdY4-js_3TqtIEXTEqcM,1350
+nulapack/doolittle.py,sha256=xGPG2UZ-OYC-ppbBElQBFr1Ir5lbVah9Jsg1Zi-7ikU,1334
+nulapack/gauss_seidel.py,sha256=kz97lJm7dl_kxg5MXYGmMdtDgkR45sk2QduEVmG-N0o,1419
+nulapack/jacobi.py,sha256=Cdy4Q8ZZMZHFSDhYdb-QBBdqbrerRVk7WZG6kEaDIl8,1403
+nulapack/thomas.py,sha256=RHi1a5ezy7McoJSJZUEPb3y9TY2uAlu2JiJW3D-TfjY,1174
+nulapack/__init__.py,sha256=q29f7X9G_onZDSRatW3iFInQtN2iA1egHl_ZdguvWzY,455
+nulapack-0.2.0.dist-info/RECORD,,

nulapack-0.1.0.dist-info/RECORD

@@ -1,12 +0,0 @@
-nulapack-0.1.0.dist-info/METADATA,sha256=bVRCi7e8mpN7gE2XmqCPR4HvkSAyJjBdQtR258959BA,45248
-nulapack-0.1.0.dist-info/WHEEL,sha256=50PeAbplA6PkI0hYOYoeacB9US1R6EguyfOnsccH0WU,85
-nulapack-0.1.0.dist-info/LICENSE,sha256=IwGE9guuL-ryRPEKi6wFPI_zOhg7zDZbTYuHbSt_SAk,35823
-_nulapack.cp312-win_amd64.pyd,sha256=c4aeWFhiLM8zuhGlOOlCCx3IKEjlJbpz5AjEYXKcPZo,170410
-_nulapack.cp312-win_amd64.dll.a,sha256=tNJcPfqGyAdLLQPl2cR8zY5Ju4yNrfCksTTkOfMwZ30,1752
-nulapack/cholesky.py,sha256=q22LfUITIW4neCoND2EVdNEKPcf5zu12XK3mLYXjZS8,1397
-nulapack/doolittle.py,sha256=xGPG2UZ-OYC-ppbBElQBFr1Ir5lbVah9Jsg1Zi-7ikU,1334
-nulapack/gauss_seidel.py,sha256=kz97lJm7dl_kxg5MXYGmMdtDgkR45sk2QduEVmG-N0o,1419
-nulapack/jacobi.py,sha256=Cdy4Q8ZZMZHFSDhYdb-QBBdqbrerRVk7WZG6kEaDIl8,1403
-nulapack/thomas.py,sha256=RHi1a5ezy7McoJSJZUEPb3y9TY2uAlu2JiJW3D-TfjY,1174
-nulapack/__init__.py,sha256=VHZWNkFUoLNkxzMMTI2GEAE2XDGoQqiBAeIa_pWg8GQ,421
-nulapack-0.1.0.dist-info/RECORD,,