added matrixExponential and all related sub functions

	matrixExponential just work with real Double for the moment

1 files changed, 34 insertions, 0 deletions

diff --git a/src/matrixOperations/infiniteNorm/cinfnorma.c b/src/matrixOperations/infiniteNorm/cinfnorma.c
new file mode 100644
index 00000000..e930d407
--- /dev/null
+++ b/src/matrixOperations/infiniteNorm/cinfnorma.c
@@ -0,0 +1,34 @@
+/*
+ * Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+ * Copyright (C) 2008-2008 - INRIA - Allan SIMON
+ *
+ * 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
+ *
+ */
+
+#include "infiniteNorm.h"
+
+float cinfnorma(floatComplex* in, int _iRows, int _iCols)
+{
+	int iIndex1 = 0, iIndex2 = 0;
+	float fltTemp = 0;
+	float fltRef = 0;
+
+		for(iIndex1 = 0 ; iIndex1 < _iRows ; iIndex1++)
+		{
+			fltTemp = 0;
+			for(iIndex2 = 0 ; iIndex2 < _iCols ; iIndex2++)
+			{
+				fltTemp += spythags( creals( in[iIndex1 + iIndex2 * _iRows]), 
+                                     cimags( in[iIndex1 + iIndex2 * _iRows]));
+			}
+			if(fltTemp > fltRef)
+				fltRef = fltTemp;
+		}
+	
+	return fltRef;
+}