From 3c1f8b710b2dde2d43038ebeb962c582ff242dae Mon Sep 17 00:00:00 2001
From: torset
Date: Fri, 20 Feb 2009 10:54:28 +0000
Subject: Update logm

---
 src/matrixOperations/logm/dlogma.c | 10 +++++++---
 1 file changed, 7 insertions(+), 3 deletions(-)

(limited to 'src/matrixOperations/logm/dlogma.c')

diff --git a/src/matrixOperations/logm/dlogma.c b/src/matrixOperations/logm/dlogma.c
index 2798e127..8df8ce3a 100644
--- a/src/matrixOperations/logm/dlogma.c
+++ b/src/matrixOperations/logm/dlogma.c
@@ -14,15 +14,19 @@
 #include "logm.h"
 
 
-void dlogma (double* in, int rows, doubleComplex* out){
-	doubleComplex *inCpx;
+void dlogma (double* in, int rows, double* out){
+	doubleComplex *inCpx,*outCpx;
 	int i;
 	
 	inCpx=malloc((uint)(rows*rows)*sizeof(doubleComplex));
+	outCpx=malloc((uint)(rows*rows)*sizeof(doubleComplex));
 	
 	for (i=0;i<rows*rows;i++) inCpx[i] = DoubleComplex(in[i],0);
 	
-	zlogma(inCpx, rows, out);
+	zlogma(inCpx, rows, outCpx);
+	
+	for (i=0;i<rows*rows;i++) out[i]=zreals(outCpx[i]);
 	
 	free(inCpx);
+	free(outCpx);
 }
-- 
cgit