1 files changed, 100 insertions, 0 deletions
diff --git a/2.3-1/src/c/matrixOperations/logm/zlogma.c b/2.3-1/src/c/matrixOperations/logm/zlogma.c
new file mode 100644
index 00000000..3fb7b274
--- /dev/null
+++ b/2.3-1/src/c/matrixOperations/logm/zlogma.c
@@ -0,0 +1,100 @@
+/*
+ *  Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+ *  Copyright (C) 2008 - INRIA - Arnaud TORSET
+ *
+ *  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 <stdlib.h>
+#include "lapack.h"
+#include "logm.h"
+#include "log.h"
+#include "matrixMultiplication.h"
+#include "matrixTranspose.h"
+#include "matrixInversion.h"
+#include "max.h"
+#include "conj.h"
+#include "spec.h"
+
+void zlogma (doubleComplex* in, int rows, doubleComplex* out){
+
+	/* Algo : */
+	/* find eigenvalues and eigenvectors */
+	/* use formula
+	*	logm = Vp * diag(log(diag(vp)) * inv(Vp) */
+
+
+	int i = 0,j = 0;
+	doubleComplex *eigenvalues, *eigenvectors, *tmp, *tmp2;
+	int hermitienne = 0;
+
+
+	/* Data initialization */
+	eigenvalues = (doubleComplex*) malloc( (unsigned int) (rows*rows)*sizeof(doubleComplex));
+	eigenvectors = (doubleComplex*) malloc( (unsigned int) (rows*rows)*sizeof(doubleComplex));
+	tmp = (doubleComplex*) malloc( (unsigned int) (rows*rows)*sizeof(doubleComplex));
+	tmp2 = (doubleComplex*) malloc( (unsigned int) (rows*rows)*sizeof(doubleComplex));
+
+
+
+
+
+	/* hermitian test */
+	for (i=0;i<rows;i++){
+		for(j=0;j<rows;j++)
+			if ( (zreals(in[i+j*rows])!=zreals(in[j+i*rows])) ||
+			     (zimags(in[i+j*rows])!=zimags(in[j+i*rows])) )  break;
+		if (j!=rows) break;
+	}
+	if ((i==rows) && (j==rows)) {
+	  hermitienne=0;
+	}
+
+	/* find eigenvalues and eigenvectors */
+	zspec2a(in,rows,eigenvalues,eigenvectors);
+
+
+	/* make operation on eigenvalues and eigenvectors */
+	/*
+	for (i=0;i < rows;i++) {
+		eigenvalues[i+i*rows] = zlogs(eigenvalues[i+i*rows]);
+	}
+	*/
+
+	for (i=0;i < rows*rows;i++) {
+	  if (zreals(eigenvalues[i]) == 0 && zimags(eigenvalues[i]) == 0) {
+	    eigenvalues[i]=DoubleComplex(-0.0, 0.0);
+	  }
+	  else {
+	    eigenvalues[i]=zlogs(eigenvalues[i]);
+	  }
+
+	}
+	/*
+	  zloga(eigenvalues, rows*rows, eigenvalues);
+	*/
+
+	zmulma(eigenvectors, rows, rows, eigenvalues, rows, rows, tmp);
+
+	if (hermitienne) {
+		/* we use eigenvalues as a temporary matrix cause it's useless now*/
+		ztransposea(eigenvectors,rows,rows,eigenvalues);
+		zconja(eigenvalues,rows*rows,tmp2);
+	}
+	else {
+	  zinverma(eigenvectors, tmp2, rows);
+	}
+
+	zmulma(tmp, rows, rows, tmp2, rows, rows, out);
+
+
+	free(eigenvalues);
+	free(eigenvectors);
+	free(tmp);
+
+}