3 files changed, 83 insertions, 2 deletions
diff --git a/2.3-1/src/c/interpolation/includes/interp1.h b/2.3-1/src/c/interpolation/includes/interp1.h
index 1c01417e..738bfe67 100644
--- a/2.3-1/src/c/interpolation/includes/interp1.h
+++ b/2.3-1/src/c/interpolation/includes/interp1.h
@@ -23,6 +23,8 @@ extern "C" {
 //void dinterp13a(double *x,double *fx,double *q,int size,double *out);
 void dinterp13a(double *x,int size1,double *fx,int size2,double *q,int size3,char *a,int size4,double *out);
 
+void sinterp13a(float *x,int size1,float *fx,int size2,float *q,int size3,char *a,int size4,float *out);
+
 #ifdef  __cplusplus
 } /* extern "C" */
 #endif
diff --git a/2.3-1/src/c/interpolation/interfaces/int_interp1.h b/2.3-1/src/c/interpolation/interfaces/int_interp1.h
index 6d579e10..07d8ece4 100644
--- a/2.3-1/src/c/interpolation/interfaces/int_interp1.h
+++ b/2.3-1/src/c/interpolation/interfaces/int_interp1.h
@@ -19,9 +19,13 @@
 #include "interp1.h"
 #include <string.h>
 
-#define d2d2d2interp1d2(x,size1,fx,size2,q,size3,out)           dinterp13a(x,size1[0]*size1[1],fx,size2[0]*size1[1],q,size3[0]*size3[1],"linear",6,out)
+#define d2d2d2interp1d2(x,size1,fx,size2,q,size3,out)           	dinterp13a(x,size1[0]*size1[1],fx,size2[0]*size2[1],q,size3[0]*size3[1],"linear",6,out)
 
-#define d2d2d2g2interp1d2(x,size1,fx,size2,q,size3,ch,size4,out)	dinterp13a(x,size1[0]*size1[1],fx,size2[0]*size1[1],q,size3[0]*size3[1],ch,size4[0]*size4[1],out)
+#define d2d2d2g2interp1d2(x,size1,fx,size2,q,size3,ch,size4,out)	dinterp13a(x,size1[0]*size1[1],fx,size2[0]*size2[1],q,size3[0]*size3[1],ch,size4[0]*size4[1],out)
+
+#define s2s2s2interp1s2(x,size1,fx,size2,q,size3,out)			sinterp13a(x,size1[0]*size1[1],fx,size2[0]*size2[1],q,size3[0]*size3[1],"linear",6,out)
+
+#define s2s2s2g2interp1s2(x,size1,fx,size2,q,size3,ch,size4,out)	sinterp13a(x,size1[0]*size1[1],fx,size2[0]*size2[1],q,size3[0]*size3[1],ch,size4[0]*size4[1],out)
 
 
 #endif
diff --git a/2.3-1/src/c/interpolation/interp1/sinterp13a.c b/2.3-1/src/c/interpolation/interp1/sinterp13a.c
new file mode 100644
index 00000000..b8cb0851
--- /dev/null
+++ b/2.3-1/src/c/interpolation/interp1/sinterp13a.c
@@ -0,0 +1,75 @@
+/* Copyright (C) 2016 - IIT Bombay - FOSSEE
+
+ This file must be used under the terms of the CeCILL.
+ This source file is licensed as described in the file COPYING, which
+ you should have received as part of this distribution.  The terms
+ are also available at
+ http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+ Author: Ukasha Noor
+ Organization: FOSSEE, IIT Bombay
+ Email: toolbox@scilab.in
+*/
+
+
+#include "interp1.h"
+#include <string.h>
+
+void sinterp13a(float *x,int size1,float *fx,int size2,float *q,int size3,char *ch,int size4,float *out)
+{
+	int i,j,k,f;
+	float a,b;
+	if(strcmp(ch,"linear")==0)
+	{
+		for(i=0;i<size3;i++)
+		{
+			f=0;
+			for(j=0;j<size1;j++)
+			{
+				if(q[i]==x[j])
+				{
+					out[i]=fx[j];
+					f=1;
+					break;
+				}
+			}
+			if(f==0)
+			{
+				j=0;
+				while(q[i]>x[j])
+				{
+					j++;
+				}
+				a=x[j-1];
+				b=x[j];
+				out[i]=fx[j-1]+(q[i]-a)*((fx[j]-fx[j-1])/(b-a));
+			}
+		}
+	}
+	else if(strcmp(ch,"nearest")==0)
+	{
+		for(i=0;i<size3;i++)
+		{
+			f=0;
+			for(j=0;j<size1;j++)
+			{
+				if(q[i]==x[j])
+				{
+					out[i]=fx[j];
+					f=1;
+					break;
+				}
+			}
+			if(f==0)
+			{
+				j=0;
+				while(q[i]>x[j])
+				{
+					j++;
+				}
+				out[i]=fx[j];
+			}
+		}
+	}
+}
+			
+