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/householder/dhouseholdera.c | 90 +++++++++++++++++++++++++
 1 file changed, 90 insertions(+)
 create mode 100644 src/c/linearAlgebra/householder/dhouseholdera.c

(limited to 'src/c/linearAlgebra/householder/dhouseholdera.c')

diff --git a/src/c/linearAlgebra/householder/dhouseholdera.c b/src/c/linearAlgebra/householder/dhouseholdera.c
new file mode 100644
index 0000000..5a98bfa
--- /dev/null
+++ b/src/c/linearAlgebra/householder/dhouseholdera.c
@@ -0,0 +1,90 @@
+/* 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
+ */
+
+/* Householder orthogonal reflexion matrix */
+
+/*
+Syntax :- 
+    //u=householder(v [,w])
+    //Description
+    //given 2 column vectors v w of same size householder(v,w) returns a unitary
+    //column vector u, such that (eye-2*u*u')*v is proportional to w.
+    //(eye-2*u*u') is the orthogonal Householder reflexion matrix
+    //
+    // w default value is eye(v). In this case vector (eye-2*u*u')*v is the
+    // vector  eye(v)*(+-norm(v))
+*/
+#include <stdio.h>
+#include <stdlib.h>
+#include "householder.h"
+#include "eye.h"
+#include "matrixTranspose.h"
+#include "matrixMultiplication.h"
+#include <math.h>
+#include "norm.h"
+
+void dhouseholdera(int ninp,double *inp1,int row,double *inp2,double *out1){
+
+	int i,j;
+	double *x;
+	x = (double *)malloc(row*sizeof(double));
+	
+	for(i=0;i<row;i++)
+		x[i] = inp1[i];	
+	
+	if(ninp < 2){
+		deyea(x,row,1);	
+	}
+	else{
+		for(i=0;i<row;i++){
+			x[i] = inp2[i];		
+		}	
+	}
+
+	/* vt transpose of inp1 */
+	
+	double *vt;
+	vt = (double *)malloc(row*sizeof(double));
+	dtransposea(inp1,row,1,vt);
+	
+	/*wt transpose of inp2 */	
+	double *wt;
+	wt = (double *)malloc(row*sizeof(double));
+	dtransposea(x,row,1,wt);
+	
+	/* vvt = inp1*vt */
+	double *vvt;
+	vvt = (double *)malloc(1*1*sizeof(double));
+	dmulma(vt,1,row,inp1,row,1,vvt);
+	
+	/* wwt = inp2*wt */
+	double *wwt;
+	wwt = (double *)malloc(1*1*sizeof(double));
+	dmulma(wt,1,row,x,row,1,wwt);
+	
+	/* a=-sqrt((v'*v)/(w'*w)) */
+	double a;
+	//a = (double *)malloc(1*1*sizeof(double));
+	a = -sqrt(vvt[0]/wwt[0]);
+
+	for(i=0;i<row;i++){
+		out1[i] = x[i]*a+inp1[i];	
+	}
+	
+	/* norm of out1 */
+	double r;
+	r = dnorma(out1,row,1,2);
+
+	for(i=0;i<row;i++){
+		out1[i]=out1[i]/r;	
+	}
+}
-- 
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/householder/dhouseholdera.c | 90 -------------------------
 1 file changed, 90 deletions(-)
 delete mode 100644 src/c/linearAlgebra/householder/dhouseholdera.c

(limited to 'src/c/linearAlgebra/householder/dhouseholdera.c')

