1 /* 2 * CDDL HEADER START 3 * 4 * The contents of this file are subject to the terms of the 5 * Common Development and Distribution License (the "License"). 6 * You may not use this file except in compliance with the License. 7 * 8 * You can obtain a copy of the license at usr/src/OPENSOLARIS.LICENSE 9 * or http://www.opensolaris.org/os/licensing. 10 * See the License for the specific language governing permissions 11 * and limitations under the License. 12 * 13 * When distributing Covered Code, include this CDDL HEADER in each 14 * file and include the License file at usr/src/OPENSOLARIS.LICENSE. 15 * If applicable, add the following below this CDDL HEADER, with the 16 * fields enclosed by brackets "[]" replaced with your own identifying 17 * information: Portions Copyright [yyyy] [name of copyright owner] 18 * 19 * CDDL HEADER END 20 */ 21 /* 22 * Copyright 2011 Nexenta Systems, Inc. All rights reserved. 23 */ 24 /* 25 * Copyright 2005 Sun Microsystems, Inc. All rights reserved. 26 * Use is subject to license terms. 27 */ 28 29 #pragma weak __erff = erff 30 #pragma weak __erfcf = erfcf 31 32 #include "libm.h" 33 34 #if defined(__i386) && !defined(__amd64) 35 extern int __swapRP(int); 36 #endif 37 38 float 39 erff(float x) { 40 int ix; 41 42 ix = *(int *)&x & ~0x80000000; 43 if (ix > 0x7f800000) /* x is NaN */ 44 return (x * x); 45 return ((float)erf((double)x)); 46 } 47 48 float 49 erfcf(float x) { 50 float f; 51 int ix; 52 #if defined(__i386) && !defined(__amd64) 53 int rp; 54 #endif 55 56 ix = *(int *)&x & ~0x80000000; 57 if (ix > 0x7f800000) /* x is NaN */ 58 return (x * x); 59 60 #if defined(__i386) && !defined(__amd64) 61 rp = __swapRP(fp_extended); 62 #endif 63 f = (float)erfc((double)x); 64 #if defined(__i386) && !defined(__amd64) 65 if (rp != fp_extended) 66 (void) __swapRP(rp); 67 #endif 68 return (f); 69 } 70