Index: /trunk/src/emx/ChangeLog.LIBC
===================================================================
--- /trunk/src/emx/ChangeLog.LIBC	(revision 2012)
+++ /trunk/src/emx/ChangeLog.LIBC	(revision 2013)
@@ -1,3 +1,18 @@
 /* $Id$ */
+
+2005-06-11: knut st. osmundsen <bird-gccos2-spam@anduin.net>
+    - libc:
+        o The smallcnv stuff is gone, waste of space since we're not really into 
+          static linking. Should've done this when switching to gdtoa.
+        o Ported the msun and some other math bits from current FreeBSD libc.
+          This mostly replaces the EMX math, except for some double long variants
+          and some EMX specific ones. This is a risky change, but it's needed
+          if we're gonna catch up on the C99 bits and make libstdc++ happy.
+          There are some new headers, complex.h and fenv.h being the importan ones.
+          float.h have gotten some updates while math.h is totally replaced. 
+          More than 100 standard functions were added!
+          NOTE: msun is using MATH_ERREXCEPT (math_errhandling) and not 
+                MATH_ERRNO for error reporting.
+          I really hope this is not gonna uncover more OS/2 387 mess.
 
 2005-06-05: knut st. osmundsen <bird-gccos2-spam@anduin.net>
Index: /trunk/src/emx/include/386/asm.h
===================================================================
--- /trunk/src/emx/include/386/asm.h	(revision 2013)
+++ /trunk/src/emx/include/386/asm.h	(revision 2013)
@@ -0,0 +1,174 @@
+/*-
+ * Copyright (c) 1990 The Regents of the University of California.
+ * All rights reserved.
+ *
+ * This code is derived from software contributed to Berkeley by
+ * William Jolitz.
+ *
+ * Redistribution and use in source and binary forms, with or without
+ * modification, are permitted provided that the following conditions
+ * are met:
+ * 1. Redistributions of source code must retain the above copyright
+ *    notice, this list of conditions and the following disclaimer.
+ * 2. Redistributions in binary form must reproduce the above copyright
+ *    notice, this list of conditions and the following disclaimer in the
+ *    documentation and/or other materials provided with the distribution.
+ * 3. All advertising materials mentioning features or use of this software
+ *    must display the following acknowledgement:
+ *	This product includes software developed by the University of
+ *	California, Berkeley and its contributors.
+ * 4. Neither the name of the University nor the names of its contributors
+ *    may be used to endorse or promote products derived from this software
+ *    without specific prior written permission.
+ *
+ * THIS SOFTWARE IS PROVIDED BY THE REGENTS AND CONTRIBUTORS ``AS IS'' AND
+ * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
+ * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
+ * ARE DISCLAIMED.  IN NO EVENT SHALL THE REGENTS OR CONTRIBUTORS BE LIABLE
+ * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
+ * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
+ * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
+ * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
+ * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
+ * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
+ * SUCH DAMAGE.
+ *
+ *	from: @(#)DEFS.h	5.1 (Berkeley) 4/23/90
+ * $FreeBSD: src/sys/i386/include/asm.h,v 1.10 2002/09/23 07:27:45 peter Exp $
+ */
+
+#ifndef _MACHINE_ASM_H_
+#define	_MACHINE_ASM_H_
+
+#include <sys/cdefs.h>
+
+#ifdef PIC
+#define	PIC_PROLOGUE	\
+	pushl	%ebx;	\
+	call	1f;	\
+1:			\
+	popl	%ebx;	\
+	addl	$_GLOBAL_OFFSET_TABLE_+[.-1b],%ebx
+#define	PIC_EPILOGUE	\
+	popl	%ebx
+#define	PIC_PLT(x)	x@PLT
+#define	PIC_GOT(x)	x@GOT(%ebx)
+#else
+#define	PIC_PROLOGUE
+#define	PIC_EPILOGUE
+#define	PIC_PLT(x)	x
+#define	PIC_GOT(x)	x
+#endif
+
+/*
+ * CNAME and HIDENAME manage the relationship between symbol names in C
+ * and the equivalent assembly language names.  CNAME is given a name as
+ * it would be used in a C program.  It expands to the equivalent assembly
+ * language name.  HIDENAME is given an assembly-language name, and expands
+ * to a possibly-modified form that will be invisible to C programs.
+ */
+#define CNAME(csym)		__CONCAT(_,csym) /* bird, we need underscore */
+#define STDCNAME(csym)		__CONCAT(__std_,csym) /* bird, same as _STD() */
+#define HIDENAME(asmsym)	__CONCAT(.,asmsym)
+
+/* XXX should use .p2align 4,0x90 for -m486. */
+/*#define _START_ENTRY	.text; .p2align 2,0x90*/
+#define _START_ENTRY	.text; .p2align 4,0xcc /* bird */
+
+#define _ENTRY(x)	_START_ENTRY; \
+			.globl CNAME(x); .type CNAME(x),@function; CNAME(x):
+
+#ifdef PROF
+#define	ALTENTRY(x)	_ENTRY(x); \
+			pushl %ebp; movl %esp,%ebp; \
+			call PIC_PLT(HIDENAME(mcount)); \
+			popl %ebp; \
+			jmp 9f
+#define	ENTRY(x)	_ENTRY(x); \
+			pushl %ebp; movl %esp,%ebp; \
+			call PIC_PLT(HIDENAME(mcount)); \
+			popl %ebp; \
+			9:
+#else
+#define	ALTENTRY(x)	_ENTRY(x)
+#define	ENTRY(x)	_ENTRY(x)
+#endif
+#define	STDENTRY(x)	ENTRY(__CONCAT(_std_,x)) /* libc */
+
+#define RCSID(x)	.text; .asciz x
+
+#undef __FBSDID
+#if !defined(lint) && !defined(STRIP_FBSDID)
+#define __FBSDID(s)	.ident s
+#else
+#define __FBSDID(s)	/* nothing */
+#endif /* not lint and not STRIP_FBSDID */
+
+#ifdef _ARCH_INDIRECT
+/*
+ * Generate code to select between the generic functions and _ARCH_INDIRECT
+ * specific ones.
+ * XXX nested __CONCATs don't work with non-ANSI cpp's.
+ */
+#define	ANAME(x)	CNAME(__CONCAT(__CONCAT(__,_ARCH_INDIRECT),x))
+#define	ASELNAME(x)	CNAME(__CONCAT(__arch_select_,x))
+#define	AVECNAME(x)	CNAME(__CONCAT(__arch_,x))
+#define	GNAME(x)	CNAME(__CONCAT(__generic_,x))
+
+/* Don't bother profiling this. */
+#ifdef PIC
+#define	ARCH_DISPATCH(x) \
+			_START_ENTRY; \
+			.globl CNAME(x); .type CNAME(x),@function; CNAME(x): ; \
+			PIC_PROLOGUE; \
+			movl PIC_GOT(AVECNAME(x)),%eax; \
+			PIC_EPILOGUE; \
+			jmpl *(%eax)
+
+#define	ARCH_SELECT(x)	_START_ENTRY; \
+			.type ASELNAME(x),@function; \
+			ASELNAME(x): \
+			PIC_PROLOGUE; \
+			call PIC_PLT(CNAME(__get_hw_float)); \
+			testl %eax,%eax; \
+			movl PIC_GOT(ANAME(x)),%eax; \
+			jne 8f; \
+			movl PIC_GOT(GNAME(x)),%eax; \
+			8: \
+			movl PIC_GOT(AVECNAME(x)),%edx; \
+			movl %eax,(%edx); \
+			PIC_EPILOGUE; \
+			jmpl *%eax
+#else /* !PIC */
+#define	ARCH_DISPATCH(x) \
+			_START_ENTRY; \
+			.globl CNAME(x); .type CNAME(x),@function; CNAME(x): ; \
+			jmpl *AVECNAME(x)
+
+#define	ARCH_SELECT(x)	_START_ENTRY; \
+			.type ASELNAME(x),@function; \
+			ASELNAME(x): \
+			call CNAME(__get_hw_float); \
+			testl %eax,%eax; \
+			movl $ANAME(x),%eax; \
+			jne 8f; \
+			movl $GNAME(x),%eax; \
+			8: \
+			movl %eax,AVECNAME(x); \
+			jmpl *%eax
+#endif /* PIC */
+
+#define	ARCH_VECTOR(x)	.data; .p2align 2; \
+			.globl AVECNAME(x); \
+			.type AVECNAME(x),@object; \
+			.size AVECNAME(x),4; \
+			AVECNAME(x): .long ASELNAME(x)
+
+#undef _ENTRY
+#define	_ENTRY(x)	ARCH_VECTOR(x); ARCH_SELECT(x); ARCH_DISPATCH(x); \
+			_START_ENTRY; \
+			.globl ANAME(x); .type ANAME(x),@function; ANAME(x):
+
+#endif /* _ARCH_INDIRECT */
+
+#endif /* !_MACHINE_ASM_H_ */
Index: /trunk/src/emx/include/386/npx.h
===================================================================
--- /trunk/src/emx/include/386/npx.h	(revision 2013)
+++ /trunk/src/emx/include/386/npx.h	(revision 2013)
@@ -0,0 +1,158 @@
+/*-
+ * Copyright (c) 1990 The Regents of the University of California.
+ * All rights reserved.
+ *
+ * This code is derived from software contributed to Berkeley by
+ * William Jolitz.
+ *
+ * Redistribution and use in source and binary forms, with or without
+ * modification, are permitted provided that the following conditions
+ * are met:
+ * 1. Redistributions of source code must retain the above copyright
+ *    notice, this list of conditions and the following disclaimer.
+ * 2. Redistributions in binary form must reproduce the above copyright
+ *    notice, this list of conditions and the following disclaimer in the
+ *    documentation and/or other materials provided with the distribution.
+ * 3. All advertising materials mentioning features or use of this software
+ *    must display the following acknowledgement:
+ *	This product includes software developed by the University of
+ *	California, Berkeley and its contributors.
+ * 4. Neither the name of the University nor the names of its contributors
+ *    may be used to endorse or promote products derived from this software
+ *    without specific prior written permission.
+ *
+ * THIS SOFTWARE IS PROVIDED BY THE REGENTS AND CONTRIBUTORS ``AS IS'' AND
+ * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
+ * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
+ * ARE DISCLAIMED.  IN NO EVENT SHALL THE REGENTS OR CONTRIBUTORS BE LIABLE
+ * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
+ * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
+ * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
+ * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
+ * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
+ * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
+ * SUCH DAMAGE.
+ *
+ *	from: @(#)npx.h	5.3 (Berkeley) 1/18/91
+ * $FreeBSD: src/sys/i386/include/npx.h,v 1.26 2002/09/23 05:55:10 peter Exp $
+ */
+
+/*
+ * 287/387 NPX Coprocessor Data Structures and Constants
+ * W. Jolitz 1/90
+ */
+
+/** @file
+ * FreeBSD 5.1
+ */
+
+#ifndef _MACHINE_NPX_H_
+#define	_MACHINE_NPX_H_
+
+/* Environment information of floating point unit */
+struct env87 {
+	long	en_cw;		/* control word (16bits) */
+	long	en_sw;		/* status word (16bits) */
+	long	en_tw;		/* tag word (16bits) */
+	long	en_fip;		/* floating point instruction pointer */
+	u_short	en_fcs;		/* floating code segment selector */
+	u_short	en_opcode;	/* opcode last executed (11 bits ) */
+	long	en_foo;		/* floating operand offset */
+	long	en_fos;		/* floating operand segment selector */
+};
+
+/* Contents of each floating point accumulator */
+struct fpacc87 {
+#ifdef dontdef /* too unportable */
+	u_long	fp_mantlo;	/* mantissa low (31:0) */
+	u_long	fp_manthi;	/* mantissa high (63:32) */
+	int	fp_exp:15;	/* exponent */
+	int	fp_sgn:1;	/* mantissa sign */
+#else
+	u_char	fp_bytes[10];
+#endif
+};
+
+/* Floating point context */
+struct save87 {
+	struct	env87 sv_env;	/* floating point control/status */
+	struct	fpacc87	sv_ac[8];	/* accumulator contents, 0-7 */
+	u_char	sv_pad0[4];	/* padding for (now unused) saved status word */
+	/*
+	 * Bogus padding for emulators.  Emulators should use their own
+	 * struct and arrange to store into this struct (ending here)
+	 * before it is inspected for ptracing or for core dumps.  Some
+	 * emulators overwrite the whole struct.  We have no good way of
+	 * knowing how much padding to leave.  Leave just enough for the
+	 * GPL emulator's i387_union (176 bytes total).
+	 */
+	u_char	sv_pad[64];	/* padding; used by emulators */
+};
+
+struct  envxmm {
+	u_int16_t	en_cw;		/* control word (16bits) */
+	u_int16_t	en_sw;		/* status word (16bits) */
+	u_int16_t	en_tw;		/* tag word (16bits) */
+	u_int16_t	en_opcode;	/* opcode last executed (11 bits ) */
+	u_int32_t	en_fip;		/* floating point instruction pointer */
+	u_int16_t	en_fcs;		/* floating code segment selector */
+	u_int16_t	en_pad0;	/* padding */
+	u_int32_t	en_foo;		/* floating operand offset */
+	u_int16_t	en_fos;		/* floating operand segment selector */
+	u_int16_t	en_pad1;	/* padding */
+	u_int32_t	en_mxcsr;	/* SSE sontorol/status register */
+	u_int32_t	en_pad2;	/* padding */
+};
+
+/* Contents of each SSE extended accumulator */
+struct  xmmacc {
+	u_char	xmm_bytes[16];
+};
+
+struct  savexmm {
+	struct	envxmm	sv_env;
+	struct {
+		struct fpacc87	fp_acc;
+		u_char		fp_pad[6];      /* padding */
+	} sv_fp[8];
+	struct xmmacc	sv_xmm[8];
+	u_char sv_pad[224];
+} __aligned(16);
+
+union	savefpu {
+	struct	save87	sv_87;
+	struct	savexmm	sv_xmm;
+};
+
+/*
+ * The hardware default control word for i387's and later coprocessors is
+ * 0x37F, giving:
+ *
+ *	round to nearest
+ *	64-bit precision
+ *	all exceptions masked.
+ *
+ * We modify the affine mode bit and precision bits in this to give:
+ *
+ *	affine mode for 287's (if they work at all) (1 in bitfield 1<<12)
+ *	53-bit precision (2 in bitfield 3<<8)
+ *
+ * 64-bit precision often gives bad results with high level languages
+ * because it makes the results of calculations depend on whether
+ * intermediate values are stored in memory or in FPU registers.
+ */
+#define	__INITIAL_NPXCW__	0x127F
+
+#ifdef _KERNEL
+int	npxdna(void);
+void	npxdrop(void);
+void	npxexit(struct thread *td);
+int	npxformat(void);
+int	npxgetregs(struct thread *td, union savefpu *addr);
+void	npxinit(u_short control);
+void	npxsave(union savefpu *addr);
+void	npxsetregs(struct thread *td, union savefpu *addr);
+int	npxtrap(void);
+#endif
+
+#endif /* !_MACHINE_NPX_H_ */
Index: /trunk/src/emx/include/complex.h
===================================================================
--- /trunk/src/emx/include/complex.h	(revision 2013)
+++ /trunk/src/emx/include/complex.h	(revision 2013)
@@ -0,0 +1,65 @@
+/*-
+ * Copyright (c) 2001 The FreeBSD Project.
+ * All rights reserved.
+ *
+ * Redistribution and use in source and binary forms, with or without
+ * modification, are permitted provided that the following conditions
+ * are met:
+ * 1. Redistributions of source code must retain the above copyright
+ *    notice, this list of conditions and the following disclaimer.
+ * 2. Redistributions in binary form must reproduce the above copyright
+ *    notice, this list of conditions and the following disclaimer in the
+ *    documentation and/or other materials provided with the distribution.
+ *
+ * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``AS IS'' AND
+ * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
+ * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
+ * ARE DISCLAIMED.  IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE
+ * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
+ * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
+ * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
+ * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
+ * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
+ * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
+ * SUCH DAMAGE.
+ *
+ * $FreeBSD: src/include/complex.h,v 1.6 2004/08/14 18:03:21 stefanf Exp $
+ */
+
+/**
+ * FreeBSD 5.3
+ */
+
+#ifndef _COMPLEX_H
+#define	_COMPLEX_H
+
+#ifdef __GNUC__
+#if __STDC_VERSION__ < 199901
+#define	_Complex	__complex__
+#endif
+#define	_Complex_I	1.0fi
+#endif
+
+#define	complex		_Complex
+#define	I		_Complex_I
+
+#include <sys/cdefs.h>
+
+__BEGIN_DECLS
+
+double		cabs(double complex);
+float		cabsf(float complex);
+double		cimag(double complex);
+float		cimagf(float complex);
+long double	cimagl(long double complex);
+double complex	conj(double complex);
+float complex	conjf(float complex);
+long double complex
+		conjl(long double complex);
+double		creal(double complex);
+float		crealf(float complex);
+long double	creall(long double complex);
+
+__END_DECLS
+
+#endif /* _COMPLEX_H */
Index: /trunk/src/emx/include/emx/asm386.h
===================================================================
--- /trunk/src/emx/include/emx/asm386.h	(revision 2012)
+++ /trunk/src/emx/include/emx/asm386.h	(revision 2013)
@@ -63,7 +63,8 @@
 /* Cf. <math.h> */
 
-#define FP_ZERO         0
-#define FP_SUBNORMAL    1
-#define FP_NORMAL       2
-#define FP_INFINITE     3
-#define FP_NAN          4
+#define	FP_INFINITE	0x01
+#define	FP_NAN		0x02
+#define	FP_NORMAL	0x04
+#define	FP_SUBNORMAL	0x08
+#define	FP_ZERO		0x10
+
Index: /trunk/src/emx/include/fenv.h
===================================================================
--- /trunk/src/emx/include/fenv.h	(revision 2013)
+++ /trunk/src/emx/include/fenv.h	(revision 2013)
@@ -0,0 +1,240 @@
+/*-
+ * Copyright (c) 2004-2005 David Schultz <das@FreeBSD.ORG>
+ * All rights reserved.
+ *
+ * Redistribution and use in source and binary forms, with or without
+ * modification, are permitted provided that the following conditions
+ * are met:
+ * 1. Redistributions of source code must retain the above copyright
+ *    notice, this list of conditions and the following disclaimer.
+ * 2. Redistributions in binary form must reproduce the above copyright
+ *    notice, this list of conditions and the following disclaimer in the
+ *    documentation and/or other materials provided with the distribution.
+ *
+ * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``AS IS'' AND
+ * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
+ * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
+ * ARE DISCLAIMED.  IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE
+ * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
+ * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
+ * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
+ * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
+ * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
+ * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
+ * SUCH DAMAGE.
+ *
+ * $FreeBSD: src/lib/msun/i387/fenv.h,v 1.4 2005/03/17 22:21:46 das Exp $
+ */
+
+#ifndef	_FENV_H_
+#define	_FENV_H_
+
+#include <sys/cdefs.h>
+#include <sys/_types.h>
+
+/*                   
+ * To preserve binary compatibility with FreeBSD 5.3, we pack the
+ * mxcsr into some reserved fields, rather than changing sizeof(fenv_t).
+ */
+typedef struct {
+	__uint16_t	__control;
+	__uint16_t      __mxcsr_hi;
+	__uint16_t	__status;
+	__uint16_t      __mxcsr_lo;
+	__uint32_t	__tag;
+	char		__other[16];
+} fenv_t;
+
+#define	__get_mxcsr(env)	(((env).__mxcsr_hi << 16) |	\
+				 ((env).__mxcsr_lo))
+#define	__set_mxcsr(env, x)	do {				\
+	(env).__mxcsr_hi = (__uint32_t)(x) >> 16;		\
+	(env).__mxcsr_lo = (__uint16_t)(x);			\
+} while (0)
+
+typedef	__uint16_t	fexcept_t;
+
+/* Exception flags */
+#define	FE_INVALID	0x01
+#define	FE_DENORMAL	0x02
+#define	FE_DIVBYZERO	0x04
+#define	FE_OVERFLOW	0x08
+#define	FE_UNDERFLOW	0x10
+#define	FE_INEXACT	0x20
+#define	FE_ALL_EXCEPT	(FE_DIVBYZERO | FE_DENORMAL | FE_INEXACT | \
+			 FE_INVALID | FE_OVERFLOW | FE_UNDERFLOW)
+
+/* Rounding modes */
+#define	FE_TONEAREST	0x0000
+#define	FE_DOWNWARD	0x0400
+#define	FE_UPWARD	0x0800
+#define	FE_TOWARDZERO	0x0c00
+#define	_ROUND_MASK	(FE_TONEAREST | FE_DOWNWARD | \
+			 FE_UPWARD | FE_TOWARDZERO)
+
+/*
+ * As compared to the x87 control word, the SSE unit's control word
+ * has the rounding control bits offset by 3 and the exception mask
+ * bits offset by 7.
+ */
+#define	_SSE_ROUND_SHIFT	3
+#define	_SSE_EMASK_SHIFT	7
+
+/* After testing for SSE support once, we cache the result in __has_sse. */
+enum __sse_support { __SSE_YES, __SSE_NO, __SSE_UNK };
+extern enum __sse_support __has_sse;
+int __test_sse(void);
+#ifdef __SSE__
+#define	__HAS_SSE()	1
+#else
+#define	__HAS_SSE()	(__has_sse == __SSE_YES ||			\
+			 (__has_sse == __SSE_UNK && __test_sse()))
+#endif
+
+__BEGIN_DECLS
+
+/* Default floating-point environment */
+extern const fenv_t	__fe_dfl_env;
+#define	FE_DFL_ENV	(&__fe_dfl_env)
+
+#define	__fldcw(__cw)		__asm __volatile("fldcw %0" : : "m" (__cw))
+#define	__fldenv(__env)		__asm __volatile("fldenv %0" : : "m" (__env))
+#define	__fnclex()		__asm __volatile("fnclex")
+#define	__fnstenv(__env)	__asm __volatile("fnstenv %0" : "=m" (*(__env)))
+#define	__fnstcw(__cw)		__asm __volatile("fnstcw %0" : "=m" (*(__cw)))
+#define	__fnstsw(__sw)		__asm __volatile("fnstsw %0" : "=am" (*(__sw)))
+#define	__fwait()		__asm __volatile("fwait")
+#define	__ldmxcsr(__csr)	__asm __volatile("ldmxcsr %0" : : "m" (__csr))
+#define	__stmxcsr(__csr)	__asm __volatile("stmxcsr %0" : "=m" (*(__csr)))
+
+static __inline int
+feclearexcept(int __excepts)
+{
+	fenv_t __env;
+	int __mxcsr;
+
+	if (__excepts == FE_ALL_EXCEPT) {
+		__fnclex();
+	} else {
+		__fnstenv(&__env);
+		__env.__status &= ~__excepts;
+		__fldenv(__env);
+	}
+	if (__HAS_SSE()) {
+		__stmxcsr(&__mxcsr);
+		__mxcsr &= ~__excepts;
+		__ldmxcsr(__mxcsr);
+	}
+	return (0);
+}
+
+static __inline int
+fegetexceptflag(fexcept_t *__flagp, int __excepts)
+{
+	int __mxcsr, __status;
+
+	__fnstsw(&__status);
+	if (__HAS_SSE())
+		__stmxcsr(&__mxcsr);
+	else
+		__mxcsr = 0;
+	*__flagp = (__mxcsr | __status) & __excepts;
+	return (0);
+}
+
+int fesetexceptflag(const fexcept_t *__flagp, int __excepts);
+int feraiseexcept(int __excepts);
+
+static __inline int
+fetestexcept(int __excepts)
+{
+	int __mxcsr, __status;
+
+	__fnstsw(&__status);
+	if (__HAS_SSE())
+		__stmxcsr(&__mxcsr);
+	else
+		__mxcsr = 0;
+	return ((__status | __mxcsr) & __excepts);
+}
+
+static __inline int
+fegetround(void)
+{
+	int __control;
+
+	/*
+	 * We assume that the x87 and the SSE unit agree on the
+	 * rounding mode.  Reading the control word on the x87 turns
+	 * out to be about 5 times faster than reading it on the SSE
+	 * unit on an Opteron 244.
+	 */
+	__fnstcw(&__control);
+	return (__control & _ROUND_MASK);
+}
+
+static __inline int
+fesetround(int __round)
+{
+	int __mxcsr, __control;
+
+	if (__round & ~_ROUND_MASK)
+		return (-1);
+
+	__fnstcw(&__control);
+	__control &= ~_ROUND_MASK;
+	__control |= __round;
+	__fldcw(__control);
+
+	if (__HAS_SSE()) {
+		__stmxcsr(&__mxcsr);
+		__mxcsr &= ~(_ROUND_MASK << _SSE_ROUND_SHIFT);
+		__mxcsr |= __round << _SSE_ROUND_SHIFT;
+		__ldmxcsr(__mxcsr);
+	}
+
+	return (0);
+}
+
+int fegetenv(fenv_t *__envp);
+int feholdexcept(fenv_t *__envp);
+
+static __inline int
+fesetenv(const fenv_t *__envp)
+{
+	fenv_t __env = *__envp;
+	int __mxcsr;
+
+	__mxcsr = __get_mxcsr(__env);
+	__set_mxcsr(__env, 0xffffffff);
+	__fldenv(__env);
+	if (__HAS_SSE())
+		__ldmxcsr(__mxcsr);
+	return (0);
+}
+
+int feupdateenv(const fenv_t *__envp);
+
+#if __BSD_VISIBLE
+
+int feenableexcept(int __mask);
+int fedisableexcept(int __mask);
+
+static __inline int
+fegetexcept(void)
+{
+	int __control;
+
+	/*
+	 * We assume that the masks for the x87 and the SSE unit are
+	 * the same.
+	 */
+	__fnstcw(&__control);
+	return (~__control & FE_ALL_EXCEPT);
+}
+
+#endif /* __BSD_VISIBLE */
+
+__END_DECLS
+
+#endif	/* !_FENV_H_ */
Index: /trunk/src/emx/include/float.h
===================================================================
--- /trunk/src/emx/include/float.h	(revision 2012)
+++ /trunk/src/emx/include/float.h	(revision 2013)
@@ -1,5 +1,5 @@
 /* float.h,v 1.5 2004/09/14 22:27:33 bird Exp */
 /** @file
- * EMX + GCC
+ * EMX + GCC + Some FreeBSD Bits.
  *
  * @remark Must be used instead of the GCC one.
@@ -11,7 +11,11 @@
 #include <sys/cdefs.h>
 
+__BEGIN_DECLS
+extern int __flt_rounds(void);
+__END_DECLS
 
 #define FLT_RADIX	2		/* b */
-#define FLT_ROUNDS	1		/* FP addition rounds to nearest */
+#define FLT_ROUNDS	__flt_rounds()
+
 #define	FLT_EVAL_METHOD	(-1)		/* i387 semantics are...interesting */
 #define	DECIMAL_DIG	21		/* max precision in decimal digits */
@@ -67,20 +71,4 @@
 #if !defined (__STRICT_ANSI__) && !defined (_POSIX_SOURCE) || defined(__USE_EMX)
 
-#define M_E             2.71828182845904523536
-#define M_LOG2E         1.44269504088896340736
-#define M_LOG10E        0.434294481903251827651
-#define M_LN2           0.693147180559945309422
-#define M_LN10          2.302585092994045684018
-
-#define M_PI            3.14159265358979323846
-#define M_PI_2          1.57079632679489661923
-#define M_PI_4          0.785398163397448309616
-#define M_1_PI          0.318309886183790671538
-#define M_2_PI          0.636619772367581343076
-#define M_2_SQRTPI      1.12837916709551257390
-
-#define M_SQRT2         1.41421356237309504880
-#define M_SQRT1_2       0.707106781186547524401
-
 #if !defined (_OMIT_387_STUFF)
 
