4 files changed, 64 insertions, 1 deletions
diff --git a/src/c/statisticsFunctions/interfaces/int_mean.h b/src/c/statisticsFunctions/interfaces/int_mean.h
index bf196539..d1d301cd 100644
--- a/src/c/statisticsFunctions/interfaces/int_mean.h
+++ b/src/c/statisticsFunctions/interfaces/int_mean.h
@@ -48,6 +48,22 @@
 
 #define z2d0meanz2(in1,size,in2,out)		(in2==1) ? zrowmeana(in1,size[0],size[1],out) : zcolumnmeana(in1,size[0],size[1],out)	
 		
+	/*'r'  and 'c' case  */	
+#define s0g2means0(in1,in2,size2)			(in2[0]=='r') ? srowmeans(in1) : scolumnmeans(in1)	
+
+#define d0g2meand0(in1,in2,size2)			(in2[0]=='r') ? drowmeans(in1) : dcolumnmeans(in1)	
+
+#define c0g2meanc0(in1,in2,size2)			(in2[0]=='r') ? crowmeans(in1) : ccolumnmeans(in1)	
+
+#define z0g2meanz0(in1,in2,size2)			(in2[0]=='r') ? zrowmeans(in1) : zcolumnmeans(in1)	
+
+#define s2g2means2(in1,size,in2,size2,out)		(in2[0]=='r') ? srowmeana(in1,size[0],size[1],out) : scolumnmeana(in1,size[0],size[1],out)	
+
+#define d2g2meand2(in1,size,in2,size2,out)		(in2[0]=='r') ? drowmeana(in1,size[0],size[1],out) : dcolumnmeana(in1,size[0],size[1],out)
+
+#define c2g2meanc2(in1,size,in2,size2,out)		(in2[0]=='r') ? crowmeana(in1,size[0],size[1],out) : ccolumnmeana(in1,size[0],size[1],out)
+
+#define z2g2meanz2(in1,size,in2,size2,out)		(in2[0]=='r') ? zrowmeana(in1,size[0],size[1],out) : zcolumnmeana(in1,size[0],size[1],out)	
 		
 	
 	
diff --git a/src/c/statisticsFunctions/interfaces/int_prod.h b/src/c/statisticsFunctions/interfaces/int_prod.h
index cad1a740..dfec8548 100644
--- a/src/c/statisticsFunctions/interfaces/int_prod.h
+++ b/src/c/statisticsFunctions/interfaces/int_prod.h
@@ -48,6 +48,22 @@
 
 #define z2d0prodz2(in1,size,in2,out)		(in2==1) ? zrowproda(in1,size[0],size[1],out) : zcolumnproda(in1,size[0],size[1],out)	
 		
+	/*'r'  and 'c' case  */	
+#define s0g2prods0(in1,in2,size2)			(in2[0]=='r') ? srowprods(in1) : scolumnprods(in1)	
+
+#define d0g2prodd0(in1,in2,size2)			(in2[0]=='r') ? drowprods(in1) : dcolumnprods(in1)	
+
+#define c0g2prodc0(in1,in2,size2)			(in2[0]=='r') ? crowprods(in1) : ccolumnprods(in1)	
+
+#define z0g2prodz0(in1,in2,size2)			(in2[0]=='r') ? zrowprods(in1) : zcolumnprods(in1)	
+
+#define s2g2prods2(in1,size,in2,size2,out)		(in2[0]=='r') ? srowproda(in1,size[0],size[1],out) : scolumnproda(in1,size[0],size[1],out)	
+
+#define d2g2prodd2(in1,size,in2,size2,out)		(in2[0]=='r') ? drowproda(in1,size[0],size[1],out) : dcolumnproda(in1,size[0],size[1],out)
+
+#define c2g2prodc2(in1,size,in2,size2,out)		(in2[0]=='r') ? crowproda(in1,size[0],size[1],out) : ccolumnproda(in1,size[0],size[1],out)
+
+#define z2g2prodz2(in1,size,in2,size2,out)		(in2[0]=='r') ? zrowproda(in1,size[0],size[1],out) : zcolumnproda(in1,size[0],size[1],out)	
 		
 	
 	
