static const double huge = 1.0e300;

double ceil(double x)

{

int32_t i0,i1,j0;

uint32_t i,j;

EXTRACT_WORDS(i0, i1, x);

// FIXME signed shift

j0 = ((i0>>20)&0x7ff) - 0x3ff;

if (j0 < 20) {

if (j0 < 0) {

/* raise inexact if x != 0 */

if (huge+x > 0.0) {

if (i0 < 0) {

i0 = 0x80000000;

i1=0;

} else if ((i0|i1) != 0) {

i0=0x3ff00000;

i1=0;

}

}

} else {

i = 0x000fffff>>j0;

if (((i0&i)|i1) == 0) /* x is integral */

return x;

/* raise inexact flag */

if (huge+x > 0.0) {

if (i0 > 0)

i0 += 0x00100000>>j0;

i0 &= ~i;

i1 = 0;

}

}

} else if (j0 > 51) {

if (j0 == 0x400)  /* inf or NaN */

return x+x;

return x;         /* x is integral */

} else {

i = (uint32_t)0xffffffff>>(j0-20);

if ((i1&i) == 0)

return x; /* x is integral */

/* raise inexact flag */

if (huge+x > 0.0) {

if (i0 > 0) {

if (j0 == 20)

i0 += 1;

else {

j = i1 + (1<<(52-j0));

if (j < i1)  /* got a carry */

i0 += 1;

i1 = j;

}

}

i1 &= ~i;

}

}

INSERT_WORDS(x, i0, i1);

return x;

}‍

------------------

------------------ 原始邮件 ------------------

Test it in qemu-system-arm, with soft fp, realview-pb-a8 platform‍

we porting musl to our project, and just write test functions to test cases, i don't kown, how to test musl in origin test cases.

 

On Tue, Oct 14, 2014 at 10:17:34AM +0800, bobodog wrote:
> HI,ALL
>          We use musl in arm process, but found some bugs about math library,  x86 process is ok!

I can't reproduce this issue. Are you using soft (arm) or hard (armhf)
float based target?

>         The test code:
> int test(int argc, char ** argv)
> {
> double n = 0.0;
> int i;
>
>
> for(i = 0; i < 20; i++)
> {
> n = n - 0.125;
> printf("ceil(%f) = %f, floor(%f) = %f, round(%f) = %f\n", n, ceil(n), n, floor(n), n, round(n));
> }
>
>
> return 0;
> }‍

Is this the whole test program? I don't see a main function or any
#include directives. Compiling it with gcc -Dtest=main results in:

foo.c:10:17: warning: incompatible implicit declaration of built-in function ‘printf’ [enabled by default]
foo.c:10:78: warning: incompatible implicit declaration of built-in function ‘ceil’ [enabled by default]
foo.c:10:90: warning: incompatible implicit declaration of built-in function ‘floor’ [enabled by default]
foo.c:10:103: warning: incompatible implicit declaration of built-in function ‘round’ [enabled by default]

These warnings should all be treated as errors; calling a function
with variadic arguments (printf) or a function that returns a floating
point type without a valid declaration/prototype is a serious error
and will not work.

Rich