Index: /trunk/src/emx/include/machine/asm.h
===================================================================
--- /trunk/src/emx/include/machine/asm.h	(revision 2013)
+++ /trunk/src/emx/include/machine/asm.h	(revision 2013)
@@ -0,0 +1,1 @@
+#include <386/asm.h>
Index: /trunk/src/emx/include/machine/npx.h
===================================================================
--- /trunk/src/emx/include/machine/npx.h	(revision 2013)
+++ /trunk/src/emx/include/machine/npx.h	(revision 2013)
@@ -0,0 +1,1 @@
+#include <386/npx.h>
Index: /trunk/src/emx/include/math.h
===================================================================
--- /trunk/src/emx/include/math.h	(revision 2012)
+++ /trunk/src/emx/include/math.h	(revision 2013)
@@ -1,203 +1,477 @@
-/* math.h,v 1.5 2004/09/14 22:27:34 bird Exp */
+/*
+ * ====================================================
+ * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
+ *
+ * Developed at SunPro, a Sun Microsystems, Inc. business.
+ * Permission to use, copy, modify, and distribute this
+ * software is freely granted, provided that this notice
+ * is preserved.
+ * ====================================================
+ */
+
+/*
+ * from: @(#)fdlibm.h 5.1 93/09/24
+ * $FreeBSD: src/lib/msun/src/math.h,v 1.61 2005/04/16 21:12:47 das Exp $
+ */
+
 /** @file
- * EMX
- * @todo This needs C99 attention!
- */
-
-#ifndef _MATH_H
-#define _MATH_H
-
-#if defined (__cplusplus)
-extern "C" {
-#endif
-
-#define HUGE_VAL        1e9999
-
-double acos (double);
-double asin (double);
-double atan (double);
-double atan2 (double, double);
-double cbrt (double);
-double ceil (double);
-double cos (double);
-double cosh (double);
-double exp (double);
-double fabs (double);
-double floor (double);
-double fmod (double, double);
-double frexp (double, int *);
-double hypot (double, double);
-double ldexp (double, int);
-double log (double);
-double log10 (double);
-double modf (double, double *);
-double pow (double, double);
-double rint (double);
-double sin (double);
-double sinh (double);
-double sqrt (double);
-double tan (double);
-double tanh (double);
-double trunc (double);
-
-#if !defined (__NO_C9X)
-
-#define HUGE_VALF       1e9999F
-#define HUGE_VALL       1e9999L
-#define INFINITY        1e9999          /* TODO: float_t */
-#define NAN             (0.0/0.0)       /* TODO: Exceptions, float_t */
-
-#define DECIMAL_DIG     21
-
-#define FP_ZERO         0
-#define FP_SUBNORMAL    1
-#define FP_NORMAL       2
-#define FP_INFINITE     3
-#define FP_NAN          4
-
-#define fpclassify(x) \
-    ((sizeof (x) == sizeof (float))  ? __fpclassifyf (x) \
-   : (sizeof (x) == sizeof (double)) ? __fpclassify (x) \
-   :                                   __fpclassifyl (x))
-
-#define signbit(x) \
-    ((sizeof (x) == sizeof (float))  ? __signbitf (x) \
-   : (sizeof (x) == sizeof (double)) ? __signbit (x) \
-   :                                   __signbitl (x))
-
-#define isfinite(x) \
-    ((sizeof (x) == sizeof (float))  ? __isfinitef (x) \
-   : (sizeof (x) == sizeof (double)) ? __isfinite (x) \
-   :                                   __isfinitel (x))
-
-#define isnormal(x) \
-    ((sizeof (x) == sizeof (float))  ? __isnormalf (x) \
-   : (sizeof (x) == sizeof (double)) ? __isnormal (x) \
-   :                                   __isnormall (x))
-
-#define isnan(x) \
-    ((sizeof (x) == sizeof (float))  ? __isnanf (x) \
-   : (sizeof (x) == sizeof (double)) ? __isnan (x) \
-   :                                   __isnanl (x))
-
-float copysignf (float, float);
-double copysign (double, double);
-long double copysignl (long double, long double);
-
-float nextafterf (float, float);
-double nextafter (double, double);
-long double nextafterl (long double, long double);
-
-long double acosl (long double);
-long double asinl (long double);
-long double atanl (long double);
-long double atan2l (long double, long double);
-long double cbrtl (long double);
-long double ceill (long double);
-long double cosl (long double);
-long double coshl (long double);
-long double expl (long double);
-long double fabsl (long double);
-long double floorl (long double);
-long double fmodl (long double, long double);
-long double frexpl (long double, int *);
-long double hypotl (long double, long double);
-long double ldexpl (long double, int);
-long double logl (long double);
-long double log10l (long double);
-long double modfl (long double, long double *);
-long double powl (long double, long double);
-long double rintl (long double);
-long double sinl (long double);
-long double sinhl (long double);
-long double sqrtl (long double);
-long double tanl (long double);
-long double tanhl (long double);
-long double truncl (long double);
-
-int __fpclassify (double);
-int __fpclassifyf (float);
-int __fpclassifyl (long double);
-
-int __signbit (double);
-int __signbitf (float);
-int __signbitl (long double);
-
-int __isfinite (double);
-int __isfinitef (float);
-int __isfinitel (long double);
-
-int __isnormal (double);
-int __isnormalf (float);
-int __isnormall (long double);
-
-int __isnan (double);
-int __isnanf (float);
-int __isnanl (long double);
-
-#endif
-
-#if !defined (__STRICT_ANSI__)
-
-#define HUGE    HUGE_VAL
-
-double atof (__const__ char *);
-
-#endif
-
-
-#if !defined (__STRICT_ANSI__) || defined (_WITH_UNDERSCORE)
-
-#define _HUGE_VALF      1e9999F
-#define _HUGE_VALL      1e9999L
-#define _INFINITY       1e9999          /* TODO: float_t */
-#define _NAN            (0.0/0.0)       /* TODO: Exceptions, float_t */
-
-#define _LHUGE_VAL      _HUGE_VALL
-#define _LHUGE          _HUGE_VALL
-
-
-float _copysignf (float, float);
-double _copysign (double, double);
-long double _copysignl (long double, long double);
-
-float _nextafterf (float, float);
-double _nextafter (double, double);
-long double _nextafterl (long double, long double);
-
-long double _acosl (long double);
-long double _asinl (long double);
-long double _atanl (long double);
-long double _atan2l (long double, long double);
-long double _atofl (__const__ char *);
-long double _cbrtl (long double);
-long double _ceill (long double);
-long double _cosl (long double);
-long double _coshl (long double);
-long double _hypotl (long double, long double);
-long double _expl (long double);
-long double _fabsl (long double);
-long double _floorl (long double);
-long double _fmodl (long double, long double);
-long double _frexpl (long double, int *);
-long double _ldexpl (long double, int);
-long double _logl (long double);
-long double _log10l (long double);
-long double _modfl (long double, long double *);
-long double _powl (long double, long double);
-long double _rintl (long double);
-long double _sinl (long double);
-long double _sinhl (long double);
-long double _sqrtl (long double);
-long double _tanl (long double);
-long double _tanhl (long double);
-long double _truncl (long double);
-
-#endif
-
-
-#if defined (__cplusplus)
-}
-#endif
-
-#endif /* not _MATH_H */
+ * FreeBSD HEAD 2005-06-xx
+ */
+
+#ifndef _MATH_H_
+#define	_MATH_H_
+
+#include <sys/cdefs.h>
+#include <sys/_types.h>
+#include <machine/_limits.h>
+
+/*
+ * ANSI/POSIX
+ */
+extern const union __infinity_un {
+	unsigned char	__uc[8];
+	double		__ud;
+} __infinity;
+
+extern const union __nan_un {
+	unsigned char	__uc[sizeof(float)];
+	float		__uf;
+} __nan;
+
+#if __GNUC_PREREQ__(3, 3) || (defined(__INTEL_COMPILER) && __INTEL_COMPILER >= 800)
+#define	__MATH_BUILTIN_CONSTANTS
+#endif
+
+#if __GNUC_PREREQ__(3, 0) && !defined(__INTEL_COMPILER)
+#define	__MATH_BUILTIN_RELOPS
+#endif
+
+#ifdef __MATH_BUILTIN_CONSTANTS
+#define	HUGE_VAL	__builtin_huge_val()
+#else
+#define	HUGE_VAL	(__infinity.__ud)
+#endif
+
+#if __ISO_C_VISIBLE >= 1999
+#define	FP_ILOGB0	(-__INT_MAX)
+#define	FP_ILOGBNAN	__INT_MAX
+
+#ifdef __MATH_BUILTIN_CONSTANTS
+#define	HUGE_VALF	__builtin_huge_valf()
+#define	HUGE_VALL	__builtin_huge_vall()
+#define	INFINITY	__builtin_inf()
+#define	NAN		__builtin_nan("")
+#else
+#define	HUGE_VALF	(float)HUGE_VAL
+#define	HUGE_VALL	(long double)HUGE_VAL
+#define	INFINITY	HUGE_VALF
+#define	NAN		(__nan.__uf)
+#endif /* __MATH_BUILTIN_CONSTANTS */
+
+#define	MATH_ERRNO	1
+#define	MATH_ERREXCEPT	2
+#define	math_errhandling	MATH_ERREXCEPT
+
+/* XXX We need a <machine/math.h>. */
+#if defined(__ia64__) || defined(__sparc64__)
+#define	FP_FAST_FMA
+#endif
+#ifdef __ia64__
+#define	FP_FAST_FMAL
+#endif
+#define	FP_FAST_FMAF
+
+/* Symbolic constants to classify floating point numbers. */
+#define	FP_INFINITE	0x01
+#define	FP_NAN		0x02
+#define	FP_NORMAL	0x04
+#define	FP_SUBNORMAL	0x08
+#define	FP_ZERO		0x10
+#define	fpclassify(x) \
+    ((sizeof (x) == sizeof (float)) ? __fpclassifyf(x) \
+    : (sizeof (x) == sizeof (double)) ? __fpclassifyd(x) \
+    : __fpclassifyl(x))
+
+#define	isfinite(x)					\
+    ((sizeof (x) == sizeof (float)) ? __isfinitef(x)	\
+    : (sizeof (x) == sizeof (double)) ? __isfinite(x)	\
+    : __isfinitel(x))
+#define	isinf(x)					\
+    ((sizeof (x) == sizeof (float)) ? __isinff(x)	\
+    : (sizeof (x) == sizeof (double)) ? isinf(x)	\
+    : __isinfl(x))
+#define	isnan(x)					\
+    ((sizeof (x) == sizeof (float)) ? isnanf(x)		\
+    : (sizeof (x) == sizeof (double)) ? isnan(x)	\
+    : __isnanl(x))
+#define	isnormal(x)					\
+    ((sizeof (x) == sizeof (float)) ? __isnormalf(x)	\
+    : (sizeof (x) == sizeof (double)) ? __isnormal(x)	\
+    : __isnormall(x))
+
+#ifdef __MATH_BUILTIN_RELOPS
+#define	isgreater(x, y)		__builtin_isgreater((x), (y))
+#define	isgreaterequal(x, y)	__builtin_isgreaterequal((x), (y))
+#define	isless(x, y)		__builtin_isless((x), (y))
+#define	islessequal(x, y)	__builtin_islessequal((x), (y))
+#define	islessgreater(x, y)	__builtin_islessgreater((x), (y))
+#define	isunordered(x, y)	__builtin_isunordered((x), (y))
+#else
+#define	isgreater(x, y)		(!isunordered((x), (y)) && (x) > (y))
+#define	isgreaterequal(x, y)	(!isunordered((x), (y)) && (x) >= (y))
+#define	isless(x, y)		(!isunordered((x), (y)) && (x) < (y))
+#define	islessequal(x, y)	(!isunordered((x), (y)) && (x) <= (y))
+#define	islessgreater(x, y)	(!isunordered((x), (y)) && \
+					((x) > (y) || (y) > (x)))
+#define	isunordered(x, y)	(isnan(x) || isnan(y))
+#endif /* __MATH_BUILTIN_RELOPS */
+
+#define	signbit(x)					\
+    ((sizeof (x) == sizeof (float)) ? __signbitf(x)	\
+    : (sizeof (x) == sizeof (double)) ? __signbit(x)	\
+    : __signbitl(x))
+
+typedef	__double_t	double_t;
+typedef	__float_t	float_t;
+#endif /* __ISO_C_VISIBLE >= 1999 */
+
+/*
+ * XOPEN/SVID
+ */
+#if __BSD_VISIBLE || __XSI_VISIBLE
+#define	M_E		2.7182818284590452354	/* e */
+#define	M_LOG2E		1.4426950408889634074	/* log 2e */
+#define	M_LOG10E	0.43429448190325182765	/* log 10e */
+#define	M_LN2		0.69314718055994530942	/* log e2 */
+#define	M_LN10		2.30258509299404568402	/* log e10 */
+#define	M_PI		3.14159265358979323846	/* pi */
+#define	M_PI_2		1.57079632679489661923	/* pi/2 */
+#define	M_PI_4		0.78539816339744830962	/* pi/4 */
+#define	M_1_PI		0.31830988618379067154	/* 1/pi */
+#define	M_2_PI		0.63661977236758134308	/* 2/pi */
+#define	M_2_SQRTPI	1.12837916709551257390	/* 2/sqrt(pi) */
+#define	M_SQRT2		1.41421356237309504880	/* sqrt(2) */
+#define	M_SQRT1_2	0.70710678118654752440	/* 1/sqrt(2) */
+
+#define	MAXFLOAT	((float)3.40282346638528860e+38)
+extern int signgam;
+#endif /* __BSD_VISIBLE || __XSI_VISIBLE */
+
+#if __BSD_VISIBLE
+#if 0
+/* Old value from 4.4BSD-Lite math.h; this is probably better. */
+#define	HUGE		HUGE_VAL
+#else
+#define	HUGE		MAXFLOAT
+#endif
+#endif /* __BSD_VISIBLE */
+
+/*
+ * Most of these functions depend on the rounding mode and have the side
+ * effect of raising floating-point exceptions, so they are not declared
+ * as __pure2.  In C99, FENV_ACCESS affects the purity of these functions.
+ */
+__BEGIN_DECLS
+/*
+ * ANSI/POSIX
+ */
+int	__fpclassifyd(double) __pure2;
+int	__fpclassifyf(float) __pure2;
+int	__fpclassifyl(long double) __pure2;
+int	__isfinitef(float) __pure2;
+int	__isfinite(double) __pure2;
+int	__isfinitel(long double) __pure2;
+int	__isinff(float) __pure2;
+int	__isinfl(long double) __pure2;
+int	__isnanl(long double) __pure2;
+int	__isnormalf(float) __pure2;
+int	__isnormal(double) __pure2;
+int	__isnormall(long double) __pure2;
+int	__signbit(double) __pure2;
+int	__signbitf(float) __pure2;
+int	__signbitl(long double) __pure2;
+
+double	acos(double);
+double	asin(double);
+double	atan(double);
+double	atan2(double, double);
+double	cos(double);
+double	sin(double);
+double	tan(double);
+
+double	cosh(double);
+double	sinh(double);
+double	tanh(double);
+
+double	exp(double);
+double	frexp(double, int *);	/* fundamentally !__pure2 */
+double	ldexp(double, int);
+double	log(double);
+double	log10(double);
+double	modf(double, double *);	/* fundamentally !__pure2 */
+
+double	pow(double, double);
+double	sqrt(double);
+
+double	ceil(double);
+double	fabs(double) __pure2;
+double	floor(double);
+double	fmod(double, double);
+
+/*
+ * These functions are not in C90.
+ */
+#if __BSD_VISIBLE || __ISO_C_VISIBLE >= 1999 || __XSI_VISIBLE
+double	acosh(double);
+double	asinh(double);
+double	atanh(double);
+double	cbrt(double);
+double	erf(double);
+double	erfc(double);
+double	exp2(double);
+double	expm1(double);
+double	fma(double, double, double);
+double	hypot(double, double);
+int	ilogb(double) __pure2;
+int	(isinf)(double) __pure2;
+int	(isnan)(double) __pure2;
+double	lgamma(double);
+long long llrint(double);
+long long llround(double);
+double	log1p(double);
+double	logb(double);
+long	lrint(double);
+long	lround(double);
+double	nextafter(double, double);
+double	remainder(double, double);
+double	remquo(double, double, int *);
+double	rint(double);
+#endif /* __BSD_VISIBLE || __ISO_C_VISIBLE >= 1999 || __XSI_VISIBLE */
+
+#if __BSD_VISIBLE || __XSI_VISIBLE
+double	j0(double);
+double	j1(double);
+double	jn(int, double);
+double	scalb(double, double);
+double	y0(double);
+double	y1(double);
+double	yn(int, double);
+
+#if __XSI_VISIBLE <= 500 || __BSD_VISIBLE
+double	gamma(double);
+#endif
+#endif /* __BSD_VISIBLE || __XSI_VISIBLE */
+
+#if __BSD_VISIBLE || __ISO_C_VISIBLE >= 1999
+double	copysign(double, double) __pure2;
+double	fdim(double, double);
+double	fmax(double, double) __pure2;
+double	fmin(double, double) __pure2;
+double	nearbyint(double);
+double	round(double);
+double	scalbln(double, long);
+double	scalbn(double, int);
+double	tgamma(double);
+double	trunc(double);
+#endif
+
+/*
+ * BSD math library entry points
+ */
+#if __BSD_VISIBLE
+double	drem(double, double);
+int	finite(double) __pure2;
+int	isnanf(float) __pure2;
+
+/*
+ * Reentrant version of gamma & lgamma; passes signgam back by reference
+ * as the second argument; user must allocate space for signgam.
+ */
+double	gamma_r(double, int *);
+double	lgamma_r(double, int *);
+
+/*
+ * IEEE Test Vector
+ */
+double	significand(double);
+#endif /* __BSD_VISIBLE */
+
+/* float versions of ANSI/POSIX functions */
+#if __ISO_C_VISIBLE >= 1999
+float	acosf(float);
+float	asinf(float);
+float	atanf(float);
+float	atan2f(float, float);
+float	cosf(float);
+float	sinf(float);
+float	tanf(float);
+
+float	coshf(float);
+float	sinhf(float);
+float	tanhf(float);
+
+float	exp2f(float);
+float	expf(float);
+float	expm1f(float);
+float	frexpf(float, int *);	/* fundamentally !__pure2 */
+int	ilogbf(float) __pure2;
+float	ldexpf(float, int);
+float	log10f(float);
+float	log1pf(float);
+float	logf(float);
+float	modff(float, float *);	/* fundamentally !__pure2 */
+
+float	powf(float, float);
+float	sqrtf(float);
+
+float	ceilf(float);
+float	fabsf(float) __pure2;
+float	floorf(float);
+float	fmodf(float, float);
+float	roundf(float);
+
+float	erff(float);
+float	erfcf(float);
+float	hypotf(float, float);
+float	lgammaf(float);
+
+float	acoshf(float);
+float	asinhf(float);
+float	atanhf(float);
+float	cbrtf(float);
+float	logbf(float);
+float	copysignf(float, float) __pure2;
+long long llrintf(float);
+long long llroundf(float);
+long	lrintf(float);
+long	lroundf(float);
+float	nearbyintf(float);
+float	nextafterf(float, float);
+float	remainderf(float, float);
+float	remquof(float, float, int *);
+float	rintf(float);
+float	scalblnf(float, long);
+float	scalbnf(float, int);
+float	truncf(float);
+
+float	fdimf(float, float);
+float	fmaf(float, float, float);
+float	fmaxf(float, float) __pure2;
+float	fminf(float, float) __pure2;
+#endif
+
+/*
+ * float versions of BSD math library entry points
+ */
+#if __BSD_VISIBLE
+float	dremf(float, float);
+int	finitef(float) __pure2;
+float	gammaf(float);
+float	j0f(float);
+float	j1f(float);
+float	jnf(int, float);
+float	scalbf(float, float);
+float	y0f(float);
+float	y1f(float);
+float	ynf(int, float);
+
+/*
+ * Float versions of reentrant version of gamma & lgamma; passes
+ * signgam back by reference as the second argument; user must
+ * allocate space for signgam.
+ */
+float	gammaf_r(float, int *);
+float	lgammaf_r(float, int *);
+
+/*
+ * float version of IEEE Test Vector
+ */
+float	significandf(float);
+#endif	/* __BSD_VISIBLE */
+
+/*
+ * long double versions of ISO/POSIX math functions
+ */
+#if __ISO_C_VISIBLE >= 1999
+#if 1  /* bird: we've got these */
+long double	acoshl(long double);
+long double	acosl(long double);
+long double	asinhl(long double);
+long double	asinl(long double);
+long double	atan2l(long double, long double);
+long double	atanhl(long double);
+long double	atanl(long double);
+long double	cbrtl(long double);
+#endif
+long double	ceill(long double);
+long double	copysignl(long double, long double) __pure2;
+#if 1 /* bird: we've got some of these */
+long double	coshl(long double);
+long double	cosl(long double);
+/*long double	erfcl(long double);
+long double	erfl(long double);
+long double	exp2l(long double);*/
+long double	expl(long double);
+/*long double	expm1l(long double);*/
+#endif
+long double	fabsl(long double) __pure2;
+long double	fdiml(long double, long double);
+long double	floorl(long double);
+long double	fmal(long double, long double, long double);
+long double	fmaxl(long double, long double) __pure2;
+long double	fminl(long double, long double) __pure2;
+#if 1 /* bird: we've got this one */
+long double	fmodl(long double, long double);
+#endif
+long double	frexpl(long double value, int *); /* fundamentally !__pure2 */
+#if 1 /* bird: we've got this one */
+long double	hypotl(long double, long double);
+#endif
+int		ilogbl(long double) __pure2;
+long double	ldexpl(long double, int);
+#if 0
+long double	lgammal(long double);
+long long	llrintl(long double);
+#endif
+long long	llroundl(long double);
+#if 1 /* bird: we've got some of these */
+long double	log10l(long double);
+/*long double	log1pl(long double);
+long double	log2l(long double);
+long double	logbl(long double);*/
+long double	logl(long double);
+/*long		lrintl(long double);*/
+#endif
+long		lroundl(long double);
+#if 1 /* bird: we've got some of these. */
+long double	modfl(long double, long double *); /* fundamentally !__pure2 */
+/*long double	nanl(const char *) __pure2;
+long double	nearbyintl(long double);*/
+#endif
+long double	nextafterl(long double, long double);
+double		nexttoward(double, long double);
+float		nexttowardf(float, long double);
+long double	nexttowardl(long double, long double);
+#if 1 /* bird: we've got powl. */
+long double	powl(long double, long double);
+/*long double	remainderl(long double, long double);
+long double	remquol(long double, long double, int *);*/
+long double	rintl(long double);
+#endif
+long double	roundl(long double);
+long double	scalblnl(long double, long);
+long double	scalbnl(long double, int);
+#if 1 /* birdL we 've got most of these. */
+long double	sinhl(long double);
+long double	sinl(long double);
+long double	sqrtl(long double);
+long double	tanhl(long double);
+long double	tanl(long double);
+/*long double	tgammal(long double); */
+#endif
+long double	truncl(long double);
+
+#endif /* __ISO_C_VISIBLE >= 1999 */
+__END_DECLS
+
+#endif /* !_MATH_H_ */
Index: /trunk/src/emx/src/lib/bsd/gen/fpclassify.c
===================================================================
--- /trunk/src/emx/src/lib/bsd/gen/fpclassify.c	(revision 2012)
+++ /trunk/src/emx/src/lib/bsd/gen/fpclassify.c	(revision 2013)
@@ -28,4 +28,5 @@
  */
 
+#include "namespace.h"
 #include <sys/endian.h>
 
Index: /trunk/src/emx/src/lib/bsd/i386/gen/fabs.S
===================================================================
--- /trunk/src/emx/src/lib/bsd/i386/gen/fabs.S	(revision 2012)
+++ /trunk/src/emx/src/lib/bsd/i386/gen/fabs.S	(revision 2013)
@@ -41,5 +41,5 @@
 __FBSDID("$FreeBSD: src/lib/libc/i386/gen/fabs.S,v 1.9 2002/03/23 02:44:18 obrien Exp $");
 
-ENTRY(fabs)
+STDENTRY(fabs)
 	fldl	4(%esp)
 	fabs
Index: /trunk/src/emx/src/lib/bsd/i386/gen/flt_rounds.c
===================================================================
--- /trunk/src/emx/src/lib/bsd/i386/gen/flt_rounds.c	(revision 2012)
+++ /trunk/src/emx/src/lib/bsd/i386/gen/flt_rounds.c	(revision 2013)
@@ -4,4 +4,5 @@
  */
 
+#include "namespace.h"
 #include <sys/cdefs.h>
 __FBSDID("$FreeBSD: src/lib/libc/i386/gen/flt_rounds.c,v 1.1 2004/07/19 08:17:24 das Exp $");
Index: /trunk/src/emx/src/lib/bsd/i386/gen/infinity.c
===================================================================
--- /trunk/src/emx/src/lib/bsd/i386/gen/infinity.c	(revision 2012)
+++ /trunk/src/emx/src/lib/bsd/i386/gen/infinity.c	(revision 2013)
@@ -3,4 +3,5 @@
  */
 
+#include "namespace.h"
 #include <sys/cdefs.h>
 __FBSDID("$FreeBSD: src/lib/libc/i386/gen/infinity.c,v 1.10 2003/02/08 20:37:52 mike Exp $");
Index: /trunk/src/emx/src/lib/bsd/i386/gen/ldexp.c
===================================================================
--- /trunk/src/emx/src/lib/bsd/i386/gen/ldexp.c	(revision 2012)
+++ /trunk/src/emx/src/lib/bsd/i386/gen/ldexp.c	(revision 2013)
@@ -35,4 +35,5 @@
  */
 
+#include "namespace.h"
 #if defined(LIBC_SCCS) && !defined(lint)
 static char sccsid[] = "@(#)ldexp.c	8.1 (Berkeley) 6/4/93";
@@ -52,5 +53,5 @@
  */
 double
-ldexp (double value, int exp)
+_STD(ldexp) (double value, int exp)
 {
 	double temp, texp, temp2;
Index: /trunk/src/emx/src/lib/bsd/i386/gen/modf.S
===================================================================
--- /trunk/src/emx/src/lib/bsd/i386/gen/modf.S	(revision 2012)
+++ /trunk/src/emx/src/lib/bsd/i386/gen/modf.S	(revision 2013)
@@ -51,5 +51,5 @@
 /* With CHOP mode on, frndint behaves as TRUNC does.  Useful. */
 
-ENTRY(modf)
+STDENTRY(modf)
 	pushl %ebp
 	movl %esp,%ebp
Index: /trunk/src/emx/src/lib/bsd/i386/i386.smak
===================================================================
--- /trunk/src/emx/src/lib/bsd/i386/i386.smak	(revision 2013)
+++ /trunk/src/emx/src/lib/bsd/i386/i386.smak	(revision 2013)
@@ -0,0 +1,1 @@
+libc.TSRC.386 += $(wildcard src/lib/bsd/i386/*/*.c) $(wildcard src/lib/bsd/i386/*/*.S)
Index: /trunk/src/emx/src/lib/bsd/include/_fpmath.h
===================================================================
--- /trunk/src/emx/src/lib/bsd/include/_fpmath.h	(revision 2012)
+++ /trunk/src/emx/src/lib/bsd/include/_fpmath.h	(revision 2013)
@@ -1,4 +1,4 @@
 /*-
- * Copyright (c) 2002, 2003 David Schultz <dschultz@uclink.Berkeley.EDU>
+ * Copyright (c) 2002, 2003 David Schultz <das@FreeBSD.ORG>
  * All rights reserved.
  *
@@ -24,5 +24,5 @@
  * SUCH DAMAGE.
  *
- * $FreeBSD: src/lib/libc/i386/_fpmath.h,v 1.3 2004/01/18 07:57:01 das Exp $
+ * $FreeBSD: src/lib/libc/i386/_fpmath.h,v 1.5 2005/03/07 04:55:22 das Exp $
  */
 
@@ -38,5 +38,6 @@
 };
 
-#define	mask_nbit_l(u)	((u).bits.manh &= 0x7fffffff)
+#define	LDBL_NBIT	0x80000000
+#define	mask_nbit_l(u)	((u).bits.manh &= ~LDBL_NBIT)
 
 #define	LDBL_MANH_SIZE	32
