 @@ -1,14 +1,14 @@
-/*	$NetBSD: n_atan2.S,v 1.8 2008/03/20 18:49:39 mhitch Exp $	*/
+/*	$NetBSD: n_atan2.S,v 1.9 2014/10/10 20:58:09 martin Exp $	*/
 /*
  * Copyright (c) 1985, 1993
  *	The Regents of the University of California.  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.
  * 3. Neither the name of the University nor the names of its contributors
 @@ -74,26 +74,27 @@
 #ifdef WEAK_ALIAS
 WEAK_ALIAS(atan2f, _atan2f)
 #endif
 ENTRY(_atan2f, 0)
 	cvtfd	4(%ap),-(%sp)
 	calls	$2,_C_LABEL(_atan2)
 	cvtdf	%r0,%r0
 	ret
 #ifdef WEAK_ALIAS
 WEAK_ALIAS(atan2, _atan2)
 WEAK_ALIAS(_atan2l, _atan2)
 #endif
 ENTRY(_atan2, 0x0fc0)
 	movq	4(%ap),%r2		# %r2 = y
 	movq	12(%ap),%r4		# %r4 = x
 	bicw3	$0x7f,%r2,%r0
 	bicw3	$0x7f,%r4,%r1
 	cmpw	%r0,$0x8000		# y is the reserved operand
 	jeql	resop
 	cmpw	%r1,$0x8000		# x is the reserved operand
 	jeql	resop
 	subl2	$8,%sp
 	bicw3	$0x7fff,%r2,-4(%fp)	# copy y sign bit to -4(%fp)

 @@ -1,14 +1,14 @@
-/*	$NetBSD: n_cabs.S,v 1.6 2008/03/20 16:41:26 mhitch Exp $	*/
+/*	$NetBSD: n_cabs.S,v 1.7 2014/10/10 20:58:09 martin Exp $	*/
 /*
  * Copyright (c) 1985, 1993
  *	The Regents of the University of California.  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.
  * 3. Neither the name of the University nor the names of its contributors
 @@ -45,26 +45,28 @@
 /*	entry for c functions cabs and hypot */
 #ifdef WEAK_ALIAS
 WEAK_ALIAS(hypotf, _hypotf)
 #endif
 ENTRY(_hypotf, 0)
 	cvtfd	4(%ap),-(%sp)
 	calls	$2,_C_LABEL(_hypot)
 	cvtdf	%r0,%r0
 	ret
 #ifdef WEAK_ALIAS
 WEAK_ALIAS(hypot, _hypot)
 WEAK_ALIAS(hypotl, _hypot)
 WEAK_ALIAS(_hypotl, _hypot)
 #endif
 ALTENTRY(cabs)
 ENTRY(_hypot, 0x8040) 		# save %r6, enable floating overflow
 	movq	4(%ap),%r0	# %r0:1 = x
 	movq	12(%ap),%r2	# %r2:3 = y
 	jbr	cabs2
 /*	entry for Fortran use, call by:   d = abs(z) */
 ENTRY(z_abs, 0x8040)		# save %r6, enable floating overflow
 	movl	4(%ap),%r2	# indirect addressing is necessary here
 	movq	(%r2)+,%r0	# %r0:1 = x
 	movq	(%r2),%r2		# %r2:3 = y

 @@ -1,14 +1,14 @@
-/*      $NetBSD: n_atan2.c,v 1.6 2003/08/07 16:44:50 agc Exp $        */
+/*      $NetBSD: n_atan2.c,v 1.7 2014/10/10 20:58:09 martin Exp $        */
 /*
  * Copyright (c) 1985, 1993
  *	The Regents of the University of California.  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.
  * 3. Neither the name of the University nor the names of its contributors
 @@ -160,26 +160,30 @@ ic(a11,    1.6438029044759730479E-2  ,
 #define	a2	vccast(a2)
 #define	a3	vccast(a3)
 #define	a4	vccast(a4)
 #define	a5	vccast(a5)
 #define	a6	vccast(a6)
 #define	a7	vccast(a7)
 #define	a8	vccast(a8)
 #define	a9	vccast(a9)
 #define	a10	vccast(a10)
 #define	a11	vccast(a11)
 #define	a12	vccast(a12)
 #endif
 #ifdef __weak_alias
 __weak_alias(_atan2l, atan2);
 #endif
 double
 atan2(double y, double x)
+{
 	static const double zero=0, one=1, small=1.0E-9, big=1.0E18;
 	double t,z,signy,signx,hi,lo;
 	int k,m;
 #if !defined(__vax__)&&!defined(tahoe)
     /* if x or y is NAN */
 	if(x!=x) return(x); if(y!=y) return(y);
 #endif	/* !defined(__vax__)&&!defined(tahoe) */
     /* copy down the sign of y and x */

 @@ -1,14 +1,14 @@
