interpn 0.2.0__tar.gz → 0.2.1__tar.gz

{interpn-0.2.0 → interpn-0.2.1}/Cargo.lock

@@ -43,9 +43,9 @@ checksum = "1e186cfbae8084e513daff4240b4797e342f988cecda4fb6c939150f96315fd8"
 
 [[package]]
 name = "interpn"
-version = "0.4.1"
+version = "0.4.2"
 source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "82ec859d62d17054a7314d5bb47d7e4eeb4648b87b7c1e0e33de1d4df21e0375"
+checksum = "f37cc48f61030a3bed5eb01e583af8ae1248f3e3bf63d05d87a18f27a98a426b"
 dependencies = [
  "num-traits",
 ]
@@ -177,6 +177,12 @@ dependencies = [
  "windows-targets",
 ]
 
+[[package]]
+name = "portable-atomic"
+version = "1.6.0"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "7170ef9988bc169ba16dd36a7fa041e5c4cbeb6a35b76d4c03daded371eae7c0"
+
 [[package]]
 name = "proc-macro2"
 version = "1.0.69"
@@ -188,15 +194,16 @@ dependencies = [
 
 [[package]]
 name = "pyo3"
-version = "0.20.0"
+version = "0.20.3"
 source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "04e8453b658fe480c3e70c8ed4e3d3ec33eb74988bd186561b0cc66b85c3bc4b"
+checksum = "53bdbb96d49157e65d45cc287af5f32ffadd5f4761438b527b055fb0d4bb8233"
 dependencies = [
  "cfg-if",
  "indoc",
  "libc",
  "memoffset",
  "parking_lot",
+ "portable-atomic",
  "pyo3-build-config",
  "pyo3-ffi",
  "pyo3-macros",
@@ -205,9 +212,9 @@ dependencies = [
 
 [[package]]
 name = "pyo3-build-config"
-version = "0.20.0"
+version = "0.20.3"
 source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "a96fe70b176a89cff78f2fa7b3c930081e163d5379b4dcdf993e3ae29ca662e5"
+checksum = "deaa5745de3f5231ce10517a1f5dd97d53e5a2fd77aa6b5842292085831d48d7"
 dependencies = [
  "once_cell",
  "target-lexicon",
@@ -215,9 +222,9 @@ dependencies = [
 
 [[package]]
 name = "pyo3-ffi"
-version = "0.20.0"
+version = "0.20.3"
 source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "214929900fd25e6604661ed9cf349727c8920d47deff196c4e28165a6ef2a96b"
+checksum = "62b42531d03e08d4ef1f6e85a2ed422eb678b8cd62b762e53891c05faf0d4afa"
 dependencies = [
  "libc",
  "pyo3-build-config",
@@ -225,9 +232,9 @@ dependencies = [
 
 [[package]]
 name = "pyo3-macros"
-version = "0.20.0"
+version = "0.20.3"
 source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "dac53072f717aa1bfa4db832b39de8c875b7c7af4f4a6fe93cdbf9264cf8383b"
+checksum = "7305c720fa01b8055ec95e484a6eca7a83c841267f0dd5280f0c8b8551d2c158"
 dependencies = [
  "proc-macro2",
  "pyo3-macros-backend",
@@ -237,12 +244,13 @@ dependencies = [
 
 [[package]]
 name = "pyo3-macros-backend"
-version = "0.20.0"
+version = "0.20.3"
 source = "registry+https://github.com/rust-lang/crates.io-index"
-checksum = "7774b5a8282bd4f25f803b1f0d945120be959a36c72e08e7cd031c792fdfd424"
+checksum = "7c7e9b68bb9c3149c5b0cade5d07f953d6d125eb4337723c4ccdb665f1f96185"
 dependencies = [
  "heck",
  "proc-macro2",
+ "pyo3-build-config",
  "quote",
  "syn",
 ]

{interpn-0.2.0 → interpn-0.2.1}/Cargo.toml

@@ -9,9 +9,9 @@ name = "_interpn"
 crate-type = ["cdylib"]
 
 [dependencies]
-pyo3 = "0.20.0"
+pyo3 = "0.20.3"
 numpy = "0.20.0"
-interpn = {version = "0.4.1", default_features = false}
+interpn = {version = "0.4.2", default_features = false}
 #interpn = { path = "../interpn/interpn" }
 
 [profile.release]

{interpn-0.2.0 → interpn-0.2.1}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.3
 Name: interpn
-Version: 0.2.0
+Version: 0.2.1
 Classifier: Programming Language :: Rust
 Classifier: Programming Language :: Python :: Implementation :: CPython
 Classifier: Programming Language :: Python :: Implementation :: PyPy
@@ -42,7 +42,7 @@ Python bindings to the `interpn` Rust library for N-dimensional interpolation an
 | Feature →<br>↓ Interpolant Method | Regular<br>Grid | Rectilinear<br>Grid | Json<br>Serialization |
 |-----------------------------------|-----------------|---------------------|-----------------------|
 | Linear                            |   ✅            |     ✅              | ✅                    |
-| Cubic                             |   ✅            |     💡              | ✅                    |
+| Cubic                             |   ✅            |     ✅              | ✅                    |
 
 The methods provided here, while more limited in scope than scipy's, are
 * significantly faster for higher dimensions (1-3 orders of magnitude under most conditions)
@@ -50,18 +50,54 @@ The methods provided here, while more limited in scope than scipy's, are
 * produce significantly improved floating-point error (by 1-2 orders of magnitude)
 * are json-serializable using Pydantic
 * can also be used easily in web and embedded applications via the Rust library
+* are permissively licensed
+
+![ND throughput 1 obs](./docs/throughput_vs_dims_1_obs.svg)
 
 See [here](https://interpnpy.readthedocs.io/en/latest/perf/) for more info about quality-of-fit, throughput, and memory usage.
 
+## Installation
+
+```bash
+pip install interpn
+```
+
+## Example: Available Methods
+
+```python
+import interpn
+import numpy as np
+
+# Build grid
+x = np.linspace(0.0, 10.0, 5)
+y = np.linspace(20.0, 30.0, 4)
+grids = [x, y]
+
+xgrid, ygrid = np.meshgrid(x, y, indexing="ij")
+zgrid = (xgrid + 2.0 * ygrid)  # Values at grid points
+
+# Grid inputs for true regular grid
+dims = [x.size, y.size]
+starts = np.array([x[0], y[0]])
+steps = np.array([x[1] - x[0], y[1] - y[0]])
+
+# Initialize different interpolators
+# Call like `linear_regular.eval([xs, ys])`
+linear_regular = interpn.MultilinearRegular.new(dims, starts, steps, zgrid)
+cubic_regular = interpn.MulticubicRegular.new(dims, starts, steps, zgrid)
+linear_rectilinear = interpn.MultilinearRectilinear.new(grids, zgrid)
+cubic_rectilinear = interpn.MulticubicRectilinear.new(grids, zgrid)
+```
 
 ## Example: Multilinear Interpolation on a Regular Grid
+
 ```python
 import interpn
 import numpy as np
 
 # Build grid
 x = np.linspace(0.0, 10.0, 5)
-y = np.linspace(20.0, 30.0, 3)
+y = np.linspace(20.0, 30.0, 4)
 
 xgrid, ygrid = np.meshgrid(x, y, indexing="ij")
 zgrid = (xgrid + 2.0 * ygrid)  # Values at grid points
@@ -75,15 +111,13 @@ steps = np.array([x[1] - x[0], y[1] - y[0]])
 obs = [xgrid.flatten(), ygrid.flatten()]
 
 # Initialize
-interpolator = interpn.MultilinearRegular.new(
-    dims, starts, steps, zgrid.flatten()
-)
+interpolator = interpn.MultilinearRegular.new(dims, starts, steps, zgrid.flatten())
 
 # Interpolate
 out = interpolator.eval(obs)
 
 # Check result
-assert np.all(out == zgrid.flatten())
+assert np.allclose(out, zgrid.flatten(), rtol=1e-13)
 
 # Serialize and deserialize
 roundtrip_interpolator = interpn.MultilinearRegular.model_validate_json(

{interpn-0.2.0 → interpn-0.2.1}/README.md

@@ -11,7 +11,7 @@ Python bindings to the `interpn` Rust library for N-dimensional interpolation an
 | Feature →<br>↓ Interpolant Method | Regular<br>Grid | Rectilinear<br>Grid | Json<br>Serialization |
 |-----------------------------------|-----------------|---------------------|-----------------------|
 | Linear                            |   ✅            |     ✅              | ✅                    |
-| Cubic                             |   ✅            |     💡              | ✅                    |
+| Cubic                             |   ✅            |     ✅              | ✅                    |
 
 The methods provided here, while more limited in scope than scipy's, are
 * significantly faster for higher dimensions (1-3 orders of magnitude under most conditions)
@@ -19,18 +19,54 @@ The methods provided here, while more limited in scope than scipy's, are
 * produce significantly improved floating-point error (by 1-2 orders of magnitude)
 * are json-serializable using Pydantic
 * can also be used easily in web and embedded applications via the Rust library
+* are permissively licensed
+
+![ND throughput 1 obs](./docs/throughput_vs_dims_1_obs.svg)
 
 See [here](https://interpnpy.readthedocs.io/en/latest/perf/) for more info about quality-of-fit, throughput, and memory usage.
 
+## Installation
+
+```bash
+pip install interpn
+```
+
+## Example: Available Methods
+
+```python
+import interpn
+import numpy as np
+
+# Build grid
+x = np.linspace(0.0, 10.0, 5)
+y = np.linspace(20.0, 30.0, 4)
+grids = [x, y]
+
+xgrid, ygrid = np.meshgrid(x, y, indexing="ij")
+zgrid = (xgrid + 2.0 * ygrid)  # Values at grid points
+
+# Grid inputs for true regular grid
+dims = [x.size, y.size]
+starts = np.array([x[0], y[0]])
+steps = np.array([x[1] - x[0], y[1] - y[0]])
+
+# Initialize different interpolators
+# Call like `linear_regular.eval([xs, ys])`
+linear_regular = interpn.MultilinearRegular.new(dims, starts, steps, zgrid)
+cubic_regular = interpn.MulticubicRegular.new(dims, starts, steps, zgrid)
+linear_rectilinear = interpn.MultilinearRectilinear.new(grids, zgrid)
+cubic_rectilinear = interpn.MulticubicRectilinear.new(grids, zgrid)
+```
 
 ## Example: Multilinear Interpolation on a Regular Grid
+
 ```python
 import interpn
 import numpy as np
 
 # Build grid
 x = np.linspace(0.0, 10.0, 5)
-y = np.linspace(20.0, 30.0, 3)
+y = np.linspace(20.0, 30.0, 4)
 
 xgrid, ygrid = np.meshgrid(x, y, indexing="ij")
 zgrid = (xgrid + 2.0 * ygrid)  # Values at grid points
@@ -44,15 +80,13 @@ steps = np.array([x[1] - x[0], y[1] - y[0]])
 obs = [xgrid.flatten(), ygrid.flatten()]
 
 # Initialize
-interpolator = interpn.MultilinearRegular.new(
-    dims, starts, steps, zgrid.flatten()
-)
+interpolator = interpn.MultilinearRegular.new(dims, starts, steps, zgrid.flatten())
 
 # Interpolate
 out = interpolator.eval(obs)
 
 # Check result
-assert np.all(out == zgrid.flatten())
+assert np.allclose(out, zgrid.flatten(), rtol=1e-13)
 
 # Serialize and deserialize
 roundtrip_interpolator = interpn.MultilinearRegular.model_validate_json(