Index: unk/src/emx/src/lib/conv/386/smallcnv.s
===================================================================
--- /trunk/src/emx/src/lib/conv/386/smallcnv.s	(revision 2012)
+++ 	(revision )
@@ -1,7 +1,0 @@
-/ smallcnv.s (emx+gcc) -- Copyright (c) 1996 by Eberhard Mattes
-
-        .stabs  "___atod", 11, 0, 0, 0
-        .stabs  "___small_atod", 1, 0, 0, 0
-
-        .stabs  "___dtoa", 11, 0, 0, 0
-        .stabs  "___small_dtoa", 1, 0, 0, 0
Index: /trunk/src/emx/src/lib/conv/conv.smak
===================================================================
--- /trunk/src/emx/src/lib/conv/conv.smak	(revision 2012)
+++ /trunk/src/emx/src/lib/conv/conv.smak	(revision 2013)
@@ -46,5 +46,5 @@
 .MODULE :=
 .TKIND	:= aout
-.TSRC.386 := src/lib/conv/386/smallcnv.s
+#.TSRC.386 := src/lib/conv/386/smallcnv.s - smallcnv is dead.
 .TSRC	:= $(.TSRC.$(CPU))
 .INSDIR	:=
Index: unk/src/emx/src/lib/conv/smalatod.c
===================================================================
--- /trunk/src/emx/src/lib/conv/smalatod.c	(revision 2012)
+++ 	(revision )
@@ -1,227 +1,0 @@
-/* smalatod.c (emx+gcc) -- Copyright (c) 1990-1996 by Eberhard Mattes */
-
-#include "libc-alias.h"
-#include <stdlib.h>
-#include <math.h>
-#include <float.h>
-#include <locale.h>
-#include <errno.h>
-#define __atod  __small_atod /* bird: let's get a prototype. */
-#include <emx/float.h>
-#include <InnoTekLIBC/locale.h>
-#include <ctype.h>
-
-
-const char * __small_atod (long double *p_result, const char *string,
-                           int min_exp, int max_exp, int bias,
-                           int mant_dig, int decimal_dig,
-                           int max_10_exp, int min_den_10_exp)
-{
-  const unsigned char *s, *end;
-  long double x, lim;
-  long scale;
-  char ok, dot, neg;
-  int fpclass;
-  union
-  {
-    long double ld;
-    struct
-    {
-      unsigned long long mant;
-      unsigned short exp;
-    } x;
-  } z;
-
-  s = end = (unsigned char *)string;
-
-  switch (mant_dig)
-    {
-    case 24:
-      lim = FLT_MAX / 10.0;
-      break;
-    case 53:
-      lim = DBL_MAX / 10.0;
-      break;
-    case 64:
-      lim = LDBL_MAX / 10.0L;
-      break;
-    default:
-      *p_result = NAN;
-      return (const char *)string;
-    }
-
-  while (isspace (*s))
-    ++s;
-  neg = 0;
-  if (*s == '+')
-    ++s;
-  else if (*s == '-')
-    {
-      neg = 1; ++s;
-    }
-
-  if ((s[0] == 'i' || s[0] == 'I')
-      && (s[1] == 'n' || s[1] == 'N')
-      && (s[2] == 'f' || s[2] == 'F'))
-    {
-      s += 3;
-      if ((s[0] == 'i' || s[0] == 'I')
-          && (s[1] == 'n' || s[1] == 'N')
-          && (s[2] == 'i' || s[2] == 'I')
-          && (s[3] == 't' || s[3] == 'T')
-          && (s[4] == 'y' || s[4] == 'Y'))
-        s += 5;
-      *p_result = neg ? -INFINITY : INFINITY;
-      return (const char *)s;
-    }
-  else if ((s[0] == 'n' || s[0] == 'N')
-           && (s[1] == 'a' || s[1] == 'A')
-           && (s[2] == 'n' || s[2] == 'N'))
-    {
-      s += 3;
-      end = s;
-      if (end[0] == '(')
-        {
-          ++end;
-          while (*end == '_' || (*end >= '0' && *end <= '9')
-                 || (*end >= 'a' && *end <= 'z')
-                 || (*end >= 'A' && *end <= 'Z'))
-            ++end;
-          if (*end == ')')
-            ++end;
-          else
-            end = s;
-        }
-      *p_result = NAN;
-      return (const char *)end;
-    }
-
-  /* Scan the entire mantissa, updating SCALE. */
-
-  ok = dot = 0;
-  x = 0.0; scale = 0;
-  for (;;)
-    {
-      if (*s == (unsigned char)__libc_gLocaleLconv.s.decimal_point[0] && !dot)
-        {
-          dot = 1;
-          ++s;
-        }
-      else if (isdigit (*s))
-        {
-          ok = 1;
-          if (x > lim)
-            ++scale;            /* Overflow */
-          else
-            x = x * 10.0L + (*s - '0');
-          ++s;
-          if (dot)
-            --scale;
-        }
-      else
-        break;
-    }
-
-  if (!ok)
-    {
-      *p_result = 0.0;
-      return (const char *)string;
-    }
-
-  /* Scan the exponent. */
-
-  end = s;
-  if (*s == 'e' || *s == 'E')
-    {
-      char eneg;
-
-      ++s; eneg = 0;
-      if (*s == '+')
-        ++s;
-      else if (*s == '-')
-        {
-          eneg = 1; ++s;
-        }
-      if (isdigit (*s))         /* Note: strtol() skips whitespace! */
-        {
-          int saved_errno;
-          long exp;
-          const unsigned char *q;
-
-          saved_errno = errno; errno = 0;
-          exp = strtol ((const char *)s, (char **)(void *)&q, 10);
-          if (errno == ERANGE)
-            {
-              /* If the exponent overflows, we need a *very* long
-                 mantissa to get a valid number.  Such a mantissa
-                 doesn't fit into memory. */
-
-              end = q;
-              if (eneg)
-                goto underflow;
-              else
-                goto overflow;
-            }
-          else if (errno == 0 && q != s)
-            {
-              end = q;
-
-              /* TODO: check for integer overflow */
-
-              if (eneg)
-                scale -= exp;
-              else
-                scale += exp;
-            }
-          errno = saved_errno;
-        }
-    }
-
-  /* Apply the scaling factor and the exponent.  Don't do this if X is
-     0.0 to be able to detect underflow. */
-
-  if (x == 0.0 || scale == 0)
-    z.ld = x;
-  else
-    z.ld = x * _powl (10.0L, (long double)scale);
-
-  switch (mant_dig)
-    {
-    case 24:
-      fpclass = __fpclassifyf ((float)z.ld);
-      break;
-    case 53:
-      fpclass = __fpclassify ((double)z.ld);
-      break;
-    case 64:
-      fpclass = __fpclassifyl (z.ld);
-      break;
-    default:
-      fpclass = FP_NAN;
-      break;
-    }
-
-  switch (fpclass)
-    {
-    case FP_ZERO:
-      if (x != 0.0)
-        errno = ERANGE;
-      break;
-    case FP_INFINITE:
-      errno = ERANGE;
-      break;
-    }
-
-  *p_result = neg ? -z.ld : z.ld;
-  return (const char *)end;
-
-overflow:
-  *p_result = neg ? -HUGE_VALL : HUGE_VALL;
-  errno = ERANGE;
-  return (const char *)end;
-
-underflow:
-  *p_result = neg ? -0.0 : 0.0;
-  errno = ERANGE;
-  return (const char *)end;
-}
Index: unk/src/emx/src/lib/conv/smaldtoa.c
===================================================================
--- /trunk/src/emx/src/lib/conv/smaldtoa.c	(revision 2012)
+++ 	(revision )
@@ -1,114 +1,0 @@
-/* smalldtoa.c (emx+gcc) -- Copyright (c) 1994-1996 by Eberhard Mattes */
-
-#include "libc-alias.h"
-#include <stdlib.h>
-#include <string.h>
-#include <math.h>
-#include <float.h>
-#define __dtoa  __small_dtoa /* bird: let's get a prototype. */
-#include <emx/float.h>
-
-
-/* Convert a positive floating point number to a string of decimal
-   digits and an exponent.  The position of the decimal point is
-   assumed to be after the first digit.  Examples:
-
-   1234567890123456789 --> 1234567890123456789, xp = 18
-   123                 --> 1230000000000000000, xp =  2
-   12                  --> 1200000000000000000, xp =  1
-   1.2                 --> 1200000000000000000, xp =  0
-   0.12                --> 1200000000000000000, xp = -1
-
-   This function assumes that
-
-     LDBL_MANT_DIG <= 8 * sizeof (unsigned long long)
-
-   which is true for the extended real format (64-bit mantissa).
-   There are algorithms which are more precise, but they're more
-   complicated. */
-
-char *__small_dtoa (char *buffer, int *p_exp,
-                    long double x, int ndigits, int fmt, int dig)
-{
-  long double scaled;
-  unsigned long long intg;
-  int xp, len, a, i, next;
-
-  xp = (int)_floorl (_log10l (x));
-
-  if (LDBL_DIG - 1 - xp < LDBL_MAX_10_EXP)
-    {
-      /* This is the normal case. */
-
-      scaled = x * _powl (10.0, LDBL_DIG - 1 - xp);
-    }
-  else
-    {
-      /* Avoid overflow in _powl() for X near LDBL_MIN. */
-
-      scaled = (_powl (10.0, LDBL_DIG - 1 - xp - (LDBL_MAX_10_EXP - 1))
-                * (x * _powl (10.0, LDBL_MAX_10_EXP - 1)));
-    }
-
-  intg = (unsigned long long)scaled;
-  buffer[0] = '0';
-  _ulltoa (intg, buffer + 1, 10);
-
-  switch (fmt)
-    {
-    case DTOA_PRINTF_E:
-      a = ndigits + 1;
-      break;
-    case DTOA_PRINTF_F:
-      if (xp >= 0)
-        a = xp + 1 + ndigits;
-      else
-        a = ndigits;
-      break;
-    case DTOA_PRINTF_G:
-    case DTOA_GCVT:
-      a = ndigits;
-      break;
-    default:
-      abort ();
-    }
-  if (a > dig)
-    a = dig;
-
-  len = strlen (buffer + 1);
-  xp += len - LDBL_DIG;
-  if (len > a)
-    {
-      /* Round. */
-
-      next = 0;
-      if (buffer[a+1] > '5')
-        next = 1;
-      else if (buffer[a+1] == '5')
-        {
-          i = a + 2;
-          while (buffer[i] == '0')
-            ++i;
-          next = (buffer[i] != 0 || (buffer[a] & 1));
-        }
-      if (next)
-        while (++buffer[a] > '9')
-          --a;
-    }
-
-  i = a + 1;
-  while (i < DECIMAL_DIG + 1)
-    buffer[i++] = '0';
-  buffer[i] = 0;
-
-  if (buffer[0] != '0')
-    {
-      *p_exp = xp + 1;
-      return buffer;
-    }
-  else
-    {
-      *p_exp = xp;
-      return buffer + 1;
-    }
-}
Index: /trunk/src/emx/src/lib/libc.def
===================================================================
--- /trunk/src/emx/src/lib/libc.def	(revision 2012)
+++ /trunk/src/emx/src/lib/libc.def	(revision 2013)
@@ -237,7 +237,7 @@
     "___fmutex_request_internal_must_complete" @250
     "___fork" @251
-    "___fpclassify" @252
-    "___fpclassifyf" @253
-    "___fpclassifyl" @254
+    ;msun - dead "___fpclassify" @252
+    ;msun - dead "___fpclassifyf" @253
+    ;msun - dead "___fpclassifyl" @254
     "___free_ovflpage" @255
     "___ftime" @256
@@ -258,6 +258,6 @@
     "___isfinitef" @271
     "___isfinitel" @272
-    "___isnan" @273
-    "___isnanf" @274
+    "_isnan" @273   ;msun changed name
+    "_isnanf" @274  ;msun changed name
     "___isnanl" @275
     "___isnormal" @276
@@ -407,6 +407,6 @@
     "___signbitl" @420
     "___sigpause_bsd" @421
-    "___small_atod" @422
-    "___small_dtoa" @423
+    ; dead "___small_atod" @422
+    ; dead "___small_dtoa" @423
     "___smutex_request_internal" @424
     "___spawni" @425
@@ -876,5 +876,5 @@
     "__std_posix_memalign" @889
     "__std_pow" @890
-    "__std_powl" @891
+    ;msun - "__std_powl" @891
     "__std_printf" @892
     "__std_psignal" @893
@@ -1554,2 +1554,157 @@
     "__std_shmdt" @1559
     "__std_shmget" @1560
+    ; sunm
+    "__std_acosf" @1561
+    "__std_acosh" @1562
+    "__std_acoshf" @1563
+    "__std_asinf" @1564
+    "__std_asinh" @1565
+    "__std_asinhf" @1566
+    "__std_atan2f" @1567
+    "__std_atanf" @1568
+    "__std_atanh" @1569
+    "__std_atanhf" @1570
+    "__std_cabs" @1571
+    "__std_cabsf" @1572
+    "__std_cbrtf" @1573
+    "__std_ceilf" @1574
+    "__std_cimag" @1575
+    "__std_cimagf" @1576
+    "__std_cimagl" @1577
+    "__std_conj" @1578
+    "__std_conjf" @1579
+    "__std_conjl" @1580
+    "__std_cosf" @1581
+    "__std_coshf" @1582
+    "__std_creal" @1583
+    "__std_crealf" @1584
+    "__std_creall" @1585
+    "__std_erf" @1586
+    "__std_erfc" @1587
+    "__std_erfcf" @1588
+    "__std_erff" @1589
+    "__std_exp2" @1590
+    "__std_exp2f" @1591
+    "__std_expf" @1592
+    "__std_expm1" @1593
+    "__std_expm1f" @1594
+    "__std_truncf" @1595
+    "__std_fabsf" @1596
+    "__std_fdim" @1597
+    "__std_fdimf" @1598
+    "__std_fdiml" @1599
+    "__std_fegetenv" @1600
+    "__std_feholdexcept" @1601
+    "__std_feraiseexcept" @1602
+    "__std_fesetexceptflag" @1603
+    "__std_feupdateenv" @1604
+    "__std_floorf" @1605
+    "__std_fma" @1606
+    "__std_fmaf" @1607
+    "__std_fmal" @1608
+    "__std_fmax" @1609
+    "__std_fmaxf" @1610
+    "__std_fmaxl" @1611
+    "__std_fmin" @1612
+    "__std_fminf" @1613
+    "__std_fminl" @1614
+    "__std_fmodf" @1615
+    "__std_frexpf" @1616
+    "__std_hypotf" @1617
+    "__std_ilogb" @1618
+    "__std_ilogbf" @1619
+    "__std_ilogbl" @1620
+    "__std_ldexpf" @1621
+    "__std_lgamma" @1622
+    "__std_lgammaf" @1623
+    "__std_llrint" @1624
+    "__std_llrintf" @1625
+    "__std_llround" @1626
+    "__std_llroundf" @1627
+    "__std_llroundl" @1628
+    "__std_log10f" @1629
+    "__std_log1p" @1630
+    "__std_log1pf" @1631
+    "__std_logb" @1632
+    "__std_logbf" @1633
+    "__std_logf" @1634
+    "__std_lrint" @1635
+    "__std_lrintf" @1636
+    "__std_lround" @1637
+    "__std_lroundf" @1638
+    "__std_lroundl" @1639
+    "__std_modff" @1640
+    "__std_nearbyint" @1641
+    "__std_nearbyintf" @1642
+    "__std_nexttoward" @1643
+    "__std_nexttowardf" @1644
+    "__std_powf" @1645
+    "__std_remainder" @1646
+    "__std_remainderf" @1647
+    "__std_remquo" @1648
+    "__std_remquof" @1649
+    "__std_rintf" @1650
+    "__std_round" @1651
+    "__std_roundf" @1652
+    "__std_roundl" @1653
+    "__std_scalbf" @1654
+    "__std_scalbln" @1655
+    "__std_scalblnf" @1656
+    "__std_scalblnl" @1657
+    "__std_scalbn" @1658
+    "__std_scalbnf" @1659
+    "__std_scalbnl" @1660
+    "__std_signgam" @1661
+    "__std_sinf" @1662
+    "__std_sinhf" @1663
+    "__std_sqrtf" @1664
+    "__std_tanf" @1665
+    "__std_tanhf" @1666
+    "__std_tgamma" @1667
+    "___fpclassifyd" @1668
+    "___fpclassifyf" @1669
+    "___fpclassifyl" @1670
+    "_drem" @1671
+    "_dremf" @1672
+    "___exp__D" @1673
+    "___fe_dfl_env" @1674
+    "___fedisableexcept" @1675
+    "___feenableexcept" @1676
+    "_finite" @1677
+    "_finitef" @1678
+    "_gamma" @1679
+    "_gammaf" @1680
+    "_gammaf_r" @1681
+    "_gamma_r" @1682
+    "___has_sse" @1683
+    "___ieee754_rem_pio2" @1684
+    "___ieee754_rem_pio2f" @1685
+    "_j0" @1686
+    "_j0f" @1687
+    "_j1" @1688
+    "_j1f" @1689
+    "_jn" @1690
+    "_jnf" @1691
+    "___kernel_cos" @1692
+    "___kernel_cosf" @1693
+    "___kernel_rem_pio2" @1694
+    "___kernel_rem_pio2f" @1695
+    "___kernel_sin" @1696
+    "___kernel_sinf" @1697
+    "___kernel_tan" @1698
+    "___kernel_tanf" @1699
+    "_lgammaf_r" @1700
+    "_lgamma_r" @1701
+    "___log__D" @1702
+    "_scalb" @1703
+    "_significand" @1704
+    "_significandf" @1705
+    "__std_powl" @1706
+    "___test_sse" @1707
+    "_y0" @1708
+    "_y0f" @1709
+    "_y1" @1710
+    "_y1f" @1711
+    "_yn" @1712
+    "_ynf" @1713
+    "_z_abs" @1714
Index: /trunk/src/emx/src/lib/libc.smak
===================================================================
--- /trunk/src/emx/src/lib/libc.smak	(revision 2012)
+++ /trunk/src/emx/src/lib/libc.smak	(revision 2013)
@@ -398,5 +398,5 @@
 # files for _STD() macros invocations. Fortunately, gawk rulez forever thus
 # this happens in just a fraction of second.
