From 1e3fb075853b2455eecfe45e5f48d3b78f8eecbe Mon Sep 17 00:00:00 2001
From: torset
Date: Tue, 6 Jan 2009 15:34:54 +0000
Subject: fix bug in sqrt

---
 src/elementaryFunctions/sqrt/csqrts.c | 2 +-
 src/elementaryFunctions/sqrt/zsqrts.c | 4 ++--
 2 files changed, 3 insertions(+), 3 deletions(-)

(limited to 'src')

diff --git a/src/elementaryFunctions/sqrt/csqrts.c b/src/elementaryFunctions/sqrt/csqrts.c
index da7000e8..f06d2dd1 100644
--- a/src/elementaryFunctions/sqrt/csqrts.c
+++ b/src/elementaryFunctions/sqrt/csqrts.c
@@ -66,7 +66,7 @@ floatComplex	csqrts(floatComplex in) {
       else
 	{
 	  RealOut	= sabss(ImgIn) / Temp;
-	  ImgOut	= _sign(0.5f, ImgIn) * Temp;
+	  ImgOut	= (_sign(0.5f, ImgIn)) * Temp;
 	}
     }
   else
diff --git a/src/elementaryFunctions/sqrt/zsqrts.c b/src/elementaryFunctions/sqrt/zsqrts.c
index 357a752e..bdfc7a49 100644
--- a/src/elementaryFunctions/sqrt/zsqrts.c
+++ b/src/elementaryFunctions/sqrt/zsqrts.c
@@ -9,7 +9,7 @@
  *  http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
  *
  */
-
+#include <stdio.h>
 #include <math.h>
 #include "sqrt.h"
 #include "lapack.h"
@@ -66,7 +66,7 @@ doubleComplex	zsqrts(doubleComplex in) {
       else
 	{
 	  RealOut	= dabss(ImgIn) / Temp;
-	  ImgOut	= _sign(0.5, ImgIn) * Temp;
+	  ImgOut	= (_sign(0.5, ImgIn)) * Temp;
 	}
     }
   else
-- 
cgit