differintC 0.0.1__cp313-cp313-win_amd64.whl → 0.0.2__cp313-cp313-win_amd64.whl

differintC/__init__.py

@@ -0,0 +1,80 @@
+import os
+import sys
+import ctypes
+from pathlib import Path
+import importlib.resources
+import warnings
+
+
+def _load_dll():
+    if sys.platform != "win32":
+        return
+
+    # Get package directory
+    package_dir = Path(__file__).parent
+
+    # 1. Try package directory (where __init__.py is)
+    dll_path = package_dir / "libfftw3-3.dll"
+    if dll_path.exists():
+        os.add_dll_directory(str(package_dir))
+        print(f"Loaded FFTW DLL from package directory: {dll_path}")
+        return
+
+    # 2. Try site-packages root directory
+    site_packages_dll = package_dir.parent / "libfftw3-3.dll"
+    if site_packages_dll.exists():
+        os.add_dll_directory(str(site_packages_dll.parent))
+        print(f"Loaded FFTW DLL from site-packages root: {site_packages_dll}")
+        return
+
+    # 3. Try using importlib.resources
+    try:
+        with importlib.resources.path("differintC", "libfftw3-3.dll") as dll_path:
+            if dll_path.exists():
+                os.add_dll_directory(str(dll_path.parent))
+                print(f"Loaded FFTW DLL via importlib: {dll_path}")
+                return
+    except Exception as e:
+        warnings.warn(f"importlib.resources failed: {e}")
+
+    # 4. Try Windows system directory
+    system_dll = Path("C:/Windows/System32") / "libfftw3-3.dll"
+    if system_dll.exists():
+        os.add_dll_directory(str(system_dll.parent))
+        print(f"Loaded FFTW DLL from system directory: {system_dll}")
+        return
+
+    # 5. Try loading directly (if in PATH)
+    try:
+        ctypes.CDLL("libfftw3-3.dll")
+        print("Loaded FFTW DLL directly from PATH")
+        return
+    except OSError:
+        pass
+
+    # If all else fails
+    raise RuntimeError(
+        "FFTW DLL not found! Checked locations:\n"
+        + f"1. Package directory: {package_dir}/libfftw3-3.dll\n"
+        + f"2. Site-packages root: {package_dir.parent}/libfftw3-3.dll\n"
+        + f"3. System directory: C:/Windows/System32/libfftw3-3.dll"
+    )
+
+
+# Load DLL before importing C++ module
+try:
+    _load_dll()
+except RuntimeError as e:
+    print(f"FFTW DLL loading failed: {e}")
+    # Try to proceed anyway - might be available in PATH
+    try:
+        ctypes.CDLL("libfftw3-3.dll")
+        print("Successfully loaded FFTW DLL after initial failure")
+    except OSError:
+        raise
+
+# Import the C++ extension
+from differintC import *  # Directly import from site-packages
+
+# Add version attribute
+__version__ = "0.0.2"

