renamed fastkde to polyexpkde

calvinmccarter · Sep 12, 2023 · 9122bc5 · 9122bc5
1 parent 990ab70
commit 9122bc5
Show file tree

Hide file tree

Showing 8 changed files with 21 additions and 29 deletions.
diff --git a/examples/basic-example.ipynb b/examples/basic-example.ipynb
diff --git a/fastkde/__init__.py b/fastkde/__init__.py
diff --git a/polyexpkde/__init__.py b/polyexpkde/__init__.py
@@ -0,0 +1,4 @@
+from polyexpkde.ksum import h_Gauss_to_K
+from polyexpkde.ksum import kernel_density
+from polyexpkde.ksum import kernel_eval
+from polyexpkde.polyexp_kde import PolyExpKernelDensity
diff --git a/fastkde/ksum.py → polyexpkde/ksum.py b/fastkde/ksum.py → polyexpkde/ksum.py
@@ -157,7 +157,7 @@ def ksum_numba(x, y, x_eval, h, betas, output, counts, coefs, Ly, Ry):
         for i in range(n_eval):
             ix = np.round(counts[i])
             if ix == 0:
-                exp_mult = np.exp((x_eval[i] - x[0]) / h)
+                exp_mult = max(np.exp((x_eval[i] - x[0]) / h), K_EPS)
                 output[i] += (
                     betas[orddo] * np.power(x[0] - x_eval[i], orddo) 
                     / denom * exp_mult * y[0]
@@ -172,7 +172,7 @@ def ksum_numba(x, y, x_eval, h, betas, output, counts, coefs, Ly, Ry):
                 for j in range(orddo + 1):
                     output[i] += betas[orddo] * coefs[j] * (
                         np.power(x_eval[i], orddo - j) * Ly[j, ix - 1] * exp_mult
-                        + np.power(-x_eval[i], orddo - j) * Ry[j, ix - 1] / max(exp_mult, 1e-300)
+                        + np.power(-x_eval[i], orddo - j) * Ry[j, ix - 1] / exp_mult
                     ) / denom
 
 

diff --git a/fastkde/polyexp_kde.py → polyexpkde/polyexp_kde.py b/fastkde/polyexp_kde.py → polyexpkde/polyexp_kde.py
@@ -9,7 +9,7 @@
     check_random_state,
 )
 
-from fastkde.ksum import (
+from polyexpkde.ksum import (
     kernel_eval,
     kernel_density,
 )

diff --git a/setup.py b/setup.py
@@ -7,8 +7,8 @@ def readme():
 
 
 configuration = {
-    "name": "fastkde",
-    "version": "0.1.1",
+    "name": "polyexpkde",
+    "version": "0.1.2",
     "description": "Fast KDE with polynomial-exponential kernel",
     "long_description": readme(),
     "classifiers": [
@@ -29,12 +29,12 @@ def readme():
         "Programming Language :: Python :: 3.8",
     ],
     "keywords": "KDE, kernel density estimation",
-    "url": "http://github.com/calvinmccarter/fastkde",
+    "url": "http://github.com/calvinmccarter/polyexpkde",
     "author": "Calvin McCarter",
     "author_email": "[email protected]",
     "maintainer": "Calvin McCarter",
     "maintainer_email": "[email protected]",
-    "packages": ["fastkde"],
+    "packages": ["polyexpkde"],
     "install_requires": [
         "numba >= 0.48",
         "numpy",

diff --git a/tests/fastkde/test_basic.py → tests/polyexpkde/test_basic.py b/tests/fastkde/test_basic.py → tests/polyexpkde/test_basic.py
@@ -2,7 +2,7 @@
 import pytest
 import scipy.stats as spst
 
-from fastkde import PolyExpKernelDensity
+from polyexpkde import PolyExpKernelDensity
 
 
 def test_gaussian_approx():

diff --git a/tests/test_imports.py b/tests/test_imports.py
@@ -1,5 +1,5 @@
 def test_import_packages():
     """Test that importing our package works."""
-    import fastkde
-    from fastkde import PolyExpKernelDensity
-    from fastkde import kernel_density
+    import polyexpkde
+    from polyexpkde import PolyExpKernelDensity
+    from polyexpkde import kernel_density