Add a test of spline

AEI-ACR · Apr 4, 2023 · 2e2f98b · 2e2f98b
1 parent fd2f77f
commit 2e2f98b
Showing 1 changed file with 13 additions and 0 deletions.
diff --git a/test/test_spline.py b/test/test_spline.py
@@ -0,0 +1,13 @@
+import numpy as np
+from pygsl_lite import spline
+
+def test_cubic_spline():
+ t_old = np.linspace(-2*np.pi,2*np.pi,50)
+ y_old = np.cos(t_old)
+ n = len(t_old)
+ intrp = spline.cspline(n)
+ intrp.init(t_old, y_old)
+ t_new = np.linspace(-2*np.pi,2*np.pi,200)
+ y_interp = intrp.eval_e_vector(t_new)
+ assert np.abs(np.sum(y_interp)-0.9878615069653094)<1e-8,"Interpolation test failed"
+