splinekit 0.2.7__tar.gz → 0.2.10__tar.gz

{splinekit-0.2.7 → splinekit-0.2.10}/.gitignore

@@ -94,6 +94,9 @@ target/
 # Jupyter Notebook
 .ipynb_checkpoints
 
+# Jupyterlite
+_output/
+
 # IPython
 profile_default/
 ipython_config.py

{splinekit-0.2.7 → splinekit-0.2.10}/PKG-INFO

@@ -1,9 +1,9 @@
 Metadata-Version: 2.4
 Name: splinekit
-Version: 0.2.7
+Version: 0.2.10
 Summary: A set of functions and classes to handle periodic 1D polynomial splines
-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>
+Author-Email: Philippe Thévenaz <philippe.thevenaz@icloud.com>
+Maintainer-Email: Philippe Thévenaz <philippe.thevenaz@icloud.com>
 License-Expression: MIT
 License-File: LICENSE
 Classifier: Programming Language :: Python :: 3
@@ -51,7 +51,7 @@ Description-Content-Type: text/markdown
 
 <!-- splinekit/README.md -->
 
-# Version 0.2.7
+# Version 0.2.10
 
 # splinekit: Spline Operations
 `splinekit` is a Python-based open-source software library aimed at the manipulation of one-dimensional periodic splines.
@@ -116,7 +116,7 @@ pip install -e .[dev]
 
 ## Release Notes
 
-### 0.1.0
+### 0.1.*
 First stable version. Pure Python.
 
 ### 0.2.1
@@ -133,3 +133,9 @@ Some demos added. Documentation minor updates.
 
 ### 0.2.7
 Upgraded to jupyterlite 0.7
+
+### 0.2.8
+Some demos added.
+
+### 0.2.9 — 0.2.10
+All demos added. Operators. Faster execution.

{splinekit-0.2.7 → splinekit-0.2.10}/README.md

@@ -1,6 +1,6 @@
 <!-- splinekit/README.md -->
 
-# Version 0.2.7
+# Version 0.2.10
 
 # splinekit: Spline Operations
 `splinekit` is a Python-based open-source software library aimed at the manipulation of one-dimensional periodic splines.
@@ -65,7 +65,7 @@ pip install -e .[dev]
 
 ## Release Notes
 
-### 0.1.0
+### 0.1.*
 First stable version. Pure Python.
 
 ### 0.2.1
@@ -82,3 +82,9 @@ Some demos added. Documentation minor updates.
 
 ### 0.2.7
 Upgraded to jupyterlite 0.7
+
+### 0.2.8
+Some demos added.
+
+### 0.2.9 — 0.2.10
+All demos added. Operators. Faster execution.

{splinekit-0.2.7 → splinekit-0.2.10/dev/Accelerated}/pyproject.toml

@@ -10,7 +10,7 @@ build-backend = "scikit_build_core.build" # Activate for C dynamic library
 # Project
 [project]
 name = "splinekit"
