splinekit 0.2.6__tar.gz → 0.2.8__tar.gz

{splinekit-0.2.6 → splinekit-0.2.8}/PKG-INFO

@@ -1,9 +1,9 @@
 Metadata-Version: 2.4
 Name: splinekit
-Version: 0.2.6
+Version: 0.2.8
 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,6 +51,8 @@ Description-Content-Type: text/markdown
 
 <!-- splinekit/README.md -->
 
+# Version 0.2.8
+
 # splinekit: Spline Operations
 `splinekit` is a Python-based open-source software library aimed at the manipulation of one-dimensional periodic splines.
 
@@ -114,7 +116,7 @@ pip install -e .[dev]
 
 ## Release Notes
 
-### 0.1.0
+### 0.1.*
 First stable version. Pure Python.
 
 ### 0.2.1
@@ -129,3 +131,8 @@ Minor debug of least-significant-bit edge cases.
 ### 0.2.6
 Some demos added. Documentation minor updates.
 
+### 0.2.7
+Upgraded to jupyterlite 0.7
+
+### 0.2.8
+Some demos added.

{splinekit-0.2.6 → splinekit-0.2.8}/README.md

@@ -1,5 +1,7 @@
 <!-- splinekit/README.md -->
 
+# Version 0.2.8
+
 # splinekit: Spline Operations
 `splinekit` is a Python-based open-source software library aimed at the manipulation of one-dimensional periodic splines.
 
@@ -63,7 +65,7 @@ pip install -e .[dev]
 
 ## Release Notes
 
-### 0.1.0
+### 0.1.*
 First stable version. Pure Python.
 
 ### 0.2.1
@@ -78,3 +80,8 @@ Minor debug of least-significant-bit edge cases.
 ### 0.2.6
 Some demos added. Documentation minor updates.
 
+### 0.2.7
+Upgraded to jupyterlite 0.7
+
+### 0.2.8
+Some demos added.

