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/rowcomp/drowcompa.c | 79 +++++++++++++++++++++++++++++++++
1 file changed, 79 insertions(+)
create mode 100644 src/c/linearAlgebra/rowcomp/drowcompa.c
(limited to 'src/c/linearAlgebra/rowcomp/drowcompa.c')
diff --git a/src/c/linearAlgebra/rowcomp/drowcompa.c b/src/c/linearAlgebra/rowcomp/drowcompa.c
new file mode 100644
index 00000000..3161a2d6
--- /dev/null
+++ b/src/c/linearAlgebra/rowcomp/drowcompa.c
@@ -0,0 +1,79 @@
+/* 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
+ */
+
+/* This function is used to find row compression, range */
+
+#include "rowcomp.h"
+#include <stdio.h>
+#include <stdlib.h>
+#include <math.h>
+#include "svd.h"
+#include "norm.h"
+#include "eye.h"
+#include "matrixTranspose.h"
+#include "qr.h"
+
+/* All variable names, are in consideration of scilab documentation. for reference please check the scilab code.*/
+
+double drowcompa(int ninp,double *A,int row,int col,char *flag,double tol,double *w){
+ double rk;
+ double *U;
+ double *S;
+ double *V;
+ double *q,*r,*e;
+ if(row == 0 || col == 0){
+ w = NULL;
+ return 0;
+ }
+
+ double nA1 = dnorma(A,row,col,1);
+ if(nA1 < sqrt(pow(2,-52))/10){
+ deyea(w,row,col);
+ return 0;
+ }
+ if(ninp == 1){
+ flag = "svd";
+ tol = sqrt(pow(2,-52))*nA1;
+ }
+ else if(ninp == 2){
+ tol = sqrt(pow(2,-52))*nA1;
+ }
+ else{
+ if(tol < 0){ /* if tolerance is negative */
+ printf(" Wrong values for input argument #: Non-negative scalar expected");
+ }
+ }
+ int M = row,N=col;
+ int minMN = min(M,N);
+ char check[3]="qr";
+ //printf(" %s ",flag);
+ if(strcmp(check,flag) == 0){
+ /* calling qr function*/
+ //printf(" * ");
+ q = (double *)malloc(M*min(M,N)*sizeof(double));
+ r = (double *)malloc(minMN*N*sizeof(double));
+ e = (double *)malloc(N*N*sizeof(double));
+ rk = dqra(2,4,A,M,N,tol,q,r,e);
+ memcpy(w,q,row*col*sizeof(double));
+ dtransposea(q,row,row,w);
+ return rk;
+ }
+ else{
+ /* svd function type */
+ U = (double *)malloc(row*row*sizeof(double));
+ S = (double *)malloc(row*col*sizeof(double));
+ V = (double *)malloc(col*col*sizeof(double));
+ rk = dsvda(tol,A,row,col,0,4,U,S,V);
+ dtransposea(U,row,row,w);
+ return rk;
+ }
+}
--
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/rowcomp/drowcompa.c | 79 ---------------------------------
1 file changed, 79 deletions(-)
delete mode 100644 src/c/linearAlgebra/rowcomp/drowcompa.c
(limited to 'src/c/linearAlgebra/rowcomp/drowcompa.c')
diff --git a/src/c/linearAlgebra/rowcomp/drowcompa.c b/src/c/linearAlgebra/rowcomp/drowcompa.c
deleted file mode 100644
index 3161a2d6..00000000
--- a/src/c/linearAlgebra/rowcomp/drowcompa.c
+++ /dev/null
@@ -1,79 +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
- */
-
-/* This function is used to find row compression, range */
-
-#include "rowcomp.h"
-#include <stdio.h>
-#include <stdlib.h>
-#include <math.h>
-#include "svd.h"
-#include "norm.h"
-#include "eye.h"
-#include "matrixTranspose.h"
-#include "qr.h"
-
-/* All variable names, are in consideration of scilab documentation. for reference please check the scilab code.*/
-
-double drowcompa(int ninp,double *A,int row,int col,char *flag,double tol,double *w){
- double rk;
- double *U;
- double *S;
- double *V;
- double *q,*r,*e;
- if(row == 0 || col == 0){
- w = NULL;
- return 0;
- }
-
- double nA1 = dnorma(A,row,col,1);
- if(nA1 < sqrt(pow(2,-52))/10){
- deyea(w,row,col);
- return 0;
- }
- if(ninp == 1){
- flag = "svd";
- tol = sqrt(pow(2,-52))*nA1;
- }
- else if(ninp == 2){
- tol = sqrt(pow(2,-52))*nA1;
- }
- else{
- if(tol < 0){ /* if tolerance is negative */
- printf(" Wrong values for input argument #: Non-negative scalar expected");
- }
- }
- int M = row,N=col;
- int minMN = min(M,N);
- char check[3]="qr";
- //printf(" %s ",flag);
- if(strcmp(check,flag) == 0){
- /* calling qr function*/
- //printf(" * ");
- q = (double *)malloc(M*min(M,N)*sizeof(double));
- r = (double *)malloc(minMN*N*sizeof(double));
- e = (double *)malloc(N*N*sizeof(double));
- rk = dqra(2,4,A,M,N,tol,q,r,e);
- memcpy(w,q,row*col*sizeof(double));
- dtransposea(q,row,row,w);
- return rk;
- }
- else{
- /* svd function type */
- U = (double *)malloc(row*row*sizeof(double));
- S = (double *)malloc(row*col*sizeof(double));
- V = (double *)malloc(col*col*sizeof(double));
- rk = dsvda(tol,A,row,col,0,4,U,S,V);
- dtransposea(U,row,row,w);
- return rk;
- }
-}
--
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/rowcomp/drowcompa.c | 79 +++++++++++++++++++++++++++++++++
1 file changed, 79 insertions(+)
create mode 100644 src/c/linearAlgebra/rowcomp/drowcompa.c
(limited to 'src/c/linearAlgebra/rowcomp/drowcompa.c')
diff --git a/src/c/linearAlgebra/rowcomp/drowcompa.c b/src/c/linearAlgebra/rowcomp/drowcompa.c
new file mode 100644
index 00000000..3161a2d6
--- /dev/null
+++ b/src/c/linearAlgebra/rowcomp/drowcompa.c
@@ -0,0 +1,79 @@
+/* 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
+ */
+
+/* This function is used to find row compression, range */
+
+#include "rowcomp.h"
+#include <stdio.h>
+#include <stdlib.h>
+#include <math.h>
+#include "svd.h"
+#include "norm.h"
+#include "eye.h"
+#include "matrixTranspose.h"
+#include "qr.h"
+
+/* All variable names, are in consideration of scilab documentation. for reference please check the scilab code.*/
+
+double drowcompa(int ninp,double *A,int row,int col,char *flag,double tol,double *w){
+ double rk;
+ double *U;
+ double *S;
+ double *V;
+ double *q,*r,*e;
+ if(row == 0 || col == 0){
+ w = NULL;
+ return 0;
+ }
+
+ double nA1 = dnorma(A,row,col,1);
+ if(nA1 < sqrt(pow(2,-52))/10){
+ deyea(w,row,col);
+ return 0;
+ }
+ if(ninp == 1){
+ flag = "svd";
+ tol = sqrt(pow(2,-52))*nA1;
+ }
+ else if(ninp == 2){
+ tol = sqrt(pow(2,-52))*nA1;
+ }
+ else{
+ if(tol < 0){ /* if tolerance is negative */
+ printf(" Wrong values for input argument #: Non-negative scalar expected");
+ }
+ }
+ int M = row,N=col;
+ int minMN = min(M,N);
+ char check[3]="qr";
+ //printf(" %s ",flag);
+ if(strcmp(check,flag) == 0){
+ /* calling qr function*/
+ //printf(" * ");
+ q = (double *)malloc(M*min(M,N)*sizeof(double));
+ r = (double *)malloc(minMN*N*sizeof(double));
+ e = (double *)malloc(N*N*sizeof(double));
+ rk = dqra(2,4,A,M,N,tol,q,r,e);
+ memcpy(w,q,row*col*sizeof(double));
+ dtransposea(q,row,row,w);
+ return rk;
+ }
+ else{
+ /* svd function type */
+ U = (double *)malloc(row*row*sizeof(double));
+ S = (double *)malloc(row*col*sizeof(double));
+ V = (double *)malloc(col*col*sizeof(double));
+ rk = dsvda(tol,A,row,col,0,4,U,S,V);
+ dtransposea(U,row,row,w);
+ return rk;
+ }
+}
--
cgit