fixup #8538 -- correct rounding modes for floor, trunc

This fixes the rounding modes used for the floor and trunc C library functions.

Author: kroening

regression/cbmc-library/floor-01/main.c

@@ -3,7 +3,10 @@
 
 int main()
 {
-  floor();
-  assert(0);
+  assert(floor(1.1)==1.0);
+  assert(floor(1.9)==1.0);
+  assert(floor(-1.1)==-2.0);
+  assert(floor(-1.9)==-2.0);
+  assert(signbit(floor(-0.0)));
   return 0;
 }

regression/cbmc-library/floor-01/test.desc

@@ -1,6 +1,6 @@
-KNOWNBUG
+CORE
 main.c
---pointer-check --bounds-check
+
 ^EXIT=0$
 ^SIGNAL=0$
 ^VERIFICATION SUCCESSFUL$

regression/cbmc-library/floorf-01/main.c

@@ -3,7 +3,10 @@
 
 int main()
 {
-  floorf();
-  assert(0);
+  assert(floorf(1.1f)==1.0f);
+  assert(floorf(1.9f)==1.0f);
+  assert(floorf(-1.1f)==-2.0f);
+  assert(floorf(-1.9f)==-2.0f);
+  assert(signbit(floorf(-0.0f)));
   return 0;
 }

regression/cbmc-library/floorf-01/test.desc

@@ -1,6 +1,6 @@
-KNOWNBUG
+CORE
 main.c
---pointer-check --bounds-check
+
 ^EXIT=0$
 ^SIGNAL=0$
 ^VERIFICATION SUCCESSFUL$

regression/cbmc-library/floorl-01/main.c

@@ -3,7 +3,10 @@
 
 int main()
 {
-  floorl();
-  assert(0);
+  assert(floorl(1.1l)==1.0l);
+  assert(floorl(1.9l)==1.0l);
+  assert(floorl(-1.1l)==-2.0l);
+  assert(floorl(-1.9l)==-2.0l);
+  assert(signbit(floorl(-0.0l)));
   return 0;
 }

regression/cbmc-library/floorl-01/test.desc

@@ -1,6 +1,6 @@
-KNOWNBUG
+CORE
 main.c
---pointer-check --bounds-check
+
 ^EXIT=0$
 ^SIGNAL=0$
 ^VERIFICATION SUCCESSFUL$

regression/cbmc-library/round-01/main.c

@@ -3,7 +3,12 @@
 
 int main()
 {
-  round();
-  assert(0);
+  assert(round(1.1)==1.0);
+  assert(round(1.5)==2.0);
+  assert(round(1.9)==2.0);
+  assert(round(-1.1)==-1.0);
+  assert(round(-1.5)==-2.0);
+  assert(round(-1.9)==-2.0);
+  assert(signbit(round(-0.0)));
   return 0;
 }

regression/cbmc-library/round-01/test.desc

@@ -1,6 +1,6 @@
-KNOWNBUG
+CORE
 main.c
---pointer-check --bounds-check
+
 ^EXIT=0$
 ^SIGNAL=0$
 ^VERIFICATION SUCCESSFUL$

regression/cbmc-library/roundf-01/main.c

@@ -3,7 +3,12 @@
 
 int main()
 {
-  roundf();
-  assert(0);
+  assert(roundf(1.1f)==1.0f);
+  assert(roundf(1.5f)==2.0f);
+  assert(roundf(1.9f)==2.0f);
+  assert(roundf(-1.1f)==-1.0f);
+  assert(roundf(-1.5f)==-2.0f);
+  assert(roundf(-1.9f)==-2.0f);
+  assert(signbit(roundf(-0.0f)));
   return 0;
 }

regression/cbmc-library/roundf-01/test.desc

@@ -1,6 +1,6 @@
-KNOWNBUG
+CORE
 main.c
---pointer-check --bounds-check
+
 ^EXIT=0$
 ^SIGNAL=0$
 ^VERIFICATION SUCCESSFUL$

regression/cbmc-library/roundl-01/main.c

@@ -3,7 +3,12 @@
 
 int main()
 {
-  roundl();
-  assert(0);
+  assert(roundl(1.1l)==1.0l);
+  assert(roundl(1.5l)==2.0l);
+  assert(roundl(1.9l)==2.0l);
+  assert(roundl(-1.1l)==-1.0l);
+  assert(roundl(-1.5l)==-2.0l);
+  assert(roundl(-1.9l)==-2.0l);
+  assert(signbit(roundl(-0.0l)));
   return 0;
 }

regression/cbmc-library/roundl-01/test.desc

@@ -1,6 +1,6 @@
-KNOWNBUG
+CORE
 main.c
---pointer-check --bounds-check
+
 ^EXIT=0$
 ^SIGNAL=0$
 ^VERIFICATION SUCCESSFUL$

regression/cbmc-library/trunc-01/main.c