-$.stmp-libc-std: $(wildcard src/lib/*/*.c src/lib/*/386/*.s src/lib/bsd/*/*.c src/lib/bsd/*/glue/*.c src/lib/lgpl/*/*.c)
+$.stmp-libc-std: $(wildcard src/lib/*/*.c src/lib/*/386/*.s src/lib/bsd/*/*.c src/lib/bsd/*/glue/*.c src/lib/lgpl/*/*.c src/lib/msun/*/*.c)
 	gawk -f src/lib/mkstd.awk $^ >__tmp__
 	$(call UPDATE,__tmp__,$.libc-std.h)
Index: /trunk/src/emx/src/lib/libc06b4.def
===================================================================
--- /trunk/src/emx/src/lib/libc06b4.def	(revision 2012)
+++ /trunk/src/emx/src/lib/libc06b4.def	(revision 2013)
@@ -235,7 +235,7 @@
     "___fmutex_request_internal_must_complete" @250
     "___fork" @251
-    "___fpclassify" @252
-    "___fpclassifyf" @253
-    "___fpclassifyl" @254
+    "___fpclassify" = ___nullstub_function @252
+    "___fpclassifyf" = ___nullstub_function @253
+    "___fpclassifyl" = ___nullstub_function @254
     "___free_ovflpage" @255
     "___ftime" @256
@@ -256,6 +256,6 @@
     "___isfinitef" @271
     "___isfinitel" @272
-    "___isnan" @273
-    "___isnanf" @274
+    "_isnan" @273
+    "_isnanf" @274
     "___isnanl" @275
     "___isnormal" @276
@@ -405,6 +405,6 @@
     "___signbitl" @420
     "___sigpause_bsd" @421
-    "___small_atod" @422
-    "___small_dtoa" @423
+    "___small_atod" = ___nullstub_function @422
+    "___small_dtoa" = ___legacy_dtoa @423
     "___smutex_request_internal" @424
     "___spawni" @425
Index: /trunk/src/emx/src/lib/math/386/acos.s
===================================================================
--- /trunk/src/emx/src/lib/math/386/acos.s	(revision 2012)
+++ /trunk/src/emx/src/lib/math/386/acos.s	(revision 2013)
@@ -1,7 +1,8 @@
 / acos.s (emx+gcc) -- Copyright (c) 1991-1996 by Eberhard Mattes
 
+#define LONG_DOUBLE
 #include <emx/asm386.h>
 
-#define FUNC    MATHSUFFIX1(acos)
+#define FUNC    _STD(acosl)
 
         .globl  FUNC
Index: /trunk/src/emx/src/lib/math/386/asin.s
===================================================================
--- /trunk/src/emx/src/lib/math/386/asin.s	(revision 2012)
+++ /trunk/src/emx/src/lib/math/386/asin.s	(revision 2013)
@@ -1,7 +1,8 @@
 / asin.s (emx+gcc) -- Copyright (c) 1991-1996 by Eberhard Mattes
 
+#define LONG_DOUBLE
 #include <emx/asm386.h>
 
-#define FUNC    MATHSUFFIX1(asin)
+#define FUNC    _STD(asinl)
 
         .globl  FUNC
Index: /trunk/src/emx/src/lib/math/386/atan.s
===================================================================
--- /trunk/src/emx/src/lib/math/386/atan.s	(revision 2012)
+++ /trunk/src/emx/src/lib/math/386/atan.s	(revision 2013)
@@ -1,7 +1,8 @@
 / atan.s (emx+gcc) -- Copyright (c) 1991-1996 by Eberhard Mattes
 
+#define LONG_DOUBLE
 #include <emx/asm386.h>
 
-#define FUNC    MATHSUFFIX1(atan)
+#define FUNC    _STD(atanl)
 
         .globl  FUNC
Index: /trunk/src/emx/src/lib/math/386/atan2.s
===================================================================
--- /trunk/src/emx/src/lib/math/386/atan2.s	(revision 2012)
+++ /trunk/src/emx/src/lib/math/386/atan2.s	(revision 2013)
@@ -1,7 +1,8 @@
 / atan2.s (emx+gcc) -- Copyright (c) 1991-1996 by Eberhard Mattes
 
+#define LONG_DOUBLE
 #include <emx/asm386.h>
 
-#define FUNC    MATHSUFFIX1(atan2)
+#define FUNC    _STD(atan2l)
 
         .globl  FUNC
Index: unk/src/emx/src/lib/math/386/ceil.s
===================================================================
--- /trunk/src/emx/src/lib/math/386/ceil.s	(revision 2012)
+++ 	(revision )
@@ -1,33 +1,0 @@
-/ ceil.s (emx+gcc) -- Copyright (c) 1991-1996 by Eberhard Mattes
-
-#include <emx/asm386.h>
-
-#define FUNC    MATHSUFFIX1(ceil)
-
-        .globl  FUNC
-
-        .text
-
-        ALIGN
-
-/ double ceil (double x)
-
-#define cw1       0(%esp)
-#define cw2       2(%esp)
-/define ret_addr  4(%esp)
-#define x         8(%esp)
-
-FUNC:
-        PROFILE_NOFRAME
-        subl    $4, %esp
-        fstcw   cw1
-        movw    cw1, %ax
-        andw    $0xf3ff, %ax
-        orw     $0x0800, %ax            /* round up towards +inf */
-        movw    %ax, cw2
-        fldcw   cw2
-        FLD     x                       /* x */
-        frndint
-        fldcw   cw1
-        addl    $4, %esp
-        EPILOGUE(FUNC)
Index: unk/src/emx/src/lib/math/386/copysig.s
===================================================================
--- /trunk/src/emx/src/lib/math/386/copysig.s	(revision 2012)
+++ 	(revision )
@@ -1,61 +1,0 @@
-/ copysig.s (emx+gcc) -- Copyright (c) 1996 by Eberhard Mattes
-
-#include <emx/asm386.h>
-
-#define FUNC    MATHSUFFIX1(copysign)
-
-        .globl  FUNC
-
-        .text
-
-        ALIGN
-
-#if defined (LONG_DOUBLE)
-
-/ long double copysignl (long double x, long double y)
-
-#define x(i)    4+i(%esp)
-#define y(i)    16+i(%esp)
-
-FUNC:
-        PROFILE_NOFRAME
-        movb    y(9), %al
-        andb    $0x7f, x(9)
-        andb    $0x80, %al
-        orb     %al, x(9)
-        fldt    x(0)
-        EPILOGUE(FUNC)
-
-#elif defined (FLOAT)
-
-/ float copysignf (float x, float y)
-
-#define x(i)    4+i(%esp)
-#define y(i)    8+i(%esp)
-
-FUNC:
-        PROFILE_NOFRAME
-        movb    y(3), %al
-        andb    $0x7f, x(3)
-        andb    $0x80, %al
-        orb     %al, x(3)
-        flds    x(0)
-        EPILOGUE(FUNC)
-
-#else
-
-/ double copysign (double x, double y)
-
-#define x(i)    4+i(%esp)
-#define y(i)    12+i(%esp)
-
-FUNC:
-        PROFILE_NOFRAME
-        movb    y(7), %al
-        andb    $0x7f, x(7)
-        andb    $0x80, %al
-        orb     %al, x(7)
-        fldl    x(0)
-        EPILOGUE(FUNC)
-
-#endif
Index: /trunk/src/emx/src/lib/math/386/cos.s
===================================================================
--- /trunk/src/emx/src/lib/math/386/cos.s	(revision 2012)
+++ /trunk/src/emx/src/lib/math/386/cos.s	(revision 2013)
@@ -2,7 +2,8 @@
 /                    Modified 1993-1996 by Eberhard Mattes
 
+#define LONG_DOUBLE
 #include <emx/asm386.h>
 
-#define FUNC    MATHSUFFIX1(cos)
+#define FUNC    _STD(cosl)
 
         .globl  FUNC
@@ -21,5 +22,5 @@
         fcos                            /* cos(x) */
         fstsw   %ax
-        testb   $0x04, %ah                              
+        testb   $0x04, %ah
         jnz     Llarge                  /* C2 != 0 ? */
 Lreturn:EPILOGUE(FUNC)
Index: /trunk/src/emx/src/lib/math/386/cosh.s
===================================================================
--- /trunk/src/emx/src/lib/math/386/cosh.s	(revision 2012)
+++ /trunk/src/emx/src/lib/math/386/cosh.s	(revision 2013)
@@ -1,7 +1,8 @@
 / cosh.s (emx+gcc) -- Copyright (c) 1992-1996 by Eberhard Mattes
 
+#define LONG_DOUBLE
 #include <emx/asm386.h>
 
-#define FUNC    MATHSUFFIX1(cosh)
+#define FUNC    _STD(coshl)
 
         .globl  FUNC
@@ -31,5 +32,5 @@
         orw     $0x0400, %ax            /* round down towards -inf */
         movw    %ax, cw2
-        fldcw   cw2      
+        fldcw   cw2
         fld     %st                     /* y, y */
         frndint                         /* int (y), y */
Index: /trunk/src/emx/src/lib/math/386/exp.s
===================================================================
--- /trunk/src/emx/src/lib/math/386/exp.s	(revision 2012)
+++ /trunk/src/emx/src/lib/math/386/exp.s	(revision 2013)
@@ -1,7 +1,8 @@
 / exp.s (emx+gcc) -- Copyright (c) 1991-2000 by Eberhard Mattes
 
+#define LONG_DOUBLE
 #include <emx/asm386.h>
 
-#define FUNC    MATHSUFFIX1(exp)
+#define FUNC    _STD(expl)
 
         .globl  FUNC
@@ -36,5 +37,5 @@
         orw     $0x0400, %ax            /* round down towards -inf */
         movw    %ax, cw2
-        fldcw   cw2      
+        fldcw   cw2
         fld     %st                     /* y, y */
         frndint                         /* int (y), y */
Index: unk/src/emx/src/lib/math/386/fabs.s
===================================================================
--- /trunk/src/emx/src/lib/math/386/fabs.s	(revision 2012)
+++ 	(revision )
@@ -1,21 +1,0 @@
-/ fabs.s (emx+gcc) -- Copyright (c) 1991-1996 by Eberhard Mattes
-
-#include <emx/asm386.h>
-
-#define FUNC    MATHSUFFIX1(fabs)
-
-        .globl  FUNC
-
-        .text
-
-        ALIGN
-
-/ double fabs (double x)
-
-#define x       4(%esp)
-
-FUNC:
-        PROFILE_NOFRAME
-        FLD     x                       /* x */
-        fabs
-        EPILOGUE(FUNC)
Index: unk/src/emx/src/lib/math/386/floor.s
===================================================================
--- /trunk/src/emx/src/lib/math/386/floor.s	(revision 2012)
+++ 	(revision )
@@ -1,33 +1,0 @@
-/ floor.s (emx+gcc) -- Copyright (c) 1991-1996 by Eberhard Mattes
-
-#include <emx/asm386.h>
-
-#define FUNC    MATHSUFFIX1(floor)
-
-        .globl  FUNC
-
-        .text
-
-        ALIGN
-
-/ double floor (double x)
-
-#define cw1       0(%esp)
-#define cw2       2(%esp)
-/define ret_addr  4(%esp)
-#define x         8(%esp)
-
-FUNC:
-        PROFILE_NOFRAME
-        subl    $4, %esp
-        fstcw   cw1
-        movw    cw1, %ax
-        andw    $0xf3ff, %ax
-        orw     $0x0400, %ax            /* round down towards -inf */
-        movw    %ax, cw2
-        fldcw   cw2
-        FLD     x                       /* x */
-        frndint
-        fldcw   cw1
-        addl    $4, %esp
-        EPILOGUE(FUNC)
Index: /trunk/src/emx/src/lib/math/386/fmod.s
===================================================================
--- /trunk/src/emx/src/lib/math/386/fmod.s	(revision 2012)
+++ /trunk/src/emx/src/lib/math/386/fmod.s	(revision 2013)
@@ -2,7 +2,8 @@
 /                     Modified 1993-1996 by Eberhard Mattes
 
+#define LONG_DOUBLE
 #include <emx/asm386.h>
 
-#define FUNC    MATHSUFFIX1(fmod)
+#define FUNC    _STD(fmodl)
 
         .globl  FUNC
Index: unk/src/emx/src/lib/math/386/fpclass.s
===================================================================
--- /trunk/src/emx/src/lib/math/386/fpclass.s	(revision 2012)
+++ 	(revision )
@@ -1,116 +1,0 @@
-/ fpclass.s (emx+gcc) -- Copyright (c) 1996 by Eberhard Mattes
-
-#include <emx/asm386.h>
-
-#define FUNC    MATHSUFFIX3(fpclassify)
-
-        .globl  FUNC
-
-        .text
-
-        ALIGN
-
-#if defined (LONG_DOUBLE)
-
-/ int __fpclassifyl (long double x)
-
-#define x       4(%esp)
-
-FUNC:
-        PROFILE_NOFRAME
-        fldt    x
-        fxam
-        fstsw   %ax
-        fstp    %st(0)
-        movb    %ah, %al
-        andl    $0x47, %eax
-        movb    Ltable(%eax), %al
-        EPILOGUE(FUNC)
-
-Ltable: .byte   FP_NAN          /* +unnormal, not supported by 387 */
-        .byte   FP_NAN          /* -unnormal, not supported by 387 */
-        .byte   FP_NAN          /* +NaN */
-        .byte   FP_NAN          /* -NaN */
-        .byte   FP_NORMAL       /* +normal */
-        .byte   FP_INFINITE     /* +infinity */
-        .byte   FP_NORMAL       /* -normal */
-        .byte   FP_INFINITE     /* -infinity */
-        .byte   FP_NAN, FP_NAN, FP_NAN, FP_NAN, FP_NAN, FP_NAN, FP_NAN, FP_NAN
-        .byte   FP_NAN, FP_NAN, FP_NAN, FP_NAN, FP_NAN, FP_NAN, FP_NAN, FP_NAN
-        .byte   FP_NAN, FP_NAN, FP_NAN, FP_NAN, FP_NAN, FP_NAN, FP_NAN, FP_NAN
-        .byte   FP_NAN, FP_NAN, FP_NAN, FP_NAN, FP_NAN, FP_NAN, FP_NAN, FP_NAN
-        .byte   FP_NAN, FP_NAN, FP_NAN, FP_NAN, FP_NAN, FP_NAN, FP_NAN, FP_NAN
-        .byte   FP_NAN, FP_NAN, FP_NAN, FP_NAN, FP_NAN, FP_NAN, FP_NAN, FP_NAN
-        .byte   FP_NAN, FP_NAN, FP_NAN, FP_NAN, FP_NAN, FP_NAN, FP_NAN, FP_NAN
-        .byte   FP_ZERO         /* +zero */
-        .byte   FP_NAN          /* empty */
-        .byte   FP_ZERO         /* -zero */
-        .byte   FP_NAN          /* empty */
-        .byte   FP_SUBNORMAL    /* +denormal */
-        .byte   FP_NAN          /* empty, not supported by 387 */
-        .byte   FP_SUBNORMAL    /* -denormal */
-        .byte   FP_NAN          /* empty, not supported by 387 */
-
-#elif defined (FLOAT)
-
-/ int __fpclassifyf (float x)
-
-#define x       4(%esp)
-
-FUNC:
-        PROFILE_NOFRAME
-        movl    x, %edx
-        movl    $FP_ZERO, %eax          /* Clear upper 24 bits */
-        andl    $0x7fffffff, %edx
-        jz      Lreturn
-        movb    $FP_SUBNORMAL, %al
-        cmpl    $0x00800000, %edx
-        jb      Lreturn
-        movb    $FP_NORMAL, %al
-        cmpl    $0x7f800000, %edx
-        jb      Lreturn
-        movb    $FP_INFINITE, %al
-        je      Lreturn
-        movb    $FP_NAN, %al
-        ALIGN
-Lreturn:
-        EPILOGUE(FUNC)
-
-#else
-
-/ int __fpclassify (double x)
-
-#define x0      4(%esp)
-#define x4      8(%esp)
-
-FUNC:
-        PROFILE_NOFRAME
-        movl    x4, %edx
-        movl    $FP_NORMAL, %eax                /* Clear upper 24 bits */
-        andl    $0x7ff00000, %edx
-        jz      Lexp0
-        cmpl    $0x7ff00000, %edx
-        je      Lexpmax
-        ALIGN
-Lreturn:EPILOGUE(FUNC)
-
-        ALIGN
-Lexp0:  movb    $FP_SUBNORMAL, %al
-        testl   $0x000fffff, x4
-        jnz     Lreturn
-        cmpl    $0, x0
-        jne     Lreturn
-        movb    $FP_ZERO, %al
-        jmp     Lreturn
-
-        ALIGN
-Lexpmax:testl   $0x000fffff, x4
-        jnz     Lnan
-        movb    $FP_INFINITE, %al
-        cmpl    $0, x0
-        je      Lreturn
-        ALIGN
-Lnan:   movb    $FP_NAN, %al
-        jmp     Lreturn
-
-#endif
Index: unk/src/emx/src/lib/math/386/frexp.s
===================================================================
--- /trunk/src/emx/src/lib/math/386/frexp.s	(revision 2012)
+++ 	(revision )
@@ -1,40 +1,0 @@
-/ frexp.s (emx+gcc) -- Copyright (c) 1992-1993 by Steffen Haecker
-/                      Modified 1993-1996 by Eberhard Mattes
-
-#include <emx/asm386.h>
-
-#define FUNC    MATHSUFFIX1(frexp)
-
-        .globl  FUNC
-
-        .text
-
-        ALIGN
-
-/ double frexp (double x, int *exp_ptr)
-
-#define x       4(%esp)
-#if defined (LONG_DOUBLE)
-#define exp_ptr 16(%esp)
-#else
-#define exp_ptr 12(%esp)
-#endif
-
-FUNC:
-        PROFILE_NOFRAME
-        movl    exp_ptr, %edx
-        movl    $0, (%edx)
-        fldl    __const_M_ONE           /* -1 */
-        FLD     x                       /* x */
-        ftst
-        fstsw   %ax
-        andw    $0x4100, %ax
-        xorb    $0x40, %ah
-        jz      1f                      /* x = 0 ? */
-        fxtract
-        fxch    %st(1)
-        fistpl  (%edx)
-        fscale
-        incl    (%edx)
-1:      fstp    %st(1)
-        EPILOGUE(FUNC)
Index: /trunk/src/emx/src/lib/math/386/fxam.s
===================================================================
--- /trunk/src/emx/src/lib/math/386/fxam.s	(revision 2012)
+++ /trunk/src/emx/src/lib/math/386/fxam.s	(revision 2013)
@@ -3,5 +3,5 @@
 #include <emx/asm386.h>
 
-#define FUNC    MATHSUFFIX2(fxam)
+#define FUNC    __fxam
 
         .globl  FUNC
Index: /trunk/src/emx/src/lib/math/386/fxaml.s
===================================================================
--- /trunk/src/emx/src/lib/math/386/fxaml.s	(revision 2013)
+++ /trunk/src/emx/src/lib/math/386/fxaml.s	(revision 2013)
@@ -0,0 +1,32 @@
+/ fxam.s (emx+gcc) -- Copyright (c) 1991-1996 by Eberhard Mattes
+
+#define LONG_DOUBLE
+#include <emx/asm386.h>
+
+#define FUNC    __fxaml
+
+        .globl  FUNC
+
+        .text
+
+        ALIGN
+
+/ int _fxam (double x)
+/ int _fxaml (long double x)
+
+#define x       4(%esp)
+
+FUNC:
+        PROFILE_NOFRAME
+        FLD     x                       /* x */
+        fxam
+        fstsw   %ax
+        fstp    %st(0)
+        movl    %eax, %edx
+        shrl    $8, %eax
+        andl    $7, %eax
+        test    $0x4000, %edx
+        jz      Lreturn
+        orl     $8, %eax
+Lreturn:
+        EPILOGUE(FUNC)
Index: /trunk/src/emx/src/lib/math/386/hypot.s
===================================================================
--- /trunk/src/emx/src/lib/math/386/hypot.s	(revision 2012)
+++ /trunk/src/emx/src/lib/math/386/hypot.s	(revision 2013)
@@ -5,7 +5,8 @@
 / However, such an algorithm should perhaps be used for _hypotl().
 
+#define LONG_DOUBLE
 #include <emx/asm386.h>
 
-#define FUNC    MATHSUFFIX1(hypot)
+#define FUNC    _STD(hypotl)
 
         .globl  FUNC
Index: unk/src/emx/src/lib/math/386/isfin.s
===================================================================
--- /trunk/src/emx/src/lib/math/386/isfin.s	(revision 2012)
+++ 	(revision )
@@ -1,85 +1,0 @@
-/ isfin.s (emx+gcc) -- Copyright (c) 1996 by Eberhard Mattes
-
-#include <emx/asm386.h>
-
-#define FUNC    MATHSUFFIX3(isfinite)
-
-        .globl  FUNC
-
-        .text
-
-        ALIGN
-
-#if defined (LONG_DOUBLE)
-
-/ int __isfinitel (long double x)
-
-#define x       4(%esp)
-
-FUNC:
-        PROFILE_NOFRAME
-        fldt    x
-        fxam
-        fstsw   %ax
-        fstp    %st(0)
-        movb    %ah, %al
-        andl    $0x47, %eax
-        movb    Ltable(%eax), %al
-        EPILOGUE(FUNC)
-
-Ltable: .byte   0               /* +unnormal, not supported by 387 */
-        .byte   0               /* -unnormal, not supported by 387 */
-        .byte   0               /* +NaN */
-        .byte   0               /* -NaN */
-        .byte   1               /* +normal */
-        .byte   0               /* +infinity */
-        .byte   1               /* -normal */
-        .byte   0               /* -infinity */
-        .byte   0, 0, 0, 0, 0, 0, 0, 0
-        .byte   0, 0, 0, 0, 0, 0, 0, 0
-        .byte   0, 0, 0, 0, 0, 0, 0, 0
-        .byte   0, 0, 0, 0, 0, 0, 0, 0
-        .byte   0, 0, 0, 0, 0, 0, 0, 0
-        .byte   0, 0, 0, 0, 0, 0, 0, 0
-        .byte   0, 0, 0, 0, 0, 0, 0, 0
-        .byte   1               /* +zero */
-        .byte   0          /* empty */
-        .byte   1               /* -zero */
-        .byte   0          /* empty */
-        .byte   1               /* +denormal */
-        .byte   0          /* empty, not supported by 387 */
-        .byte   1               /* -denormal */
-        .byte   0          /* empty, not supported by 387 */
-
-#elif defined (FLOAT)
-
-/ int __isfinitef (float x)
-
-#define x       4(%esp)
-
-FUNC:
-        PROFILE_NOFRAME
-        movl    x, %edx
-        xorl    %eax, %eax
-        andl    $0x7fffffff, %edx
-        cmpl    $0x7f800000, %edx
-        setb    %al
-        EPILOGUE(FUNC)
-
-#else
-
-/ int __isfinite (double x)
-
-#define x0      4(%esp)
-#define x4      8(%esp)
-
-FUNC:
-        PROFILE_NOFRAME
-        movl    x4, %edx
-        xorl    %eax, %eax
-        andl    $0x7ff00000, %edx
-        cmpl    $0x7ff00000, %edx
-        setne   %al
-        EPILOGUE(FUNC)
-
-#endif
Index: unk/src/emx/src/lib/math/386/isnan.s
===================================================================
--- /trunk/src/emx/src/lib/math/386/isnan.s	(revision 2012)
+++ 	(revision )
@@ -1,92 +1,0 @@
-/ isnan.s (emx+gcc) -- Copyright (c) 1996 by Eberhard Mattes
-
-#include <emx/asm386.h>
-
-#define FUNC    MATHSUFFIX3(isnan)
-
-        .globl  FUNC
-
-        .text
-
-        ALIGN
-
-#if defined (LONG_DOUBLE)
-
-/ int __isnanl (long double x)
-
-#define x       4(%esp)
-
-FUNC:
-        PROFILE_NOFRAME
-        fldt    x
-        fxam
-        fstsw   %ax
-        fstp    %st(0)
-        movb    %ah, %al
-        andl    $0x47, %eax
-        movb    Ltable(%eax), %al
-        EPILOGUE(FUNC)
-
-Ltable: .byte   1               /* +unnormal, not supported by 387 */
-        .byte   1               /* -unnormal, not supported by 387 */
-        .byte   1               /* +NaN */
-        .byte   1               /* -NaN */
-        .byte   0               /* +normal */
-        .byte   0               /* +infinity */
-        .byte   0               /* -normal */
-        .byte   0               /* -infinity */
-        .byte   1, 1, 1, 1, 1, 1, 1, 1
-        .byte   1, 1, 1, 1, 1, 1, 1, 1
-        .byte   1, 1, 1, 1, 1, 1, 1, 1
-        .byte   1, 1, 1, 1, 1, 1, 1, 1
-        .byte   1, 1, 1, 1, 1, 1, 1, 1
-        .byte   1, 1, 1, 1, 1, 1, 1, 1
-        .byte   1, 1, 1, 1, 1, 1, 1, 1
-        .byte   0               /* +zero */
-        .byte   1               /* empty */
-        .byte   0               /* -zero */
-        .byte   1               /* empty */
-        .byte   0               /* +denormal */
-        .byte   1               /* empty, not supported by 387 */
-        .byte   0               /* -denormal */
-        .byte   1               /* empty, not supported by 387 */
-
-#elif defined (FLOAT)
-
-/ int __isnanf (float x)
-
-#define x       4(%esp)
-
-FUNC:
-        PROFILE_NOFRAME
-        movl    x, %edx
-        xorl    %eax, %eax
-        andl    $0x7fffffff, %edx
-        cmpl    $0x7f800000, %edx
-        seta    %al
-        EPILOGUE(FUNC)
-
-#else
-
-/ int __isnan (double x)
-
-#define x0      4(%esp)
-#define x4      8(%esp)
-
-FUNC:
-        PROFILE_NOFRAME
-        movl    x4, %edx
-        xorl    %eax, %eax
-        andl    $0x7fffffff, %edx
-        cmpl    $0x7ff00000, %edx
-        je      Linf
-        seta    %al
-        ALIGN
-Lreturn:EPILOGUE(FUNC)
-
-        ALIGN
-Linf:   cmpl    $0, x0
-        setne   %al
-        jmp     Lreturn
-
-#endif
Index: unk/src/emx/src/lib/math/386/isnorm.s
===================================================================
--- /trunk/src/emx/src/lib/math/386/isnorm.s	(revision 2012)
+++ 	(revision )
@@ -1,64 +1,0 @@
-/ isnorm.s (emx+gcc) -- Copyright (c) 1996 by Eberhard Mattes
-
-#include <emx/asm386.h>
-
-#define FUNC    MATHSUFFIX3(isnormal)
-
-        .globl  FUNC
-
-        .text
-
-        ALIGN
-
-#if defined (LONG_DOUBLE)
-
-/ int __isnormall (long double x)
-
-#define x       4(%esp)
-
-FUNC:
-        PROFILE_NOFRAME
-        fldt    x
-        fxam
-        fstsw   %ax
-        fstp    %st(0)
-        movb    %ah, %al
-        andl    $0x45, %eax
-        cmpb    $0x04, %al
-        sete    %al
-        EPILOGUE(FUNC)
-
-#elif defined (FLOAT)
-
-/ int __isnormalf (float x)
-
-#define x       4(%esp)
-
-FUNC:
-        PROFILE_NOFRAME
-        movl    x, %edx
-        xorl    %eax, %eax
-        andl    $0x7fffffff, %edx
-        subl    $0x00800000, %edx
-        cmpl    $0x7f000000, %edx
-        setb    %al
-        EPILOGUE(FUNC)
-
-#else
-
-/ int __isnormal (double x)
-
-#define x0      4(%esp)
-#define x4      8(%esp)
-
-FUNC:
-        PROFILE_NOFRAME
-        movl    x4, %edx
-        xorl    %eax, %eax
-        andl    $0x7ff00000, %edx
-        subl    $0x00100000, %edx
-        cmpl    $0x7fe00000, %edx
-        setb    %al
-        EPILOGUE(FUNC)
-
-#endif
Index: unk/src/emx/src/lib/math/386/ldexp.s
===================================================================
--- /trunk/src/emx/src/lib/math/386/ldexp.s	(revision 2012)
+++ 	(revision )
@@ -1,35 +1,0 @@
-/ ldexp.s (emx+gcc) -- Copyright (c) 1991-1996 by Eberhard Mattes
-
-#include <emx/asm386.h>
-
-#define FUNC    MATHSUFFIX1(ldexp)
-
-        .globl  FUNC
-
-        .text
-
-        ALIGN
-
-/ double ldexp (double x, int exp)
-
-#define x        4(%esp)
-#if defined (LONG_DOUBLE)
-#define exp     16(%esp)
-#else
-#define exp     12(%esp)
-#endif
-
-FUNC:
-        PROFILE_NOFRAME
-        fildl   exp                     /* exp */
-        FLD     x                       /* x */
-        fscale
-        fstp    %st(1)
-        CONV(x)                         /* convert to double */
-        _xam
-        j_inf   Lerror
-Lreturn:EPILOGUE(FUNC)
-
-        ALIGN
-Lerror: SET_ERRNO_CONST($ERANGE)
-        jmp     Lreturn
Index: /trunk/src/emx/src/lib/math/386/log.s
===================================================================
--- /trunk/src/emx/src/lib/math/386/log.s	(revision 2012)
+++ /trunk/src/emx/src/lib/math/386/log.s	(revision 2013)
@@ -1,7 +1,8 @@
 / log.s (emx+gcc) -- Copyright (c) 1991-2000 by Eberhard Mattes
 
+#define LONG_DOUBLE
 #include <emx/asm386.h>
 
-#define FUNC    MATHSUFFIX1(log)
+#define FUNC    _STD(logl)
 
         .globl  FUNC
Index: /trunk/src/emx/src/lib/math/386/log10.s
===================================================================
--- /trunk/src/emx/src/lib/math/386/log10.s	(revision 2012)
+++ /trunk/src/emx/src/lib/math/386/log10.s	(revision 2013)
@@ -1,7 +1,8 @@
 / log10.s (emx+gcc) -- Copyright (c) 1991-2000 by Eberhard Mattes
 
+#define LONG_DOUBLE
 #include <emx/asm386.h>
 
-#define FUNC    MATHSUFFIX1(log10)
+#define FUNC    _STD(log10l)
 
         .globl  FUNC
Index: /trunk/src/emx/src/lib/math/386/modf.s
===================================================================
--- /trunk/src/emx/src/lib/math/386/modf.s	(revision 2012)
+++ /trunk/src/emx/src/lib/math/386/modf.s	(revision 2013)
@@ -1,7 +1,8 @@
 / modf.s (emx+gcc) -- Copyright (c) 1991-1996 by Eberhard Mattes
 
+#define LONG_DOUBLE
 #include <emx/asm386.h>
 
-#define FUNC    MATHSUFFIX1(modf)
+#define FUNC    _STD(modfl)
 
         .globl  FUNC
Index: unk/src/emx/src/lib/math/386/nextaft.s
===================================================================
--- /trunk/src/emx/src/lib/math/386/nextaft.s	(revision 2012)
+++ 	(revision )
@@ -1,323 +1,0 @@
-/ nextaft.s (emx+gcc) -- Copyright (c) 1996-1998 by Eberhard Mattes
-
-#include <emx/asm386.h>
-
-#define FUNC    MATHSUFFIX1(nextafter)
-
-        .globl  FUNC
-
-        .text
-
-        ALIGN
-
-#if defined (LONG_DOUBLE)
-
-/ long double nextafterl (long double x, long double y)
-
-#define x(i)    4+i(%esp)
-#define y(i)    16+i(%esp)
-
-FUNC:
-        PROFILE_NOFRAME
-        fldt    x(0)
-        fldt    y(0)
-        fucompp
-        fstsw   %ax
-        testb   $0x04, %ah
-        jnz     Lunordered              /* Unordered (at least one NaN) */
-        testb   $0x40, %ah
-        jnz     Leq                     /* x==y */
-        movb    x(9), %al
-        shrb    $7, %al                 /* signbit(x) -> bit 0 of %al */
-        xorb    %ah, %al                /* Note: bit 0 of %ah is 1 iff x>y */
-        testb   $0x01, %al
-        movw    x(8), %ax               /* Fetch exponent */
-        jnz     Ldown                   /* Decrement mantissa */
-
-/* Increment mantissa. */
-
-Lup:    andl    $0x7fff, %eax           /* Clear sign bit */
-        cmpl    $0x7fff, %eax           /* Infinity? */
-        je      Lreturn_x               /* Yes => return x */
-        addl    $1, x(0)                /* Increment lower word */
-        adcl    $0, x(4)                /* Propagate carry into upper word */
-        jc      Lup_of                  /* Overflow */
-        testl   %eax, %eax              /* Denormal? */
-        jnz     Lreturn_x               /* No => return x */
-        testb   $0x80, x(7)             /* Pseudo-denormal? */
-        jz      Lreturn_x_exception     /* No => return x, raise exceptions */
-        incw    x(8)                    /* Turn it into a normal */
-        jmp     Lreturn_x
-
-/* Return x and raise appropriate exceptions (inexact and
-   underflow/overflow). */
-
-        ALIGN
-Lreturn_x_exception:
-        fldt    x(0)
-        fstl    y(0)                    /* This raises the exceptions */
-        jmp     Lreturn                 /* (we store as 64-bit number!) */
-
-/* Mantissa overflow.  Increment exponent. */
-
-        ALIGN
-Lup_of: cmpl    $0x7ffe, %eax           /* Exponent overflow? */
-        je      Lup_make_inf            /* Yes => return infinity */
-        movb    $0x80, x(7)             /* Set mantissa to 100...000 */
-        incw    x(8)                    /* Increment exponent */
-        jmp     Lreturn_x               /* Return x */
-
-/* x is LDBL_MAX, raise overflow and inexact exceptions,
-   return infinity. */
-
-        ALIGN
-Lup_make_inf:
-        movl    $0xffffffff, x(0)       /* Restore mantissa of x */
-        movl    $0xffffffff, x(4)
-        fldt    x(0)                    /* Compute x+x */
-        fadd    %st
-        jmp     Lreturn
-
-/* Decrement mantissa. */
-
-        ALIGN
-Ldown:  andl    $0x7fff, %eax           /* Clear sign bit */
-        cmpl    $0x7fff, %eax           /* Infinity? */
-        je      Ldown_is_inf            /* Yes => return biggest normal */
-        subl    $1, x(0)                /* Decrement lower word */
-        sbbl    $0, x(4)                /* Propagate carry into upper word */
-        jc      Ldown_uf0               /* Underflow (integer bit was 0) */
-        jo      Ldown_uf1               /* Underflow (integer bit was 1) */
-        test    %eax, %eax              /* Is the exponent zero? */
-        jnz     Lreturn_x               /* No => return x */
-        jmp     Lreturn_x_exception     /* Yes => maybe raise exceptions */
-
-/* Underflow (mantissa was 100...000 and is now 011...111).  Decrement
-   the exponent.  Return the biggest denormal if the exponent becomes 0. */
-
-        ALIGN
-Ldown_uf1:
-        decw    x(8)                    /* Decrement exponent */
-        cmp     $1, %eax                /* Exponent now zero? */
-        je      Lreturn_x_exception     /* Yes => return x (mantissa OK) */
-
-/* It's still a normal.  Set the mantissa to 111...111. */
-
-        orb     $0x80, x(7)
-        jmp     Lreturn_x
-
-/* Underflow (mantissa was 000...000 and is now 111...111.  Return
-   the smallest denormal with y's sign. */
-
-        ALIGN
-Ldown_uf0:
-        movl    $1, y(0)
-        movl    $0, y(4)
-        andw    $0x8000, y(8)
-        fldt    y(0)
-        fstl    y(0)                    /* This raises the exceptions */
-        jmp     Lreturn
-
-/* x is infinite.  Return the biggest normal with x's sign. */
-
-        ALIGN
-Ldown_is_inf:
-        movl    $0xffffffff, x(0)
-        movl    $0xffffffff, x(4)
-        andw    $0x8000, x(8)
-        orw     $0x7ffe, x(8)
-        jmp     Lreturn_x
-
-/* Unordered: return one of the NaNs. */
-
-        ALIGN
-Lunordered:
-        fldt    x(0)
-        _xam
-        j_nan   Lreturn                 /* x is a Nan => return x */
-        fstp    %st(0)                  /* Remove x */
-
-/* Fall through to return y */
-
-/* x==y: return y. */
-
-        ALIGN
-Leq:    fldt    y(0)
-        jmp     Lreturn
-
-        ALIGN
-Lreturn_x:
-        fldt    x(0)
-
-        ALIGN
-Lreturn:
-        EPILOGUE(FUNC)
-
-#elif defined (FLOAT)
-
-/ float nextafterf (float x, float y)
-
-#define x       4(%esp)
-#define y       8(%esp)
-
-FUNC:
-        PROFILE_NOFRAME
-        flds    x
-        flds    y
-        fucompp
-        fstsw   %ax
-        testb   $0x04, %ah
-        jnz     Lunordered              /* Unordered (at least one NaN) */
-        testb   $0x40, %ah
-        jnz     Leq                     /* x==y */
-        movl    x, %edx
-        shrl    $31, %edx               /* signbit(x) -> bit 0 of %dl */
-        xorb    %ah, %dl                /* Note: bit 0 of %ah is 1 iff x>y */
-        testb   $0x01, %dl
-        movl    $1, %ecx                /* Increment */
-        jz      Lcommon
-        movl    $-1, %ecx               /* TODO: Do this without jump */
-        ALIGN
-Lcommon:movl    x, %eax
-        andl    $0x7fffffff, %eax       /* Clear sign bit */
-        jz      Lx_zero                 /* x = 0.0 => smallest denormal */
-        addl    %ecx, x
-        movl    x, %edx
-        andl    $0x7fffffff, %edx
-        jz      Lreturn_x               /* 0.0 => no underflow */
-        cmpl    $0x7f800000, %edx
-        je      Loverflow
-        cmpl    $0x00800000, %edx
-        jb      Lunderflow
-        ALIGN
-Lreturn_x:
-        flds    x
-        ALIGN
-Lreturn:
-        EPILOGUE(FUNC)
-
-        ALIGN
-Lx_zero:
-        movl    y, %eax
-        andl    $0x80000000, %eax       /* Extract sign bit of y */
-        orl     $0x00000001, %eax       /* Smallest denormal */
-        movl    %eax, x
-        ALIGN
-Lunderflow:
-        fldl    Lsmall
-        fstps   y
-        jmp     Lreturn_x
-
-        ALIGN
-Loverflow:
-        fldl    Lbig
-        fstps   y                       /* Raise overflow and inexact */
-        jmp     Lreturn_x
-
-        ALIGN
-Lunordered:
-        movl    x, %eax
-        andl    $0x7fffffff, %eax
-        cmpl    $0x7f800000, %eax
-        ja      Lreturn_x               /* x is a NaN => return x */
-
-/* Fall through to return y (which is a NaN) */
-
-        ALIGN
-Leq:    flds    y
-        jmp     Lreturn
-
-Lsmall: .long   0x00000001, 0x00000000
-Lbig:   .long   0xffffffff, 0x7fefffff
-
-#else
-
-/ double nextafter (double x, double y)
-
-#define x(i)    4+i(%esp)
-#define y(i)    12+i(%esp)
-
-FUNC:
-        PROFILE_NOFRAME
-        fldl    x(0)
-        fldl    y(0)
-        fucompp
-        fstsw   %ax
-        testb   $0x04, %ah
-        jnz     Lunordered              /* Unordered (at least one NaN) */
-        testb   $0x40, %ah
-        jnz     Leq                     /* x==y */
-
-        movb    x(7), %dl
-        shrb    $7, %dl                 /* signbit(x) -> bit 0 of %dl */
-        xorb    %ah, %dl                /* Note: bit 0 of %ah is 1 iff x>y */
-
-        movl    x(4), %eax
-        andl    $0x7fffffff, %eax       /* Clear sign bit */
-        orl     x(0), %eax              /* x == 0.0? */
-        jz      Lx_zero                 /* Yes => smallest denormal */
-        testb   $0x01, %dl
-        jnz     Ldec
-        addl    $1, x(0)
-        adcl    $0, x(4)
-        ALIGN
-Lcommon:movl    x(4), %eax
-        andl    $0x7fffffff, %eax       /* Clear sign bit */
-        orl     x(0), %eax
-        jz      Lreturn_x               /* 0.0 => no underflow */
-        movl    x(4), %eax
-        andl    $0x7ff00000, %eax       /* Extract exponent */
-        cmpl    $0x7ff00000, %eax
-        je      Loverflow
-        cmpl    $0x00100000, %eax
-        jb      Lunderflow
-        ALIGN
-Lreturn_x:
-        fldl    x(0)
-        ALIGN
-Lreturn:
-        EPILOGUE(FUNC)
-
-        ALIGN
-Ldec:   subl    $1, x(0)
-        sbbl    $0, x(4)
-        jmp     Lcommon
-
-        ALIGN
-Lx_zero:
-        andl    $0x80000000, y(4)       /* Keep sign bit of y */
-        movl    $0x00000001, y(0)       /* Smallest denormal */
-        fldl    Lsmall
-        fstps   x(0)
-        fldl    y(0)
-        jmp     Lreturn
-
-        ALIGN
-Lunderflow:
-        fldl    Lsmall
-        fstps   y(0)
-        jmp     Lreturn_x
-
-        ALIGN
-Loverflow:
-        fldl    Lbig
-        fstps   y(0)                    /* Raise overflow and inexact */
-        jmp     Lreturn_x
-
-        ALIGN
-Lunordered:
-        movl    x(4), %eax
-        andl    $0x7fffffff, %eax
-        cmpl    $0x7ff80000, %eax
-        jae     Lreturn_x               /* x is a NaN => return x */
-
-/* Fall through to return y (which is a NaN) */
-
-        ALIGN
-Leq:    fldl    y(0)
-        jmp     Lreturn
-
-Lsmall: .long   0x00000001, 0x00000000
-Lbig:   .long   0xffffffff, 0x7fefffff
-
-#endif
Index: /trunk/src/emx/src/lib/math/386/pow.s
===================================================================
--- /trunk/src/emx/src/lib/math/386/pow.s	(revision 2012)
+++ /trunk/src/emx/src/lib/math/386/pow.s	(revision 2013)
@@ -2,7 +2,8 @@
 /                    Modified 1992-1996 by Eberhard Mattes
 
+#define LONG_DOUBLE
 #include <emx/asm386.h>
 
-#define FUNC    MATHSUFFIX1(pow)
+#define FUNC    _STD(powl)
 
         .globl  FUNC
@@ -116,5 +117,5 @@
         jmp     Ldone
 
-/ x = +INF, y non-zero and finite 
+/ x = +INF, y non-zero and finite
         ALIGN
 Lx_pinf:fstp    %st(1)                  /* pop x, st(0):=x */
@@ -125,5 +126,5 @@
         jmp     Ldone
 
-/ x = -INF, y non-zero and finite 
+/ x = -INF, y non-zero and finite
         ALIGN
 Lx_minf:testb   $2, %cl                 /* y < 0? */
Index: /trunk/src/emx/src/lib/math/386/rint.s
===================================================================
--- /trunk/src/emx/src/lib/math/386/rint.s	(revision 2012)
+++ /trunk/src/emx/src/lib/math/386/rint.s	(revision 2013)
@@ -1,7 +1,8 @@
 / rint.s (emx+gcc) -- Copyright (c) 1992-1996 by Eberhard Mattes
 
+#define LONG_DOUBLE
 #include <emx/asm386.h>
 
-#define FUNC    MATHSUFFIX1(rint)
+#define FUNC    _STD(rintl)
 
         .globl  FUNC
Index: unk/src/emx/src/lib/math/386/signbit.s
===================================================================
--- /trunk/src/emx/src/lib/math/386/signbit.s	(revision 2012)
+++ 	(revision )
@@ -1,55 +1,0 @@
-/ signbit.s (emx+gcc) -- Copyright (c) 1996 by Eberhard Mattes
-
-#include <emx/asm386.h>
-
-#define FUNC    MATHSUFFIX3(signbit)
-
-        .globl  FUNC
-
-        .text
-
-        ALIGN
-
-#if defined (LONG_DOUBLE)
-
-/ int __signbitl (long double x)
-
-#define x0      4(%esp)
-#define x9      13(%esp)
-
-FUNC:
-        PROFILE_NOFRAME
-        xorl    %eax, %eax
-        testb   $0x80, x9
-        setnz   %al
-        EPILOGUE(FUNC)
-
-#elif defined (FLOAT)
-
-/ int __signbitf (float x)
-
-#define x0      4(%esp)
-#define x3      7(%esp)
-
-FUNC:
-        PROFILE_NOFRAME
-        xorl    %eax, %eax
-        testb   $0x80, x3
-        setnz   %al
-        EPILOGUE(FUNC)
-
-#else
-
-/ int __signbit (double x)
-
-#define x0      4(%esp)
-#define x7      11(%esp)
-
-FUNC:
-        PROFILE_NOFRAME
-        xorl    %eax, %eax
-        testb   $0x80, x7
-        setnz   %al
-        EPILOGUE(FUNC)
-
-#endif
Index: /trunk/src/emx/src/lib/math/386/sin.s
===================================================================
--- /trunk/src/emx/src/lib/math/386/sin.s	(revision 2012)
+++ /trunk/src/emx/src/lib/math/386/sin.s	(revision 2013)
@@ -2,7 +2,8 @@
 /                    Modified 1993-1996 by Eberhard Mattes
 
+#define LONG_DOUBLE
 #include <emx/asm386.h>
 
-#define FUNC    MATHSUFFIX1(sin)
+#define FUNC    _STD(sinl)
 
         .globl  FUNC
@@ -21,5 +22,5 @@
         fsin                            /* sin(x) */
         fstsw   %ax
-        testb   $0x04, %ah                              
+        testb   $0x04, %ah
         jnz     Llarge                  /* C2 != 0 ? */
 Lreturn:EPILOGUE(FUNC)
Index: /trunk/src/emx/src/lib/math/386/sinh.s
===================================================================
--- /trunk/src/emx/src/lib/math/386/sinh.s	(revision 2012)
+++ /trunk/src/emx/src/lib/math/386/sinh.s	(revision 2013)
@@ -1,7 +1,8 @@
 / sinh.s (emx+gcc) -- Copyright (c) 1992-1996 by Eberhard Mattes
 
+#define LONG_DOUBLE
 #include <emx/asm386.h>
 
-#define FUNC    MATHSUFFIX1(sinh)
+#define FUNC    _STD(sinhl)
 
         .globl  FUNC
@@ -31,5 +32,5 @@
         orw     $0x0400, %ax            /* round down towards -inf */
         movw    %ax, cw2
-        fldcw   cw2      
+        fldcw   cw2
         fld     %st                     /* y, y */
         frndint                         /* int (y), y */
Index: /trunk/src/emx/src/lib/math/386/sqrt.s
===================================================================
--- /trunk/src/emx/src/lib/math/386/sqrt.s	(revision 2012)
+++ /trunk/src/emx/src/lib/math/386/sqrt.s	(revision 2013)
@@ -1,7 +1,8 @@
 / sqrt.s (emx+gcc) -- Copyright (c) 1991-1996 by Eberhard Mattes
 
+#define LONG_DOUBLE
 #include <emx/asm386.h>
 
-#define FUNC    MATHSUFFIX1(sqrt)
+#define FUNC    _STD(sqrtl)
 
         .globl  FUNC
Index: /trunk/src/emx/src/lib/math/386/tan.s
===================================================================
--- /trunk/src/emx/src/lib/math/386/tan.s	(revision 2012)
+++ /trunk/src/emx/src/lib/math/386/tan.s	(revision 2013)
@@ -2,7 +2,8 @@
 /                    Modified 1993-1996 by Eberhard Mattes
 
+#define LONG_DOUBLE
 #include <emx/asm386.h>
 
-#define FUNC    MATHSUFFIX1(tan)
+#define FUNC    _STD(tanl)
 
         .globl  FUNC
Index: /trunk/src/emx/src/lib/math/386/tanh.s
===================================================================
--- /trunk/src/emx/src/lib/math/386/tanh.s	(revision 2012)
+++ /trunk/src/emx/src/lib/math/386/tanh.s	(revision 2013)
@@ -1,7 +1,8 @@
 / tanh.s (emx+gcc) -- Copyright (c) 1992-1996 by Eberhard Mattes
 
+#define LONG_DOUBLE
 #include <emx/asm386.h>
 
-#define FUNC    MATHSUFFIX1(tanh)
+#define FUNC    _STD(tanhl)
 
         .globl  FUNC
@@ -32,5 +33,5 @@
         orw     $0x0400, %ax            /* round down towards -inf */
         movw    %ax, cw2
-        fldcw   cw2      
+        fldcw   cw2
         fld     %st                     /* y, y */
         frndint                         /* int (y), y */
Index: unk/src/emx/src/lib/math/386/trunc.s
===================================================================
--- /trunk/src/emx/src/lib/math/386/trunc.s	(revision 2012)
+++ 	(revision )
@@ -1,33 +1,0 @@
-/ trunc.s (emx+gcc) -- Copyright (c) 1992-1996 by Eberhard Mattes
-
-#include <emx/asm386.h>
-
-#define FUNC    MATHSUFFIX1(trunc)
-
-        .globl  FUNC
-
-        .text
-
-        ALIGN
-
-/ double _trunc (double x)
-
-#define cw1       0(%esp)
-#define cw2       2(%esp)
-/define ret_addr  4(%esp)
-#define x         8(%esp)
-
-FUNC:
-        PROFILE_NOFRAME
-        subl    $4, %esp
-        fstcw   cw1
-        movw    cw1, %ax
-        andw    $0xf3ff, %ax
-        orw     $0x0c00, %ax            /* chop by truncating toward 0 */
-        movw    %ax, cw2
-        fldcw   cw2
-        FLD     x                       /* x */
-        frndint
-        fldcw   cw1
-        addl    $4, %esp
-        EPILOGUE(FUNC)
Index: unk/src/emx/src/lib/math/cbrt.c
===================================================================
--- /trunk/src/emx/src/lib/math/cbrt.c	(revision 2012)
+++ 	(revision )
@@ -1,12 +1,0 @@
-/* cbrt.c (emx+gcc) -- Copyright (c) 1992-1995 by Eberhard Mattes */
-
-#include "libc-alias.h"
-#include <math.h>
-
-double _STD(cbrt) (double x)
-{
-  if (x >= 0)
-    return pow (x, 1.0 / 3.0);
-  else
-    return -pow (-x, 1.0 / 3.0);
-}
Index: /trunk/src/emx/src/lib/math/cbrtl.c
===================================================================
--- /trunk/src/emx/src/lib/math/cbrtl.c	(revision 2012)
+++ /trunk/src/emx/src/lib/math/cbrtl.c	(revision 2013)
@@ -7,6 +7,6 @@
 {
   if (x >= 0)
-    return _powl (x, 1.0 / 3.0);
+    return powl (x, 1.0 / 3.0);
   else
-    return -_powl (-x, 1.0 / 3.0);
+    return -powl (-x, 1.0 / 3.0);
 }
Index: /trunk/src/emx/src/lib/math/math.smak
===================================================================
--- /trunk/src/emx/src/lib/math/math.smak	(revision 2012)
+++ /trunk/src/emx/src/lib/math/math.smak	(revision 2013)
@@ -2,31 +2,2 @@
 libc.TSRC += $(wildcard src/lib/math/*.c)
 
-# Virtual targets, e.g. compiled from same input files but with different opts
-libc.TDEP.386 += $(addprefix @O@src/lib/math/386/, copysigf.o copysigl.o \
-  fpclassf.o fpclassl.o isfinf.o isfinl.o isnanf.o isnanl.o isnormf.o \
-  isnorml.o nextaftf.o nextaftl.o signbitf.o signbitl.o acosl.o asinl.o \
-  atanl.o atan2l.o ceill.o cosl.o coshl.o expl.o fabsl.o floorl.o fmodl.o \
-  frexpl.o fxaml.o hypotl.o ldexpl.o logl.o log10l.o modfl.o powl.o rintl.o \
-  sinl.o sinhl.o sqrtl.o tanl.o tanhl.o truncl.o)
-
-# bird: this doesn't work any better now. sorry. :/
-#$(call .MVER,math/386/%f.o): src/lib/math/386/%.s
-#	$(call DO.COMPILE.s, -DFLOAT)
-#$(call .MVER,math/386/%l.o): src/lib/math/386/%.s
-#	$(call DO.COMPILE.s, -DLONG_DOUBLE)
-	
-$.aout-prof/src/lib/math/386/%f.o: src/lib/math/386/%.s
-	$(call DO.COMPILE.s, -DFLOAT)
-$.aout-prof/src/lib/math/386/%l.o: src/lib/math/386/%.s
-	$(call DO.COMPILE.s, -DLONG_DOUBLE)
-
-$.aout-log/src/lib/math/386/%f.o: src/lib/math/386/%.s
-	$(call DO.COMPILE.s, -DFLOAT)
-$.aout-log/src/lib/math/386/%l.o: src/lib/math/386/%.s
-	$(call DO.COMPILE.s, -DLONG_DOUBLE)
-	
-$.aout/src/lib/math/386/%f.o: src/lib/math/386/%.s
-	$(call DO.COMPILE.s, -DFLOAT)
-$.aout/src/lib/math/386/%l.o: src/lib/math/386/%.s
-	$(call DO.COMPILE.s, -DLONG_DOUBLE)
-	
Index: /trunk/src/emx/src/lib/mkstd.awk
===================================================================
--- /trunk/src/emx/src/lib/mkstd.awk	(revision 2012)
+++ /trunk/src/emx/src/lib/mkstd.awk	(revision 2013)
@@ -54,4 +54,12 @@
 }
 
+/STDENTRY\( *[_0-9A-Za-z]+ *\)/{
+  while (match($0, "STDENTRY\\( *[_0-9A-Za-z ]+ *\\)"))
+  {
+    std_fun[gensub("STDENTRY\\( *([_0-9A-Za-z]+) *\\).*", "\\1", "", substr($0,RSTART,RLENGTH))]=1;
+    $0 = substr($0, 1, RSTART-1) substr($0, RSTART+RLENGTH);
+  }
+}
+
 /__weak_reference\( *[_0-9A-Za-z]+ *, *[_0-9A-Za-z]+ *\)/{
   while (match($0, "__weak_reference\\( *[_0-9A-Za-z ]+ *, *[_0-9A-Za-z ]+ *\\)"))
Index: /trunk/src/emx/src/lib/msun/bsdsrc/b_exp.c
===================================================================
--- /trunk/src/emx/src/lib/msun/bsdsrc/b_exp.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/bsdsrc/b_exp.c	(revision 2013)
@@ -32,4 +32,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char sccsid[] = "@(#)exp.c	8.1 (Berkeley) 6/4/93";
Index: /trunk/src/emx/src/lib/msun/bsdsrc/b_log.c
===================================================================
--- /trunk/src/emx/src/lib/msun/bsdsrc/b_log.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/bsdsrc/b_log.c	(revision 2013)
@@ -32,4 +32,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char sccsid[] = "@(#)log.c	8.2 (Berkeley) 11/30/93";
Index: /trunk/src/emx/src/lib/msun/bsdsrc/b_tgamma.c
===================================================================
--- /trunk/src/emx/src/lib/msun/bsdsrc/b_tgamma.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/bsdsrc/b_tgamma.c	(revision 2013)
@@ -32,4 +32,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char sccsid[] = "@(#)gamma.c	8.1 (Berkeley) 6/4/93";
Index: /trunk/src/emx/src/lib/msun/msun.smak
===================================================================
--- /trunk/src/emx/src/lib/msun/msun.smak	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/msun.smak	(revision 2013)
@@ -59,5 +59,5 @@
 
 # Add to libc modules.
-libc.TSRC.386 += $(addprefix src/lib/msun/i386/,$(ARCH_SRCS))
+libc.TSRC.386 += $(addprefix src/lib/msun/i387/,$(ARCH_SRCS))
 libc.TSRC += $(MSUN_SRCS)
 
Index: /trunk/src/emx/src/lib/msun/src/e_acos.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_acos.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_acos.c	(revision 2013)
@@ -7,9 +7,10 @@
  * Developed at SunSoft, a Sun Microsystems, Inc. business.
  * Permission to use, copy, modify, and distribute this
- * software is freely granted, provided that this notice 
+ * software is freely granted, provided that this notice
  * is preserved.
  * ====================================================
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_acos.c,v 1.10 2005/02/04 18:26:05 das Exp $";
@@ -17,5 +18,5 @@
 
 /* __ieee754_acos(x)
- * Method :                  
+ * Method :
  *	acos(x)  = pi/2 - asin(x)
  *	acos(-x) = pi/2 + asin(x)
@@ -24,5 +25,5 @@
  * For x>0.5
  * 	acos(x) = pi/2 - (pi/2 - 2asin(sqrt((1-x)/2)))
- *		= 2asin(sqrt((1-x)/2))  
+ *		= 2asin(sqrt((1-x)/2))
  *		= 2s + 2s*z*R(z) 	...z=(1-x)/2, s=sqrt(z)
  *		= 2f + (2c + 2s*z*R(z))
Index: /trunk/src/emx/src/lib/msun/src/e_acosf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_acosf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_acosf.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_acosf.c,v 1.7 2002/05/28 17:03:12 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/e_acosh.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_acosh.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_acosh.c	(revision 2013)
@@ -7,5 +7,5 @@
  * Developed at SunSoft, a Sun Microsystems, Inc. business.
  * Permission to use, copy, modify, and distribute this
- * software is freely granted, provided that this notice 
+ * software is freely granted, provided that this notice
  * is preserved.
  * ====================================================
@@ -13,4 +13,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_acosh.c,v 1.8 2005/02/04 18:26:05 das Exp $";
@@ -19,5 +20,5 @@
 /* __ieee754_acosh(x)
  * Method :
- *	Based on 
+ *	Based on
  *		acosh(x) = log [ x + sqrt(x*x-1) ]
  *	we have
@@ -50,5 +51,5 @@
 	    if(hx >=0x7ff00000) {	/* x is inf of NaN */
 	        return x+x;
-	    } else 
+	    } else
 		return __ieee754_log(x)+ln2;	/* acosh(huge)=log(2x) */
 	} else if(((hx-0x3ff00000)|lx)==0) {
Index: /trunk/src/emx/src/lib/msun/src/e_acoshf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_acoshf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_acoshf.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_acoshf.c,v 1.7 2002/05/28 17:03:12 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/e_asin.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_asin.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_asin.c	(revision 2013)
@@ -7,9 +7,10 @@
  * Developed at SunSoft, a Sun Microsystems, Inc. business.
  * Permission to use, copy, modify, and distribute this
- * software is freely granted, provided that this notice 
+ * software is freely granted, provided that this notice
  * is preserved.
  * ====================================================
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_asin.c,v 1.11 2005/02/04 18:26:05 das Exp $";
@@ -17,10 +18,10 @@
 
 /* __ieee754_asin(x)
- * Method :                  
+ * Method :
  *	Since  asin(x) = x + x^3/6 + x^5*3/40 + x^7*15/336 + ...
  *	we approximate asin(x) on [0,0.5] by
  *		asin(x) = x + x*x^2*R(x^2)
  *	where
- *		R(x^2) is a rational approximation of (asin(x)-x)/x^3 
+ *		R(x^2) is a rational approximation of (asin(x)-x)/x^3
  *	and its remez error is bounded by
  *		|(asin(x)-x)/x^3 - R(x^2)| < 2^(-58.75)
@@ -81,9 +82,9 @@
 		    /* asin(1)=+-pi/2 with inexact */
 		return x*pio2_hi+x*pio2_lo;	
-	    return (x-x)/(x-x);		/* asin(|x|>1) is NaN */   
+	    return (x-x)/(x-x);		/* asin(|x|>1) is NaN */
 	} else if (ix<0x3fe00000) {	/* |x|<0.5 */
 	    if(ix<0x3e400000) {		/* if |x| < 2**-27 */
 		if(huge+x>one) return x;/* return x with inexact if x!=0*/
-	    } else 
+	    } else
 		t = x*x;
 		p = t*(pS0+t*(pS1+t*(pS2+t*(pS3+t*(pS4+t*pS5)))));
