From 06337f0dc8114c70fd0c7767083971a0d091750a Mon Sep 17 00:00:00 2001
From: Sandeep Gupta
Date: Wed, 5 Jul 2017 12:41:25 +0530
Subject: LinearAlgebra and MatrixOperation Update

---
 src/c/linearAlgebra/sqroot/dsqroota.c | 130 ++++++++++++++++++++++++++++++++++
 1 file changed, 130 insertions(+)
 create mode 100644 src/c/linearAlgebra/sqroot/dsqroota.c

(limited to 'src/c/linearAlgebra/sqroot/dsqroota.c')

diff --git a/src/c/linearAlgebra/sqroot/dsqroota.c b/src/c/linearAlgebra/sqroot/dsqroota.c
new file mode 100644
index 0000000..038afc9
--- /dev/null
+++ b/src/c/linearAlgebra/sqroot/dsqroota.c
@@ -0,0 +1,130 @@
+/* Copyright (C) 2017 - 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: Sandeep Gupta
+ Organization: FOSSEE, IIT Bombay
+ Email: toolbox@scilab.in
+ */
+
+/* Function - sqroot of scilab, W*W' hermitian factorization */
+
+#include <stdio.h> 
+#include "stdlib.h"
+#include "string.h"
+#include <math.h>
+#include "matrixTranspose.h"
+#include "svd.h"
+#include "matrixMultiplication.h"
+
+#define eps 2.22044604925e-16
+
+/*It would be good, if you are doing conversoins for only valid inputs before
+ running the program, check all the conditions before hand.
+*/
+
+void dsqroota(double *inp,int row,int col, double *out){	
+	if(row != col){
+		printf("Enter valid inputs only - matrix should be symetric\n");
+		return;
+	}
+	int i,j;
+	double *U,*S,*V;
+	double *A,*B;
+	int rk;
+	
+	U = (double *)malloc((double)row*row*sizeof(double));
+	S = (double *)malloc((double)Min(row,col)*Min(row,col)*sizeof(double));		
+	V = (double *)malloc((double)col*col*sizeof(double));
+	A = (double *)malloc(rk*rk*sizeof(double));
+	B = (double *)malloc(rk*row*sizeof(double));
+	
+	double *Q1;	
+	Q1 = (double *)malloc(row*col*sizeof(double));	
+	
+	dtransposea(inp,row,col,Q1);
+	
+	double *Q2;
+	Q2 = (double *)malloc(col*row*sizeof(double));
+
+	/* Q2 = (inp+inp1')/2; */
+	for(i=0;i<row;i++){
+		for(j=0;j<row;j++){
+			Q2[i+j*row] = ((inp[i+j*row]+Q1[i+j*row])/2);	
+		}	
+	}
+	
+	/*  norm(Q2-Q1,1) - finding the max value from sum of columns */ 
+	double sum = 0;
+	double maxi=0;
+	for(i=0;i<row;i++){
+		sum = 0;			
+		for(j=0;j<col;j++) sum+=(Q2[i*row+j]-inp[i*row+j]);	
+		if(maxi < sum){
+			maxi = sum;	
+		}
+	}		
+	
+	/* if norm(Q1-Q,1) > 100*%eps then */
+	if(maxi > 100*eps){
+		printf("Warning: Wrong size for input argument and Symmetric expected\n");	
+	}
+	maxi = 0;
+	for(i=0;i<col;i++){
+		sum=0;
+		for(j=0;j<row;j++) sum+=Q1[i*row+j];
+		if(maxi < sum){
+			maxi = sum;		
+		}
+	}
+
+	/*if norm(Q,1) < sqrt(%eps) then S=[];return;end*/
+
+	if(maxi < sqrt(eps)){
+		out = NULL;	
+	}
+	else{		
+		rk = dsvda(0,inp,row,col,0,4,U,S,V);
+		
+		/*Will be used in complex numbers*/		
+		//C = (double *)malloc(rk*row*sizeof(double));
+		
+		/*for(i=0;i<row;i++){
+			for(j=0;j<row;j++){
+				printf("%lf ",S[i*row+j]);			
+			}	
+			printf("\n");	
+		}*/
+		/*sqrt of S*/
+		printf("%d ",rk);
+		for(i=0;i<rk;i++){
+			for(j=0;j<rk;j++){
+				A[i*rk+j] = sqrt(S[i*row+j]);
+				//printf("%lf ",A[i*rk+j]);
+			}
+		}
+		for(i=0;i<col*rk;i++){
+			B[i] = V[i];
+			printf("%lf ",B[i]);	
+		}
+		printf("\n");
+		/*for(i=0;i<rk;i++){
+			for(j=0;j<rk;j++){
+				//A[i*rk+j] = sqrt(S[i*row+j]);
+				printf("%lf ",A[i*rk+j]);
+			}
+			printf("\n");
+		}*/
+		/*for(i=0;i<col;i++){
+			for(j=0;j<rk;j++){
+				//B[i*col+j] = VT[i*col+j];	
+				printf("%lf ",B[i*col+j]);	
+			}	
+			printf("\n");
+		}*/
+		dmulma(B,col,rk,A,rk,rk,out);
+	}
+}
-- 
cgit 