@@ -3,7 +3,10 @@
 
 int main()
 {
-  trunc();
-  assert(0);
+  assert(trunc(1.1)==1.0);
+  assert(trunc(1.9)==1.0);
+  assert(trunc(-1.1)==-1.0);
+  assert(trunc(-1.9)==-1.0);
+  assert(signbit(trunc(-0.0)));
   return 0;
 }

regression/cbmc-library/trunc-01/test.desc

@@ -1,6 +1,6 @@
-KNOWNBUG
+CORE
 main.c
---pointer-check --bounds-check
+
 ^EXIT=0$
 ^SIGNAL=0$
 ^VERIFICATION SUCCESSFUL$

regression/cbmc-library/truncf-01/main.c

@@ -3,7 +3,10 @@
 
 int main()
 {
-  truncf();
-  assert(0);
+  assert(truncf(1.1f)==1.0f);
+  assert(truncf(1.9f)==1.0f);
+  assert(truncf(-1.1f)==-1.0f);
+  assert(truncf(-1.9f)==-1.0f);
+  assert(signbit(trunc(-0.0f)));
   return 0;
 }

regression/cbmc-library/truncf-01/test.desc

@@ -1,6 +1,6 @@
-KNOWNBUG
+CORE
 main.c
---pointer-check --bounds-check
+
 ^EXIT=0$
 ^SIGNAL=0$
 ^VERIFICATION SUCCESSFUL$

regression/cbmc-library/truncl-01/main.c

@@ -3,7 +3,10 @@
 
 int main()
 {
-  truncl();
-  assert(0);
+  assert(truncl(1.1l)==1.0l);
+  assert(truncl(1.9l)==1.0l);
+  assert(truncl(-1.1l)==-1.0l);
+  assert(truncl(-1.9l)==-1.0l);
+  assert(signbit(truncl(-0.0l)));
   return 0;
 }

regression/cbmc-library/truncl-01/test.desc

@@ -1,6 +1,6 @@
-KNOWNBUG
+CORE
 main.c
---pointer-check --bounds-check
+
 ^EXIT=0$
 ^SIGNAL=0$
 ^VERIFICATION SUCCESSFUL$

src/ansi-c/library/math.c

@@ -1322,7 +1322,7 @@ long double ceill(long double x)
 
 double floor(double x)
 {
-  return __CPROVER_round_to_integrald(x, 3); // FE_DOWNWARD
+  return __CPROVER_round_to_integrald(x, 1); // FE_DOWNWARD
 }
 
 /* FUNCTION: floorf */
@@ -1339,7 +1339,7 @@ double floor(double x)
 
 float floorf(float x)
 {
-  return __CPROVER_round_to_integralf(x, 3); // FE_DOWNWARD
+  return __CPROVER_round_to_integralf(x, 1); // FE_DOWNWARD
 }
 
 
@@ -1357,7 +1357,7 @@ float floorf(float x)
 
 long double floorl(long double x)
 {
-  return __CPROVER_round_to_integralld(x, 3); // FE_DOWNWARD
+  return __CPROVER_round_to_integralld(x, 1); // FE_DOWNWARD
 }
 
 
@@ -1381,7 +1381,7 @@ long double floorl(long double x)
 
 double trunc(double x)
 {
-  return __CPROVER_round_to_integrald(x, 0); // FE_TOWARDZERO
+  return __CPROVER_round_to_integrald(x, 3); // FE_TOWARDZERO
 }
 
 /* FUNCTION: truncf */
@@ -1398,7 +1398,7 @@ double trunc(double x)
 
 float truncf(float x)
 {
-  return __CPROVER_round_to_integralf(x, 0); // FE_TOWARDZERO
+  return __CPROVER_round_to_integralf(x, 3); // FE_TOWARDZERO
 }
 
 
@@ -1416,7 +1416,7 @@ float truncf(float x)
 
 long double truncl(long double x)
 {
-  return __CPROVER_round_to_integralld(x, 0); // FE_TOWARDZERO
+  return __CPROVER_round_to_integralld(x, 3); // FE_TOWARDZERO
 }
 
 
@@ -1889,7 +1889,7 @@ long long int llroundl(long double x)
 
 double modf(double x, double *iptr)
 {
-  *iptr = __CPROVER_round_to_integrald(x, 0); // FE_TOWARDZERO
+  *iptr = __CPROVER_round_to_integrald(x, 3); // FE_TOWARDZERO
   return (x - *iptr);
 }
 
@@ -1907,7 +1907,7 @@ double modf(double x, double *iptr)
 
 float modff(float x, float *iptr)
 {
-  *iptr = __CPROVER_round_to_integralf(x, 0); // FE_TOWARDZERO
+  *iptr = __CPROVER_round_to_integralf(x, 3); // FE_TOWARDZERO
   return (x - *iptr);
 }
 
@@ -1926,7 +1926,7 @@ float modff(float x, float *iptr)
 
 long double modfl(long double x, long double *iptr)
 {
-  *iptr = __CPROVER_round_to_integralld(x, 0); // FE_TOWARDZERO
+  *iptr = __CPROVER_round_to_integralld(x, 3); // FE_TOWARDZERO
   return (x - *iptr);
 }