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/matrixOperations/norm/dnorma.c | 12 ++++++++++--
 1 file changed, 10 insertions(+), 2 deletions(-)

(limited to 'src/c/matrixOperations/norm')

diff --git a/src/c/matrixOperations/norm/dnorma.c b/src/c/matrixOperations/norm/dnorma.c
index c912f85..2bb9b8b 100644
--- a/src/c/matrixOperations/norm/dnorma.c
+++ b/src/c/matrixOperations/norm/dnorma.c
@@ -6,6 +6,8 @@
  are also available at
  http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
  Author: Siddhesh Wani
+
+ Edited by: Sandeep Gupta
  Organization: FOSSEE, IIT Bombay
  Email: toolbox@scilab.in
  */
@@ -14,13 +16,17 @@
  /*Acceptable norms are: 1, 2, 'inf', 'fro', */
 
 #include "norm.h"
+#include "svd.h"
 #include <math.h>
+#include <stdio.h>
+#include<stdlib.h>
 
 double dnorma (double *in, int row, int col, int norm)
 {
 	double res = 0, sum = 0;
 	int col_count, row_count;
-
+	double *S;
+	S = (double *)malloc(min(row,col)*sizeof(double));
 	switch (norm)
 	{
 		case 1:  /*largest column sum*/
@@ -37,6 +43,8 @@ double dnorma (double *in, int row, int col, int norm)
 			break;
 
 		case 2:  /*Largest singular value of the matrix*/
+			dsvda(0,in,row,col,0,1,S,NULL,NULL);			
+			res = S[0];
 			break;
 
 		case 3: /*inf: largest row sum*/
@@ -69,4 +77,4 @@ double dnorma (double *in, int row, int col, int norm)
 
 
 	return res;
-}
\ No newline at end of file
+}
-- 
cgit