From cb1d99232e521c34e9f0c271a6c4176cc7b9cbe4 Mon Sep 17 00:00:00 2001
From: Sandeep Gupta
Date: Thu, 6 Jul 2017 15:16:16 +0530
Subject: Changes

---
 src/c/linearAlgebra/sqroot/dsqroota.c | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

(limited to 'src/c/linearAlgebra/sqroot/dsqroota.c')

diff --git a/src/c/linearAlgebra/sqroot/dsqroota.c b/src/c/linearAlgebra/sqroot/dsqroota.c
index 038afc9..a9062e6 100644
--- a/src/c/linearAlgebra/sqroot/dsqroota.c
+++ b/src/c/linearAlgebra/sqroot/dsqroota.c
@@ -108,9 +108,9 @@ void dsqroota(double *inp,int row,int col, double *out){
 		}
 		for(i=0;i<col*rk;i++){
 			B[i] = V[i];
-			printf("%lf ",B[i]);	
+			//printf("%lf ",B[i]);	
 		}
-		printf("\n");
+		//printf("\n");
 		/*for(i=0;i<rk;i++){
 			for(j=0;j<rk;j++){
 				//A[i*rk+j] = sqrt(S[i*row+j]);
-- 
cgit 


From ea958d3c401761dcc24865d9639b2fab31038db8 Mon Sep 17 00:00:00 2001
From: Brijeshcr
Date: Thu, 6 Jul 2017 15:48:47 +0530
Subject: Revert "LinearAlgebra Function Added"

---
 src/c/linearAlgebra/sqroot/dsqroota.c | 130 ----------------------------------
 1 file changed, 130 deletions(-)
 delete mode 100644 src/c/linearAlgebra/sqroot/dsqroota.c

(limited to 'src/c/linearAlgebra/sqroot/dsqroota.c')

diff --git a/src/c/linearAlgebra/sqroot/dsqroota.c b/src/c/linearAlgebra/sqroot/dsqroota.c
deleted file mode 100644
index a9062e6..0000000
--- a/src/c/linearAlgebra/sqroot/dsqroota.c
+++ /dev/null
@@ -1,130 +0,0 @@
-/* Copyright (C) 2017 - 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: Sandeep Gupta
- Organization: FOSSEE, IIT Bombay
- Email: toolbox@scilab.in
- */
-
-/* Function - sqroot of scilab, W*W' hermitian factorization */
-
-#include <stdio.h> 
-#include "stdlib.h"
-#include "string.h"
-#include <math.h>
-#include "matrixTranspose.h"
-#include "svd.h"
-#include "matrixMultiplication.h"
-
-#define eps 2.22044604925e-16
-
-/*It would be good, if you are doing conversoins for only valid inputs before
- running the program, check all the conditions before hand.
-*/
-
-void dsqroota(double *inp,int row,int col, double *out){	
-	if(row != col){
-		printf("Enter valid inputs only - matrix should be symetric\n");
-		return;
-	}
-	int i,j;
-	double *U,*S,*V;
-	double *A,*B;
-	int rk;
-	
-	U = (double *)malloc((double)row*row*sizeof(double));
-	S = (double *)malloc((double)Min(row,col)*Min(row,col)*sizeof(double));		
-	V = (double *)malloc((double)col*col*sizeof(double));
-	A = (double *)malloc(rk*rk*sizeof(double));
-	B = (double *)malloc(rk*row*sizeof(double));
-	
-	double *Q1;	
-	Q1 = (double *)malloc(row*col*sizeof(double));	
-	
-	dtransposea(inp,row,col,Q1);
-	
-	double *Q2;
-	Q2 = (double *)malloc(col*row*sizeof(double));
-
-	/* Q2 = (inp+inp1')/2; */
-	for(i=0;i<row;i++){
-		for(j=0;j<row;j++){
-			Q2[i+j*row] = ((inp[i+j*row]+Q1[i+j*row])/2);	
-		}	
-	}
-	
-	/*  norm(Q2-Q1,1) - finding the max value from sum of columns */ 
-	double sum = 0;
-	double maxi=0;
-	for(i=0;i<row;i++){
-		sum = 0;			
-		for(j=0;j<col;j++) sum+=(Q2[i*row+j]-inp[i*row+j]);	
-		if(maxi < sum){
-			maxi = sum;	
-		}
-	}		
-	
-	/* if norm(Q1-Q,1) > 100*%eps then */
-	if(maxi > 100*eps){
-		printf("Warning: Wrong size for input argument and Symmetric expected\n");	
-	}
-	maxi = 0;
-	for(i=0;i<col;i++){
-		sum=0;
-		for(j=0;j<row;j++) sum+=Q1[i*row+j];
-		if(maxi < sum){
-			maxi = sum;		
-		}
-	}
-
-	/*if norm(Q,1) < sqrt(%eps) then S=[];return;end*/
-
-	if(maxi < sqrt(eps)){
-		out = NULL;	
-	}
-	else{		
-		rk = dsvda(0,inp,row,col,0,4,U,S,V);
-		
-		/*Will be used in complex numbers*/		
-		//C = (double *)malloc(rk*row*sizeof(double));
-		
-		/*for(i=0;i<row;i++){
-			for(j=0;j<row;j++){
-				printf("%lf ",S[i*row+j]);			
-			}	
-			printf("\n");	
-		}*/
-		/*sqrt of S*/
-		printf("%d ",rk);
-		for(i=0;i<rk;i++){
-			for(j=0;j<rk;j++){
-				A[i*rk+j] = sqrt(S[i*row+j]);
-				//printf("%lf ",A[i*rk+j]);
-			}
-		}
-		for(i=0;i<col*rk;i++){
-			B[i] = V[i];
-			//printf("%lf ",B[i]);	
-		}
-		//printf("\n");
-		/*for(i=0;i<rk;i++){
-			for(j=0;j<rk;j++){
-				//A[i*rk+j] = sqrt(S[i*row+j]);
-				printf("%lf ",A[i*rk+j]);
-			}
-			printf("\n");
-		}*/
-		/*for(i=0;i<col;i++){
-			for(j=0;j<rk;j++){
-				//B[i*col+j] = VT[i*col+j];	
-				printf("%lf ",B[i*col+j]);	
-			}	
-			printf("\n");
-		}*/
-		dmulma(B,col,rk,A,rk,rk,out);
-	}
-}
-- 
cgit 