{splinekit-0.2.6 → splinekit-0.2.8/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.6"
+version = "0.2.8"
 dependencies = [
     "numpy >=1.26",
     "scipy >=1.15",

splinekit-0.2.8/dev/Jupyter/bspline/members/bspline_member.ipynb

@@ -0,0 +1,109 @@
+{
+ "cells": [
+  {
+   "cell_type": "markdown",
+   "id": "0c723be4-e5f6-4b5a-a6a2-11ece85e5f1b",
+   "metadata": {},
+   "source": [
+    "# Polynomial B-Spline\n",
+    "Here is the plot of a B-spline of a given degree, over precisely its support. The hollow (blue or red) circles indicate the samples at the integers. The function is made of polynomial pieces that lie between so-called knots. The small filled (black) discs indicate the location of the knots of the spline. They lie at the half integers for even degrees and at the integers for odd degrees."
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "65fb414e-a36a-4bff-ae5a-a03e6498deb8",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "# Load the required libraries\n",
+    "from ipywidgets import interactive\n",
+    "import matplotlib.pyplot as plt\n",
+    "\n",
+    "import splinekit as sk # This library\n",
+    "\n",
+    "# Define the plot function\n",
+    "def b_spline_plot (\n",
+    "    degree = 3\n",
+    "):\n",
+    "    # Support of a B-spline\n",
+    "    supp = sk.b_spline_support(degree)\n",
+    "    # Avoid overlap and leave a margin\n",
+    "    # The margin serves a purpose only for the B-spline of degree 0\n",
+    "    period = supp.diameter + 1\n",
+    "    # Construct a spline from a non-overlapping periodic B-spline\n",
+    "    s = sk.PeriodicSpline1D.periodized_b_spline(\n",
+    "        period = int(period),\n",
+    "        degree = degree\n",
+    "    )\n",
+    "    # Plot the B-spline\n",
+    "    s.plot(\n",
+    "        plt.subplots(),\n",
+    "        plotdomain = sk.interval.Closed((-0.5 * supp.diameter, 0.5 * supp.diameter)),\n",
+    "        plotrange = sk.interval.Closed((-0.05, 1.05)),\n",
+    "        plotpoints = 200 + 1\n",
+    "    )\n",
+    "\n",
+    "# Interact with the degree\n",
+    "interactive(b_spline_plot, degree = (0, 30))\n"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "b80ced66-8818-4b26-9ea0-effa40d17b3e",
+   "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.6/dev/Jupyter → splinekit-0.2.8/dev/Jupyter/bspline/stability}/bspline_numeric_stability.ipynb

@@ -172,7 +172,8 @@
     "#                 #\n",
     "###################\n",
     "\n",
-    "# Duration of the computations on a desktop computer of year 2021: ~20s\n",
+    "# Duration of the native computations on a desktop computer of year 2021: ~20s\n",
+    "# Duration of the WebAssembly computations: about twice\n",
     "\n",
     "\n"
    ]
@@ -321,7 +322,8 @@
     "#                 #\n",
     "###################\n",
     "\n",
-    "# Duration of the computations on a desktop computer of year 2021: ~20s\n",
+    "# Duration of the native computations on a desktop computer of year 2021: ~20s\n",
+    "# Duration of the WebAssembly computations: about twice\n",
     "\n",
     "\n"
    ]
@@ -447,7 +449,8 @@
     "#                 #\n",
     "###################\n",
     "\n",
-    "# Duration of the computations on a desktop computer of year 2021: ~5s\n",
+    "# Duration of the native computations on a desktop computer of year 2021: ~5s\n",
+    "# Duration of the WebAssembly computations: about twice\n",
     "\n",
     "\n"
    ]

splinekit-0.2.8/dev/Jupyter/padding/change-basis/padding_bases.ipynb

@@ -0,0 +1,222 @@
+{
+ "cells": [
+  {
+   "cell_type": "markdown",
+   "id": "02240a05-3526-4530-ae12-d9522465b594",
+   "metadata": {},
+   "source": [
+    "# Change of Basis for Periodic Splines\n",
+    "The data samples at the integers are 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 thin colored curves are the weighted and shifted basis components whose sum is the curve depicted in thick blue."
+   ]
+  },
+  {
+   "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_degree = 9 # Maximal spline degree\n",
+    "max_samples = 15 # Maximal support\n",
+    "\n",
+    "# Initialize the generator of random numbers\n",
+    "rng = np.random.default_rng()\n",
+    "\n",
+    "# Initial random periodic cubic spline\n",
+    "s0 = sk.PeriodicSpline1D.from_spline_coeff(rng.standard_normal(6), degree = 3)\n",
+    "\n",
+    "# Plot\n",
+    "def update_plot (\n",
+    "    degree = 3,\n",
+    "    samples = 6,\n",
+    "    basis = 0\n",
+    "):\n",
+    "    global s0\n",
+    "\n",
+    "    # Update of the spline\n",
+    "    f = s0.get_samples(0, support_length = s0.period)\n",
+    "    if s0.period < samples:\n",
+    "        f = np.append(f, rng.standard_normal(samples - len(f)))\n",
+    "    else:\n",
+    "        f = f[ : samples]\n",
+    "    s0 = sk.PeriodicSpline1D.from_samples(f, degree = degree)\n",
+    "\n",
+    "    # Decomposition according to the basis\n",
+    "    weights = np.copy(s0.spline_coeff)\n",
+    "    if 0 == basis: # BASIC\n",
+    "        sk.change_basis_p(\n",
+    "            weights,\n",
+    "            degree = degree,\n",
+    "            source_basis = sk.Bases.BASIC,\n",
+    "            target_basis = sk.Bases.BASIC\n",
+    "        )\n",
+    "        bases = [\n",
+    "            sk.PeriodicSpline1D.periodized_b_spline(\n",
+    "                period = s0.period,\n",
+    "                degree = degree\n",
+    "            ).delayed_by(k).times(wk)\n",
+    "            for (k, wk) in enumerate(weights)\n",
+    "        ]\n",
+    "    elif 1 == basis: # CARDINAL\n",
+    "        sk.change_basis_p(\n",
+    "            weights,\n",
+    "            degree = degree,\n",
+    "            source_basis = sk.Bases.BASIC,\n",
+    "            target_basis = sk.Bases.CARDINAL\n",
+    "        )\n",
+    "        bases = [\n",
+    "            sk.PeriodicSpline1D.periodized_cardinal_b_spline(\n",
+    "                period = s0.period,\n",
+    "                degree = degree\n",
+    "            ).delayed_by(k).times(wk)\n",
+    "            for (k, wk) in enumerate(weights)\n",
+    "        ]\n",
+    "    elif 2 == basis: # DUAL\n",
+    "        sk.change_basis_p(\n",
+    "            weights,\n",
+    "            degree = degree,\n",
+    "            source_basis = sk.Bases.BASIC,\n",
+    "            target_basis = sk.Bases.DUAL\n",
+    "        )\n",
+    "        bases = [\n",
+    "            sk.PeriodicSpline1D.periodized_dual_b_spline(\n",
+    "                period = s0.period,\n",
+    "                dual_degree = degree,\n",
+    "                primal_degree = degree\n",
+    "            ).delayed_by(k).times(wk)\n",
+    "            for (k, wk) in enumerate(weights)\n",
+    "        ]\n",
+    "    elif 3 == basis: # ORTHONORMAL\n",
+    "        sk.change_basis_p(\n",
+    "            weights,\n",
+    "            degree = degree,\n",
+    "            source_basis = sk.Bases.BASIC,\n",
+    "            target_basis = sk.Bases.ORTHONORMAL\n",
+    "        )\n",
+    "        bases = [\n",
+    "            sk.PeriodicSpline1D.periodized_orthonormal_b_spline(\n",
+    "                period = s0.period,\n",
+    "                degree = degree\n",
+    "            ).delayed_by(k).times(wk)\n",
+    "            for (k, wk) in enumerate(weights)\n",
+    "        ]\n",
+    "\n",
+    "    # Dynamic range\n",
+    "    image = {b.image() for b in bases}\n",
+    "    image.add(s0.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",
+    "    s0.plot(\n",
+    "        subplot,\n",
+    "        plotrange = plotrange,\n",
+    "        plotpoints = 200 + 1,\n",
+    "        curve_fmt = \"-C0\",\n",
+    "        knot_marker = \"\"\n",
+    "    )\n",
+    "\n",
+    "    # Independent plot of each weighted basis component\n",
+    "    for (k, b) in enumerate(bases):\n",
+    "        b.plot(\n",
+    "            subplot,\n",
+    "            plotpoints = 200 + 1,\n",
+    "            curve_fmt = \"-C\" + str(1 + k % 9),\n",
+    "            curve_lw = 0.25,\n",
+    "            curve_markerfmt = \"\",\n",
+    "            curvestem_linefmt = \"None\",\n",
+    "            knot_marker = \"\",\n",
+    "            periodbound_markerfmt = \"\",\n",
+    "            periodboundstem_linefmt = \"None\"\n",
+    "        )\n",
+    "    plt.show()\n",
+    "\n",
+    "widgets.interactive(\n",
+    "    update_plot,\n",
+    "    degree = (0, max_degree),\n",
+    "    samples = (1, max_samples),\n",
+    "    basis = widgets.RadioButtons(\n",
+    "        options = [\n",
+    "            (\"BASIC\", 0),\n",
+    "            (\"CARDINAL\", 1),\n",
+    "            (\"DUAL\", 2),\n",
+    "            (\"ORTHONORMAL\", 3)\n",
+    "        ],\n",
+    "        value = 0,\n",
+    "        description = \"Basis:\",\n",
+    "        disabled = False\n",
+    "    )\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
+}

{splinekit-0.2.6/dev/Jupyter → splinekit-0.2.8/dev/Jupyter/padding/padding}/padding_data.ipynb

@@ -26,9 +26,13 @@
     "# Setup\n",
     "max_degree = 5 # Maximal spline degree\n",
     "max_samples = 12 # Maximal support of the observed data samples\n",
+    "\n",
+    "# Initialize the generator of random numbers\n",
+    "rng = np.random.default_rng()\n",
+    "\n",
     "# Random periodic cubic spline\n",
     "k0 = 6 # Initial number of observed samples\n",
-    "s0 = sk.PeriodicSpline1D.from_spline_coeff(np.random.standard_normal(k0), degree = 3)\n",
+    "s0 = sk.PeriodicSpline1D.from_spline_coeff(rng.standard_normal(k0), degree = 3)\n",
     "\n",
     "# Plot\n",
     "def update_plot (\n",
@@ -47,7 +51,7 @@
     "    # Update of the number of samples\n",
     "    f0 = s0.get_samples(0, support_length = k0)\n",
     "    if k0 < samples:\n",
-    "        f0 = np.append(f0, np.random.standard_normal(samples - len(f0)))\n",
+    "        f0 = np.append(f0, rng.standard_normal(samples - len(f0)))\n",
     "    else:\n",
     "        f0 = f0[ : samples]\n",
     "    k0 = len(f0) # Number of observed samples\n",

{splinekit-0.2.6/dev/Jupyter → splinekit-0.2.8/dev/Jupyter/padding/padding}/padding_speed.ipynb

@@ -159,8 +159,9 @@
     "#                 #\n",
     "###################\n",
     "\n",
+    "# Duration of the native computations on a desktop computer of year 2021: ~25s\n",
+    "# Duration of the WebAssembly computations: about twice\n",
     "\n",
-    "# Duration of the computations on a desktop computer of year 2021: ~25s\n",
     "\n"
    ]
   },

{splinekit-0.2.6/dev/Jupyter → splinekit-0.2.8/dev/Jupyter/periodic-spline/bounds}/spline_bounds.ipynb

@@ -28,6 +28,9 @@
     "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",
     "# Persistent spline\n",
     "f = sk.PeriodicSpline1D()\n",
     "\n",
@@ -46,7 +49,7 @@
     "    # Update of the period\n",
     "    c = f.spline_coeff\n",
     "    if f.period < period:\n",
-    "        c = np.append(c, np.random.standard_normal(period - f.period))\n",
+    "        c = np.append(c, rng.standard_normal(period - f.period))\n",
     "    else:\n",
     "        c = c[ : period]\n",
     "    # Spline being bounded\n",

splinekit-0.2.8/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
+}