ecos: packages/language/c/libm/current/src/double/ieee754-api/w

annotate packages/language/c/libm/current/src/double/ieee754-api/w_pow.c @ 66:bf00f99aec69 ecos-sw-2000-02-02

Merge from eCos master repository on 2000-02-02-19:16:44-GMT

author	jlarmour
date	Wed, 02 Feb 2000 19:57:02 +0000
parents	c38311975d4f
children	e0c0827131d1

rev	line source
0 3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	1 //========================================================================
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	2 //
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	3 // w_pow.c
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	4 //
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	5 // Part of the standard mathematical function library
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	6 //
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	7 //========================================================================
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	8 //####COPYRIGHTBEGIN####
64 c38311975d4f Merge from eCos master repository on 2000-01-28-04:28:11-GMT jlarmour parents: 2 diff changeset	9 //
c38311975d4f Merge from eCos master repository on 2000-01-28-04:28:11-GMT jlarmour parents: 2 diff changeset	10 // -------------------------------------------
c38311975d4f Merge from eCos master repository on 2000-01-28-04:28:11-GMT jlarmour parents: 2 diff changeset	11 // The contents of this file are subject to the Red Hat eCos Public License
66 bf00f99aec69 Merge from eCos master repository on 2000-02-02-19:16:44-GMT jlarmour parents: 64 diff changeset	12 // Version 1.1 (the "License"); you may not use this file except in
64 c38311975d4f Merge from eCos master repository on 2000-01-28-04:28:11-GMT jlarmour parents: 2 diff changeset	13 // compliance with the License. You may obtain a copy of the License at
66 bf00f99aec69 Merge from eCos master repository on 2000-02-02-19:16:44-GMT jlarmour parents: 64 diff changeset	14 // http://www.redhat.com/
64 c38311975d4f Merge from eCos master repository on 2000-01-28-04:28:11-GMT jlarmour parents: 2 diff changeset	15 //
66 bf00f99aec69 Merge from eCos master repository on 2000-02-02-19:16:44-GMT jlarmour parents: 64 diff changeset	16 // Software distributed under the License is distributed on an "AS IS"
64 c38311975d4f Merge from eCos master repository on 2000-01-28-04:28:11-GMT jlarmour parents: 2 diff changeset	17 // basis, WITHOUT WARRANTY OF ANY KIND, either express or implied. See the
c38311975d4f Merge from eCos master repository on 2000-01-28-04:28:11-GMT jlarmour parents: 2 diff changeset	18 // License for the specific language governing rights and limitations under
c38311975d4f Merge from eCos master repository on 2000-01-28-04:28:11-GMT jlarmour parents: 2 diff changeset	19 // the License.
c38311975d4f Merge from eCos master repository on 2000-01-28-04:28:11-GMT jlarmour parents: 2 diff changeset	20 //
c38311975d4f Merge from eCos master repository on 2000-01-28-04:28:11-GMT jlarmour parents: 2 diff changeset	21 // The Original Code is eCos - Embedded Configurable Operating System,
c38311975d4f Merge from eCos master repository on 2000-01-28-04:28:11-GMT jlarmour parents: 2 diff changeset	22 // released September 30, 1998.
c38311975d4f Merge from eCos master repository on 2000-01-28-04:28:11-GMT jlarmour parents: 2 diff changeset	23 //
c38311975d4f Merge from eCos master repository on 2000-01-28-04:28:11-GMT jlarmour parents: 2 diff changeset	24 // The Initial Developer of the Original Code is Red Hat.
c38311975d4f Merge from eCos master repository on 2000-01-28-04:28:11-GMT jlarmour parents: 2 diff changeset	25 // Portions created by Red Hat are
c38311975d4f Merge from eCos master repository on 2000-01-28-04:28:11-GMT jlarmour parents: 2 diff changeset	26 // Copyright (C) 1998, 1999, 2000 Red Hat, Inc.
c38311975d4f Merge from eCos master repository on 2000-01-28-04:28:11-GMT jlarmour parents: 2 diff changeset	27 // All Rights Reserved.
c38311975d4f Merge from eCos master repository on 2000-01-28-04:28:11-GMT jlarmour parents: 2 diff changeset	28 // -------------------------------------------
c38311975d4f Merge from eCos master repository on 2000-01-28-04:28:11-GMT jlarmour parents: 2 diff changeset	29 //
0 3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	30 //####COPYRIGHTEND####
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	31 //========================================================================
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	32 //#####DESCRIPTIONBEGIN####
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	33 //
2 443894e2e912 Block commit of eCos version 1.2.1 jlarmour parents: 0 diff changeset	34 // Author(s): jlarmour
443894e2e912 Block commit of eCos version 1.2.1 jlarmour parents: 0 diff changeset	35 // Contributors: jlarmour
0 3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	36 // Date: 1998-02-13
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	37 // Purpose:
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	38 // Description:
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	39 // Usage:
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	40 //
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	41 //####DESCRIPTIONEND####
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	42 //
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	43 //========================================================================
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	44
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	45 // CONFIGURATION
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	46
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	47 #include <pkgconf/libm.h> // Configuration header
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	48
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	49 // Include the Math library?
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	50 #ifdef CYGPKG_LIBM
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	51
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	52 // Derived from code with the following copyright
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	53
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	54
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	55
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	56 /* @(#)w_pow.c 1.3 95/01/18 */
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	57 /*
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	58 * ====================================================
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	59 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	60 *
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	61 * Developed at SunSoft, a Sun Microsystems, Inc. business.
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	62 * Permission to use, copy, modify, and distribute this
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	63 * software is freely granted, provided that this notice
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	64 * is preserved.
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	65 * ====================================================
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	66 */
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	67
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	68 /*
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	69 * wrapper pow(x,y) return x**y
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	70 */
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	71
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	72 #include "mathincl/fdlibm.h"
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	73
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	74
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	75 double pow(double x, double y) /* wrapper pow */
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	76 {
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	77 #ifdef CYGSEM_LIBM_COMPAT_IEEE_ONLY
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	78 return __ieee754_pow(x,y);
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	79 #else
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	80 double z;
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	81 z=__ieee754_pow(x,y);
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	82 if(cyg_libm_get_compat_mode() == CYGNUM_LIBM_COMPAT_IEEE\|\|
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	83 isnan(y))
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	84 return z;
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	85 if(isnan(x)) {
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	86 if(y==0.0)
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	87 return __kernel_standard(x,y,42); /* pow(NaN,0.0) */
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	88 else
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	89 return z;
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	90 }
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	91 if(x==0.0){
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	92 if(y==0.0)
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	93 return __kernel_standard(x,y,20); /* pow(0.0,0.0) */
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	94 if(finite(y)&&y<0.0)
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	95 return __kernel_standard(x,y,23); /* pow(0.0,negative) */
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	96 return z;
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	97 }
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	98 if(!finite(z)) {
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	99 if(finite(x)&&finite(y)) {
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	100 if(isnan(z))
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	101 return __kernel_standard(x,y,24);/* pow neg*non-int /
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	102 else
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	103 return __kernel_standard(x,y,21); /* pow overflow */
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	104 }
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	105 }
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	106 if(z==0.0&&finite(x)&&finite(y))
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	107 return __kernel_standard(x,y,22); /* pow underflow */
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	108 return z;
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	109 #endif
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	110 }
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	111
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	112 #endif // ifdef CYGPKG_LIBM
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	113
3111d98ba7b3 Initial commit of eCos version 1.1 jlarmour parents: diff changeset	114 // EOF w_pow.c

Mercurial > ecos

annotate packages/language/c/libm/current/src/double/ieee754-api/w_pow.c @ 66:bf00f99aec69 ecos-sw-2000-02-02