-version = "0.2.7"
+version = "0.2.10"
 dependencies = [
     "numpy >=1.26",
     "scipy >=1.15",

splinekit-0.2.10/dev/Jupyter/periodic-spline/addition/spline_add_prjct.ipynb

@@ -0,0 +1,173 @@
+{
+ "cells": [
+  {
+   "cell_type": "markdown",
+   "id": "006cebf8-f66b-48aa-9f8b-0840d37ca07d",
+   "metadata": {},
+   "source": [
+    "# Addition of Splines of Identical Period\n",
+    "We display the augend in the <span style=\"color:#1f77b4\">**azure**</span> color; it is the realization of a random periodic polynomial spline of a specified degree $n_{1}$ and delay $\\delta x_{1}.$ For some other (independent) degree $n_{2}$ and delay $\\delta x_{2},$ we display the addend in the <span style=\"color:#17becf\">**aqua**</span> color; it is another realization of a random periodic polynomial spline. We then plot as a thick curve in the <span style=\"color:#8c564b\">**sepia**</span> color the spline of specified degree $n$ and delay $\\delta x$ that best approximates the addition of the addend to the augend. The augend, the addend, and the approximation of their sum have the same period, but may not share the location of their knots."
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "b3185aea-65df-4ca8-8011-efe4d523876c",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "# Load the required libraries\n",
+    "import ipywidgets as widgets\n",
+    "import matplotlib.pyplot as plt\n",
+    "import numpy as np\n",
+    "\n",
+    "import splinekit as sk # This library\n",
+    "\n",
+    "# Setup\n",
+    "max_period = 15 # Maximal period\n",
+    "max_degree = 9 # Maximal spline degree\n",
+    "max_delay = 3.0 # Maximal absolute delay\n",
+    "\n",
+    "# Initialize the generator of random numbers\n",
+    "rng = np.random.default_rng()\n",
+    "\n",
+    "# Initial random periodic splines\n",
+    "f1 = sk.PeriodicSpline1D.from_spline_coeff(rng.standard_normal(6), degree = 2)\n",
+    "f2 = sk.PeriodicSpline1D.from_spline_coeff(rng.standard_normal(6), degree = 3)\n",
+    "\n",
+    "# Plot\n",
+    "def update_plot (\n",
+    "    period = 6,\n",
+    "    degree_1 = 2,\n",
+    "    delay_1 = 0.0,\n",
+    "    degree_2 = 3,\n",
+    "    delay_2 = 0.0,\n",
+    "    degree = 1,\n",
+    "    delay = 0.0\n",
+    "):\n",
+    "    global f1\n",
+    "    global f2\n",
+    "\n",
+    "    # Update of the splines\n",
+    "    if f1.period != period:\n",
+    "        f1 = sk.PeriodicSpline1D.from_spline_coeff(\n",
+    "            rng.standard_normal(period),\n",
+    "            degree = f1.degree\n",
+    "        )\n",
+    "    f1.degree = degree_1\n",
+    "    f1.delay = delay_1\n",
+    "    if f2.period != period:\n",
+    "        f2 = sk.PeriodicSpline1D.from_spline_coeff(\n",
+    "            rng.standard_normal(period),\n",
+    "            degree = f2.degree\n",
+    "        )\n",
+    "    f2.degree = degree_2\n",
+    "    f2.delay = delay_2\n",
+    "\n",
+    "    # Addition\n",
+    "    f = sk.PeriodicSpline1D.add_projected(f1, f2, degree = degree, delay = delay)\n",
+    "\n",
+    "    # Dynamic range\n",
+    "    image = {f1.image(), f2.image(), f.image()}\n",
+    "    plotrange = sk.interval.Interval.enclosure(image)\n",
+    "    plotrange = sk.interval.Closed((\n",
+    "        plotrange.midpoint - 0.55 * plotrange.diameter,\n",
+    "        plotrange.midpoint + 0.55 * plotrange.diameter\n",
+    "    ))\n",
+    "\n",
+    "    # Plots\n",
+    "    subplot = plt.subplots()\n",
+    "    # Augend\n",
+    "    f1.plot(subplot, plotrange = plotrange, plotpoints = 200 + 1)\n",
+    "    # Addend\n",
+    "    f2.plot(\n",
+    "        subplot,\n",
+    "        plotrange = plotrange,\n",
+    "        plotpoints = 200 + 1,\n",
+    "        curve_fmt = \"-C9\",\n",
+    "        curve_markerfmt = \"oC9\",\n",
+    "    )\n",
+    "    # Sum\n",
+    "    f.plot(\n",
+    "        subplot,\n",
+    "        plotrange = plotrange,\n",
+    "        plotpoints = 200 + 1,\n",
+    "        curve_fmt = \"-C5\",\n",
+    "        curve_markerfmt = \"oC5\",\n",
+    "        curve_lw = 2.0\n",
+    "    )\n",
+    "\n",
+    "    # Final display\n",
+    "    plt.show()\n",
+    "\n",
+    "widgets.interactive(\n",
+    "    update_plot,\n",
+    "    period = (1, max_period),\n",
+    "    degree_1 = (0, max_degree),\n",
+    "    delay_1 = (-max_delay, max_delay),\n",
+    "    degree_2 = (0, max_degree),\n",
+    "    delay_2 = (-max_delay, max_delay),\n",
+    "    degree = (0, max_degree),\n",
+    "    delay = (-max_delay, max_delay)\n",
+    ")\n"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "b5490a8b-70d3-4b94-a4ad-ea140f8cd9de",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n"
+   ]
+  }
+ ],
+ "metadata": {
+  "kernelspec": {
+   "display_name": "Python 3 (ipykernel)",
+   "language": "python",
+   "name": "python3"
+  },
+  "language_info": {
+   "codemirror_mode": {
+    "name": "ipython",
+    "version": 3
+   },
+   "file_extension": ".py",
+   "mimetype": "text/x-python",
+   "name": "python",
+   "nbconvert_exporter": "python",
+   "pygments_lexer": "ipython3",
+   "version": "3.12.6"
+  },
+  "widgets": {
+   "application/vnd.jupyter.widget-state+json": {
+    "state": {},
+    "version_major": 2,
+    "version_minor": 0
+   }
+  }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 5
+}

splinekit-0.2.10/dev/Jupyter/periodic-spline/addition/spline_addition.ipynb

@@ -0,0 +1,165 @@
+{
+ "cells": [
+  {
+   "cell_type": "markdown",
+   "id": "006cebf8-f66b-48aa-9f8b-0840d37ca07d",
+   "metadata": {},
+   "source": [
+    "# Addition of Two Compatible Splines\n",
+    "We display the augend in the <span style=\"color:#1f77b4\">**azure**</span> color; it is the realization of a random periodic polynomial spline of a specified period, degree, and delay. For the same period, degree, and delay, we display the addend in the <span style=\"color:#17becf\">**aqua**</span> color; it is another realization of a random periodic polynomial spline. We then plot as a thick curve in the <span style=\"color:#8c564b\">**sepia**</span> color the sum resulting from the addition of the addend to the augend."
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "b3185aea-65df-4ca8-8011-efe4d523876c",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "# Load the required libraries\n",
+    "import ipywidgets as widgets\n",
+    "import matplotlib.pyplot as plt\n",
+    "import numpy as np\n",
+    "\n",
+    "import splinekit as sk # This library\n",
+    "\n",
+    "# Setup\n",
+    "max_period = 15 # Maximal period\n",
+    "max_degree = 9 # Maximal spline degree\n",
+    "max_delay = 3.0 # Maximal absolute delay\n",
+    "\n",
+    "# Initialize the generator of random numbers\n",
+    "rng = np.random.default_rng()\n",
+    "\n",
+    "# Initial random periodic cubic splines\n",
+    "f1 = sk.PeriodicSpline1D.from_spline_coeff(rng.standard_normal(6), degree = 3)\n",
+    "f2 = sk.PeriodicSpline1D.from_spline_coeff(rng.standard_normal(6), degree = 3)\n",
+    "\n",
+    "# Plot\n",
+    "def update_plot (\n",
+    "    period = 6,\n",
+    "    degree = 3,\n",
+    "    delay = 0.0\n",
+    "):\n",
+    "    global f1\n",
+    "    global f2\n",
+    "\n",
+    "    # Update of the splines\n",
+    "    if f1.period != period:\n",
+    "        f1 = sk.PeriodicSpline1D.from_spline_coeff(\n",
+    "            rng.standard_normal(period),\n",
+    "            degree = f1.degree\n",
+    "        )\n",
+    "    f1.degree = degree\n",
+    "    f1.delay = delay\n",
+    "    if f2.period != period:\n",
+    "        f2 = sk.PeriodicSpline1D.from_spline_coeff(\n",
+    "            rng.standard_normal(period),\n",
+    "            degree = f2.degree\n",
+    "        )\n",
+    "    f2.degree = degree\n",
+    "    f2.delay = delay\n",
+    "\n",
+    "    # Addition\n",
+    "    f = sk.PeriodicSpline1D.add(f1, f2)\n",
+    "\n",
+    "    # Dynamic range\n",
+    "    image = {f1.image(), f2.image(), f.image()}\n",
+    "    plotrange = sk.interval.Interval.enclosure(image)\n",
+    "    plotrange = sk.interval.Closed((\n",
+    "        plotrange.midpoint - 0.55 * plotrange.diameter,\n",
+    "        plotrange.midpoint + 0.55 * plotrange.diameter\n",
+    "    ))\n",
+    "\n",
+    "    # Plots\n",
+    "    subplot = plt.subplots()\n",
+    "    # Augend\n",
+    "    f1.plot(subplot, plotrange = plotrange, plotpoints = 200 + 1)\n",
+    "    # Addend\n",
+    "    f2.plot(\n",
+    "        subplot,\n",
+    "        plotrange = plotrange,\n",
+    "        plotpoints = 200 + 1,\n",
+    "        curve_fmt = \"-C9\",\n",
+    "        curve_markerfmt = \"oC9\",\n",
+    "    )\n",
+    "    # Sum\n",
+    "    f.plot(\n",
+    "        subplot,\n",
+    "        plotrange = plotrange,\n",
+    "        plotpoints = 200 + 1,\n",
+    "        curve_fmt = \"-C5\",\n",
+    "        curve_markerfmt = \"oC5\",\n",
+    "        curve_lw = 2.0\n",
+    "    )\n",
+    "\n",
+    "    # Final display\n",
+    "    plt.show()\n",
+    "\n",
+    "widgets.interactive(\n",
+    "    update_plot,\n",
+    "    period = (1, max_period),\n",
+    "    degree = (0, max_degree),\n",
+    "    delay = (-max_delay, max_delay)\n",
+    ")\n"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "b5490a8b-70d3-4b94-a4ad-ea140f8cd9de",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n"
+   ]
+  }
+ ],
+ "metadata": {
+  "kernelspec": {
+   "display_name": "Python 3 (ipykernel)",
+   "language": "python",
+   "name": "python3"
+  },
+  "language_info": {
+   "codemirror_mode": {
+    "name": "ipython",
+    "version": 3
+   },
+   "file_extension": ".py",
+   "mimetype": "text/x-python",
+   "name": "python",
+   "nbconvert_exporter": "python",
+   "pygments_lexer": "ipython3",
+   "version": "3.12.6"
+  },
+  "widgets": {
+   "application/vnd.jupyter.widget-state+json": {
+    "state": {},
+    "version_major": 2,
+    "version_minor": 0
+   }
+  }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 5
+}

splinekit-0.2.10/dev/Jupyter/periodic-spline/calculus/spline_antigrad.ipynb

@@ -0,0 +1,146 @@
+{
+ "cells": [
+  {
+   "cell_type": "markdown",
+   "id": "02240a05-3526-4530-ae12-d9522465b594",
+   "metadata": {},
+   "source": [
+    "# Anti-Gradient of Periodic Splines\n",
+    "The spline $f$ whose anti-gradient is sought is shown as a <span style=\"color:#c20078\">**fuchsia**</span> curve, with data samples at the integers represented with circles and stem lines. The sample at the origin, as well as its periodized replicates, is indicated by a red circle and stem line. The small black disks depict the knots of the spline. The anti-gradient is shown in the <span style=\"color:#1f77b4\">**azure**</span> color, while the gradient of the anti-gradient is shown in <span style=\"color:#2ca02c\">**green**</span>. Up to some value offset that depends on $E\\{f\\},$ the green and the fuchsia curves are identical.\n",
+    "\n",
+    "It is instructive to realize that the knots of all three splines share their abscissa. This indicates that the delay associated to the anti-gradient is a half-integer away from the delay associated to the spline $f$ because the degree of one of these two splines is even when the degree of the other one is odd, and reciprocally."
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "230ca8a8-c806-4207-b416-9954d55f1377",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "# Load the required libraries\n",
+    "import ipywidgets as widgets\n",
+    "import matplotlib.pyplot as plt\n",
+    "import numpy as np\n",
+    "\n",
+    "import splinekit as sk # This library\n",
+    "\n",
+    "# Setup\n",
+    "max_period = 15 # Maximal period\n",
+    "max_degree = 5 # Maximal spline degree\n",
+    "max_delay = 8.0 # Maximal absolute delay\n",
+    "\n",
+    "# Initialize the generator of random numbers\n",
+    "rng = np.random.default_rng()\n",
+    "\n",
+    "# Initial random periodic cubic splines with normal Gaussian coefficients\n",
+    "f0 = sk.PeriodicSpline1D.from_spline_coeff(rng.standard_normal(6), degree = 3)\n",
+    "f0 = f0.plus(rng.choice([-1.0, 1.0]) * rng.uniform(0.05, 0.15) - f0.mean())\n",
+    "\n",
+    "# Plot\n",
+    "def update_plot (\n",
+    "    period = 6,\n",
+    "    degree = 1,\n",
+    "    delay = 0.0\n",
+    "):\n",
+    "    global f0\n",
+    "\n",
+    "    # Update of the spline\n",
+    "    if f0.period != period:\n",
+    "        f0 = sk.PeriodicSpline1D.from_spline_coeff(\n",
+    "            rng.standard_normal(period),\n",
+    "            degree = f0.degree\n",
+    "        )\n",
+    "        # Slight departure from zero-mean\n",
+    "        f0 = f0.plus(rng.choice([-1.0, 1.0]) * rng.uniform(0.05, 0.15) - f0.mean())\n",
+    "    f0.degree = degree\n",
+    "    f0.delay = delay\n",
+    "\n",
+    "    # Anti-gradient\n",
+    "    g0 = f0.anti_grad()\n",
+    "    # Gradient of the anti-gradient\n",
+    "    f = g0.gradient()\n",
+    "\n",
+    "    # Dynamic range\n",
+    "    image = {f0.image(), g0.image(), f.image()}\n",
+    "    plotrange = sk.interval.Interval.enclosure(image)\n",
+    "    plotrange = sk.interval.Closed((\n",
+    "        plotrange.midpoint - 0.55 * plotrange.diameter,\n",
+    "        plotrange.midpoint + 0.55 * plotrange.diameter\n",
+    "    ))\n",
+    "\n",
+    "    # Plot of the spline\n",
+    "    subplot = plt.subplots()\n",
+    "    f0.plot(subplot, plotrange = plotrange, plotpoints = 200 + 1, curve_fmt = \"-m\")\n",
+    "    g0.plot(subplot, plotrange = plotrange, plotpoints = 200 + 1, curve_fmt = \"-C0\")\n",
+    "    f.plot(subplot, plotrange = plotrange, plotpoints = 200 + 1, curve_fmt = \"-C2\")\n",
+    "    plt.show()\n",
+    "\n",
+    "# Interaction\n",
+    "widgets.interactive(\n",
+    "    update_plot,\n",
+    "    period = (1, max_period),\n",
+    "    degree = (0, max_degree),\n",
+    "    delay = (-max_delay, max_delay)\n",
+    ")\n"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "c719d9b2-273a-431b-9812-11db4edcb69c",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n",
+    "\n"
+   ]
+  }
+ ],
+ "metadata": {
+  "kernelspec": {
+   "display_name": "Python 3 (ipykernel)",
+   "language": "python",
+   "name": "python3"
+  },
+  "language_info": {
+   "codemirror_mode": {
+    "name": "ipython",
+    "version": 3
+   },
+   "file_extension": ".py",
+   "mimetype": "text/x-python",
+   "name": "python",
+   "nbconvert_exporter": "python",
+   "pygments_lexer": "ipython3",
+   "version": "3.12.6"
+  },
+  "widgets": {
+   "application/vnd.jupyter.widget-state+json": {
+    "state": {},
+    "version_major": 2,
+    "version_minor": 0
+   }
+  }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 5
+}