00001 /* 00002 * Copyright 2004 PathScale, Inc. All Rights Reserved. 00003 */ 00004 00005 /* 00006 00007 Copyright (C) 2000, 2001 Silicon Graphics, Inc. All Rights Reserved. 00008 00009 This program is free software; you can redistribute it and/or modify it 00010 under the terms of version 2.1 of the GNU Lesser General Public License 00011 as published by the Free Software Foundation. 00012 00013 This program is distributed in the hope that it would be useful, but 00014 WITHOUT ANY WARRANTY; without even the implied warranty of 00015 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. 00016 00017 Further, this software is distributed without any warranty that it is 00018 free of the rightful claim of any third person regarding infringement 00019 or the like. Any license provided herein, whether implied or 00020 otherwise, applies only to this software file. Patent licenses, if 00021 any, provided herein do not apply to combinations of this program with 00022 other software, or any other product whatsoever. 00023 00024 You should have received a copy of the GNU Lesser General Public 00025 License along with this program; if not, write the Free Software 00026 Foundation, Inc., 59 Temple Place - Suite 330, Boston MA 02111-1307, 00027 USA. 00028 00029 Contact information: Silicon Graphics, Inc., 1600 Amphitheatre Pky, 00030 Mountain View, CA 94043, or: 00031 00032 http://www.sgi.com 00033 00034 For further information regarding this notice, see: 00035 00036 http://oss.sgi.com/projects/GenInfo/NoticeExplan 00037 00038 */ 00039 00040 00041 #include <math.h> 00042 #include "moremath.h" 00043 #include "cmplx.h" 00044 #include "defalias.h" 00045 00046 complex __csqrt(float zreal, float zimag) 00047 { 00048 register float mag; 00049 complex result; 00050 00051 if ((mag = hypotf(zreal, zimag)) == 0.0F) { 00052 result.real = result.imag = 0.0F; 00053 } else if (zreal > 0.0F) { 00054 result.real = sqrtf(0.5F * (mag + zreal)); 00055 result.imag = zimag/result.real/2.0F; 00056 } else { 00057 result.imag = sqrtf(0.5F * (mag - zreal)); 00058 if(zimag < 0.0F) 00059 result.imag = - result.imag; 00060 result.real = zimag/result.imag/2.0F; 00061 } 00062 return result; 00063 } 00064 00065 complex c_sqrt_(complex *z) 00066 { 00067 return __csqrt(z->real, z->imag); 00068 } 00069 00070 #if defined(BUILD_OS_DARWIN) 00071 /* Mach-O doesn't support aliases */ 00072 complex c_sqrt(complex *z) { return c_sqrt_(z); } 00073 #else /* defined(BUILD_OS_DARWIN) */ 00074 defalias(c_sqrt_, c_sqrt); 00075 #endif /* defined(BUILD_OS_DARWIN) */
1.5.6