diff --git a/src/c/statisticsFunctions/interfaces/int_sum.h b/src/c/statisticsFunctions/interfaces/int_sum.h
index 4a6240ea..08f01fee 100644
--- a/src/c/statisticsFunctions/interfaces/int_sum.h
+++ b/src/c/statisticsFunctions/interfaces/int_sum.h
@@ -49,7 +49,22 @@
 #define z2d0sumz2(in1,size,in2,out)		(in2==1) ? zrowsuma(in1,size[0],size[1],out) : zcolumnsuma(in1,size[0],size[1],out)	
 		
 		
-	
+/*'r'  and 'c' case  */	
+#define s0g2sums0(in1,in2,size2)			(in2[0]=='r') ? srowsums(in1) : scolumnsums(in1)	
+
+#define d0g2sumd0(in1,in2,size2)			(in2[0]=='r') ? drowsums(in1) : dcolumnsums(in1)	
+
+#define c0g2sumc0(in1,in2,size2)			(in2[0]=='r') ? crowsums(in1) : ccolumnsums(in1)	
+
+#define z0g2sumz0(in1,in2,size2)			(in2[0]=='r') ? zrowsums(in1) : zcolumnsums(in1)	
+
+#define s2g2sums2(in1,size,in2,size2,out)		(in2[0]=='r') ? srowsuma(in1,size[0],size[1],out) : scolumnsuma(in1,size[0],size[1],out)	
+
+#define d2g2sumd2(in1,size,in2,size2,out)		(in2[0]=='r') ? drowsuma(in1,size[0],size[1],out) : dcolumnsuma(in1,size[0],size[1],out)
+
+#define c2g2sumc2(in1,size,in2,size2,out)		(in2[0]=='r') ? crowsuma(in1,size[0],size[1],out) : ccolumnsuma(in1,size[0],size[1],out)
+
+#define z2g2sumz2(in1,size,in2,size2,out)		(in2[0]=='r') ? zrowsuma(in1,size[0],size[1],out) : zcolumnsuma(in1,size[0],size[1],out)	
 	
 
 #endif /* !__INT_SUM_H__ */
diff --git a/src/c/statisticsFunctions/interfaces/int_variance.h b/src/c/statisticsFunctions/interfaces/int_variance.h
index 5b30d5bb..dfc191ca 100644
--- a/src/c/statisticsFunctions/interfaces/int_variance.h
+++ b/src/c/statisticsFunctions/interfaces/int_variance.h
@@ -49,6 +49,22 @@
 #define z2d0variancez2(in1,size,in2,out)		(in2==1) ? zrowvariancea(in1,size[0],size[1],out) : zcolumnvariancea(in1,size[0],size[1],out)	
 		
 		
+	/*'r'  and 'c' case  */	
+#define s0g2variances0(in1,in2,size2)			(in2[0]=='r') ? srowvariances(in1) : scolumnvariances(in1)	
+
+#define d0g2varianced0(in1,in2,size2)			(in2[0]=='r') ? drowvariances(in1) : dcolumnvariances(in1)	
+
+#define c0g2variancec0(in1,in2,size2)			(in2[0]=='r') ? crowvariances(in1) : ccolumnvariances(in1)	
+
+#define z0g2variancez0(in1,in2,size2)			(in2[0]=='r') ? zrowvariances(in1) : zcolumnvariances(in1)	
+
+#define s2g2variances2(in1,size,in2,size2,out)		(in2[0]=='r') ? srowvariancea(in1,size[0],size[1],out) : scolumnvariancea(in1,size[0],size[1],out)	
+
+#define d2g2varianced2(in1,size,in2,size2,out)		(in2[0]=='r') ? drowvariancea(in1,size[0],size[1],out) : dcolumnvariancea(in1,size[0],size[1],out)
+
+#define c2g2variancec2(in1,size,in2,size2,out)		(in2[0]=='r') ? crowvariancea(in1,size[0],size[1],out) : ccolumnvariancea(in1,size[0],size[1],out)
+
+#define z2g2variancez2(in1,size,in2,size2,out)		(in2[0]=='r') ? zrowvariancea(in1,size[0],size[1],out) : zcolumnvariancea(in1,size[0],size[1],out)