Aren/haiku
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

adding POSIX-functions llround*(), which Maxime reported missing

Browse Source 
* copied implementations for llround(), llroundf() and llroundl() from
  glibc-2.3.2
* added corresponding declarations to math.h


git-svn-id: file:///srv/svn/repos/haiku/haiku/trunk@33863 a95241bf-73f2-0310-859d-f6bbb57e9c96
This commit is contained in:
Oliver Tappe 2009-11-02 19:13:15 +00:00
parent dd9740dad4
commit 7e0d60f5b7
5 changed files with 229 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
headers/posix/math.h
View File

@ -98,6 +98,7 @@ extern float hypotf(float x, float y);
extern float ldexpf(float x, int exponent);
extern float lgammaf(float x);
extern long long llrintf(float x);
extern long long llroundf(float x);
extern float log10f(float x);
extern float log1pf(float x);
extern float logbf(float x);
@ -135,6 +136,7 @@ extern double gamma(double x);
extern double ldexp(double x, int exponent);
extern double lgamma(double x);
extern long long llrint(double x);
extern long long llround(double x);
extern double log(double x);
extern double log10(double x);
extern long lrint(double x);
@ -161,6 +163,7 @@ extern long double atan2l(long double y, long double x);
extern long double lgammal(long double x);
extern long double nearbyintl(long double x);
extern long long llrintl(long double x);
extern long long llroundl(long double x);
extern long lrintl(long double x);
extern long lroundl(long double x);
extern long double remquol(long double x, long double y, int *quo);

81
src/system/libroot/posix/glibc/arch/generic/s_llround.c Normal file
View File

@ -0,0 +1,81 @@
/* Round double value to long long int.
Copyright (C) 1997 Free Software Foundation, Inc.
This file is part of the GNU C Library.
Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
The GNU C Library is free software; you can redistribute it and/or
modify it under the terms of the GNU Lesser General Public
License as published by the Free Software Foundation; either
version 2.1 of the License, or (at your option) any later version.
The GNU C Library is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
Lesser General Public License for more details.
You should have received a copy of the GNU Lesser General Public
License along with the GNU C Library; if not, write to the Free
Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
02111-1307 USA. */
#include <math.h>
#include "math_private.h"
long long int
__llround (double x)
{
int32_t j0;
u_int32_t i1, i0;
long long int result;
int sign;
EXTRACT_WORDS (i0, i1, x);
j0 = ((i0 >> 20) & 0x7ff) - 0x3ff;
sign = (i0 & 0x80000000) != 0 ? -1 : 1;
i0 &= 0xfffff;
i0 |= 0x100000;
if (j0 < 20)
{
if (j0 < 0)
return j0 < -1 ? 0 : sign;
else
{
i0 += 0x80000 >> j0;
result = i0 >> (20 - j0);
}
}
else if (j0 < (int32_t) (8 * sizeof (long long int)) - 1)
{
if (j0 >= 52)
result = (((long long int) i0 << 32) | i1) << (j0 - 52);
else
{
u_int32_t j = i1 + (0x80000000 >> (j0 - 20));
if (j < i1)
++i0;
if (j0 == 20)
result = (long long int) i0;
else
result = ((long long int) i0 << (j0 - 20)) | (j >> (52 - j0));
}
}
else
{
/* The number is too large. It is left implementation defined
what happens. */
return (long long int) x;
}
return sign * result;
}
weak_alias (__llround, llround)
#ifdef NO_LONG_DOUBLE
strong_alias (__llround, __llroundl)
weak_alias (__llround, llroundl)
#endif

63
src/system/libroot/posix/glibc/arch/generic/s_llroundf.c Normal file
View File

@ -0,0 +1,63 @@
/* Round float value to long long int.
Copyright (C) 1997 Free Software Foundation, Inc.
This file is part of the GNU C Library.
Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
The GNU C Library is free software; you can redistribute it and/or
modify it under the terms of the GNU Lesser General Public
License as published by the Free Software Foundation; either
version 2.1 of the License, or (at your option) any later version.
The GNU C Library is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
Lesser General Public License for more details.
You should have received a copy of the GNU Lesser General Public
License along with the GNU C Library; if not, write to the Free
Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
02111-1307 USA. */
#include <math.h>
#include "math_private.h"
long long int
__llroundf (float x)
{
int32_t j0;
u_int32_t i;
long long int result;
int sign;
GET_FLOAT_WORD (i, x);
j0 = ((i >> 23) & 0xff) - 0x7f;
sign = (i & 0x80000000) != 0 ? -1 : 1;
i &= 0x7fffff;
i |= 0x800000;
if (j0 < (int32_t) (8 * sizeof (long long int)) - 1)
{
if (j0 < 0)
return j0 < -1 ? 0 : sign;
else if (j0 >= 23)
result = (long long int) i << (j0 - 23);
else
{
i += 0x400000 >> j0;
result = i >> (23 - j0);
}
}
else
{
/* The number is too large. It is left implementation defined
what happens. */
return (long long int) x;
}
return sign * result;
}
weak_alias (__llroundf, llroundf)

81
src/system/libroot/posix/glibc/arch/generic/s_llroundl.c Normal file
View File

@ -0,0 +1,81 @@
/* Round long double value to long long int.
Copyright (C) 1997, 2001 Free Software Foundation, Inc.
This file is part of the GNU C Library.
Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
The GNU C Library is free software; you can redistribute it and/or
modify it under the terms of the GNU Lesser General Public
License as published by the Free Software Foundation; either
version 2.1 of the License, or (at your option) any later version.
The GNU C Library is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
Lesser General Public License for more details.
You should have received a copy of the GNU Lesser General Public
License along with the GNU C Library; if not, write to the Free
Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
02111-1307 USA. */
#include <math.h>
#include "math_private.h"
long long int
__llroundl (long double x)
{
int32_t j0;
u_int32_t se, i1, i0;
long long int result;
int sign;
GET_LDOUBLE_WORDS (se, i0, i1, x);
j0 = (se & 0x7fff) - 0x3fff;
sign = (se & 0x8000) != 0 ? -1 : 1;
if (j0 < 31)
{
if (j0 < 0)
return j0 < -1 ? 0 : sign;
else
{
u_int32_t j = i0 + (0x40000000 >> j0);
if (j < i0)
{
j >>= 1;
j |= 0x80000000;
++j0;
}
result = j >> (31 - j0);
}
}
else if (j0 < (int32_t) (8 * sizeof (long long int)) - 1)
{
if (j0 >= 63)
result = (((long long int) i0 << 32) | i1) << (j0 - 63);
else
{
u_int32_t j = i1 + (0x80000000 >> (j0 - 31));
result = (long long int) i0;
if (j < i1)
++result;
if (j0 > 31)
result = (result << (j0 - 31)) | (j >> (63 - j0));
}
}
else
{
/* The number is too large. It is left implementation defined
what happens. */
return (long long int) x;
}
return sign * result;
}
weak_alias (__llroundl, llroundl)

1
src/system/libroot/posix/glibc/arch/x86/Jamfile
View File

@ -58,6 +58,7 @@ local genericSources =
s_isinf.c s_isinff.c
s_isnan.c s_isnanf.c
s_ldexp.c s_ldexpf.c #s_ldexpl.c
s_llround.c s_llroundf.c s_llroundl.c
s_lround.c s_lroundf.c s_lroundl.c
s_matherr.c
s_modf.c s_modff.c # s_modfl.c