wrapping up

Author: SwayamInSync

quaddtype/README.md

@@ -27,7 +27,7 @@ np.array([1,2,3], dtype=QuadPrecDType("longdouble"))
 
 The code needs the quad precision pieces of the sleef library, which
 is not available on most systems by default, so we have to generate
-that first.  The below assumes one has the required pieces to build
+that first. The below assumes one has the required pieces to build
 sleef (cmake and libmpfr-dev), and that one is in the package
 directory locally.
 
@@ -40,6 +40,7 @@ cd ..
 ```
 
 Building the `numpy-quaddtype` package from locally installed sleef:
+
 ```bash
 export SLEEF_DIR=$PWD/sleef/build
 export LIBRARY_PATH=$SLEEF_DIR/lib
@@ -57,10 +58,13 @@ export LDFLAGS="-Wl,-rpath,$SLEEF_DIR/lib -fopenmp -latomic -lpthread"
 export CFLAGS="-fPIC"
 export CXXFLAGS="-fPIC"
 
+# To build without QBLAS (default for MSVC)
+# export CFLAGS="-fPIC -DDISABLE_QUADBLAS"
+# export CXXFLAGS="-fPIC -DDISABLE_QUADBLAS"
+
 python -m pip install . -v --no-build-isolation -Cbuilddir=build -C'compile-args=-v'
 
 # Run the tests
 cd ..
 python -m pytest
 ```
-

quaddtype/numpy_quaddtype/src/umath/matmul.h

@@ -1,35 +1,12 @@
 #ifndef _QUADDTYPE_MATMUL_H
 #define _QUADDTYPE_MATMUL_H
 
-/**
- * Quad Precision Matrix Multiplication for NumPy
- *
- * This module implements matrix multiplication functionality for the QuadPrecDType
- * by registering custom loops with numpy's matmul generalized ufunc.
- *
- * Supports all matmul operation types:
- * - Vector-vector (dot product): (n,) @ (n,) -> scalar
- * - Matrix-vector: (m,n) @ (n,) -> (m,)
- * - Vector-matrix: (n,) @ (n,p) -> (p,)
- * - Matrix-matrix: (m,n) @ (n,p) -> (m,p)
- *
- * Uses naive algorithms optimized for correctness rather than performance.
- * For production use, consider integration with QBLAS optimized routines.
- */
-
 #include <Python.h>
 
 #ifdef __cplusplus
 extern "C" {
 #endif
 
-/**
- * Initialize the matmul operations for the quad precision dtype.
- * This function registers the matmul generalized ufunc with numpy.
- *
- * @param numpy The numpy module object
- * @return 0 on success, -1 on failure
- */
 int
 init_matmul_ops(PyObject *numpy);
 

quaddtype/release_tracker.md

@@ -91,5 +91,3 @@
 | floor         | ✅    |
 | ceil          | ✅    |
 | trunc         | ✅    |
-
-- Fixing QBLAS integration to work unaligned arrays without or recovering from bad allocation fallback