@@ -109,5 +110,5 @@
 	    q  = pio4_hi-2.0*w;
 	    t  = pio4_hi-(p-q);
-	}    
-	if(hx>0) return t; else return -t;    
+	}
+	if(hx>0) return t; else return -t;
 }
Index: /trunk/src/emx/src/lib/msun/src/e_asinf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_asinf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_asinf.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_asinf.c,v 1.8 2002/05/28 17:03:12 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/e_atan2.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_atan2.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_atan2.c	(revision 2013)
@@ -7,5 +7,5 @@
  * Developed at SunSoft, a Sun Microsystems, Inc. business.
  * Permission to use, copy, modify, and distribute this
- * software is freely granted, provided that this notice 
+ * software is freely granted, provided that this notice
  * is preserved.
  * ====================================================
@@ -13,4 +13,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_atan2.c,v 1.10 2005/02/04 18:26:05 das Exp $";
@@ -20,5 +21,5 @@
  * Method :
  *	1. Reduce y to positive by atan2(y,x)=-atan2(-y,x).
- *	2. Reduce x to positive by (if x and y are unexceptional): 
+ *	2. Reduce x to positive by (if x and y are unexceptional):
  *		ARG (x+iy) = arctan(y/x)   	   ... if x > 0,
  *		ARG (x+iy) = pi - arctan[y/(-x)]   ... if x < 0,
@@ -38,7 +39,7 @@
  *
  * Constants:
- * The hexadecimal values are the intended ones for the following 
- * constants. The decimal values may be used, provided that the 
- * compiler will convert from decimal to binary accurately enough 
+ * The hexadecimal values are the intended ones for the following
+ * constants. The decimal values may be used, provided that the
+ * compiler will convert from decimal to binary accurately enough
  * to produce the hexadecimal values shown.
  */