-/*	$NetBSD: n_sincos.c,v 1.6 2003/08/07 16:44:52 agc Exp $	*/
+/*	$NetBSD: n_sincos.c,v 1.7 2014/10/10 20:58:09 martin Exp $	*/
 /*
  * Copyright (c) 1987, 1993
  *	The Regents of the University of California.  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.
  * 3. Neither the name of the University nor the names of its contributors
 @@ -28,26 +28,31 @@
  * SUCH DAMAGE.
  */
 #ifndef lint
 #if 0
 static char sccsid[] = "@(#)sincos.c	8.1 (Berkeley) 6/4/93";
 #endif
 #endif /* not lint */
 #define _LIBM_DECLARE
 #include "mathimpl.h"
 #include "trig.h"
 #ifdef __weak_alias
 __weak_alias(_sinl, sin);
 __weak_alias(_cosl, cos);
 #endif
 double
 sin(double x)
+{
 	double a,c,z;
         if(!finite(x))		/* sin(NaN) and sin(INF) must be NaN */
 		return x-x;
 	x=drem(x,PI2);		/* reduce x into [-PI,PI] */
 	a=copysign(x,__one);
 	if (a >= PIo4) {
 		if(a >= PI3o4)		/* ... in [3PI/4,PI] */
 			x = copysign((a = PI-a),x);
 		else {			/* ... in [PI/4,3PI/4]  */

 @@ -1,14 +1,14 @@
-/*      $NetBSD: n_cosh.c,v 1.8 2008/03/20 16:41:26 mhitch Exp $ */
+/*      $NetBSD: n_cosh.c,v 1.9 2014/10/10 20:58:09 martin Exp $ */
 /*
  * Copyright (c) 1985, 1993
  *	The Regents of the University of California.  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.
  * 3. Neither the name of the University nor the names of its contributors
 @@ -77,26 +77,27 @@ static char sccsid[] = "@(#)cosh.c	8.1 (
  * 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 to produce the hexadecimal values
  * shown.
  */
 #define _LIBM_STATIC
 #include "../src/namespace.h"
 #include "mathimpl.h"
 #ifdef __weak_alias
 __weak_alias(cosh, _cosh);
 __weak_alias(_coshl, _cosh);
 __weak_alias(coshf, _coshf);
 #endif
 vc(mln2hi, 8.8029691931113054792E1   ,0f33,43b0,2bdb,c7e2,   7, .B00F33C7E22BDB)
 vc(mln2lo,-4.9650192275318476525E-16 ,1b60,a70f,582a,279e, -50,-.8F1B60279E582A)
 vc(lnovfl, 8.8029691931113053016E1   ,0f33,43b0,2bda,c7e2,   7, .B00F33C7E22BDA)
 ic(mln2hi, 7.0978271289338397310E2,    10, 1.62E42FEFA39EF)
 ic(mln2lo, 2.3747039373786107478E-14, -45, 1.ABC9E3B39803F)
 ic(lnovfl, 7.0978271289338397310E2,     9, 1.62E42FEFA39EF)
 #ifdef vccast
 #define   mln2hi    vccast(mln2hi)

 @@ -1,14 +1,14 @@
-/*      $NetBSD: n_exp.c,v 1.8 2008/03/20 16:41:26 mhitch Exp $ */
+/*      $NetBSD: n_exp.c,v 1.9 2014/10/10 20:58:09 martin Exp $ */
 /*
  * Copyright (c) 1985, 1993
  *	The Regents of the University of California.  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.
  * 3. Neither the name of the University nor the names of its contributors
 @@ -72,26 +72,27 @@ static char sccsid[] = "@(#)exp.c	8.1 (B
  * 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 to produce the hexadecimal values
  * shown.
  */
 #define _LIBM_STATIC
 #include "../src/namespace.h"
 #include "mathimpl.h"
 #ifdef __weak_alias
 __weak_alias(exp, _exp);
 __weak_alias(_expl, _exp);
 __weak_alias(expf, _expf);
 #endif
 vc(ln2hi,  6.9314718055829871446E-1  ,7217,4031,0000,f7d0,   0, .B17217F7D00000)
 vc(ln2lo,  1.6465949582897081279E-12 ,bcd5,2ce7,d9cc,e4f1, -39, .E7BCD5E4F1D9CC)
 vc(lnhuge, 9.4961163736712506989E1   ,ec1d,43bd,9010,a73e,   7, .BDEC1DA73E9010)
 vc(lntiny,-9.5654310917272452386E1   ,4f01,c3bf,33af,d72e,   7,-.BF4F01D72E33AF)
 vc(invln2, 1.4426950408889634148E0   ,aa3b,40b8,17f1,295c,   1, .B8AA3B295C17F1)
 vc(p1,     1.6666666666666602251E-1  ,aaaa,3f2a,a9f1,aaaa,  -2, .AAAAAAAAAAA9F1)
 vc(p2,    -2.7777777777015591216E-3  ,0b60,bc36,ec94,b5f5,  -8,-.B60B60B5F5EC94)
 vc(p3,     6.6137563214379341918E-5  ,b355,398a,f15f,792e, -13, .8AB355792EF15F)
 vc(p4,    -1.6533902205465250480E-6  ,ea0e,b6dd,5f84,2e93, -19,-.DDEA0E2E935F84)
 vc(p5,     4.1381367970572387085E-8  ,bb4b,3431,2683,95f5, -24, .B1BB4B95F52683)

 @@ -1,14 +1,14 @@
-/*      $NetBSD: n_log.c,v 1.7 2008/03/20 16:41:26 mhitch Exp $ */
+/*      $NetBSD: n_log.c,v 1.8 2014/10/10 20:58:09 martin Exp $ */
 /*
  * Copyright (c) 1992, 1993
  *	The Regents of the University of California.  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.
  * 3. Neither the name of the University nor the names of its contributors
 @@ -33,26 +33,27 @@
 static char sccsid[] = "@(#)log.c	8.2 (Berkeley) 11/30/93";
 #endif
 #endif /* not lint */
 #include "../src/namespace.h"
 #include <math.h>
 #include <errno.h>
 #include "mathimpl.h"
 #ifdef __weak_alias
 __weak_alias(log, _log);
 __weak_alias(_logl, _log);
 __weak_alias(logf, _logf);
 #endif
 /* Table-driven natural logarithm.
+ *
  * This code was derived, with minor modifications, from:
  *	Peter Tang, "Table-Driven Implementation of the
  *	Logarithm in IEEE Floating-Point arithmetic." ACM Trans.
  *	Math Software, vol 16. no 4, pp 378-400, Dec 1990).
+ *
  * Calculates log(2^m*F*(1+f/F)), |f/j| <= 1/256,
  * where F = j/128 for j an integer in [0, 128].
+ *

 @@ -1,14 +1,14 @@
-/*      $NetBSD: n_sinh.c,v 1.7 2008/03/20 16:41:26 mhitch Exp $ */
+/*      $NetBSD: n_sinh.c,v 1.8 2014/10/10 20:58:09 martin Exp $ */
 /*
  * Copyright (c) 1985, 1993
  *	The Regents of the University of California.  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.
  * 3. Neither the name of the University nor the names of its contributors
 @@ -70,26 +70,27 @@ static char sccsid[] = "@(#)sinh.c	8.1 (
  * 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 to produce the hexadecimal values
  * shown.
  */
 #define _LIBM_STATIC
 #include "../src/namespace.h"
 #include "mathimpl.h"
 #ifdef __weak_alias
 __weak_alias(sinh, _sinh);
 __weak_alias(_sinhl, _sinh);
 __weak_alias(sinhf, _sinhf);
 #endif
 vc(mln2hi, 8.8029691931113054792E1   ,0f33,43b0,2bdb,c7e2,   7, .B00F33C7E22BDB)
 vc(mln2lo,-4.9650192275318476525E-16 ,1b60,a70f,582a,279e, -50,-.8F1B60279E582A)
 vc(lnovfl, 8.8029691931113053016E1   ,0f33,43b0,2bda,c7e2,   7, .B00F33C7E22BDA)
 ic(mln2hi, 7.0978271289338397310E2,    10, 1.62E42FEFA39EF)
 ic(mln2lo, 2.3747039373786107478E-14, -45, 1.ABC9E3B39803F)
 ic(lnovfl, 7.0978271289338397310E2,     9, 1.62E42FEFA39EF)
 #ifdef vccast
 #define	mln2hi	vccast(mln2hi)

 @@ -1,14 +1,14 @@
-/*      $NetBSD: n_pow.c,v 1.9 2013/11/24 14:46:18 martin Exp $ */
+/*      $NetBSD: n_pow.c,v 1.10 2014/10/10 20:58:09 martin Exp $ */
 /*
  * Copyright (c) 1985, 1993
  *	The Regents of the University of California.  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.
  * 3. Neither the name of the University nor the names of its contributors
 @@ -112,38 +112,38 @@ static char sccsid[] = "@(#)pow.c	8.1 (B
 #define TRUNC(x)	x = (double) (float) x
 #define _IEEE		0
 #else
 #define _IEEE		1
 #define endian		(((*(int *) &one)) ? 1 : 0)
 #define TRUNC(x) 	*(((int *) &x)+endian) &= 0xf8000000
 #define infnan(x)	0.0
 #endif		/* __vax__ or tahoe */
 static const double zero=0.0, one=1.0, two=2.0, negone= -1.0;
 static double pow_P (double, double);
 #ifdef __weak_alias
 __weak_alias(_powf, powf);
 __weak_alias(_pow, pow);
 __weak_alias(_powl, pow);
 #endif
 float
 powf(float x, float y)
+{
    return pow((double) x, (double) (y));
+}
 long double
 powl(long double x, long double y)
    return pow((double) x, (double) (y));
 double
 pow(double x, double y)
+{
 	double t;
 	if (y==zero)
 		return (one);
 	else if (y==one || (_IEEE && x != x))
 		return (x);		/* if x is NaN or y=1 */
 	else if (_IEEE && y!=y)		/* if y is NaN */
 		return (y);
 	else if (!finite(y))		/* if y is INF */
 		if ((t=fabs(x))==one)	/* +-1 ** +-INF is NaN */
 			return (y - y);