e_sqrt.c (97409) | e_sqrt.c (97413) |
---|---|
1/* @(#)e_sqrt.c 5.1 93/09/24 */ 2/* 3 * ==================================================== 4 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved. 5 * 6 * Developed at SunPro, a Sun Microsystems, Inc. business. 7 * Permission to use, copy, modify, and distribute this 8 * software is freely granted, provided that this notice 9 * is preserved. 10 * ==================================================== 11 */ 12 13#ifndef lint | 1/* @(#)e_sqrt.c 5.1 93/09/24 */ 2/* 3 * ==================================================== 4 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved. 5 * 6 * Developed at SunPro, a Sun Microsystems, Inc. business. 7 * Permission to use, copy, modify, and distribute this 8 * software is freely granted, provided that this notice 9 * is preserved. 10 * ==================================================== 11 */ 12 13#ifndef lint |
14static char rcsid[] = "$FreeBSD: head/lib/msun/src/e_sqrt.c 97409 2002-05-28 17:51:46Z alfred $"; | 14static char rcsid[] = "$FreeBSD: head/lib/msun/src/e_sqrt.c 97413 2002-05-28 18:15:04Z alfred $"; |
15#endif 16 17/* __ieee754_sqrt(x) 18 * Return correctly rounded sqrt. 19 * ------------------------------------------ 20 * | Use the hardware sqrt if you have one | 21 * ------------------------------------------ 22 * Method: --- 61 unchanged lines hidden (view full) --- 84 *--------------- 85 */ 86 87#include "math.h" 88#include "math_private.h" 89 90static const double one = 1.0, tiny=1.0e-300; 91 | 15#endif 16 17/* __ieee754_sqrt(x) 18 * Return correctly rounded sqrt. 19 * ------------------------------------------ 20 * | Use the hardware sqrt if you have one | 21 * ------------------------------------------ 22 * Method: --- 61 unchanged lines hidden (view full) --- 84 *--------------- 85 */ 86 87#include "math.h" 88#include "math_private.h" 89 90static const double one = 1.0, tiny=1.0e-300; 91 |
92 double __generic___ieee754_sqrt(double x) | 92double 93__generic___ieee754_sqrt(double x) |
93{ 94 double z; 95 int32_t sign = (int)0x80000000; 96 int32_t ix0,s0,q,m,t,i; 97 u_int32_t r,t1,s1,ix1,q1; 98 99 EXTRACT_WORDS(ix0,ix1,x); 100 --- 344 unchanged lines hidden --- | 94{ 95 double z; 96 int32_t sign = (int)0x80000000; 97 int32_t ix0,s0,q,m,t,i; 98 u_int32_t r,t1,s1,ix1,q1; 99 100 EXTRACT_WORDS(ix0,ix1,x); 101 --- 344 unchanged lines hidden --- |