@@ -75,5 +76,5 @@
 	if((iy|ly)==0) {
 	    switch(m) {
-		case 0: 
+		case 0:
 		case 1: return y; 	/* atan(+-0,+anything)=+-0 */
 		case 2: return  pi+tiny;/* atan(+0,-anything) = pi */
@@ -83,5 +84,5 @@
     /* when x = 0 */
 	if((ix|lx)==0) return (hy<0)?  -pi_o_2-tiny: pi_o_2+tiny;
-	    
+	
     /* when x is INF */
 	if(ix==0x7ff00000) {
Index: /trunk/src/emx/src/lib/msun/src/e_atan2f.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_atan2f.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_atan2f.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_atan2f.c,v 1.7 2004/06/02 17:09:05 bde Exp $";
Index: /trunk/src/emx/src/lib/msun/src/e_atanh.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_atanh.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_atanh.c	(revision 2013)
@@ -7,5 +7,5 @@
  * Developed at SunSoft, a Sun Microsystems, Inc. business.
  * Permission to use, copy, modify, and distribute this
- * software is freely granted, provided that this notice 
+ * software is freely granted, provided that this notice
  * is preserved.
  * ====================================================
@@ -13,4 +13,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_atanh.c,v 1.7 2005/02/04 18:26:05 das Exp $";
@@ -51,5 +52,5 @@
 	if ((ix|((lx|(-lx))>>31))>0x3ff00000) /* |x|>1 */
 	    return (x-x)/(x-x);
-	if(ix==0x3ff00000) 
+	if(ix==0x3ff00000)
 	    return x/zero;
 	if(ix<0x3e300000&&(huge+x)>zero) return x;	/* x<2**-28 */
@@ -58,5 +59,5 @@
 	    t = x+x;
 	    t = 0.5*log1p(t+t*x/(one-x));
-	} else 
+	} else
 	    t = 0.5*log1p((x+x)/(one-x));
 	if(hx>=0) return t; else return -t;
Index: /trunk/src/emx/src/lib/msun/src/e_atanhf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_atanhf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_atanhf.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_atanhf.c,v 1.6 2002/05/28 17:03:12 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/e_cosh.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_cosh.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_cosh.c	(revision 2013)
@@ -7,9 +7,10 @@
  * Developed at SunSoft, a Sun Microsystems, Inc. business.
  * Permission to use, copy, modify, and distribute this
- * software is freely granted, provided that this notice 
+ * software is freely granted, provided that this notice
  * is preserved.
  * ====================================================
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_cosh.c,v 1.8 2005/02/04 18:26:05 das Exp $";
@@ -17,9 +18,9 @@
 
 /* __ieee754_cosh(x)
- * Method : 
+ * Method :
  * mathematically cosh(x) if defined to be (exp(x)+exp(-x))/2
- *	1. Replace x by |x| (cosh(x) = cosh(-x)). 
- *	2. 
- *		                                        [ exp(x) - 1 ]^2 
+ *	1. Replace x by |x| (cosh(x) = cosh(-x)).
+ *	2.
+ *		                                        [ exp(x) - 1 ]^2
  *	    0        <= x <= ln2/2  :  cosh(x) := 1 + -------------------
  *			       			           2*exp(x)
@@ -28,5 +29,5 @@
  *	    ln2/2    <= x <= 22     :  cosh(x) := -------------------
  *			       			          2
- *	    22       <= x <= lnovft :  cosh(x) := exp(x)/2 
+ *	    22       <= x <= lnovft :  cosh(x) := exp(x)/2
  *	    lnovft   <= x <= ln2ovft:  cosh(x) := exp(x/2)/2 * exp(x/2)
  *	    ln2ovft  <  x	    :  cosh(x) := huge*huge (overflow)
Index: /trunk/src/emx/src/lib/msun/src/e_coshf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_coshf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_coshf.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_coshf.c,v 1.6 2002/05/28 17:03:12 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/e_exp.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_exp.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_exp.c	(revision 2013)
@@ -6,9 +6,10 @@
  *
  * Permission to use, copy, modify, and distribute this
- * software is freely granted, provided that this notice 
+ * software is freely granted, provided that this notice
  * is preserved.
  * ====================================================
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_exp.c,v 1.10 2005/02/04 18:26:05 das Exp $";
@@ -23,7 +24,7 @@
  *	Given x, find r and integer k such that
  *
- *               x = k*ln2 + r,  |r| <= 0.5*ln2.  
+ *               x = k*ln2 + r,  |r| <= 0.5*ln2.
  *
- *      Here r will be represented as r = hi-lo for better 
+ *      Here r will be represented as r = hi-lo for better
  *	accuracy.
  *
@@ -32,6 +33,6 @@
  *	Write
  *	    R(r**2) = r*(exp(r)+1)/(exp(r)-1) = 2 + r*r/6 - r**4/360 + ...
- *      We use a special Remes algorithm on [0,0.34658] to generate 
- * 	a polynomial of degree 5 to approximate R. The maximum error 
+ *      We use a special Remes algorithm on [0,0.34658] to generate
+ * 	a polynomial of degree 5 to approximate R. The maximum error
  *	of this polynomial approximation is bounded by 2**-59. In
  *	other words,
@@ -40,5 +41,5 @@
  *	and
  *	    |                  5          |     -59
- *	    | 2.0+P1*z+...+P5*z   -  R(z) | <= 2 
+ *	    | 2.0+P1*z+...+P5*z   -  R(z) | <= 2
  *	    |                             |
  *	The computation of exp(r) thus becomes
@@ -67,11 +68,11 @@
  *
  * Misc. info.
- *	For IEEE double 
+ *	For IEEE double
  *	    if x >  7.09782712893383973096e+02 then exp(x) overflow
  *	    if x < -7.45133219101941108420e+02 then exp(x) underflow
  *
  * Constants:
- * The hexadecimal values are the intended ones for the following 
- * constants. The decimal values may be used, provided that the 
+ * The hexadecimal values are the intended ones for the following
+ * constants. The decimal values may be used, provided that the
  * compiler will convert from decimal to binary accurately enough
  * to produce the hexadecimal values shown.
@@ -125,5 +126,5 @@
 
     /* argument reduction */
-	if(hx > 0x3fd62e42) {		/* if  |x| > 0.5 ln2 */ 
+	if(hx > 0x3fd62e42) {		/* if  |x| > 0.5 ln2 */
 	    if(hx < 0x3FF0A2B2) {	/* and |x| < 1.5 ln2 */
 		hi = x-ln2HI[xsb]; lo=ln2LO[xsb]; k = 1-xsb-xsb;
@@ -135,5 +136,5 @@
 	    }
 	    x  = hi - lo;
-	} 
+	}
 	else if(hx < 0x3e300000)  {	/* when |x|<2**-28 */
 	    if(huge+x>one) return one+x;/* trigger inexact */
@@ -144,5 +145,5 @@
 	t  = x*x;
 	c  = x - t*(P1+t*(P2+t*(P3+t*(P4+t*P5))));
-	if(k==0) 	return one-((x*c)/(c-2.0)-x); 
+	if(k==0) 	return one-((x*c)/(c-2.0)-x);
 	else 		y = one-((lo-(x*c)/(2.0-c))-hi);
 	if(k >= -1021) {
Index: /trunk/src/emx/src/lib/msun/src/e_expf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_expf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_expf.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_expf.c,v 1.9 2005/02/24 06:32:13 das Exp $";
Index: /trunk/src/emx/src/lib/msun/src/e_fmod.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_fmod.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_fmod.c	(revision 2013)
@@ -7,14 +7,15 @@
  * Developed at SunSoft, a Sun Microsystems, Inc. business.
  * Permission to use, copy, modify, and distribute this
- * software is freely granted, provided that this notice 
+ * software is freely granted, provided that this notice
  * is preserved.
  * ====================================================
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_fmod.c,v 1.9 2005/02/04 18:26:05 das Exp $";
 #endif
 
-/* 
+/*
  * __ieee754_fmod(x,y)
  * Return x mod y in exact arithmetic
@@ -45,5 +46,5 @@
 	if(hx<=hy) {
 	    if((hx<hy)||(lx<ly)) return x;	/* |x|<|y| return x */
-	    if(lx==ly) 
+	    if(lx==ly)
 		return Zero[(u_int32_t)sx>>31];	/* |x|=|y| return x*0*/
 	}
@@ -68,5 +69,5 @@
 
     /* set up {hx,lx}, {hy,ly} and align y to x */
-	if(ix >= -1022) 
+	if(ix >= -1022)
 	    hx = 0x00100000|(0x000fffff&hx);
 	else {		/* subnormal x, shift x to normal */
@@ -80,5 +81,5 @@
 	    }
 	}
-	if(iy >= -1022) 
+	if(iy >= -1022)
 	    hy = 0x00100000|(0x000fffff&hy);
 	else {		/* subnormal y, shift y to normal */
Index: /trunk/src/emx/src/lib/msun/src/e_fmodf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_fmodf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_fmodf.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_fmodf.c,v 1.6 2002/05/28 17:03:12 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/e_gamma.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_gamma.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_gamma.c	(revision 2013)
@@ -7,5 +7,5 @@
  * Developed at SunSoft, a Sun Microsystems, Inc. business.
  * Permission to use, copy, modify, and distribute this
- * software is freely granted, provided that this notice 
+ * software is freely granted, provided that this notice
  * is preserved.
  * ====================================================
@@ -13,4 +13,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_gamma.c,v 1.7 2005/02/04 18:26:05 das Exp $";
Index: /trunk/src/emx/src/lib/msun/src/e_gamma_r.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_gamma_r.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_gamma_r.c	(revision 2013)
@@ -7,5 +7,5 @@
  * Developed at SunSoft, a Sun Microsystems, Inc. business.
  * Permission to use, copy, modify, and distribute this
- * software is freely granted, provided that this notice 
+ * software is freely granted, provided that this notice
  * is preserved.
  * ====================================================
@@ -13,4 +13,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_gamma_r.c,v 1.7 2005/02/04 18:26:05 das Exp $";
@@ -18,6 +19,6 @@
 
 /* __ieee754_gamma_r(x, signgamp)
- * Reentrant version of the logarithm of the Gamma function 
- * with user provide pointer for the sign of Gamma(x). 
+ * Reentrant version of the logarithm of the Gamma function
+ * with user provide pointer for the sign of Gamma(x).
  *
  * Method: See __ieee754_lgamma_r
Index: /trunk/src/emx/src/lib/msun/src/e_gammaf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_gammaf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_gammaf.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_gammaf.c,v 1.6 2002/05/28 17:03:12 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/e_gammaf_r.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_gammaf_r.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_gammaf_r.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_gammaf_r.c,v 1.7 2002/05/28 18:15:03 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/e_hypot.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_hypot.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_hypot.c	(revision 2013)
@@ -7,9 +7,10 @@
  * Developed at SunSoft, a Sun Microsystems, Inc. business.
  * Permission to use, copy, modify, and distribute this
- * software is freely granted, provided that this notice 
+ * software is freely granted, provided that this notice
  * is preserved.
  * ====================================================
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_hypot.c,v 1.9 2005/02/04 18:26:05 das Exp $";
@@ -18,10 +19,10 @@
 /* __ieee754_hypot(x,y)
  *
- * Method :                  
- *	If (assume round-to-nearest) z=x*x+y*y 
- *	has error less than sqrt(2)/2 ulp, than 
+ * Method :
+ *	If (assume round-to-nearest) z=x*x+y*y
+ *	has error less than sqrt(2)/2 ulp, than
  *	sqrt(z) has error less than 1 ulp (exercise).
  *
- *	So, compute sqrt(x*x+y*y) with some care as 
+ *	So, compute sqrt(x*x+y*y) with some care as
  *	follows to get the error below 1 ulp:
  *
@@ -33,8 +34,8 @@
  *	2. if x <= 2y use
  *		t1*y1+((x-y)*(x-y)+(t1*y2+t2*y))
- *	where t1 = 2x with lower 32 bits cleared, t2 = 2x-t1, 
+ *	where t1 = 2x with lower 32 bits cleared, t2 = 2x-t1,
  *	y1= y with lower 32 bits chopped, y2 = y-y1.
  *		
- *	NOTE: scaling may be necessary if some argument is too 
+ *	NOTE: scaling may be necessary if some argument is too
  *	      large or too tiny
  *
@@ -44,6 +45,6 @@
  *
  * Accuracy:
- * 	hypot(x,y) returns sqrt(x^2+y^2) with error less 
- * 	than 1 ulps (units in the last place) 
+ * 	hypot(x,y) returns sqrt(x^2+y^2) with error less
+ * 	than 1 ulps (units in the last place)
  */
 
Index: /trunk/src/emx/src/lib/msun/src/e_hypotf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_hypotf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_hypotf.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_hypotf.c,v 1.9 2002/05/28 18:15:03 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/e_j0.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_j0.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_j0.c	(revision 2013)
@@ -7,9 +7,10 @@
  * Developed at SunSoft, a Sun Microsystems, Inc. business.
  * Permission to use, copy, modify, and distribute this
- * software is freely granted, provided that this notice 
+ * software is freely granted, provided that this notice
  * is preserved.
  * ====================================================
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_j0.c,v 1.8 2005/02/04 18:26:05 das Exp $";
@@ -35,5 +36,5 @@
  *		sin(x) +- cos(x) = -cos(2x)/(sin(x) -+ cos(x))
  * 	    to compute the worse one.)
- *	   
+ *	
  *	3 Special cases
  *		j0(nan)= nan
@@ -43,10 +44,10 @@
  * Method -- y0(x):
  *	1. For x<2.
- *	   Since 
+ *	   Since
  *		y0(x) = 2/pi*(j0(x)*(ln(x/2)+Euler) + x^2/4 - ...)
  *	   therefore y0(x)-2/pi*j0(x)*ln(x) is an even function.
  *	   We use the following function to approximate y0,
  *		y0(x) = U(z)/V(z) + (2/pi)*(j0(x)*ln(x)), z= x^2
- *	   where 
+ *	   where
  *		U(z) = u00 + u01*z + ... + u06*z^6
  *		V(z) = 1  + v01*z + ... + v04*z^4
@@ -153,5 +154,5 @@
         ix = 0x7fffffff&hx;
     /* Y0(NaN) is NaN, y0(-inf) is Nan, y0(inf) is 0  */
-	if(ix>=0x7ff00000) return  one/(x+x*x); 
+	if(ix>=0x7ff00000) return  one/(x+x*x);
         if((ix|lx)==0) return -one/zero;
         if(hx<0) return zero/zero;
Index: /trunk/src/emx/src/lib/msun/src/e_j0f.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_j0f.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_j0f.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_j0f.c,v 1.7 2002/05/28 18:15:03 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/e_j1.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_j1.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_j1.c	(revision 2013)
@@ -7,9 +7,10 @@
  * Developed at SunSoft, a Sun Microsystems, Inc. business.
  * Permission to use, copy, modify, and distribute this
- * software is freely granted, provided that this notice 
+ * software is freely granted, provided that this notice
  * is preserved.
  * ====================================================
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_j1.c,v 1.8 2005/02/04 18:26:06 das Exp $";
@@ -36,5 +37,5 @@
  *		sin(x) +- cos(x) = -cos(2x)/(sin(x) -+ cos(x))
  * 	    to compute the worse one.)
- *	   
+ *	
  *	3 Special cases
  *		j1(nan)= nan
@@ -43,7 +44,7 @@
  *		
  * Method -- y1(x):
- *	1. screen out x<=0 cases: y1(0)=-inf, y1(x<0)=NaN 
+ *	1. screen out x<=0 cases: y1(0)=-inf, y1(x<0)=NaN
  *	2. For x<2.
- *	   Since 
+ *	   Since
  *		y1(x) = 2/pi*(j1(x)*(ln(x/2)+Euler)-1/x-x/2+5/64*x^3-...)
  *	   therefore y1(x)-2/pi*j1(x)*ln(x)-1/x is an odd function.
@@ -150,5 +151,5 @@
         ix = 0x7fffffff&hx;
     /* if Y1(NaN) is NaN, Y1(-inf) is NaN, Y1(inf) is 0 */
-	if(ix>=0x7ff00000) return  one/(x+x*x); 
+	if(ix>=0x7ff00000) return  one/(x+x*x);
         if((ix|lx)==0) return -one/zero;
         if(hx<0) return zero/zero;
@@ -180,8 +181,8 @@
                 }
                 return z;
-        } 
+        }
         if(ix<=0x3c900000) {    /* x < 2**-54 */
             return(-tpi/x);
-        } 
+        }
         z = x*x;
         u = U0[0]+z*(U0[1]+z*(U0[2]+z*(U0[3]+z*U0[4])));
Index: /trunk/src/emx/src/lib/msun/src/e_j1f.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_j1f.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_j1f.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_j1f.c,v 1.7 2002/05/28 18:15:04 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/e_jn.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_jn.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_jn.c	(revision 2013)
@@ -7,9 +7,10 @@
  * Developed at SunSoft, a Sun Microsystems, Inc. business.
  * Permission to use, copy, modify, and distribute this
- * software is freely granted, provided that this notice 
+ * software is freely granted, provided that this notice
  * is preserved.
  * ====================================================
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_jn.c,v 1.9 2005/02/04 18:26:06 das Exp $";
@@ -20,5 +21,5 @@
  * floating point Bessel's function of the 1st and 2nd kind
  * of order n
- *          
+ *
  * Special cases:
  *	y0(0)=y1(0)=yn(n,0) = -inf with division by zero signal;
@@ -77,11 +78,11 @@
 	if((ix|lx)==0||ix>=0x7ff00000) 	/* if x is 0 or inf */
 	    b = zero;
-	else if((double)n<=x) {   
+	else if((double)n<=x) {
 		/* Safe to use J(n+1,x)=2n/x *J(n,x)-J(n-1,x) */
 	    if(ix>=0x52D00000) { /* x > 2**302 */
-    /* (x >> n**2) 
+    /* (x >> n**2)
      *	    Jn(x) = cos(x-(2n+1)*pi/4)*sqrt(2/x*pi)
      *	    Yn(x) = sin(x-(2n+1)*pi/4)*sqrt(2/x*pi)
-     *	    Let s=sin(x), c=cos(x), 
+     *	    Let s=sin(x), c=cos(x),
      *		xn=x-(2n+1)*pi/4, sqt2 = sqrt(2),then
      *
@@ -111,5 +112,5 @@
 	} else {
 	    if(ix<0x3e100000) {	/* x < 2**-29 */
-    /* x is tiny, return the first Taylor expansion of J(n,x) 
+    /* x is tiny, return the first Taylor expansion of J(n,x)
      * J(n,x) = 1/n!*(x/2)^n  - ...
      */
@@ -126,12 +127,12 @@
 	    } else {
 		/* use backward recurrence */
-		/* 			x      x^2      x^2       
+		/* 			x      x^2      x^2
 		 *  J(n,x)/J(n-1,x) =  ----   ------   ------   .....
 		 *			2n  - 2(n+1) - 2(n+2)
 		 *
-		 * 			1      1        1       
+		 * 			1      1        1
 		 *  (for large x)   =  ----  ------   ------   .....
 		 *			2n   2(n+1)   2(n+2)
-		 *			-- - ------ - ------ - 
+		 *			-- - ------ - ------ -
 		 *			 x     x         x
 		 *
@@ -149,7 +150,7 @@
 		 * Q(0) = w, Q(1) = w(w+h) - 1,
 		 * Q(k) = (w+k*h)*Q(k-1) - Q(k-2),
-		 * When Q(k) > 1e4	good for single 
-		 * When Q(k) > 1e9	good for double 
-		 * When Q(k) > 1e17	good for quadruple 
+		 * When Q(k) > 1e4	good for single
+		 * When Q(k) > 1e9	good for double
+		 * When Q(k) > 1e17	good for quadruple
 		 */
 	    /* determine k */
@@ -230,8 +231,8 @@
 	if(ix==0x7ff00000) return zero;
 	if(ix>=0x52D00000) { /* x > 2**302 */
-    /* (x >> n**2) 
+    /* (x >> n**2)
      *	    Jn(x) = cos(x-(2n+1)*pi/4)*sqrt(2/x*pi)
      *	    Yn(x) = sin(x-(2n+1)*pi/4)*sqrt(2/x*pi)
-     *	    Let s=sin(x), c=cos(x), 
+     *	    Let s=sin(x), c=cos(x),
      *		xn=x-(2n+1)*pi/4, sqt2 = sqrt(2),then
      *
Index: /trunk/src/emx/src/lib/msun/src/e_jnf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_jnf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_jnf.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_jnf.c,v 1.8 2002/05/28 18:15:04 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/e_lgamma.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_lgamma.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_lgamma.c	(revision 2013)
@@ -7,5 +7,5 @@
  * Developed at SunSoft, a Sun Microsystems, Inc. business.
  * Permission to use, copy, modify, and distribute this
- * software is freely granted, provided that this notice 
+ * software is freely granted, provided that this notice
  * is preserved.
  * ====================================================
@@ -13,4 +13,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_lgamma.c,v 1.8 2005/02/04 18:26:06 das Exp $";
Index: /trunk/src/emx/src/lib/msun/src/e_lgamma_r.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_lgamma_r.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_lgamma_r.c	(revision 2013)
@@ -7,5 +7,5 @@
  * Developed at SunSoft, a Sun Microsystems, Inc. business.
  * Permission to use, copy, modify, and distribute this
- * software is freely granted, provided that this notice 
+ * software is freely granted, provided that this notice
  * is preserved.
  * ====================================================
@@ -13,4 +13,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_lgamma_r.c,v 1.8 2005/02/04 18:26:06 das Exp $";
@@ -18,10 +19,10 @@
 
 /* __ieee754_lgamma_r(x, signgamp)
- * Reentrant version of the logarithm of the Gamma function 
- * with user provide pointer for the sign of Gamma(x). 
+ * Reentrant version of the logarithm of the Gamma function
+ * with user provide pointer for the sign of Gamma(x).
  *
  * Method:
  *   1. Argument Reduction for 0 < x <= 8
- * 	Since gamma(1+s)=s*gamma(s), for x in [0,8], we may 
+ * 	Since gamma(1+s)=s*gamma(s), for x in [0,8], we may
  * 	reduce x to a number in [1.5,2.5] by
  * 		lgamma(1+s) = log(s) + lgamma(s)
@@ -61,5 +62,5 @@
  *	  			    3       5             11
  *		w = w0 + w1*z + w2*z  + w3*z  + ... + w6*z
- *	where 
+ *	where
  *		|w - f(z)| < 2**-58.74
  *		
@@ -69,8 +70,8 @@
  * 		G(x) = pi/(sin(pi*x)*(-x)*G(-x))
  *	since G(-x) is positive, sign(G(x)) = sign(sin(pi*x)) for x<0
- *	Hence, for x<0, signgam = sign(sin(pi*x)) and 
+ *	Hence, for x<0, signgam = sign(sin(pi*x)) and
  *		lgamma(x) = log(|Gamma(x)|)
  *			  = log(pi/(|x*sin(pi*x)|)) - lgamma(-x);
- *	Note: one should avoid compute pi*(-x) directly in the 
+ *	Note: one should avoid compute pi*(-x) directly in the
  *	      computation of sin(pi*(-x)).
  *		
@@ -190,7 +191,7 @@
 	switch (n) {
 	    case 0:   y =  __kernel_sin(pi*y,zero,0); break;
-	    case 1:   
+	    case 1:
 	    case 2:   y =  __kernel_cos(pi*(0.5-y),zero); break;
-	    case 3:  
+	    case 3:
 	    case 4:   y =  __kernel_sin(pi*(one-y),zero,0); break;
 	    case 5:
@@ -290,5 +291,5 @@
 	    w = w0+z*(w1+y*(w2+y*(w3+y*(w4+y*(w5+y*w6)))));
 	    r = (x-half)*(t-one)+w;
-	} else 
+	} else
     /* 2**58 <= x <= inf */
 	    r =  x*(__ieee754_log(x)-one);
Index: /trunk/src/emx/src/lib/msun/src/e_lgammaf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_lgammaf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_lgammaf.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_lgammaf.c,v 1.7 2002/05/28 18:15:04 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/e_lgammaf_r.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_lgammaf_r.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_lgammaf_r.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_lgammaf_r.c,v 1.7 2002/05/28 18:15:04 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/e_log.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_log.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_log.c	(revision 2013)
@@ -7,9 +7,10 @@
  * Developed at SunSoft, a Sun Microsystems, Inc. business.
  * Permission to use, copy, modify, and distribute this
- * software is freely granted, provided that this notice 
+ * software is freely granted, provided that this notice
  * is preserved.
  * ====================================================
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_log.c,v 1.10 2005/02/04 18:26:06 das Exp $";
@@ -19,7 +20,7 @@
  * Return the logrithm of x
  *
- * Method :                  
- *   1. Argument Reduction: find k and f such that 
- *			x = 2^k * (1+f), 
+ * Method :
+ *   1. Argument Reduction: find k and f such that
+ *			x = 2^k * (1+f),
  *	   where  sqrt(2)/2 < 1+f < sqrt(2) .
  *
@@ -28,6 +29,6 @@
  *		 = 2s + 2/3 s**3 + 2/5 s**5 + .....,
  *	     	 = 2s + s*R
- *      We use a special Reme algorithm on [0,0.1716] to generate 
- * 	a polynomial of degree 14 to approximate R The maximum error 
+ *      We use a special Reme algorithm on [0,0.1716] to generate
+ * 	a polynomial of degree 14 to approximate R The maximum error
  *	of this polynomial approximation is bounded by 2**-58.45. In
  *	other words,
@@ -37,5 +38,5 @@
  *	and
  *	    |      2          14          |     -58.45
- *	    | Lg1*s +...+Lg7*s    -  R(z) | <= 2 
+ *	    | Lg1*s +...+Lg7*s    -  R(z) | <= 2
  *	    |                             |
  *	Note that 2s = f - s*f = f - hfsq + s*hfsq, where hfsq = f*f/2.
@@ -45,12 +46,12 @@
  *		log(1+f) = f - (hfsq - s*(hfsq+R)).	(better accuracy)
  *	
- *	3. Finally,  log(x) = k*ln2 + log(1+f).  
+ *	3. Finally,  log(x) = k*ln2 + log(1+f).
  *			    = k*ln2_hi+(f-(hfsq-(s*(hfsq+R)+k*ln2_lo)))
- *	   Here ln2 is split into two floating point number: 
+ *	   Here ln2 is split into two floating point number:
  *			ln2_hi + ln2_lo,
  *	   where n*ln2_hi is always exact for |n| < 2000.
  *
  * Special cases:
- *	log(x) is NaN with signal if x < 0 (including -INF) ; 
+ *	log(x) is NaN with signal if x < 0 (including -INF) ;
  *	log(+INF) is +INF; log(0) is -INF with signal;
  *	log(NaN) is that NaN with no signal.
@@ -61,7 +62,7 @@
  *
  * Constants:
- * The hexadecimal values are the intended ones for the following 
- * constants. The decimal values may be used, provided that the 
- * compiler will convert from decimal to binary accurately enough 
+ * The hexadecimal values are the intended ones for the following
+ * constants. The decimal values may be used, provided that the
+ * compiler will convert from decimal to binary accurately enough
  * to produce the hexadecimal values shown.
  */
@@ -95,10 +96,10 @@
 	k=0;
 	if (hx < 0x00100000) {			/* x < 2**-1022  */
-	    if (((hx&0x7fffffff)|lx)==0) 
+	    if (((hx&0x7fffffff)|lx)==0)
 		return -two54/zero;		/* log(+-0)=-inf */
 	    if (hx<0) return (x-x)/zero;	/* log(-#) = NaN */
 	    k -= 54; x *= two54; /* subnormal number, scale up x */
 	    GET_HIGH_WORD(hx,x);
-	} 
+	}
 	if (hx >= 0x7ff00000) return x+x;
 	k += (hx>>20)-1023;
@@ -115,5 +116,5 @@
 	    	     return dk*ln2_hi-((R-dk*ln2_lo)-f);}
 	}
- 	s = f/(2.0+f); 
+ 	s = f/(2.0+f);
 	dk = (double)k;
 	z = s*s;
@@ -121,6 +122,6 @@
 	w = z*z;
 	j = 0x6b851-hx;
-	t1= w*(Lg2+w*(Lg4+w*Lg6)); 
-	t2= z*(Lg1+w*(Lg3+w*(Lg5+w*Lg7))); 
+	t1= w*(Lg2+w*(Lg4+w*Lg6));
+	t2= z*(Lg1+w*(Lg3+w*(Lg5+w*Lg7)));
 	i |= j;
 	R = t2+t1;
Index: /trunk/src/emx/src/lib/msun/src/e_log10.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_log10.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_log10.c	(revision 2013)
@@ -7,9 +7,10 @@
  * Developed at SunSoft, a Sun Microsystems, Inc. business.
  * Permission to use, copy, modify, and distribute this
- * software is freely granted, provided that this notice 
+ * software is freely granted, provided that this notice
  * is preserved.
  * ====================================================
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_log10.c,v 1.11 2005/02/04 18:26:06 das Exp $";
@@ -18,5 +19,5 @@
 /* __ieee754_log10(x)
  * Return the base 10 logarithm of x
- * 
+ *
  * Method :
  *	Let log10_2hi = leading 40 bits of log10(2) and
@@ -24,5 +25,5 @@
  *	    ivln10   = 1/log(10) rounded.
  *	Then
- *		n = ilogb(x), 
+ *		n = ilogb(x),
  *		if(n<0)  n = n+1;
  *		x = scalbn(x,-n);
@@ -30,5 +31,5 @@
  *
  * Note 1:
- *	To guarantee log10(10**n)=n, where 10**n is normal, the rounding 
+ *	To guarantee log10(10**n)=n, where 10**n is normal, the rounding
  *	mode must set to Round-to-Nearest.
  * Note 2:
@@ -37,5 +38,5 @@
  *
  * Special cases:
- *	log10(x) is NaN with signal if x < 0; 
+ *	log10(x) is NaN with signal if x < 0;
  *	log10(+INF) is +INF with no signal; log10(0) is -INF with signal;
  *	log10(NaN) is that NaN with no signal;
Index: /trunk/src/emx/src/lib/msun/src/e_log10f.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_log10f.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_log10f.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_log10f.c,v 1.8 2002/05/28 18:15:04 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/e_logf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_logf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_logf.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_logf.c,v 1.7 2002/05/28 18:15:04 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/e_pow.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_pow.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_pow.c	(revision 2013)
@@ -5,9 +5,10 @@
  *
  * Permission to use, copy, modify, and distribute this
- * software is freely granted, provided that this notice 
+ * software is freely granted, provided that this notice
  * is preserved.
  * ====================================================
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_pow.c,v 1.11 2005/02/04 18:26:06 das Exp $";
@@ -21,5 +22,5 @@
  *		log2(x) = w1 + w2,
  *	   where w1 has 53-24 = 29 bit trailing zeros.
- *	2. Perform y*log2(x) = n+y' by simulating muti-precision 
+ *	2. Perform y*log2(x) = n+y' by simulating muti-precision
  *	   arithmetic, where |y'|<=0.5.
  *	3. Return x**y = 2**n*exp(y'*log2)
@@ -49,11 +50,11 @@
  *	pow(x,y) returns x**y nearly rounded. In particular
  *			pow(integer,integer)
- *	always returns the correct integer provided it is 
+ *	always returns the correct integer provided it is
  *	representable.
  *
  * Constants :
- * The hexadecimal values are the intended ones for the following 
- * constants. The decimal values may be used, provided that the 
- * compiler will convert from decimal to binary accurately enough 
+ * The hexadecimal values are the intended ones for the following
+ * constants. The decimal values may be used, provided that the
+ * compiler will convert from decimal to binary accurately enough
  * to produce the hexadecimal values shown.
  */
@@ -113,5 +114,5 @@
     /* +-NaN return x+y */
 	if(ix > 0x7ff00000 || ((ix==0x7ff00000)&&(lx!=0)) ||
-	   iy > 0x7ff00000 || ((iy==0x7ff00000)&&(ly!=0))) 
+	   iy > 0x7ff00000 || ((iy==0x7ff00000)&&(ly!=0)))
 		return x+y;	
 