From 7a7f685a8436b456b246c49baf76bb8af930b214 Mon Sep 17 00:00:00 2001
From: Sandeep Gupta
Date: Fri, 7 Jul 2017 00:29:35 +0530
Subject: NORM

---
 src/c/linearAlgebra/sqroot/dsqroota.c | 130 ++++++++++++++++++++++++++++++++++
 1 file changed, 130 insertions(+)
 create mode 100644 src/c/linearAlgebra/sqroot/dsqroota.c

(limited to 'src/c/linearAlgebra/sqroot/dsqroota.c')

diff --git a/src/c/linearAlgebra/sqroot/dsqroota.c b/src/c/linearAlgebra/sqroot/dsqroota.c
new file mode 100644
index 0000000..a9062e6
--- /dev/null
+++ b/src/c/linearAlgebra/sqroot/dsqroota.c
@@ -0,0 +1,130 @@
+/* Copyright (C) 2017 - 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: Sandeep Gupta
+ Organization: FOSSEE, IIT Bombay
+ Email: toolbox@scilab.in
+ */
+
+/* Function - sqroot of scilab, W*W' hermitian factorization */
+
+#include <stdio.h> 
+#include "stdlib.h"
+#include "string.h"
+#include <math.h>
+#include "matrixTranspose.h"
+#include "svd.h"
+#include "matrixMultiplication.h"
+
+#define eps 2.22044604925e-16
+
+/*It would be good, if you are doing conversoins for only valid inputs before
+ running the program, check all the conditions before hand.
+*/
+
+void dsqroota(double *inp,int row,int col, double *out){	
+	if(row != col){
+		printf("Enter valid inputs only - matrix should be symetric\n");
+		return;
+	}
+	int i,j;
+	double *U,*S,*V;
+	double *A,*B;
+	int rk;
+	
+	U = (double *)malloc((double)row*row*sizeof(double));
+	S = (double *)malloc((double)Min(row,col)*Min(row,col)*sizeof(double));		
+	V = (double *)malloc((double)col*col*sizeof(double));
+	A = (double *)malloc(rk*rk*sizeof(double));
+	B = (double *)malloc(rk*row*sizeof(double));
+	
+	double *Q1;	
+	Q1 = (double *)malloc(row*col*sizeof(double));	
+	
+	dtransposea(inp,row,col,Q1);
+	
+	double *Q2;
+	Q2 = (double *)malloc(col*row*sizeof(double));
+
+	/* Q2 = (inp+inp1')/2; */
+	for(i=0;i<row;i++){
+		for(j=0;j<row;j++){
+			Q2[i+j*row] = ((inp[i+j*row]+Q1[i+j*row])/2);	
+		}	
+	}
+	
+	/*  norm(Q2-Q1,1) - finding the max value from sum of columns */ 
+	double sum = 0;
+	double maxi=0;
+	for(i=0;i<row;i++){
+		sum = 0;			
+		for(j=0;j<col;j++) sum+=(Q2[i*row+j]-inp[i*row+j]);	
+		if(maxi < sum){
+			maxi = sum;	
+		}
+	}		
+	
+	/* if norm(Q1-Q,1) > 100*%eps then */
+	if(maxi > 100*eps){
+		printf("Warning: Wrong size for input argument and Symmetric expected\n");	
+	}
+	maxi = 0;
+	for(i=0;i<col;i++){
+		sum=0;
+		for(j=0;j<row;j++) sum+=Q1[i*row+j];
+		if(maxi < sum){
+			maxi = sum;		
+		}
+	}
+
+	/*if norm(Q,1) < sqrt(%eps) then S=[];return;end*/
+
+	if(maxi < sqrt(eps)){
+		out = NULL;	
+	}
+	else{		
+		rk = dsvda(0,inp,row,col,0,4,U,S,V);
+		
+		/*Will be used in complex numbers*/		
+		//C = (double *)malloc(rk*row*sizeof(double));
+		
+		/*for(i=0;i<row;i++){
+			for(j=0;j<row;j++){
+				printf("%lf ",S[i*row+j]);			
+			}	
+			printf("\n");	
+		}*/
+		/*sqrt of S*/
+		printf("%d ",rk);
+		for(i=0;i<rk;i++){
+			for(j=0;j<rk;j++){
+				A[i*rk+j] = sqrt(S[i*row+j]);
+				//printf("%lf ",A[i*rk+j]);
+			}
+		}
+		for(i=0;i<col*rk;i++){
+			B[i] = V[i];
+			//printf("%lf ",B[i]);	
+		}
+		//printf("\n");
+		/*for(i=0;i<rk;i++){
+			for(j=0;j<rk;j++){
+				//A[i*rk+j] = sqrt(S[i*row+j]);
+				printf("%lf ",A[i*rk+j]);
+			}
+			printf("\n");
+		}*/
+		/*for(i=0;i<col;i++){
+			for(j=0;j<rk;j++){
+				//B[i*col+j] = VT[i*col+j];	
+				printf("%lf ",B[i*col+j]);	
+			}	
+			printf("\n");
+		}*/
+		dmulma(B,col,rk,A,rk,rk,out);
+	}
+}
-- 
cgit 


From 39e059c0f4f7b517718bcc803c1e24873e25b164 Mon Sep 17 00:00:00 2001
From: Sandeep Gupta
Date: Fri, 7 Jul 2017 02:44:54 +0530
Subject: SQROOT

---
 src/c/linearAlgebra/sqroot/dsqroota.c | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

(limited to 'src/c/linearAlgebra/sqroot/dsqroota.c')

diff --git a/src/c/linearAlgebra/sqroot/dsqroota.c b/src/c/linearAlgebra/sqroot/dsqroota.c
index a9062e6..1ba9726 100644
--- a/src/c/linearAlgebra/sqroot/dsqroota.c
+++ b/src/c/linearAlgebra/sqroot/dsqroota.c
@@ -99,7 +99,7 @@ void dsqroota(double *inp,int row,int col, double *out){
 			printf("\n");	
 		}*/
 		/*sqrt of S*/
-		printf("%d ",rk);
+		//printf("%d ",rk);
 		for(i=0;i<rk;i++){
 			for(j=0;j<rk;j++){
 				A[i*rk+j] = sqrt(S[i*row+j]);
-- 
cgit