sigma arrays do need to be length N

for the interface, even though the last value is not used, to my understanding
noaa-oar-arl · Jan 24, 2024 · 2459506 · 2459506
1 parent 75d7574
commit 2459506
Showing 1 changed file with 4 additions and 2 deletions.
diff --git a/pytspack/__init__.py b/pytspack/__init__.py
@@ -43,7 +43,8 @@ def hval(xp, x, y, yp, sigma):
  where HP is the derivative of H.
  sigma : array
  Tension factors, for each interval in the original X points
- (as such, length N - 1).
+ (element I corresponds to the interval (I,I+1);
+ the last value in the array is not used).
 
  Returns
  -------
@@ -75,7 +76,8 @@ def hpval(xp, x, y, yp, sigma):
  First derivatives at original X points. HP(X(I)) = YP(I).
  sigma : array_like
  Tension factors, for each interval in the original X points
- (as such, length N - 1).
+ (element I corresponds to the interval (I,I+1);
+ the last value in the array is not used).
 
  Returns
  -------