@@ -134,5 +135,5 @@
 		}
 	    }		
-	} 
+	}
 
     /* special value of y */
@@ -145,5 +146,5 @@
 	        else			/* (|x|<1)**-,+inf = inf,0 */
 		    return (hy<0)?-y: zero;
-	    } 
+	    }
 	    if(iy==0x3ff00000) {	/* y is  +-1 */
 		if(hy<0) return one/x; else return x;
@@ -165,5 +166,5 @@
 		    if(((ix-0x3ff00000)|yisint)==0) {
 			z = (z-z)/(z-z); /* (-1)**non-int is NaN */
-		    } else if(yisint==1) 
+		    } else if(yisint==1)
 			z = -z;		/* (x<0)**odd = -(|x|**odd) */
 		}
@@ -171,5 +172,5 @@
 	    }
 	}
-    
+
     /* CYGNUS LOCAL + fdlibm-5.3 fix: This used to be
 	n = (hx>>31)+1;
@@ -193,5 +194,5 @@
 	    if(ix<0x3fefffff) return (hy<0)? s*huge*huge:s*tiny*tiny;
 	    if(ix>0x3ff00000) return (hy>0)? s*huge*huge:s*tiny*tiny;
-	/* now |1-x| is tiny <= 2**-20, suffice to compute 
+	/* now |1-x| is tiny <= 2**-20, suffice to compute
 	   log(x) by x-x^2/2+x^3/3-x^4/4 */
 	    t = ax-one;		/* t has 20 trailing zeros */
@@ -286,5 +287,5 @@
 	    if(j<0) n = -n;
 	    p_h -= t;
-	} 
+	}
 	t = p_l+p_h;
 	SET_LOW_WORD(t,0);
Index: /trunk/src/emx/src/lib/msun/src/e_powf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_powf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_powf.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_powf.c,v 1.12 2004/06/01 19:33:30 bde Exp $";
Index: /trunk/src/emx/src/lib/msun/src/e_rem_pio2.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_rem_pio2.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_rem_pio2.c	(revision 2013)
@@ -7,5 +7,5 @@
  * Developed at SunSoft, a Sun Microsystems, Inc. business.
  * Permission to use, copy, modify, and distribute this
- * software is freely granted, provided that this notice 
+ * software is freely granted, provided that this notice
  * is preserved.
  * ====================================================
@@ -13,4 +13,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_rem_pio2.c,v 1.8 2005/02/04 18:26:06 das Exp $";
@@ -18,6 +19,6 @@
 
 /* __ieee754_rem_pio2(x,y)
- * 
- * return the remainder of x rem pi/2 in y[0]+y[1] 
+ *
+ * return the remainder of x rem pi/2 in y[0]+y[1]
  * use __kernel_rem_pio2()
  */
@@ -27,18 +28,18 @@
 
 /*
- * Table of constants for 2/pi, 396 Hex digits (476 decimal) of 2/pi 
+ * Table of constants for 2/pi, 396 Hex digits (476 decimal) of 2/pi
  */
 static const int32_t two_over_pi[] = {
-0xA2F983, 0x6E4E44, 0x1529FC, 0x2757D1, 0xF534DD, 0xC0DB62, 
-0x95993C, 0x439041, 0xFE5163, 0xABDEBB, 0xC561B7, 0x246E3A, 
-0x424DD2, 0xE00649, 0x2EEA09, 0xD1921C, 0xFE1DEB, 0x1CB129, 
-0xA73EE8, 0x8235F5, 0x2EBB44, 0x84E99C, 0x7026B4, 0x5F7E41, 
-0x3991D6, 0x398353, 0x39F49C, 0x845F8B, 0xBDF928, 0x3B1FF8, 
-0x97FFDE, 0x05980F, 0xEF2F11, 0x8B5A0A, 0x6D1F6D, 0x367ECF, 
-0x27CB09, 0xB74F46, 0x3F669E, 0x5FEA2D, 0x7527BA, 0xC7EBE5, 
-0xF17B3D, 0x0739F7, 0x8A5292, 0xEA6BFB, 0x5FB11F, 0x8D5D08, 
-0x560330, 0x46FC7B, 0x6BABF0, 0xCFBC20, 0x9AF436, 0x1DA9E3, 
-0x91615E, 0xE61B08, 0x659985, 0x5F14A0, 0x68408D, 0xFFD880, 
-0x4D7327, 0x310606, 0x1556CA, 0x73A8C9, 0x60E27B, 0xC08C6B, 
+0xA2F983, 0x6E4E44, 0x1529FC, 0x2757D1, 0xF534DD, 0xC0DB62,
+0x95993C, 0x439041, 0xFE5163, 0xABDEBB, 0xC561B7, 0x246E3A,
+0x424DD2, 0xE00649, 0x2EEA09, 0xD1921C, 0xFE1DEB, 0x1CB129,
+0xA73EE8, 0x8235F5, 0x2EBB44, 0x84E99C, 0x7026B4, 0x5F7E41,
+0x3991D6, 0x398353, 0x39F49C, 0x845F8B, 0xBDF928, 0x3B1FF8,
+0x97FFDE, 0x05980F, 0xEF2F11, 0x8B5A0A, 0x6D1F6D, 0x367ECF,
+0x27CB09, 0xB74F46, 0x3F669E, 0x5FEA2D, 0x7527BA, 0xC7EBE5,
+0xF17B3D, 0x0739F7, 0x8A5292, 0xEA6BFB, 0x5FB11F, 0x8D5D08,
+0x560330, 0x46FC7B, 0x6BABF0, 0xCFBC20, 0x9AF436, 0x1DA9E3,
+0x91615E, 0xE61B08, 0x659985, 0x5F14A0, 0x68408D, 0xFFD880,
+0x4D7327, 0x310606, 0x1556CA, 0x73A8C9, 0x60E27B, 0xC08C6B,
 };
 
