tsInterp1Extrap.m: offsetting by default negative values, for which i…

…t is impossible to compute logarithm
menta78 · Aug 30, 2018 · 612f06a · 612f06a
1 parent eff6bf5
commit 612f06a
Showing 1 changed file with 11 additions and 1 deletion.
diff --git a/tsInterp1Extrap.m b/tsInterp1Extrap.m
@@ -36,6 +36,16 @@
 
 end
 
+minX = min(min(Xq(:)), min(X(:)));
+if minX < 1
+  xoffset = -minX + 1;
+else
+  xoffset = 0;
+end
+
+X = X + xoffset;
+Xq = Xq + xoffset;
+
 % sort x independent variable
 [X,ii]=sort(X);
 
@@ -78,4 +88,4 @@
 
     end
 
-end
+end