newlib/libm/test/test_ieee.c

   1
   2 #include "test.h"
   3 #include <ieeefp.h>
   4
   5
   6 /* Test fp getround and fp setround */
   7
   8 void
   9 test_getround (void)
  10 {
  11
  12   newfunc("fpgetround/fpsetround");
  13   line(1);
  14   fpsetround(FP_RN);
  15   test_iok(fpgetround(), FP_RN);
  16   line(2);
  17   fpsetround(FP_RM);
  18   test_iok(fpgetround(), FP_RM);
  19   line(3);
  20   fpsetround(FP_RP);
  21   test_iok(fpgetround(), FP_RP);
  22   line(4);
  23   fpsetround(FP_RZ);
  24   test_iok(fpgetround(), FP_RZ);
  25 }
  26
  27 /* And fpset/fpgetmask */
  28 void
  29 test_getmask (void)
  30 {
  31   newfunc("fpsetmask/fpgetmask");
  32   line(1);
  33   fpsetmask(FP_X_INV);
  34   test_iok(fpgetmask(),FP_X_INV);
  35   line(2);
  36   fpsetmask(FP_X_DX);
  37   test_iok(fpgetmask(),FP_X_DX);
  38   line(3);
  39   fpsetmask(FP_X_OFL );
  40   test_iok(fpgetmask(),FP_X_OFL);
  41   line(4);
  42   fpsetmask(FP_X_UFL);
  43   test_iok(fpgetmask(),FP_X_UFL);
  44   line(5);
  45   fpsetmask(FP_X_IMP);
  46   test_iok(fpgetmask(),FP_X_IMP);
  47 }
  48
  49 void
  50 test_getsticky (void)
  51 {
  52   newfunc("fpsetsticky/fpgetsticky");
  53   line(1);
  54   fpsetsticky(FP_X_INV);
  55   test_iok(fpgetsticky(),FP_X_INV);
  56   line(2);
  57   fpsetsticky(FP_X_DX);
  58   test_iok(fpgetsticky(),FP_X_DX);
  59   line(3);
  60   fpsetsticky(FP_X_OFL );
  61   test_iok(fpgetsticky(),FP_X_OFL);
  62   line(4);
  63   fpsetsticky(FP_X_UFL);
  64   test_iok(fpgetsticky(),FP_X_UFL);
  65   line(5);
  66   fpsetsticky(FP_X_IMP);
  67   test_iok(fpgetsticky(),FP_X_IMP);
  68 }
  69
  70 void
  71 test_getroundtoi (void)
  72 {
  73   newfunc("fpsetroundtoi/fpgetroundtoi");
  74   line(1);
  75   fpsetroundtoi(FP_RDI_TOZ);
  76   test_iok(fpgetroundtoi(),FP_RDI_TOZ);
  77
  78   line(2);
  79   fpsetroundtoi(FP_RDI_RD);
  80   test_iok(fpgetroundtoi(),FP_RDI_RD);
  81
  82 }
  83
  84 double
  85  dnumber (int msw,
  86         int lsw)
  87 {
  88
  89   __ieee_double_shape_type v;
  90   v.parts.lsw = lsw;
  91   v.parts.msw = msw;
  92   return v.value;
  93 }
  94
  95   /* Lets see if changing the rounding alters the arithmetic.
  96      Test by creating numbers which will have to be rounded when
  97      added, and seeing what happens to them */
  98  /* Keep them out here to stop  the compiler from folding the results */
  99 double n;
 100 double m;
 101 double add_rounded_up;
 102 double add_rounded_down;
 103 double sub_rounded_down ;
 104 double sub_rounded_up ;
 105   double r1,r2,r3,r4;
 106 void
 107 test_round (void)
 108 {
 109   n =                dnumber(0x40000000, 0x00000008); /* near 2 */
 110   m =                dnumber(0x40400000, 0x00000003); /* near 3.4 */
 111
 112   add_rounded_up   = dnumber(0x40410000, 0x00000004); /* For RN, RP */
 113   add_rounded_down = dnumber(0x40410000, 0x00000003); /* For RM, RZ */
 114   sub_rounded_down = dnumber(0xc0410000, 0x00000004); /* for RN, RM */
 115   sub_rounded_up   = dnumber(0xc0410000, 0x00000003); /* for RP, RZ */
 116
 117   newfunc("fpsetround");
 118
 119   line(1);
 120
 121   fpsetround(FP_RN);
 122   r1 = n + m;
 123   test_mok(r1, add_rounded_up, 64);
 124
 125   line(2);
 126   fpsetround(FP_RM);
 127   r2 = n + m;
 128   test_mok(r2, add_rounded_down, 64);
 129
 130   fpsetround(FP_RP);
 131   line(3);
 132   r3 = n + m;
 133   test_mok(r3,add_rounded_up, 64);
 134
 135   fpsetround(FP_RZ);
 136   line(4);
 137   r4 = n + m;
 138   test_mok(r4,add_rounded_down,64);
 139
 140
 141   fpsetround(FP_RN);
 142   r1 = - n - m;
 143   line(5);
 144   test_mok(r1,sub_rounded_down,64);
 145
 146   fpsetround(FP_RM);
 147   r2 = - n - m;
 148   line(6);
 149   test_mok(r2,sub_rounded_down,64);
 150
 151
 152   fpsetround(FP_RP);
 153   r3 = - n - m;
 154   line(7);
 155   test_mok(r3,sub_rounded_up,64);
 156
 157   fpsetround(FP_RZ);
 158   r4 = - n - m;
 159   line(8);
 160   test_mok(r4,sub_rounded_up,64);
 161 }
 162
 163
 164 void
 165 test_ieee (void)
 166 {
 167   fp_rnd old = fpgetround();
 168   test_getround();
 169   test_getmask();
 170   test_getsticky();
 171   test_getroundtoi();
 172
 173   test_round();
 174   fpsetround(old);
 175
 176
 177 }
 178
 179