From f5b9993f4356224e478f91ddd57ea013292b4d84 Mon Sep 17 00:00:00 2001
From: Evgeny Prilepin <esp.home@gmail.com>
Date: Thu, 17 Apr 2025 21:35:52 +0100
Subject: [PATCH 1/2] docs: add extrapolation section to tutorial documentation

---
 csaps/_sspumv.py  |  8 +++++++-
 docs/tutorial.rst | 40 ++++++++++++++++++++++++++++++++++++++--
 2 files changed, 45 insertions(+), 3 deletions(-)

diff --git a/csaps/_sspumv.py b/csaps/_sspumv.py
index b343de2..83abc12 100644
--- a/csaps/_sspumv.py
+++ b/csaps/_sspumv.py
@@ -78,7 +78,13 @@ def __repr__(self) -> str:  # pragma: no cover
 
 
 class CubicSmoothingSpline(
-    ISmoothingSpline[SplinePPForm, float, UnivariateDataType, int, Union[bool, Literal['periodic']]],
+    ISmoothingSpline[
+        SplinePPForm,
+        float,
+        UnivariateDataType,
+        int,
+        Union[bool, Literal['periodic']],
+    ]
 ):
     """Cubic smoothing spline
 
diff --git a/docs/tutorial.rst b/docs/tutorial.rst
index 0bf9d83..386f4b4 100644
--- a/docs/tutorial.rst
+++ b/docs/tutorial.rst
@@ -368,7 +368,7 @@ Let's show it on a simple example.
     yi1_n = csaps(x1, y, xi1, smooth=0.8, normalizedsmooth=True)
     yi2_n = csaps(x2, y, xi2, smooth=0.8, normalizedsmooth=True)
 
-    f, ((ax1, ax2), (ax3, ax4)) = plt.subplots(2, 2, figsize=(8, 6))
+    f, ((ax1, ax2), (ax3, ax4)) = plt.subplots(2, 2, figsize=(10, 6))
     ax1.plot(x1, y, 'o', xi1, yi1, '-')
     ax2.plot(x2, y, 'o', xi2, yi2, '-')
     ax3.plot(x1, y, 'o', xi1, yi1_n, '-')
@@ -406,7 +406,7 @@ The example for univariate data:
     yi1 = spline(xi1)
     yi2 = spline(xi2)
 
-    f, (ax1, ax2) = plt.subplots(2, 1, figsize=(5, 6))
+    f, (ax1, ax2) = plt.subplots(1, 2, figsize=(10, 4))
     ax1.plot(x, y, 'o', xi1, yi1, '.-')
     ax2.plot(x, y, 'o', xi2, yi2, '.-')
 
@@ -414,6 +414,42 @@ The example for univariate data:
     ax2.set_title('50 evaluated points')
 
 
+.. _tutorial-extrapolation:
+
+Extrapolation
+~~~~~~~~~~~~~
+
+Spline values can be evaluated out-of-bounds input X-values (the input grid) based on first and last intervals.
+These values will be extrapolated. The ``extrapolate`` parameter in the spline evaluation method is used for this.
+Extrapolation can be used for all evaluated splines: univariate, multivariate and N-D grid.
+
+Here is an example for univariate data:
+
+.. plot::
+
+    x, y = univariate_data()
+    xi = np.linspace(x[0] - 2.0, x[-1] + 2.0, 50)
+
+    s1 = csaps(x, y, smooth=0.85)
+    s2 = csaps(x, y, smooth=0.85)
+
+    yi1 = s1(xi, extrapolate=True)
+    yi2 = s2(xi, extrapolate='periodic')
+
+    _, (ax1, ax2) = plt.subplots(1, 2, figsize=(10, 4))
+    ax1.plot(x, y, 'o:', xi, yi1, '.-')
+    ax2.plot(x, y, 'o:', xi, yi2, '.-')
+
+    ax1.set_title('extrapolate=True')
+    ax2.set_title('extrapolate="periodic"')
+
+.. attention::
+
+    How we can see ``'periodic'`` extrapolation method works incorrectly with the spline. There is a discontinuity.
+    This can be fixed in the future. Please see `the issue <https://github.com/espdev/csaps/issues/46>`_ for details.
+    So don't use this method for now if you need extrapolation.
+
+
 .. _tutorial-analysis:
 
 Analysis

From 1bef7adb3b9078d4e7269ced63a8bfeb4cd4bdff Mon Sep 17 00:00:00 2001
From: Evgeny Prilepin <esp.home@gmail.com>
Date: Thu, 17 Apr 2025 21:37:32 +0100
Subject: [PATCH 2/2] build: bump version to dev0

---
 pyproject.toml | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/pyproject.toml b/pyproject.toml
index 6d902db..e8a474f 100644
--- a/pyproject.toml
+++ b/pyproject.toml
@@ -1,6 +1,6 @@
 [tool.poetry]
 name = "csaps"
-version = "1.3.2"
+version = "1.3.3.dev0"
 description = "Cubic spline approximation (smoothing)"
 authors = ["Evgeny Prilepin <esp.home@gmail.com>"]
 license = "MIT"