From 602b6a9a6598265a59bf2c23da1f08ac13e3252e Mon Sep 17 00:00:00 2001
From: jofret
Date: Thu, 27 Nov 2008 15:10:42 +0000
Subject: Correct asin BUG

---
 src/elementaryFunctions/asin/casins.c | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

(limited to 'src/elementaryFunctions/asin/casins.c')

diff --git a/src/elementaryFunctions/asin/casins.c b/src/elementaryFunctions/asin/casins.c
index 957dd8b4..9c15ce53 100644
--- a/src/elementaryFunctions/asin/casins.c
+++ b/src/elementaryFunctions/asin/casins.c
@@ -79,12 +79,12 @@ floatComplex		casins(floatComplex z) {
 
 	  if(dblAbsReal < 1)
 	    /* Am1 = 0.5d0*((y**2)/(R+(x+1.d0))+(y**2)/(S+(1.d0-x))) */
-	    dblImg1 = 0.5f * (dblAbsImg * dblAbsImg / (dblR + (dblAbsReal + 1)) + dblAbsImg * dblAbsImg / (dblS + (dblAbsReal - 1)));
+	    dblImg1 = 0.5f * (dblAbsImg * dblAbsImg / (dblR + (dblAbsReal + 1)) + dblAbsImg * dblAbsImg / (dblS + (1 - dblAbsReal)));
 	  else
 	    /* Am1 = 0.5d0*((y**2)/(R+(x+1.d0))+(S+(x-1.d0))) */
 	    dblImg1 = 0.5f * (dblAbsImg * dblAbsImg / (dblR + (dblAbsReal + 1)) + (dblS + (dblAbsReal - 1)));
 	  /* ai = logp1(Am1 + sqrt(Am1*(A+1.d0))) */
-	  dblTemp = dblImg1 + ssqrts(dblImg1 * dblA + 1);
+	  dblTemp = dblImg1 + ssqrts(dblImg1 * (dblA + 1));
 	  _pdblImg = slog1ps(dblTemp);
 	}
       else
-- 
cgit