sleef based quadprecision

Author: SwayamInSync

quaddtype2/meson.build

@@ -0,0 +1,45 @@
+project(
+  'quaddtype2',
+  'c',
+)
+
+py_mod = import('python')
+py = py_mod.find_installation('/home/rootacess/numpy-user-dtypes/quaddtype2/.venv/bin/python3')
+
+# incdir_numpy = run_command(py,
+#   [
+#     '-c',
+#     'import numpy; print(numpy.get_include())'
+#   ],
+#   check: true
+# ).stdout().strip()
+
+includes = include_directories(
+  [
+    '.venv/lib/python3.12/site-packages/numpy/_core/include',
+    'quaddtype/src'
+  ]
+)
+
+srcs = [
+  'quaddtype/src/dtype.c',
+  'quaddtype/src/quaddtype_main.c',
+]
+
+py.install_sources(
+  [
+    'quaddtype/__init__.py',
+    'quaddtype/quadscalar.py'
+  ],
+  subdir: 'quaddtype',
+  pure: false
+)
+
+py.extension_module(
+  '_quaddtype_main',
+  srcs,
+  c_args: ['-g', '-O0' ,'-lsleef', '-lsleefquad'],
+  install: true,
+  subdir: 'quaddtype',
+  include_directories: includes
+)

quaddtype2/quaddtype/__init__.py

@@ -1,4 +1,4 @@
-from .scalar import QuadScalar
+from .quadscalar import QuadScalar
 from ._quaddtype_main import QuadDType
 
 __all__ = ["QuadScalar", "QuadDType"]

quaddtype2/quaddtype/scalar.py renamed to quaddtype2/quaddtype/quadscalar.py

@@ -1,4 +1,6 @@
 #include <Python.h>
+#include<sleef.h>
+#include<sleefquad.h>
 
 #define PY_ARRAY_UNIQUE_SYMBOL quaddtype_ARRAY_API
 #define PY_UFUNC_UNIQUE_SYMBOL quaddtype_UFUNC_API
@@ -23,8 +25,8 @@ QuadDTypeObject *new_quaddtype_instance(void)
         return NULL;
     }
 
-    new->base.elsize = sizeof(__float128);
-    new->base.alignment = _Alignof(__float128);
+    new->base.elsize = sizeof(Sleef_quad);
+    new->base.alignment = _Alignof(Sleef_quad);
     return new;
 }
 
@@ -40,7 +42,7 @@ static void quaddtype_dealloc(QuadDTypeObject *self)
 
 static PyObject *quaddtype_repr(QuadDTypeObject *self)
 {
-    PyObject *res = PyUnicode_FromString("This is a quad (128-bit float) dtype.");
+    PyObject *res = PyUnicode_FromString("This is a Sleef based quad (128-bit float) dtype.");
     return res;
 }
 

quaddtype2/quaddtype/src/dtype.c

@@ -0,0 +1,14 @@
+#!/bin/bash
+set -xeuo pipefail
+IFS=$'\n\t'
+
+if [ -d "build/" ]
+then
+    rm -r build
+fi
+
+#meson setup build -Db_sanitize=address,undefined
+meson setup build
+python -m pip uninstall -y quaddtype
+python -m pip install . -v --no-build-isolation -Cbuilddir=build -C'compile-args=-v' -Csetup-args="-Dbuildtype=debug"
+#python -m pip install . -v --no-build-isolation -Cbuilddir=build -C'compile-args=-v'