From 19ee28d3164afe2038eaab82c0805e52fb462050 Mon Sep 17 00:00:00 2001
From: torset
Date: Thu, 29 Jan 2009 13:42:57 +0000
Subject: fix the bug with conjugate in zrdivma.c

---
 src/matrixOperations/division/zrdivma.c | 14 +++++++++++++-
 1 file changed, 13 insertions(+), 1 deletion(-)

(limited to 'src/matrixOperations/division/zrdivma.c')

diff --git a/src/matrixOperations/division/zrdivma.c b/src/matrixOperations/division/zrdivma.c
index 6d48988b..16fb0398 100644
--- a/src/matrixOperations/division/zrdivma.c
+++ b/src/matrixOperations/division/zrdivma.c
@@ -98,7 +98,19 @@ void	zrdivma(	doubleComplex* in1, int lines1, int columns1 ,
 			{
 				cNorm	= 'N';
 				C2F(zgetrs)(&cNorm, &columns2, &lines1, poAf, &columns2, pIpiv, poBt, &columns2, &iInfo);
-				ztransposea(poBt, columns2, lines2, out);
+				{
+				int i,j,ij,ji;
+				for(j = 0 ; j < lines2 ; j++)
+				{
+					for(i = 0 ; i < lines1 ; i++)
+					{
+						ij = i + j * lines1;
+						ji = j + i * max(lines2, columns2);
+						out[ij] 	=  DoubleComplex (  zreals( poBt[ji]) , -zimags ( poBt[ji]));
+					    printf ( "\n\t\t\t<debug>%e + %e\n " ,  zreals( poBt[ji]) , -zimags ( poBt[ji]));
+                   		}
+                   	}
+                   	}
 				
 				iExit = 1;
 			}
-- 
cgit