@@ -86,5 +87,5 @@
 	    {y[0] = x; y[1] = 0; return 0;}
 	if(ix<0x4002d97c) {  /* |x| < 3pi/4, special case with n=+-1 */
-	    if(hx>0) { 
+	    if(hx>0) {
 		z = x - pio2_1;
 		if(ix!=0x3ff921fb) { 	/* 33+53 bit pi is good enough */
@@ -121,5 +122,5 @@
 	        u_int32_t high;
 	        j  = ix>>20;
-	        y[0] = r-w; 
+	        y[0] = r-w;
 		GET_HIGH_WORD(high,y[0]);
 	        i = j-((high>>20)&0x7ff);
@@ -145,5 +146,5 @@
 	    else	 return n;
 	}
-    /* 
+    /*
      * all other (large) arguments
      */
Index: /trunk/src/emx/src/lib/msun/src/e_rem_pio2f.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_rem_pio2f.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_rem_pio2f.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_rem_pio2f.c,v 1.7 2002/05/28 17:51:46 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/e_remainder.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_remainder.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_remainder.c	(revision 2013)
@@ -7,9 +7,10 @@
  * Developed at SunSoft, a Sun Microsystems, Inc. business.
  * Permission to use, copy, modify, and distribute this
- * software is freely granted, provided that this notice 
+ * software is freely granted, provided that this notice
  * is preserved.
  * ====================================================
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_remainder.c,v 1.10 2005/02/04 18:26:06 das Exp $";
@@ -17,9 +18,9 @@
 
 /* __ieee754_remainder(x,p)
- * Return :                  
- * 	returns  x REM p  =  x - [x/p]*p as if in infinite 
- * 	precise arithmetic, where [x/p] is the (infinite bit) 
+ * Return :
+ * 	returns  x REM p  =  x - [x/p]*p as if in infinite
+ * 	precise arithmetic, where [x/p] is the (infinite bit)
  *	integer nearest x/p (in half way case choose the even one).
- * Method : 
+ * Method :
  *	Based on fmod() return x-[x/p]chopped*p exactlp.
  */
Index: /trunk/src/emx/src/lib/msun/src/e_remainderf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_remainderf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_remainderf.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_remainderf.c,v 1.7 2002/05/28 18:15:04 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/e_scalb.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_scalb.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_scalb.c	(revision 2013)
@@ -7,9 +7,10 @@
  * Developed at SunSoft, a Sun Microsystems, Inc. business.
  * Permission to use, copy, modify, and distribute this
- * software is freely granted, provided that this notice 
+ * software is freely granted, provided that this notice
  * is preserved.
  * ====================================================
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_scalb.c,v 1.12 2005/02/04 18:26:06 das Exp $";
@@ -18,5 +19,5 @@
 /*
  * __ieee754_scalb(x, fn) is provide for
- * passing various standard test suite. One 
+ * passing various standard test suite. One
  * should use scalbn() instead.
  */
Index: /trunk/src/emx/src/lib/msun/src/e_scalbf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_scalbf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_scalbf.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_scalbf.c,v 1.9 2004/07/09 10:01:10 das Exp $";
Index: /trunk/src/emx/src/lib/msun/src/e_sinh.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_sinh.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_sinh.c	(revision 2013)
@@ -7,9 +7,10 @@
  * Developed at SunSoft, a Sun Microsystems, Inc. business.
  * Permission to use, copy, modify, and distribute this
- * software is freely granted, provided that this notice 
+ * software is freely granted, provided that this notice
  * is preserved.
  * ====================================================
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_sinh.c,v 1.9 2005/02/04 18:26:06 das Exp $";
@@ -17,13 +18,13 @@
 
 /* __ieee754_sinh(x)
- * Method : 
+ * Method :
  * mathematically sinh(x) if defined to be (exp(x)-exp(-x))/2
- *	1. Replace x by |x| (sinh(-x) = -sinh(x)). 
- *	2. 
+ *	1. Replace x by |x| (sinh(-x) = -sinh(x)).
+ *	2.
  *		                                    E + E/(E+1)
  *	    0        <= x <= 22     :  sinh(x) := --------------, E=expm1(x)
  *			       			        2
  *
- *	    22       <= x <= lnovft :  sinh(x) := exp(x)/2 
+ *	    22       <= x <= lnovft :  sinh(x) := exp(x)/2
  *	    lnovft   <= x <= ln2ovft:  sinh(x) := exp(x/2)/2 * exp(x/2)
  *	    ln2ovft  <  x	    :  sinh(x) := x*shuge (overflow)
Index: /trunk/src/emx/src/lib/msun/src/e_sinhf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_sinhf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_sinhf.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_sinhf.c,v 1.7 2002/05/28 18:15:04 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/e_sqrt.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_sqrt.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_sqrt.c	(revision 2013)
@@ -7,9 +7,10 @@
  * Developed at SunSoft, a Sun Microsystems, Inc. business.
  * Permission to use, copy, modify, and distribute this
- * software is freely granted, provided that this notice 
+ * software is freely granted, provided that this notice
  * is preserved.
  * ====================================================
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_sqrt.c,v 1.10 2005/02/04 18:26:06 das Exp $";
@@ -21,8 +22,8 @@
  *	     |  Use the hardware sqrt if you have one |
  *           ------------------------------------------
- * Method: 
- *   Bit by bit method using integer arithmetic. (Slow, but portable) 
+ * Method:
+ *   Bit by bit method using integer arithmetic. (Slow, but portable)
  *   1. Normalization
- *	Scale x to y in [1,4) with even powers of 2: 
+ *	Scale x to y in [1,4) with even powers of 2:
  *	find an integer k such that  1 <= (y=x*2^(2k)) < 4, then
  *		sqrt(x) = 2^k * sqrt(y)
@@ -33,7 +34,7 @@
  *	    s  = 2*q , and	y  =  2   * ( y - q  ).		(1)
  *	     i      i            i                 i
- *                                                        
- *	To compute q    from q , one checks whether 
- *		    i+1       i                       
+ *
+ *	To compute q    from q , one checks whether
+ *		    i+1       i
  *
  *			      -(i+1) 2
@@ -45,10 +46,10 @@
  *
  *	With some algebric manipulation, it is not difficult to see
- *	that (2) is equivalent to 
+ *	that (2) is equivalent to
  *                             -(i+1)
  *			s  +  2       <= y			(3)
  *			 i                i
  *
- *	The advantage of (3) is that s  and y  can be computed by 
+ *	The advantage of (3) is that s  and y  can be computed by
  *				      i      i
  *	the following recurrence formula:
@@ -63,7 +64,7 @@
  *           i+1      i          i+1    i     i
  *				
- *	One may easily use induction to prove (4) and (5). 
+ *	One may easily use induction to prove (4) and (5).
  *	Note. Since the left hand side of (3) contain only i+2 bits,
- *	      it does not necessary to do a full (53-bit) comparison 
+ *	      it does not necessary to do a full (53-bit) comparison
  *	      in (3).
  *   3. Final rounding
@@ -105,5 +106,5 @@
 	    return x*x+x;		/* sqrt(NaN)=NaN, sqrt(+inf)=+inf
 					   sqrt(-inf)=sNaN */
-	} 
+	}
     /* take care of zero */
 	if(ix0<=0) {
@@ -139,10 +140,10 @@
 
 	while(r!=0) {
-	    t = s0+r; 
-	    if(t<=ix0) { 
-		s0   = t+r; 
-		ix0 -= t; 
-		q   += r; 
-	    } 
+	    t = s0+r;
+	    if(t<=ix0) {
+		s0   = t+r;
+		ix0 -= t;
+		q   += r;
+	    }
 	    ix0 += ix0 + ((ix1&sign)>>31);
 	    ix1 += ix1;
@@ -152,7 +153,7 @@
 	r = sign;
 	while(r!=0) {
-	    t1 = s1+r; 
+	    t1 = s1+r;
 	    t  = s0;
-	    if((t<ix0)||((t==ix0)&&(t1<=ix1))) { 
+	    if((t<ix0)||((t==ix0)&&(t1<=ix1))) {
 		s1  = t1+r;
 		if(((t1&sign)==sign)&&(s1&sign)==0) s0 += 1;
@@ -175,5 +176,5 @@
 		else if (z>one) {
 		    if (q1==(u_int32_t)0xfffffffe) q+=1;
-		    q1+=2; 
+		    q1+=2;
 		} else
 	            q1 += (q1&1);
@@ -191,8 +192,8 @@
 Other methods  (use floating-point arithmetic)
 -------------
-(This is a copy of a drafted paper by Prof W. Kahan 
+(This is a copy of a drafted paper by Prof W. Kahan
 and K.C. Ng, written in May, 1986)
 
-	Two algorithms are given here to implement sqrt(x) 
+	Two algorithms are given here to implement sqrt(x)
 	(IEEE double precision arithmetic) in software.
 	Both supply sqrt(x) correctly rounded. The first algorithm (in
@@ -200,7 +201,7 @@
 	The second one uses reciproot iterations to avoid division, but
 	requires more multiplications. Both algorithms need the ability
-	to chop results of arithmetic operations instead of round them, 
+	to chop results of arithmetic operations instead of round them,
 	and the INEXACT flag to indicate when an arithmetic operation
-	is executed exactly with no roundoff error, all part of the 
+	is executed exactly with no roundoff error, all part of the
 	standard (IEEE 754-1985). The ability to perform shift, add,
 	subtract and logical AND operations upon 32-bit words is needed
@@ -212,5 +213,5 @@
 
 	Let x0 and x1 be the leading and the trailing 32-bit words of
-	a floating point number x (in IEEE double format) respectively 
+	a floating point number x (in IEEE double format) respectively
 
 	    1    11		     52				  ...widths
@@ -220,5 +221,5 @@
 	      msb    lsb  msb				      lsb ...order
 
- 
+
 	     ------------------------  	     ------------------------
 	x0:  |s|   e    |    f1     |	 x1: |          f2           |
@@ -245,5 +246,5 @@
     (2)	Iterative refinement
 
-	Apply Heron's rule three times to y, we have y approximates 
+	Apply Heron's rule three times to y, we have y approximates
 	sqrt(x) to within 1 ulp (Unit in the Last Place):
 
@@ -270,10 +271,10 @@
 	scaled in advance to avoid spurious overflow in evaluating the
 	expression 3y*y+x. Hence it is not recommended uless division
-	is slow. If division is very slow, then one should use the 
+	is slow. If division is very slow, then one should use the
 	reciproot algorithm given in section B.
 
     (3) Final adjustment
 
-	By twiddling y's last bit it is possible to force y to be 
+	By twiddling y's last bit it is possible to force y to be
 	correctly rounded according to the prevailing rounding mode
 	as follows. Let r and i be copies of the rounding mode and
@@ -306,5 +307,5 @@
 	        R := r;  		... restore rounded mode
 	        return sqrt(x):=y.
-		    
+		
     (4)	Special cases
 
@@ -325,5 +326,5 @@
 	    y0:= k - T2[63&(k>>14)].	... y ~ 1/sqrt(x) to 7.8 bits
 
-	Here k is a 32-bit integer and T2[] is an integer array 
+	Here k is a 32-bit integer and T2[] is an integer array
 	containing correction terms. Now magically the floating
 	value of y (y's leading 32-bit word is y0, the value of
@@ -346,5 +347,5 @@
 	Apply Reciproot iteration three times to y and multiply the
 	result by x to get an approximation z that matches sqrt(x)
-	to about 1 ulp. To be exact, we will have 
+	to about 1 ulp. To be exact, we will have
 		-1ulp < sqrt(x)-z<1.0625ulp.
 	
@@ -357,5 +358,5 @@
 
 	Remark 2. The constant 1.5-2^-30 is chosen to bias the error so that
-	(a) the term z*y in the final iteration is always less than 1; 
+	(a) the term z*y in the final iteration is always less than 1;
 	(b) the error in the final result is biased upward so that
 		-1 ulp < sqrt(x) - z < 1.0625 ulp
@@ -364,5 +365,5 @@
     (3)	Final adjustment
 
-	By twiddling y's last bit it is possible to force y to be 
+	By twiddling y's last bit it is possible to force y to be
 	correctly rounded according to the prevailing rounding mode
 	as follows. Let r and i be copies of the rounding mode and
@@ -404,5 +405,5 @@
 	else {
 	    j := 1 - [(x0>>20)&1]	... j = logb(x) mod 2
-	    k := z1 >> 26;		... get z's 25-th and 26-th 
+	    k := z1 >> 26;		... get z's 25-th and 26-th
 					    fraction bits
 	    I := i or (k&j) or ((k&(j+j+1))!=(x1&3));
@@ -411,5 +412,5 @@
 	return sqrt(x):=z.
 
-	If multiplication is cheaper then the foregoing red tape, the 
+	If multiplication is cheaper then the foregoing red tape, the
 	Inexact flag can be evaluated by
 
@@ -417,5 +418,5 @@
 	    I := (z*z!=x) or I.
 
-	Note that z*z can overwrite I; this value must be sensed if it is 
+	Note that z*z can overwrite I; this value must be sensed if it is
 	True.
 
@@ -424,5 +425,5 @@
 
 		    --------------------
-		z1: |        f2        | 
+		z1: |        f2        |
 		    --------------------
 		bit 31		   bit 0
@@ -442,5 +443,5 @@
 
     (4)	Special cases (see (4) of Section A).	
- 
+
  */
- 
+
Index: /trunk/src/emx/src/lib/msun/src/e_sqrtf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/e_sqrtf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/e_sqrtf.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/e_sqrtf.c,v 1.7 2002/05/28 18:15:04 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/k_cos.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/k_cos.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/k_cos.c	(revision 2013)
@@ -7,9 +7,10 @@
  * Developed at SunSoft, a Sun Microsystems, Inc. business.
  * Permission to use, copy, modify, and distribute this
- * software is freely granted, provided that this notice 
+ * software is freely granted, provided that this notice
  * is preserved.
  * ====================================================
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/k_cos.c,v 1.8 2005/02/04 18:26:06 das Exp $";
@@ -20,5 +21,5 @@
  * kernel cos function on [-pi/4, pi/4], pi/4 ~ 0.785398164
  * Input x is assumed to be bounded by ~pi/4 in magnitude.
- * Input y is the tail of x. 
+ * Input y is the tail of x.
  *
  * Algorithm
@@ -33,10 +34,10 @@
  * 	|              2     4     6     8     10    12     14 |     -58
  * 	|cos(x)-(1-.5*x +C1*x +C2*x +C3*x +C4*x +C5*x  +C6*x  )| <= 2
- * 	|    					               | 
- * 
- * 	               4     6     8     10    12     14 
+ * 	|    					               |
+ *
+ * 	               4     6     8     10    12     14
  *	4. let r = C1*x +C2*x +C3*x +C4*x +C5*x  +C6*x  , then
  *	       cos(x) = 1 - x*x/2 + r
- *	   since cos(x+y) ~ cos(x) - sin(x)*y 
+ *	   since cos(x+y) ~ cos(x) - sin(x)*y
  *			  ~ cos(x) - x*y,
  *	   a correction term is necessary in cos(x) and hence
@@ -75,5 +76,5 @@
 	z  = x*x;
 	r  = z*(C1+z*(C2+z*(C3+z*(C4+z*(C5+z*C6)))));
-	if(ix < 0x3FD33333) 			/* if |x| < 0.3 */ 
+	if(ix < 0x3FD33333) 			/* if |x| < 0.3 */
 	    return one - (0.5*z - (z*r - x*y));
 	else {
Index: /trunk/src/emx/src/lib/msun/src/k_cosf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/k_cosf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/k_cosf.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/k_cosf.c,v 1.7 2002/05/28 18:15:04 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/k_rem_pio2.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/k_rem_pio2.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/k_rem_pio2.c	(revision 2013)
@@ -7,9 +7,10 @@
  * Developed at SunSoft, a Sun Microsystems, Inc. business.
  * Permission to use, copy, modify, and distribute this
- * software is freely granted, provided that this notice 
+ * software is freely granted, provided that this notice
  * is preserved.
  * ====================================================
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/k_rem_pio2.c,v 1.7 2005/02/04 18:26:06 das Exp $";
@@ -19,10 +20,10 @@
  * __kernel_rem_pio2(x,y,e0,nx,prec,ipio2)
  * double x[],y[]; int e0,nx,prec; int ipio2[];
- * 
- * __kernel_rem_pio2 return the last three digits of N with 
+ *
+ * __kernel_rem_pio2 return the last three digits of N with
  *		y = x - N*pi/2
  * so that |y| < pi/2.
  *
- * The method is to compute the integer (mod 8) and fraction parts of 
+ * The method is to compute the integer (mod 8) and fraction parts of
  * (2/pi)*x without doing the full multiplication. In general we
  * skip the part of the product that are known to be a huge integer (
@@ -33,8 +34,8 @@
  *
  * Input parameters:
- * 	x[]	The input value (must be positive) is broken into nx 
+ * 	x[]	The input value (must be positive) is broken into nx
  *		pieces of 24-bit integers in double precision format.
- *		x[i] will be the i-th 24 bit of x. The scaled exponent 
- *		of x[0] is given in input parameter e0 (i.e., x[0]*2^e0 
+ *		x[i] will be the i-th 24 bit of x. The scaled exponent
+ *		of x[0] is given in input parameter e0 (i.e., x[0]*2^e0
  *		match x's up to 24 bits.
  *
@@ -73,6 +74,6 @@
  *
  *	ipio2[]
- *		integer array, contains the (24*i)-th to (24*i+23)-th 
- *		bit of 2/pi after binary point. The corresponding 
+ *		integer array, contains the (24*i)-th to (24*i+23)-th
+ *		bit of 2/pi after binary point. The corresponding
  *		floating value is
  *
@@ -89,6 +90,6 @@
  *		6 for single, double, extended,and quad.
  *
- * 	jz	local integer variable indicating the number of 
- *		terms of ipio2[] used. 
+ * 	jz	local integer variable indicating the number of
+ *		terms of ipio2[] used.
  *
  *	jx	nx - 1
@@ -110,7 +111,7 @@
  *
  *	PIo2[]	double precision array, obtained by cutting pi/2
- *		into 24 bits chunks. 
- *
- *	f[]	ipio2[] in floating point 
+ *		into 24 bits chunks.
+ *
+ *	f[]	ipio2[] in floating point
  *
  *	iq[]	integer array by breaking up q[] in 24-bits chunk.
@@ -126,7 +127,7 @@
 /*
  * Constants:
- * The hexadecimal values are the intended ones for the following 
- * constants. The decimal values may be used, provided that the 
- * compiler will convert from decimal to binary accurately enough 
+ * The hexadecimal values are the intended ones for the following
+ * constants. The decimal values may be used, provided that the
+ * compiler will convert from decimal to binary accurately enough
  * to produce the hexadecimal values shown.
  */
@@ -196,5 +197,5 @@
 	    iq[jz-1] -= i<<(24-q0);
 	    ih = iq[jz-1]>>(23-q0);
-	} 
+	}
 	else if(q0==0) ih = iq[jz-1]>>23;
 	else if(z>=0.5) ih=2;
@@ -247,5 +248,5 @@
 	} else { /* break z into 24-bit if necessary */
 	    z = scalbn(z,-q0);
-	    if(z>=two24) { 
+	    if(z>=two24) {
 		fw = (double)((int32_t)(twon24*z));
 		iq[jz] = (int32_t)(z-two24*fw);
@@ -272,27 +273,27 @@
 		fw = 0.0;
 		for (i=jz;i>=0;i--) fw += fq[i];
-		y[0] = (ih==0)? fw: -fw; 
+		y[0] = (ih==0)? fw: -fw;
 		break;
 	    case 1:
 	    case 2:
 		fw = 0.0;
-		for (i=jz;i>=0;i--) fw += fq[i]; 
-		y[0] = (ih==0)? fw: -fw; 
+		for (i=jz;i>=0;i--) fw += fq[i];
+		y[0] = (ih==0)? fw: -fw;
 		fw = fq[0]-fw;
 		for (i=1;i<=jz;i++) fw += fq[i];
-		y[1] = (ih==0)? fw: -fw; 
+		y[1] = (ih==0)? fw: -fw;
 		break;
 	    case 3:	/* painful */
 		for (i=jz;i>0;i--) {
-		    fw      = fq[i-1]+fq[i]; 
+		    fw      = fq[i-1]+fq[i];
 		    fq[i]  += fq[i-1]-fw;
 		    fq[i-1] = fw;
 		}
 		for (i=jz;i>1;i--) {
-		    fw      = fq[i-1]+fq[i]; 
+		    fw      = fq[i-1]+fq[i];
 		    fq[i]  += fq[i-1]-fw;
 		    fq[i-1] = fw;
 		}
-		for (fw=0.0,i=jz;i>=2;i--) fw += fq[i]; 
+		for (fw=0.0,i=jz;i>=2;i--) fw += fq[i];
 		if(ih==0) {
 		    y[0] =  fq[0]; y[1] =  fq[1]; y[2] =  fw;
Index: /trunk/src/emx/src/lib/msun/src/k_rem_pio2f.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/k_rem_pio2f.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/k_rem_pio2f.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/k_rem_pio2f.c,v 1.7 2005/04/05 23:27:47 das Exp $";
Index: /trunk/src/emx/src/lib/msun/src/k_sin.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/k_sin.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/k_sin.c	(revision 2013)
@@ -7,9 +7,10 @@
  * Developed at SunSoft, a Sun Microsystems, Inc. business.
  * Permission to use, copy, modify, and distribute this
- * software is freely granted, provided that this notice 
+ * software is freely granted, provided that this notice
  * is preserved.
  * ====================================================
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/k_sin.c,v 1.8 2005/02/04 18:26:06 das Exp $";
@@ -20,8 +21,8 @@
  * Input x is assumed to be bounded by ~pi/4 in magnitude.
  * Input y is the tail of x.
- * Input iy indicates whether y is 0. (if iy=0, y assume to be 0). 
+ * Input iy indicates whether y is 0. (if iy=0, y assume to be 0).
  *
  * Algorithm
- *	1. Since sin(-x) = -sin(x), we need only to consider positive x. 
+ *	1. Since sin(-x) = -sin(x), we need only to consider positive x.
  *	2. if x < 2^-27 (hx<0x3e400000 0), return x with inexact if x!=0.
  *	3. sin(x) is approximated by a polynomial of degree 13 on
@@ -33,9 +34,9 @@
  * 	|sin(x)         2     4     6     8     10     12  |     -58
  * 	|----- - (1+S1*x +S2*x +S3*x +S4*x +S5*x  +S6*x   )| <= 2
- * 	|  x 					           | 
- * 
+ * 	|  x 					           |
+ *
  *	4. sin(x+y) = sin(x) + sin'(x')*y
  *		    ~ sin(x) + (1-x*x/2)*y
- *	   For better accuracy, let 
+ *	   For better accuracy, let
  *		     3      2      2      2      2
  *		r = x *(S2+x *(S3+x *(S4+x *(S5+x *S6))))
Index: /trunk/src/emx/src/lib/msun/src/k_sinf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/k_sinf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/k_sinf.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/k_sinf.c,v 1.7 2002/05/28 18:15:04 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/k_tan.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/k_tan.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/k_tan.c	(revision 2013)
@@ -12,4 +12,5 @@
 
 /* INDENT OFF */
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/k_tan.c,v 1.10 2005/02/04 18:26:06 das Exp $";
Index: /trunk/src/emx/src/lib/msun/src/k_tanf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/k_tanf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/k_tanf.c	(revision 2013)
@@ -13,4 +13,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/k_tanf.c,v 1.8 2004/06/02 04:39:44 das Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_asinh.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_asinh.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_asinh.c	(revision 2013)
@@ -11,4 +11,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_asinh.c,v 1.8 2002/05/28 18:15:04 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_asinhf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_asinhf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_asinhf.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_asinhf.c,v 1.8 2002/05/28 18:15:04 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_atan.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_atan.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_atan.c	(revision 2013)
@@ -11,4 +11,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_atan.c,v 1.9 2003/07/23 04:53:46 peter Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_atanf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_atanf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_atanf.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_atanf.c,v 1.7 2002/05/28 18:15:04 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_cbrt.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_cbrt.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_cbrt.c	(revision 2013)
@@ -11,4 +11,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_cbrt.c,v 1.7 2002/05/28 18:15:04 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_cbrtf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_cbrtf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_cbrtf.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_cbrtf.c,v 1.7 2002/05/28 18:15:04 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_ceil.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_ceil.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_ceil.c	(revision 2013)
@@ -11,4 +11,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_ceil.c,v 1.9 2003/07/23 04:53:46 peter Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_ceilf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_ceilf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_ceilf.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_ceilf.c,v 1.7 2002/05/28 18:15:04 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_ceill.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_ceill.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_ceill.c	(revision 2013)
@@ -12,4 +12,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_ceill.c,v 1.3 2005/04/22 09:57:55 stefanf Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_cimag.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_cimag.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_cimag.c	(revision 2013)
@@ -27,4 +27,5 @@
  */
 
+#include "namespace.h"
 #include <complex.h>
 
Index: /trunk/src/emx/src/lib/msun/src/s_cimagf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_cimagf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_cimagf.c	(revision 2013)
@@ -27,4 +27,5 @@
  */
 
+#include "namespace.h"
 #include <complex.h>
 
Index: /trunk/src/emx/src/lib/msun/src/s_cimagl.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_cimagl.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_cimagl.c	(revision 2013)
@@ -27,4 +27,5 @@
  */
 
+#include "namespace.h"
 #include <complex.h>
 
Index: /trunk/src/emx/src/lib/msun/src/s_conj.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_conj.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_conj.c	(revision 2013)
@@ -27,4 +27,5 @@
  */
 
+#include "namespace.h"
 #include <complex.h>
 
Index: /trunk/src/emx/src/lib/msun/src/s_conjf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_conjf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_conjf.c	(revision 2013)
@@ -27,4 +27,5 @@
  */
 
+#include "namespace.h"
 #include <complex.h>
 
Index: /trunk/src/emx/src/lib/msun/src/s_conjl.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_conjl.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_conjl.c	(revision 2013)
@@ -27,4 +27,5 @@
  */
 
+#include "namespace.h"
 #include <complex.h>
 
Index: /trunk/src/emx/src/lib/msun/src/s_copysign.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_copysign.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_copysign.c	(revision 2013)
@@ -11,4 +11,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_copysign.c,v 1.9 2003/07/23 04:53:46 peter Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_copysignf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_copysignf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_copysignf.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_copysignf.c,v 1.9 2003/07/23 04:53:46 peter Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_copysignl.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_copysignl.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_copysignl.c	(revision 2013)
@@ -27,4 +27,5 @@
  */
 
+#include "namespace.h"
 #include <math.h>
 
Index: /trunk/src/emx/src/lib/msun/src/s_cos.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_cos.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_cos.c	(revision 2013)
@@ -11,4 +11,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_cos.c,v 1.9 2003/07/23 04:53:46 peter Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_cosf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_cosf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_cosf.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_cosf.c,v 1.7 2002/05/28 18:15:04 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_creal.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_creal.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_creal.c	(revision 2013)
@@ -27,4 +27,5 @@
  */
 
+#include "namespace.h"
 #include <complex.h>
 
Index: /trunk/src/emx/src/lib/msun/src/s_crealf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_crealf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_crealf.c	(revision 2013)
@@ -27,4 +27,5 @@
  */
 
+#include "namespace.h"
 #include <complex.h>
 
Index: /trunk/src/emx/src/lib/msun/src/s_creall.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_creall.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_creall.c	(revision 2013)
@@ -27,4 +27,5 @@
  */
 
+#include "namespace.h"
 #include <complex.h>
 
Index: /trunk/src/emx/src/lib/msun/src/s_erf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_erf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_erf.c	(revision 2013)
@@ -11,4 +11,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_erf.c,v 1.7 2002/05/28 18:15:04 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_erff.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_erff.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_erff.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_erff.c,v 1.7 2002/05/28 18:15:04 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_exp2.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_exp2.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_exp2.c	(revision 2013)
@@ -25,4 +25,5 @@
  */
 
+#include "namespace.h"
 #include <sys/cdefs.h>
 __FBSDID("$FreeBSD: src/lib/msun/src/s_exp2.c,v 1.1 2005/04/05 02:57:15 das Exp $");
@@ -350,5 +351,5 @@
 			if(((ix & 0xfffff) | lx) != 0 || (hx & 0x80000000) == 0)
 				return (x); 	/* x is NaN or +Inf */
-			else 
+			else
 				return (0.0);	/* x is -Inf */
 		}
Index: /trunk/src/emx/src/lib/msun/src/s_exp2f.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_exp2f.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_exp2f.c	(revision 2013)
@@ -25,4 +25,5 @@
  */
 
+#include "namespace.h"
 #include <sys/cdefs.h>
 __FBSDID("$FreeBSD: src/lib/msun/src/s_exp2f.c,v 1.1 2005/04/05 02:57:15 das Exp $");
@@ -103,5 +104,5 @@
 			if ((ix & 0x7fffff) != 0 || (hx & 0x80000000) == 0)
 				return (x); 	/* x is NaN or +Inf */
-			else 
+			else
 				return (0.0);	/* x is -Inf */
 		}
Index: /trunk/src/emx/src/lib/msun/src/s_expm1.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_expm1.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_expm1.c	(revision 2013)
@@ -11,4 +11,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_expm1.c,v 1.7 2002/05/28 18:15:04 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_expm1f.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_expm1f.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_expm1f.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_expm1f.c,v 1.7 2002/05/28 18:15:04 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_fabs.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_fabs.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_fabs.c	(revision 2013)
@@ -11,4 +11,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_fabs.c,v 1.7 2002/05/28 18:15:04 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_fabsf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_fabsf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_fabsf.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_fabsf.c,v 1.7 2002/05/28 18:15:04 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_fabsl.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_fabsl.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_fabsl.c	(revision 2013)
@@ -29,4 +29,5 @@
  */
 
+#include "namespace.h"
 #include <math.h>
 
Index: /trunk/src/emx/src/lib/msun/src/s_fdim.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_fdim.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_fdim.c	(revision 2013)
@@ -25,4 +25,5 @@
  */
 
+#include "namespace.h"
 #include <sys/cdefs.h>
 __FBSDID("$FreeBSD: src/lib/msun/src/s_fdim.c,v 1.1 2004/06/30 07:04:01 das Exp $");
Index: /trunk/src/emx/src/lib/msun/src/s_finite.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_finite.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_finite.c	(revision 2013)
@@ -11,4 +11,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_finite.c,v 1.8 2003/07/23 04:53:46 peter Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_finitef.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_finitef.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_finitef.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_finitef.c,v 1.6 2002/05/28 17:51:46 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_floor.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_floor.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_floor.c	(revision 2013)
@@ -11,4 +11,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_floor.c,v 1.9 2003/07/23 04:53:46 peter Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_floorf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_floorf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_floorf.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_floorf.c,v 1.7 2002/05/28 18:15:04 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_floorl.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_floorl.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_floorl.c	(revision 2013)
@@ -12,4 +12,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_floorl.c,v 1.3 2005/04/22 09:57:55 stefanf Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_fma.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_fma.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_fma.c	(revision 2013)
@@ -25,4 +25,5 @@
  */
 
+#include "namespace.h"
 #include <sys/cdefs.h>
 __FBSDID("$FreeBSD: src/lib/msun/src/s_fma.c,v 1.4 2005/03/18 02:27:59 das Exp $");
Index: /trunk/src/emx/src/lib/msun/src/s_fmaf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_fmaf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_fmaf.c	(revision 2013)
@@ -25,4 +25,5 @@
  */
 
+#include "namespace.h"
 #include <sys/cdefs.h>
 __FBSDID("$FreeBSD: src/lib/msun/src/s_fmaf.c,v 1.1 2005/01/22 09:53:18 das Exp $");
Index: /trunk/src/emx/src/lib/msun/src/s_fmal.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_fmal.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_fmal.c	(revision 2013)
@@ -25,4 +25,5 @@
  */
 
+#include "namespace.h"
 #include <sys/cdefs.h>
 __FBSDID("$FreeBSD: src/lib/msun/src/s_fmal.c,v 1.2 2005/03/18 02:27:59 das Exp $");
Index: /trunk/src/emx/src/lib/msun/src/s_fmax.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_fmax.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_fmax.c	(revision 2013)
@@ -25,4 +25,5 @@
  */
 
+#include "namespace.h"
 #include <sys/cdefs.h>
 __FBSDID("$FreeBSD: src/lib/msun/src/s_fmax.c,v 1.1 2004/06/30 07:04:01 das Exp $");
Index: /trunk/src/emx/src/lib/msun/src/s_fmaxf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_fmaxf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_fmaxf.c	(revision 2013)
@@ -25,4 +25,5 @@
  */
 
+#include "namespace.h"
 #include <sys/cdefs.h>
 __FBSDID("$FreeBSD: src/lib/msun/src/s_fmaxf.c,v 1.1 2004/06/30 07:04:01 das Exp $");
Index: /trunk/src/emx/src/lib/msun/src/s_fmaxl.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_fmaxl.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_fmaxl.c	(revision 2013)
@@ -25,4 +25,5 @@
  */
 
+#include "namespace.h"
 #include <sys/cdefs.h>
 __FBSDID("$FreeBSD: src/lib/msun/src/s_fmaxl.c,v 1.1 2004/06/30 07:04:01 das Exp $");
Index: /trunk/src/emx/src/lib/msun/src/s_fmin.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_fmin.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_fmin.c	(revision 2013)
@@ -25,4 +25,5 @@
  */
 
+#include "namespace.h"
 #include <sys/cdefs.h>
 __FBSDID("$FreeBSD: src/lib/msun/src/s_fmin.c,v 1.1 2004/06/30 07:04:01 das Exp $");
Index: /trunk/src/emx/src/lib/msun/src/s_fminf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_fminf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_fminf.c	(revision 2013)
@@ -25,4 +25,5 @@
  */
 
+#include "namespace.h"
 #include <sys/cdefs.h>
 __FBSDID("$FreeBSD: src/lib/msun/src/s_fminf.c,v 1.1 2004/06/30 07:04:01 das Exp $");
Index: /trunk/src/emx/src/lib/msun/src/s_fminl.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_fminl.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_fminl.c	(revision 2013)
@@ -25,4 +25,5 @@
  */
 
+#include "namespace.h"
 #include <sys/cdefs.h>
 __FBSDID("$FreeBSD: src/lib/msun/src/s_fminl.c,v 1.1 2004/06/30 07:04:01 das Exp $");
Index: /trunk/src/emx/src/lib/msun/src/s_frexp.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_frexp.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_frexp.c	(revision 2013)
@@ -11,4 +11,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_frexp.c,v 1.10 2005/03/07 21:27:37 das Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_frexpf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_frexpf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_frexpf.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_frexpf.c,v 1.8 2002/05/28 18:15:04 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_frexpl.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_frexpl.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_frexpl.c	(revision 2013)
@@ -27,4 +27,5 @@
  */
 
+#include "namespace.h"
 #include <float.h>
 #include <math.h>
Index: /trunk/src/emx/src/lib/msun/src/s_ilogb.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_ilogb.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_ilogb.c	(revision 2013)
@@ -11,4 +11,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_ilogb.c,v 1.9 2004/10/09 17:14:28 stefanf Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_ilogbf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_ilogbf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_ilogbf.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_ilogbf.c,v 1.7 2004/10/09 17:14:28 stefanf Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_ilogbl.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_ilogbl.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_ilogbl.c	(revision 2013)
@@ -11,4 +11,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_ilogbl.c,v 1.1 2004/10/11 18:13:52 stefanf Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_isfinite.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_isfinite.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_isfinite.c	(revision 2013)
@@ -27,4 +27,5 @@
  */
 
+#include "namespace.h"
 #include <math.h>
 
Index: /trunk/src/emx/src/lib/msun/src/s_isnan.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_isnan.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_isnan.c	(revision 2013)
@@ -27,4 +27,5 @@
  */
 
+#include "namespace.h"
 #include <math.h>
 
@@ -32,7 +33,7 @@
 
 /* Provided by libc */
-#if 0
+/*#if 0 - not by my libc */
 int
-isnan(double d)
+(isnan)(double d)
 {
 	union IEEEd2bits u;
@@ -41,8 +42,8 @@
 	return (u.bits.exp == 2047 && (u.bits.manl != 0 || u.bits.manh != 0));
 }
-#endif
+/*#endif*/
 
 int
-isnanf(float f)
+(isnanf)(float f)
 {
 	union IEEEf2bits u;
Index: /trunk/src/emx/src/lib/msun/src/s_isnormal.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_isnormal.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_isnormal.c	(revision 2013)
@@ -27,4 +27,5 @@
  */
 
+#include "namespace.h"
 #include <math.h>
 
Index: /trunk/src/emx/src/lib/msun/src/s_llrint.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_llrint.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_llrint.c	(revision 2013)
@@ -1,2 +1,3 @@
+#include "namespace.h"
 #include <sys/cdefs.h>
 __FBSDID("$FreeBSD: src/lib/msun/src/s_llrint.c,v 1.1 2005/01/11 23:12:55 das Exp $");
Index: /trunk/src/emx/src/lib/msun/src/s_llrintf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_llrintf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_llrintf.c	(revision 2013)
@@ -1,2 +1,3 @@
+#include "namespace.h"
 #include <sys/cdefs.h>
 __FBSDID("$FreeBSD: src/lib/msun/src/s_llrintf.c,v 1.1 2005/01/11 23:12:55 das Exp $");
Index: /trunk/src/emx/src/lib/msun/src/s_llround.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_llround.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_llround.c	(revision 2013)
@@ -1,2 +1,3 @@
+#include "namespace.h"
 #include <sys/cdefs.h>
 __FBSDID("$FreeBSD: src/lib/msun/src/s_llround.c,v 1.2 2005/04/08 00:52:27 das Exp $");
Index: /trunk/src/emx/src/lib/msun/src/s_llroundf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_llroundf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_llroundf.c	(revision 2013)
@@ -1,2 +1,3 @@
+#include "namespace.h"
 #include <sys/cdefs.h>
 __FBSDID("$FreeBSD: src/lib/msun/src/s_llroundf.c,v 1.2 2005/04/08 00:52:27 das Exp $");
Index: /trunk/src/emx/src/lib/msun/src/s_llroundl.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_llroundl.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_llroundl.c	(revision 2013)
@@ -1,2 +1,3 @@
+#include "namespace.h"
 #include <sys/cdefs.h>
 __FBSDID("$FreeBSD: src/lib/msun/src/s_llroundl.c,v 1.1 2005/04/08 01:24:08 das Exp $");
Index: /trunk/src/emx/src/lib/msun/src/s_log1p.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_log1p.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_log1p.c	(revision 2013)
@@ -11,4 +11,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_log1p.c,v 1.7 2002/05/28 18:15:04 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_log1pf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_log1pf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_log1pf.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_log1pf.c,v 1.7 2002/05/28 18:15:04 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_logb.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_logb.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_logb.c	(revision 2013)
@@ -11,4 +11,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_logb.c,v 1.9 2003/07/23 04:53:46 peter Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_logbf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_logbf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_logbf.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_logbf.c,v 1.7 2002/05/28 18:15:04 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_lrint.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_lrint.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_lrint.c	(revision 2013)
@@ -25,4 +25,5 @@
  */
 
+#include "namespace.h"
 #include <sys/cdefs.h>
 #include <fenv.h>
Index: /trunk/src/emx/src/lib/msun/src/s_lrintf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_lrintf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_lrintf.c	(revision 2013)
@@ -1,2 +1,3 @@
+#include "namespace.h"
 #include <sys/cdefs.h>
 __FBSDID("$FreeBSD: src/lib/msun/src/s_lrintf.c,v 1.1 2005/01/11 23:12:55 das Exp $");
Index: /trunk/src/emx/src/lib/msun/src/s_lround.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_lround.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_lround.c	(revision 2013)
@@ -25,4 +25,5 @@
  */
 
+#include "namespace.h"
 #include <sys/cdefs.h>
 #include <sys/limits.h>
@@ -65,2 +66,3 @@
 	}
 }
+
Index: /trunk/src/emx/src/lib/msun/src/s_lroundf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_lroundf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_lroundf.c	(revision 2013)
@@ -1,2 +1,3 @@
+#include "namespace.h"
 #include <sys/cdefs.h>
 __FBSDID("$FreeBSD: src/lib/msun/src/s_lroundf.c,v 1.2 2005/04/08 00:52:27 das Exp $");
Index: /trunk/src/emx/src/lib/msun/src/s_lroundl.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_lroundl.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_lroundl.c	(revision 2013)
@@ -1,2 +1,3 @@
+#include "namespace.h"
 #include <sys/cdefs.h>
 __FBSDID("$FreeBSD: src/lib/msun/src/s_lroundl.c,v 1.1 2005/04/08 01:24:08 das Exp $");
Index: /trunk/src/emx/src/lib/msun/src/s_modf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_modf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_modf.c	(revision 2013)
@@ -11,4 +11,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_modf.c,v 1.7 2002/05/28 18:15:04 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_modff.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_modff.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_modff.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_modff.c,v 1.7 2002/05/28 18:15:04 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_nearbyint.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_nearbyint.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_nearbyint.c	(revision 2013)
@@ -25,4 +25,5 @@
  */
 
+#include "namespace.h"
 #include <sys/cdefs.h>
 __FBSDID("$FreeBSD: src/lib/msun/src/s_nearbyint.c,v 1.1 2004/07/06 04:46:08 das Exp $");
Index: /trunk/src/emx/src/lib/msun/src/s_nextafter.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_nextafter.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_nextafter.c	(revision 2013)
@@ -11,4 +11,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_nextafter.c,v 1.11 2005/03/07 21:27:37 das Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_nextafterf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_nextafterf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_nextafterf.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_nextafterf.c,v 1.10 2005/03/07 04:55:58 das Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_nextafterl.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_nextafterl.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_nextafterl.c	(revision 2013)
@@ -11,4 +11,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_nextafterl.c,v 1.1 2005/03/07 04:56:46 das Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_nexttoward.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_nexttoward.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_nexttoward.c	(revision 2013)
@@ -11,4 +11,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_nexttoward.c,v 1.1 2005/03/07 04:56:46 das Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_nexttowardf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_nexttowardf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_nexttowardf.c	(revision 2013)
@@ -10,4 +10,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_nexttowardf.c,v 1.1 2005/03/07 04:57:38 das Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_remquo.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_remquo.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_remquo.c	(revision 2013)
@@ -6,9 +6,10 @@
  * Developed at SunSoft, a Sun Microsystems, Inc. business.
  * Permission to use, copy, modify, and distribute this
- * software is freely granted, provided that this notice 
+ * software is freely granted, provided that this notice
  * is preserved.
  * ====================================================
  */
 
+#include "namespace.h"
 #include <sys/cdefs.h>
 __FBSDID("$FreeBSD: src/lib/msun/src/s_remquo.c,v 1.1 2005/03/25 04:40:44 das Exp $");
@@ -74,5 +75,5 @@
 
     /* set up {hx,lx}, {hy,ly} and align y to x */
-	if(ix >= -1022) 
+	if(ix >= -1022)
 	    hx = 0x00100000|(0x000fffff&hx);
 	else {		/* subnormal x, shift x to normal */
@@ -86,5 +87,5 @@
 	    }
 	}
-	if(iy >= -1022) 
+	if(iy >= -1022)
 	    hy = 0x00100000|(0x000fffff&hy);
 	else {		/* subnormal y, shift y to normal */
Index: /trunk/src/emx/src/lib/msun/src/s_remquof.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_remquof.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_remquof.c	(revision 2013)
@@ -6,9 +6,10 @@
  * Developed at SunSoft, a Sun Microsystems, Inc. business.
  * Permission to use, copy, modify, and distribute this
- * software is freely granted, provided that this notice 
+ * software is freely granted, provided that this notice
  * is preserved.
  * ====================================================
  */
 
+#include "namespace.h"
 #include <sys/cdefs.h>
 __FBSDID("$FreeBSD: src/lib/msun/src/s_remquof.c,v 1.1 2005/03/25 04:40:44 das Exp $");
Index: /trunk/src/emx/src/lib/msun/src/s_rint.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_rint.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_rint.c	(revision 2013)
@@ -11,4 +11,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_rint.c,v 1.11 2005/03/07 20:28:26 stefanf Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_rintf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_rintf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_rintf.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_rintf.c,v 1.9 2004/06/09 21:24:52 das Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_round.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_round.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_round.c	(revision 2013)
@@ -25,4 +25,5 @@
  */
 
+#include "namespace.h"
 #include <sys/cdefs.h>
 __FBSDID("$FreeBSD: src/lib/msun/src/s_round.c,v 1.3 2005/01/13 15:43:41 das Exp $");
Index: /trunk/src/emx/src/lib/msun/src/s_roundf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_roundf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_roundf.c	(revision 2013)
@@ -25,4 +25,5 @@
  */
 
+#include "namespace.h"
 #include <sys/cdefs.h>
 __FBSDID("$FreeBSD: src/lib/msun/src/s_roundf.c,v 1.3 2005/01/13 15:43:41 das Exp $");
Index: /trunk/src/emx/src/lib/msun/src/s_roundl.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_roundl.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_roundl.c	(revision 2013)
@@ -25,4 +25,5 @@
  */
 
+#include "namespace.h"
 #include <sys/cdefs.h>
 __FBSDID("$FreeBSD: src/lib/msun/src/s_roundl.c,v 1.1 2005/04/08 01:24:08 das Exp $");
Index: /trunk/src/emx/src/lib/msun/src/s_scalbln.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_scalbln.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_scalbln.c	(revision 2013)
@@ -25,4 +25,5 @@
  */
 
+#include "namespace.h"
 #include <sys/cdefs.h>
 __FBSDID("$FreeBSD: src/lib/msun/src/s_scalbln.c,v 1.2 2005/03/07 04:57:50 das Exp $");
Index: /trunk/src/emx/src/lib/msun/src/s_scalbn.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_scalbn.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_scalbn.c	(revision 2013)
@@ -11,4 +11,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_scalbn.c,v 1.11 2005/03/07 21:27:37 das Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_scalbnf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_scalbnf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_scalbnf.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_scalbnf.c,v 1.8 2005/03/07 04:52:43 das Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_scalbnl.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_scalbnl.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_scalbnl.c	(revision 2013)
@@ -11,4 +11,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_scalbnl.c,v 1.1 2005/03/07 04:52:58 das Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_signbit.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_signbit.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_signbit.c	(revision 2013)
@@ -27,4 +27,5 @@
  */
 
+#include "namespace.h"
 #include <math.h>
 
Index: /trunk/src/emx/src/lib/msun/src/s_signgam.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_signgam.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_signgam.c	(revision 2013)
@@ -1,2 +1,3 @@
+#include "namespace.h"
 #include "math.h"
 #include "math_private.h"
Index: /trunk/src/emx/src/lib/msun/src/s_significand.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_significand.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_significand.c	(revision 2013)
@@ -11,4 +11,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_significand.c,v 1.9 2003/07/23 04:53:47 peter Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_significandf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_significandf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_significandf.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_significandf.c,v 1.7 2002/05/28 18:15:04 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_sin.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_sin.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_sin.c	(revision 2013)
@@ -11,4 +11,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_sin.c,v 1.9 2003/07/23 04:53:47 peter Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_sinf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_sinf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_sinf.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_sinf.c,v 1.7 2002/05/28 18:15:04 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_tan.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_tan.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_tan.c	(revision 2013)
@@ -11,4 +11,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_tan.c,v 1.9 2003/07/23 04:53:47 peter Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_tanf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_tanf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_tanf.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_tanf.c,v 1.7 2002/05/28 18:15:04 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_tanh.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_tanh.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_tanh.c	(revision 2013)
@@ -11,4 +11,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_tanh.c,v 1.7 2002/05/28 18:15:04 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_tanhf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_tanhf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_tanhf.c	(revision 2013)
@@ -14,4 +14,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_tanhf.c,v 1.7 2002/05/28 18:15:04 alfred Exp $";
Index: /trunk/src/emx/src/lib/msun/src/s_trunc.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_trunc.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_trunc.c	(revision 2013)
@@ -11,4 +11,5 @@
  */
 
+#include "namespace.h"
 #include <sys/cdefs.h>
 __FBSDID("$FreeBSD: src/lib/msun/src/s_trunc.c,v 1.1 2004/06/20 09:25:43 das Exp $");
Index: /trunk/src/emx/src/lib/msun/src/s_truncf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_truncf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_truncf.c	(revision 2013)
@@ -11,4 +11,5 @@
  */
 
+#include "namespace.h"
 #include <sys/cdefs.h>
 __FBSDID("$FreeBSD: src/lib/msun/src/s_truncf.c,v 1.1 2004/06/20 09:25:43 das Exp $");
Index: /trunk/src/emx/src/lib/msun/src/s_truncl.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/s_truncl.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/s_truncl.c	(revision 2013)
@@ -12,4 +12,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static char rcsid[] = "$FreeBSD: src/lib/msun/src/s_truncl.c,v 1.3 2005/04/22 09:57:55 stefanf Exp $";
Index: /trunk/src/emx/src/lib/msun/src/w_cabs.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/w_cabs.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/w_cabs.c	(revision 2013)
@@ -6,4 +6,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static const char rcsid[] =
Index: /trunk/src/emx/src/lib/msun/src/w_cabsf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/w_cabsf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/w_cabsf.c	(revision 2013)
@@ -6,4 +6,5 @@
  */
 
+#include "namespace.h"
 #ifndef lint
 static const char rcsid[] =
Index: /trunk/src/emx/src/lib/msun/src/w_drem.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/w_drem.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/w_drem.c	(revision 2013)
@@ -6,4 +6,5 @@
  */
 
+#include "namespace.h"
 #include <math.h>
 
Index: /trunk/src/emx/src/lib/msun/src/w_dremf.c
===================================================================
--- /trunk/src/emx/src/lib/msun/src/w_dremf.c	(revision 2012)
+++ /trunk/src/emx/src/lib/msun/src/w_dremf.c	(revision 2013)
@@ -7,4 +7,5 @@
 /* $FreeBSD: src/lib/msun/src/w_dremf.c,v 1.3 2004/07/28 05:53:18 kan Exp $ */
 
+#include "namespace.h"
 #include "math.h"
 #include "math_private.h"
