splinekit 0.1.0__tar.gz → 0.2.2__tar.gz

splinekit-0.2.2/.github/workflows/release.yml

@@ -0,0 +1,140 @@
+# splinekit/.github/workflows/release.yml
+
+### Pure Python
+# name: Build & publish to PyPI
+# 
+# on:
+#   push:
+#     tags:
+#       - "v*"
+# 
+# permissions:
+#   contents: read
+#   id-token: write  # OIDC for PyPI
+# 
+# jobs:
+#   build-and-publish:
+#     runs-on: ubuntu-latest
+# 
+#     steps:
+#       - name: Checkout repository
+#         uses: actions/checkout@v4
+# 
+#       - name: Set up Python
+#         uses: actions/setup-python@v5
+#         with:
+#           python-version: "3.12"
+# 
+#       - name: Install build tooling
+#         run: |
+#           python -m pip install -U pip build
+# 
+#       - name: Build sdist and wheel (pure Python)
+#         run: python -m build  # creates dist/*.tar.gz and dist/*.whl
+# 
+#       - name: Publish to PyPI
+#         uses: pypa/gh-action-pypi-publish@release/v1
+
+### Activate for C dynamic library
+name: Build & publish to PyPI
+
+on:
+  push:
+    tags:
+      - "v*"
+
+permissions:
+  contents: read
+  id-token: write  # OIDC for PyPI
+
+jobs:
+  build_wheels:
+    name: Build wheels on ${{ matrix.os }}
+    runs-on: ${{ matrix.os }}
+
+    strategy:
+      fail-fast: false
+      matrix:
+        include:
+          # Linux + Windows
+          - os: ubuntu-latest
+            mac_arch: x86_64   # unused on Linux, but keeps env simple
+          - os: windows-latest
+            mac_arch: x86_64   # unused on Windows, but keeps env simple
+
+          # macOS Intel (x86_64) and Apple Silicon (arm64)
+          - os: macos-15-intel
+            mac_arch: x86_64
+          - os: macos-15   # or macos-latest
+            mac_arch: arm64
+
+    steps:
+      - uses: actions/checkout@v4
+
+      - uses: actions/setup-python@v5
+        with:
+          python-version: "3.12"
+
+      - name: Install build tools
+        run: python -m pip install -U pip cibuildwheel build
+
+      - name: Build wheels
+        run: python -m cibuildwheel --output-dir dist
+        env:
+          # CPython only, for versions we support
+          CIBW_BUILD: "cp311-* cp312-*"
+
+          # Avoid unsupported/undesired targets
+          CIBW_SKIP: "*musllinux*"
+
+          # Architectures
+          CIBW_ARCHS_LINUX: "x86_64"
+          CIBW_ARCHS_WINDOWS: "AMD64"
+          CIBW_ARCHS_MACOS: ${{ matrix.mac_arch }}
+
+          # macOS: ensure wheels target a realistic minimum for C++17
+          # (Intel: 10.13+ matches cp312 minimum; arm64 will still be 11.0+)
+          CIBW_ENVIRONMENT_MACOS: MACOSX_DEPLOYMENT_TARGET=10.13
+
+          # Simple smoke test: import and check native module presence
+          CIBW_TEST_COMMAND: >
+            python -c "import importlib.util, numpy as np, splinekit;
+            print('ok', importlib.util.find_spec('splinekit._padding') is not None)"
+
+      - name: Build sdist (once on Ubuntu)
+        if: matrix.os == 'ubuntu-latest'
+        run: python -m build -s -o dist
+
+      - name: Upload wheel artifacts
+        uses: actions/upload-artifact@v4
+        with:
+          name: wheels-${{ matrix.os }}-${{ matrix.mac_arch }}
+          path: dist/*.whl
+          if-no-files-found: error
+
+      - name: Upload sdist artifact
+        if: matrix.os == 'ubuntu-latest'
+        uses: actions/upload-artifact@v4
+        with:
+          name: sdist
+          path: dist/*.tar.gz
+          if-no-files-found: error
+
+  publish:
+    needs: [build_wheels]
+    runs-on: ubuntu-latest
+
+    steps:
+      - uses: actions/download-artifact@v4
+        with:
+          pattern: wheels-*
+          merge-multiple: true
+          path: dist
+
+      - uses: actions/download-artifact@v4
+        with:
+          name: sdist
+          path: dist
+
+      - name: Publish to PyPI
+        uses: pypa/gh-action-pypi-publish@release/v1

splinekit-0.2.2/CMakeLists.txt

@@ -0,0 +1,36 @@
+cmake_minimum_required(VERSION 3.23)
+project(padding LANGUAGES C)
+
+set(CMAKE_C_STANDARD 99)
+set(CMAKE_C_STANDARD_REQUIRED ON)
+set(CMAKE_C_EXTENSIONS OFF)
+set(CMAKE_POSITION_INDEPENDENT_CODE ON)
+
+if (MSVC)
+  add_compile_options(/W4 /permissive- /Zc:preprocessor /Zc:__cplusplus)
+else()
+  add_compile_options(-Wall -Wextra -Wpedantic)
+endif()
+
+# ---- Absolute source paths (robust on Windows) ----
+set(PADDING_SRC_DIR "${CMAKE_CURRENT_SOURCE_DIR}/ansi_c/spline_padding/src")
+
+add_library(padding SHARED)
+
+target_sources(padding
+  PRIVATE
+    "${PADDING_SRC_DIR}/spline_padding.c"
+  PUBLIC
+    FILE_SET spline_padding
+    TYPE HEADERS
+    FILES
+      "${PADDING_SRC_DIR}/spline_padding.h"
+)
+target_compile_features(padding PUBLIC c_std_99)
+
+install(
+  TARGETS padding
+  LIBRARY DESTINATION splinekit
+  RUNTIME DESTINATION splinekit
+  ARCHIVE DESTINATION splinekit
+)

{splinekit-0.1.0 → splinekit-0.2.2}/PKG-INFO

@@ -1,54 +1,52 @@
 Metadata-Version: 2.4
 Name: splinekit
-Version: 0.1.0
+Version: 0.2.2
 Summary: A set of functions and classes to handle periodic 1D polynomial splines
-Project-URL: Homepage, https://github.com/splinekit/splinekit
-Project-URL: Documentation, https://splinekit.github.io/
-Project-URL: Repository, https://github.com/splinekit/splinekit.git
-Project-URL: Issues, https://github.com/splinekit/splinekit/issues
-Project-URL: Source, https://github.com/splinekit/splinekit
-Project-URL: PyPI, https://pypi.org/project/splinekit/
-Author-email: Philippe Thévenaz <philippe.thevenaz@icloud.com>
-Maintainer-email: Philippe Thévenaz <philippe.thevenaz@icloud.com>
+Author-Email: =?utf-8?q?Philippe_Th=C3=A9venaz?= <philippe.thevenaz@icloud.com>
+Maintainer-Email: =?utf-8?q?Philippe_Th=C3=A9venaz?= <philippe.thevenaz@icloud.com>
 License-Expression: MIT
 License-File: LICENSE
+Classifier: Programming Language :: Python :: 3
+Classifier: Operating System :: OS Independent
 Classifier: Intended Audience :: Developers
 Classifier: Intended Audience :: Education
 Classifier: Intended Audience :: Science/Research
-Classifier: Operating System :: OS Independent
-Classifier: Programming Language :: Python :: 3
 Classifier: Topic :: Scientific/Engineering
 Classifier: Topic :: Scientific/Engineering :: Mathematics
 Classifier: Topic :: Software Development :: Libraries :: Python Modules
+Project-URL: Homepage, https://github.com/splinekit/splinekit
+Project-URL: Documentation, https://splinekit.github.io/
+Project-URL: Repository, https://github.com/splinekit/splinekit.git
+Project-URL: Issues, https://github.com/splinekit/splinekit/issues
+Project-URL: Source, https://github.com/splinekit/splinekit
+Project-URL: PyPI, https://pypi.org/project/splinekit/
 Requires-Python: >=3.11
-Requires-Dist: matplotlib>=3.8
 Requires-Dist: numpy>=1.26
 Requires-Dist: scipy>=1.15
 Requires-Dist: sympy>=1.12
-Provides-Extra: dev
-Requires-Dist: build; extra == 'dev'
-Requires-Dist: ipykernel>=7.1.0; extra == 'dev'
-Requires-Dist: ipython>=8.26; extra == 'dev'
-Requires-Dist: ipywidgets>=8; extra == 'dev'
-Requires-Dist: mypy>=1.8; extra == 'dev'
-Requires-Dist: pytest>=6.0; extra == 'dev'
-Requires-Dist: tox>=4.13; extra == 'dev'
+Requires-Dist: matplotlib>=3.8
+Provides-Extra: visualization
+Requires-Dist: IPython>=8.26; extra == "visualization"
 Provides-Extra: docs
-Requires-Dist: ipython>=8.26; extra == 'docs'
-Requires-Dist: myst-parser>=3.0; extra == 'docs'
-Requires-Dist: nbsphinx>=0.9; extra == 'docs'
-Requires-Dist: sphinx-copybutton>=0.5; extra == 'docs'
-Requires-Dist: sphinx-design>=0.6; extra == 'docs'
-Requires-Dist: sphinx-gallery>=0.17; extra == 'docs'
-Requires-Dist: sphinx-prompt>=1.8; extra == 'docs'
-Requires-Dist: sphinx-remove-toctrees>=1.0; extra == 'docs'
-Requires-Dist: sphinx-rtd-theme>=3.0.2; extra == 'docs'
-Requires-Dist: sphinx>=7.3; extra == 'docs'
+Requires-Dist: splinekit[visualization]; extra == "docs"
+Requires-Dist: sphinx>=7.3; extra == "docs"
+Requires-Dist: sphinx-gallery>=0.17; extra == "docs"
+Requires-Dist: sphinx-prompt>=1.8; extra == "docs"
+Requires-Dist: sphinx-copybutton>=0.5; extra == "docs"
+Requires-Dist: sphinx-remove-toctrees>=1.0; extra == "docs"
+Requires-Dist: sphinx-design>=0.6; extra == "docs"
+Requires-Dist: sphinx-rtd-theme>=3.0.2; extra == "docs"
+Requires-Dist: myst-parser>=3.0; extra == "docs"
+Requires-Dist: nbsphinx>=0.9; extra == "docs"
 Provides-Extra: tests
-Requires-Dist: pytest>=6.0; extra == 'tests'
-Requires-Dist: tox>=4.13; extra == 'tests'
-Provides-Extra: visualization
-Requires-Dist: ipython>=8.26; extra == 'visualization'
+Requires-Dist: tox>=4.13; extra == "tests"
+Requires-Dist: pytest>=6.0; extra == "tests"
+Provides-Extra: dev
+Requires-Dist: splinekit[tests,visualization]; extra == "dev"
+Requires-Dist: mypy>=1.8; extra == "dev"
+Requires-Dist: ipykernel>=7.1.0; extra == "dev"
+Requires-Dist: ipywidgets>=8; extra == "dev"
+Requires-Dist: build; extra == "dev"
 Description-Content-Type: text/markdown
 
 <!-- splinekit/README.md -->
@@ -111,3 +109,14 @@ Install `splinekit` development environment in editable mode
 ```shell
 pip install -e .[dev]
 ```
+
+## Release Notes
+
+### 0.1.0
+First stable version. Pure Python.
+
+### 0.2.1
+Same API as 0.1.0. Ansi-C dynamic library.
+
+### 0.2.2
+Some demos added

{splinekit-0.1.0 → splinekit-0.2.2}/README.md

@@ -58,3 +58,14 @@ Install `splinekit` development environment in editable mode
 ```shell
 pip install -e .[dev]
 ```
+
+## Release Notes
+
+### 0.1.0
+First stable version. Pure Python.
+
+### 0.2.1
+Same API as 0.1.0. Ansi-C dynamic library.
+
+### 0.2.2
+Some demos added

splinekit-0.2.2/ansi_c/spline_padding/src/spline_padding.c

@@ -0,0 +1,275 @@
+#include "spline_padding.h"
+
+extern void get_samples_to_coeff_p (
+	double data[],
+	int32_t dataCount,
+	double pole[],
+	int32_t poleCount
+) {
+	if ((NULL == data) || (0 >= dataCount)) {
+		return;
+	}
+	if ((NULL == pole) || (0 >= poleCount)) {
+		return;
+	}
+	if ((0 != poleCount) && (1 != dataCount)) {
+		for (int q = 0; q < poleCount; q++) {
+			const double z = pole[q];
+			double sigma = data[0];
+			double zeta = z;
+			int k = 1;
+			while ((k < dataCount) && (0.0 != zeta)) {
+				sigma += zeta * data[dataCount - k];
+				zeta *= z;
+				k++;
+			}
+			data[0] = sigma / (1.0 - zeta);
+			for (int k = 1; k < dataCount; k++) {
+				data[k] += z * data[k - 1];
+			}
+			sigma = data[dataCount - 1];
+			zeta = z;
+			k = 0;
+			while ((k < (dataCount - 1)) && (0.0 != zeta)) {
+				sigma += zeta * data[k];
+				zeta *= z;
+				k++;
+			}
+			const double z12 = (1.0 - z) * (1.0 - z);
+			data[dataCount - 1] = z12 * sigma / (1.0 - zeta);
+			for (int k = 1; k < dataCount; k++) {
+				data[dataCount - 1 - k] = z * data[dataCount - k] + z12 * data[dataCount - 1 - k];
+			}
+		}
+	}
+} /* end get_samples_to_coeff_p */
+
+extern void get_samples_to_coeff_n (
+	double data[],
+	int32_t dataCount,
+	double pole[],
+	int32_t poleCount
+) {
+	if ((NULL == data) || (0 >= dataCount)) {
+		return;
+	}
+	if ((NULL == pole) || (0 >= poleCount)) {
+		return;
+	}
+	if ((0 != poleCount) && (1 != dataCount)) {
+		for (int q = 0; q < poleCount; q++) {
+			const double z = pole[q];
+			double sigma1 = data[0];
+			double sigma2 = data[dataCount - 1];
+			double zeta = z;
+			int k = 1;
+			while ((k < (dataCount - 1)) && (0.0 != zeta)) {
+				sigma1 += zeta * data[k];
+				sigma2 += zeta * data[dataCount - 1 - k];
+				zeta *= z;
+				k++;
+			}
+			data[0] = (sigma1 + zeta * sigma2) / (1.0 - zeta * zeta);
+			for (int k = 1; k < dataCount; k++) {
+				data[k] += z * data[k - 1];
+			}
+			const double z12 = (1.0 - z) * (1.0 - z);
+			data[dataCount - 1] = z12 * (z * data[dataCount - 2] + data[dataCount - 1]) / (1.0 - z * z);
+			for (int k = 1; k < dataCount; k++) {
+				data[dataCount - 1 - k] = z * data[dataCount - k] + z12 * data[dataCount - 1 - k];
+			}
+		}
+	}
+} /* end get_samples_to_coeff_n */
+
+extern void get_samples_to_coeff_w (
+	double data[],
+	int32_t dataCount,
+	double pole[],
+	int32_t poleCount
+) {
+	if ((NULL == data) || (0 >= dataCount)) {
+		return;
+	}
+	if ((NULL == pole) || (0 >= poleCount)) {
+		return;
+	}
+	if ((0 != poleCount) && (1 != dataCount)) {
+		for (int q = 0; q < poleCount; q++) {
+			const double z = pole[q];
+			double sigma1 = 0.0;
+			double sigma2 = 0.0;
+			double zeta = 1.0;
+			int k = 0;
+			while ((k < dataCount) && (0.0 != zeta)) {
+				sigma1 += zeta * data[k];
+				sigma2 += zeta * data[dataCount - 1 - k];
+				zeta *= z;
+				k++;
+			}
+			data[0] += z * (sigma1 + zeta * sigma2) / (1.0 - zeta * zeta);
+			for (int k = 1; k < dataCount; k++) {
+				data[k] += z * data[k - 1];
+			}
+			const double z12 = (1.0 - z) * (1.0 - z);
+			data[dataCount - 1] *= 1.0 - z;
+			for (int k = 1; k < dataCount; k++) {
+				data[dataCount - 1 - k] = z * data[dataCount - k] + z12 * data[dataCount - 1 - k];
+			}
+		}
+	}
+} /* end get_samples_to_coeff_w */
+
+extern void get_samples_to_coeff_a (
+	double data[],
+	int32_t dataCount,
+	double pole[],
+	int32_t poleCount
+) {
+	if ((NULL == data) || (0 >= dataCount)) {
+		return;
+	}
+	if ((NULL == pole) || (0 >= poleCount)) {
+		return;
+	}
+	if ((0 != poleCount) && (1 != dataCount)) {
+		for (int q = 0; q < poleCount; q++) {
+			const double z = pole[q];
+			double sigma1 = 0.0;
+			double sigma2 = 0.0;
+			double zeta = z;
+			int k = 1;
+			while ((k < (dataCount - 1)) && (0.0 != zeta)) {
+				sigma1 += zeta * data[k];
+				sigma2 += zeta * data[dataCount - 1 - k];
+				zeta *= z;
+				k++;
+			}
+			data[0] = ((data[0] - zeta * data[dataCount - 1]) * (1.0 + z) / (1.0 - z) - sigma1 + zeta * sigma2) / (1.0 - zeta * zeta);
+			for (int k = 1; k < dataCount; k++) {
+				data[k] += z * data[k - 1];
+			}
+			const double z12 = (1.0 - z) * (1.0 - z);
+			data[dataCount - 1] -= z * data[dataCount - 2];
+			for (int k = 1; k < dataCount; k++) {
+				data[dataCount - 1 - k] = z * data[dataCount - k] + z12 * data[dataCount - 1 - k];
+			}
+		}
+	}
+} /* end get_samples_to_coeff_a */
+
+extern void get_samples_to_coeff_np (
+	double data[],
+	int32_t dataCount,
+	double pole[],
+	int32_t poleCount
+) {
+	if ((NULL == data) || (0 >= dataCount)) {
+		return;
+	}
+	if ((NULL == pole) || (0 >= poleCount)) {
+		return;
+	}
+	for (int q = 0; q < poleCount; q++) {
+		const double z = pole[q];
+		double sigma1 = data[0];
+		double sigma2 = data[dataCount - 1];
+		double zeta = z;
+		int k = 1;
+		while ((k < dataCount) && (0.0 != zeta)) {
+			sigma1 += zeta * data[k];
+			sigma2 += zeta * data[dataCount - 1 - k];
+			zeta *= z;
+			k++;
+		}
+		const double Z0 = z / (1.0 + zeta);
+		sigma1 *= Z0;
+		sigma2 *= Z0;
+		data[0] -= sigma2;
+		for (int k = 1; k < dataCount; k++) {
+			data[k] += z * data[k - 1];
+		}
+		zeta *= zeta;
+		const double z12 = (1.0 - z) * (1.0 - z);
+		const double Z1 = (1.0 - z) / (1.0 + z);
+		data[dataCount - 1] *= (1.0 + zeta) * Z1;
+		data[dataCount - 1] -= (sigma2 * zeta / z + sigma1) * Z1;
+		for (int k = 1; k < dataCount; k++) {
+			data[dataCount - 1 - k] = z * data[dataCount - k] + z12 * data[dataCount - 1 - k];
+		}
+	}
+} /* end get_samples_to_coeff_np */
+
+extern void get_samples_to_coeff_nn (
+	double data[],
+	int32_t dataCount,
+	double pole[],
+	int32_t poleCount
+) {
+	if ((NULL == data) || (0 >= dataCount)) {
+		return;
+	}
+	if ((NULL == pole) || (0 >= poleCount)) {
+		return;
+	}
+	for (int q = 0; q < poleCount; q++) {
+		const double z = pole[q];
+		double sigma1 = 0.0;
+		double sigma2 = 0.0;
+		double zeta = 1.0;
+		int k = 0;
+		while ((k < dataCount) && (0.0 != zeta)) {
+			sigma1 += zeta * data[k];
+			sigma2 += zeta * data[dataCount - 1 - k];
+			zeta *= z;
+			k++;
+		}
+		zeta *= z;
+		data[0] -= (sigma1 - zeta * sigma2) * z * z / (1.0 - zeta * zeta);
+		for (int k = 1; k < dataCount; k++) {
+			data[k] += z * data[k - 1];
+		}
+		const double z12 = (1.0 - z) * (1.0 - z);
+		data[dataCount - 1] *= z12;
+		for (int k = 1; k < dataCount; k++) {
+			data[dataCount - 1 - k] = z * data[dataCount - k] + z12 * data[dataCount - 1 - k];
+		}
+	}
+} /* end get_samples_to_coeff_nn */
+
+extern void get_samples_to_coeff_nw (
+	double data[],
+	int32_t dataCount,
+	double pole[],
+	int32_t poleCount
+) {
+	if ((NULL == data) || (0 >= dataCount)) {
+		return;
+	}
+	if ((NULL == pole) || (0 >= poleCount)) {
+		return;
+	}
+	for (int q = 0; q < poleCount; q++) {
+		const double z = pole[q];
+		double sigma1 = 0.0;
+		double sigma2 = 0.0;
+		double zeta = 1.0;
+		int k = 0;
+		while ((k < dataCount) && (0.0 != zeta)) {
+			sigma1 += zeta * data[k];
+			sigma2 += zeta * data[dataCount - 1 - k];
+			zeta *= z;
+			k++;
+		}
+		data[0] -= (sigma1 - zeta * sigma2) * z / (1.0 - zeta * zeta);
+		for (int k = 1; k < dataCount; k++) {
+			data[k] += z * data[k - 1];
+		}
+		const double z12 = (1.0 - z) * (1.0 - z);
+		data[dataCount - 1] *= z12 / (1.0 + z);
+		for (int k = 1; k < dataCount; k++) {
+			data[dataCount - 1 - k] = z * data[dataCount - k] + z12 * data[dataCount - 1 - k];
+		}
+	}
+} /* end get_samples_to_coeff_nw */
+