diff --git a/src/c/linearAlgebra/householder/dhouseholdera.c b/src/c/linearAlgebra/householder/dhouseholdera.c
deleted file mode 100644
index 5a98bfa..0000000
--- a/src/c/linearAlgebra/householder/dhouseholdera.c
+++ /dev/null
@@ -1,90 +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
- */
-
-/* Householder orthogonal reflexion matrix */
-
-/*
-Syntax :- 
-    //u=householder(v [,w])
-    //Description
-    //given 2 column vectors v w of same size householder(v,w) returns a unitary
-    //column vector u, such that (eye-2*u*u')*v is proportional to w.
-    //(eye-2*u*u') is the orthogonal Householder reflexion matrix
-    //
-    // w default value is eye(v). In this case vector (eye-2*u*u')*v is the
-    // vector  eye(v)*(+-norm(v))
-*/
-#include <stdio.h>
-#include <stdlib.h>
-#include "householder.h"
-#include "eye.h"
-#include "matrixTranspose.h"
-#include "matrixMultiplication.h"
-#include <math.h>
-#include "norm.h"
-
-void dhouseholdera(int ninp,double *inp1,int row,double *inp2,double *out1){
-
-	int i,j;
-	double *x;
-	x = (double *)malloc(row*sizeof(double));
-	
-	for(i=0;i<row;i++)
-		x[i] = inp1[i];	
-	
-	if(ninp < 2){
-		deyea(x,row,1);	
-	}
-	else{
-		for(i=0;i<row;i++){
-			x[i] = inp2[i];		
-		}	
-	}
-
-	/* vt transpose of inp1 */
-	
-	double *vt;
-	vt = (double *)malloc(row*sizeof(double));
-	dtransposea(inp1,row,1,vt);
-	
-	/*wt transpose of inp2 */	
-	double *wt;
-	wt = (double *)malloc(row*sizeof(double));
-	dtransposea(x,row,1,wt);
-	
-	/* vvt = inp1*vt */
-	double *vvt;
-	vvt = (double *)malloc(1*1*sizeof(double));
-	dmulma(vt,1,row,inp1,row,1,vvt);
-	
-	/* wwt = inp2*wt */
-	double *wwt;
-	wwt = (double *)malloc(1*1*sizeof(double));
-	dmulma(wt,1,row,x,row,1,wwt);
-	
-	/* a=-sqrt((v'*v)/(w'*w)) */
-	double a;
-	//a = (double *)malloc(1*1*sizeof(double));
-	a = -sqrt(vvt[0]/wwt[0]);
-
-	for(i=0;i<row;i++){
-		out1[i] = x[i]*a+inp1[i];	
-	}
-	
-	/* norm of out1 */
-	double r;
-	r = dnorma(out1,row,1,2);
-
-	for(i=0;i<row;i++){
-		out1[i]=out1[i]/r;	
-	}
-}
-- 
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/householder/dhouseholdera.c | 90 +++++++++++++++++++++++++
 1 file changed, 90 insertions(+)
 create mode 100644 src/c/linearAlgebra/householder/dhouseholdera.c

(limited to 'src/c/linearAlgebra/householder/dhouseholdera.c')

diff --git a/src/c/linearAlgebra/householder/dhouseholdera.c b/src/c/linearAlgebra/householder/dhouseholdera.c
new file mode 100644
index 0000000..5a98bfa
--- /dev/null
+++ b/src/c/linearAlgebra/householder/dhouseholdera.c
@@ -0,0 +1,90 @@
+/* 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
+ */
+
+/* Householder orthogonal reflexion matrix */
+
+/*
+Syntax :- 
+    //u=householder(v [,w])
+    //Description
+    //given 2 column vectors v w of same size householder(v,w) returns a unitary
+    //column vector u, such that (eye-2*u*u')*v is proportional to w.
+    //(eye-2*u*u') is the orthogonal Householder reflexion matrix
+    //
+    // w default value is eye(v). In this case vector (eye-2*u*u')*v is the
+    // vector  eye(v)*(+-norm(v))
+*/
+#include <stdio.h>
+#include <stdlib.h>
+#include "householder.h"
+#include "eye.h"
+#include "matrixTranspose.h"
+#include "matrixMultiplication.h"
+#include <math.h>
+#include "norm.h"
+
+void dhouseholdera(int ninp,double *inp1,int row,double *inp2,double *out1){
+
+	int i,j;
+	double *x;
+	x = (double *)malloc(row*sizeof(double));
+	
+	for(i=0;i<row;i++)
+		x[i] = inp1[i];	
+	
+	if(ninp < 2){
+		deyea(x,row,1);	
+	}
+	else{
+		for(i=0;i<row;i++){
+			x[i] = inp2[i];		
+		}	
+	}
+
+	/* vt transpose of inp1 */
+	
+	double *vt;
+	vt = (double *)malloc(row*sizeof(double));
+	dtransposea(inp1,row,1,vt);
+	
+	/*wt transpose of inp2 */	
+	double *wt;
+	wt = (double *)malloc(row*sizeof(double));
+	dtransposea(x,row,1,wt);
+	
+	/* vvt = inp1*vt */
+	double *vvt;
+	vvt = (double *)malloc(1*1*sizeof(double));
+	dmulma(vt,1,row,inp1,row,1,vvt);
+	
+	/* wwt = inp2*wt */
+	double *wwt;
+	wwt = (double *)malloc(1*1*sizeof(double));
+	dmulma(wt,1,row,x,row,1,wwt);
+	
+	/* a=-sqrt((v'*v)/(w'*w)) */
+	double a;
+	//a = (double *)malloc(1*1*sizeof(double));
+	a = -sqrt(vvt[0]/wwt[0]);
+
+	for(i=0;i<row;i++){
+		out1[i] = x[i]*a+inp1[i];	
+	}
+	
+	/* norm of out1 */
+	double r;
+	r = dnorma(out1,row,1,2);
+
+	for(i=0;i<row;i++){
+		out1[i]=out1[i]/r;	
+	}
+}
-- 
cgit