{differintc-0.0.1.dist-info → differintc-0.0.2.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.2
 Name: differintC
-Version: 0.0.1
+Version: 0.0.2
 Summary: Fast C++ implementation of fractional calculus operators via pybind11
 Author: Parsa Roshanak
 License: MIT

differintc-0.0.2.dist-info/RECORD

@@ -0,0 +1,7 @@
+differintC/__init__.py,sha256=3S0KJat-cNbhGuzPDJ--skatjDg05e8FyX8MwURX-_U,2521
+differintC/libfftw3-3.dll,sha256=ABqDXS8lrAZhWA1n21YiFXwoSBwGp_PvRa2TmDm_VKE,2712765
+differintC/libfftw3-3.lib,sha256=O1n7vwNRUjJuvlQ91RWn7xDXMD8PYDM0dHXXETX_ruo,241990
+differintc-0.0.2.dist-info/METADATA,sha256=cZMvoN1HL9Gwg0-YLTIieH_vuNj64rW8Xy7ftu69Cwg,3285
+differintc-0.0.2.dist-info/WHEEL,sha256=_PdgJ-W7uMAsWwnGnEdiicDU_PKcobRdffEM23gmN6g,106
+differintc-0.0.2.dist-info/licenses/LICENSE,sha256=jqB4xKqXwsqHxw3nhNfX51kOCtYhxNDQFqCIl5Qj8h4,1091
+differintc-0.0.2.dist-info/RECORD,,

differintc-0.0.1.dist-info/RECORD

@@ -1,6 +0,0 @@
-python/CMakeLists.txt,sha256=YH4PMGwp8NqUb5ZsKmOjtoRmwdjqAr9T8bcRf59sJcE,401
-python/module.cpp,sha256=2ptoCWNYUXnhYlJfxMm_q6NM66goo1lB9WV6mbysTF0,5663
-differintc-0.0.1.dist-info/METADATA,sha256=TPcJPxQPP4AAdefMQfHPLoj-i-IjFkDpx0ixQbe0rh8,3285
-differintc-0.0.1.dist-info/WHEEL,sha256=_PdgJ-W7uMAsWwnGnEdiicDU_PKcobRdffEM23gmN6g,106
-differintc-0.0.1.dist-info/licenses/LICENSE,sha256=jqB4xKqXwsqHxw3nhNfX51kOCtYhxNDQFqCIl5Qj8h4,1091
-differintc-0.0.1.dist-info/RECORD,,

python/CMakeLists.txt

@@ -1,17 +0,0 @@
-pybind11_add_module(differintC module.cpp ../src/differint.cpp)
-
-target_include_directories(differintC
-    PRIVATE ${CMAKE_SOURCE_DIR}/include
-)
-
-target_link_libraries(differintC
-    PRIVATE differint_core
-)
-
-set_target_properties(differintC PROPERTIES
-    PREFIX ""
-    SUFFIX ".pyd"
-)
-
-install(TARGETS differintC
-        LIBRARY DESTINATION ${CMAKE_INSTALL_PREFIX}/Lib/site-packages)

python/module.cpp

@@ -1,110 +0,0 @@
-#include <pybind11/stl.h>
-#include <pybind11/pybind11.h>
-#include <pybind11/numpy.h>
-#include <vector>
-#include <stdexcept>
-
-namespace py = pybind11;
-
-namespace differint {
-    std::vector<double> RL(double, const std::vector<double>&, double, double, size_t);
-    double RLpoint(double, const std::vector<double>&, double, double, size_t);
-    std::vector<double> GL(double, const std::vector<double>&, double, double, size_t);
-    double GLpoint(double, const std::vector<double>&, double, double, size_t);
-}
-using namespace differint;
-
-
-
-// Helper to sample a Python callable into a std::vector<double>
-template <typename Func>
-std::vector<double> call_and_sample(Func&& f, double a, double b, size_t n) {
-    std::vector<double> vals(n);
-    double step = (b - a) / (n - 1);
-    for (size_t i = 0; i < n; ++i)
-        vals[i] = f(a + step * i);
-    return vals;
-}
-
-// Overload 1: Input is a NumPy array or list (accept py::array_t<double>)
-std::vector<double> prepare_fvals(py::array_t<double> arr, size_t expected_size) {
-    if (arr.size() != expected_size)
-        throw std::runtime_error("Input array length must equal num_points");
-    // Copy data from NumPy array (you could do zero-copy with caution, but copying is safer here)
-    std::vector<double> v(arr.size());
-    std::memcpy(v.data(), arr.data(), sizeof(double) * arr.size());
-    return v;
-}
-
-// Overload 2: Input is a Python callable (py::function)
-std::vector<double> prepare_fvals(py::function func, double a, double b, size_t n) {
-    return call_and_sample([&func](double x) {
-        return func(x).cast<double>();
-    }, a, b, n);
-}
-
-
-
-
-
-PYBIND11_MODULE(differintC, m) {
-    m.doc() = "Fast fractional calculus operators in C++ with Python bindings";
-
-    // RL (whole array)
-    m.def("RL", [](double alpha, py::object f, double domain_start, double domain_end, size_t num_points) {
-        // Dispatch based on type
-        if (py::isinstance<py::array>(f)) {
-            return RL(alpha, prepare_fvals(f.cast<py::array_t<double>>(), num_points), domain_start, domain_end, num_points);
-        } else if (py::isinstance<py::function>(f)) {
-            return RL(alpha, prepare_fvals(f.cast<py::function>(), domain_start, domain_end, num_points), domain_start, domain_end, num_points);
-        } else if (py::isinstance<py::list>(f)) {
-            // treat list as array
-            return RL(alpha, prepare_fvals(f.cast<py::array_t<double>>(), num_points), domain_start, domain_end, num_points);
-        } else {
-            throw std::runtime_error("Unsupported input type for function f");
-        }
-    }, py::arg("alpha"), py::arg("f"), py::arg("domain_start") = 0.0, py::arg("domain_end") = 1.0, py::arg("num_points") = 100);
-
-    // RLpoint (single point)
-    m.def("RLpoint", [](double alpha, py::object f, double domain_start, double domain_end, size_t num_points) {
-        if (py::isinstance<py::array>(f)) {
-            return RLpoint(alpha, prepare_fvals(f.cast<py::array_t<double>>(), num_points), domain_start, domain_end, num_points);
-        } else if (py::isinstance<py::function>(f)) {
-            return RLpoint(alpha, prepare_fvals(f.cast<py::function>(), domain_start, domain_end, num_points), domain_start, domain_end, num_points);
-        } else if (py::isinstance<py::list>(f)) {
-            return RLpoint(alpha, prepare_fvals(f.cast<py::array_t<double>>(), num_points), domain_start, domain_end, num_points);
-        } else {
-            throw std::runtime_error("Unsupported input type for function f");
-        }
-    }, py::arg("alpha"), py::arg("f"), py::arg("domain_start") = 0.0, py::arg("domain_end") = 1.0, py::arg("num_points") = 100);
-
-
-
-    m.def("GL", [](double alpha, py::object f, double domain_start, double domain_end, size_t num_points) {
-        // Dispatch based on type
-        if (py::isinstance<py::array>(f)) {
-            return GL(alpha, prepare_fvals(f.cast<py::array_t<double>>(), num_points), domain_start, domain_end, num_points);
-        } else if (py::isinstance<py::function>(f)) {
-            return GL(alpha, prepare_fvals(f.cast<py::function>(), domain_start, domain_end, num_points), domain_start, domain_end, num_points);
-        } else if (py::isinstance<py::list>(f)) {
-            // treat list as array
-            return GL(alpha, prepare_fvals(f.cast<py::array_t<double>>(), num_points), domain_start, domain_end, num_points);
-        } else {
-            throw std::runtime_error("Unsupported input type for function f");
-        }
-    }, py::arg("alpha"), py::arg("f"), py::arg("domain_start") = 0.0, py::arg("domain_end") = 1.0, py::arg("num_points") = 100);
-
-    // GLpoint (single point)
-    m.def("GLpoint", [](double alpha, py::object f, double domain_start, double domain_end, size_t num_points) {
-        if (py::isinstance<py::array>(f)) {
-            return GLpoint(alpha, prepare_fvals(f.cast<py::array_t<double>>(), num_points), domain_start, domain_end, num_points);
-        } else if (py::isinstance<py::function>(f)) {
-            return GLpoint(alpha, prepare_fvals(f.cast<py::function>(), domain_start, domain_end, num_points), domain_start, domain_end, num_points);
-        } else if (py::isinstance<py::list>(f)) {
-            return GLpoint(alpha, prepare_fvals(f.cast<py::array_t<double>>(), num_points), domain_start, domain_end, num_points);
-        } else {
-            throw std::runtime_error("Unsupported input type for function f");
-        }
-    }, py::arg("alpha"), py::arg("f"), py::arg("domain_start") = 0.0, py::arg("domain_end") = 1.0, py::arg("num_points") = 100);
-
-}