splinekit-0.2.2/ansi_c/spline_padding/src/spline_padding.h

@@ -0,0 +1,57 @@
+#ifndef spline_padding_h
+#define spline_padding_h
+
+#include <stdint.h>
+#include <stdio.h>
+#include <stdlib.h>
+
+extern void get_samples_to_coeff_p (
+	double *data,
+	int32_t dataCount,
+	double *pole,
+	int32_t poleCount
+);
+
+extern void get_samples_to_coeff_n (
+	double *data,
+	int32_t dataCount,
+	double *pole,
+	int32_t poleCount
+);
+
+extern void get_samples_to_coeff_w (
+	double *data,
+	int32_t dataCount,
+	double *pole,
+	int32_t poleCount
+);
+
+extern void get_samples_to_coeff_a (
+	double *data,
+	int32_t dataCount,
+	double *pole,
+	int32_t poleCount
+);
+
+extern void get_samples_to_coeff_np (
+	double *data,
+	int32_t dataCount,
+	double *pole,
+	int32_t poleCount
+);
+
+extern void get_samples_to_coeff_nn (
+	double *data,
+	int32_t dataCount,
+	double *pole,
+	int32_t poleCount
+);
+
+extern void get_samples_to_coeff_nw (
+	double *data,
+	int32_t dataCount,
+	double *pole,
+	int32_t poleCount
+);
+
+#endif /* spline_padding_h */

splinekit-0.2.2/dev/ANSI-C/spline_padding_cmake/CMakeLists.txt

@@ -0,0 +1,15 @@
+cmake_minimum_required(VERSION 3.23)
+
+project(Padding)
+
+add_library(Padding SHARED)
+
+target_sources(Padding
+  PRIVATE
+    spline_padding.c
+  PUBLIC
+    FILE_SET spline_padding
+    TYPE HEADERS
+    FILES
+      spline_padding.h
+)