l4/pkg/valgrind/src/valgrind-3.6.0-svn/none/tests/ppc32/power5+_round.c

   1 /*  Copyright (C) 2007 IBM
   2
   3     Author: Pete Eberlein  eberlein@us.ibm.com
   4
   5     This program is free software; you can redistribute it and/or
   6     modify it under the terms of the GNU General Public License as
   7     published by the Free Software Foundation; either version 2 of the
   8     License, or (at your option) any later version.
   9
  10     This program is distributed in the hope that it will be useful, but
  11     WITHOUT ANY WARRANTY; without even the implied warranty of
  12     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  13     General Public License for more details.
  14
  15     You should have received a copy of the GNU General Public License
  16     along with this program; if not, write to the Free Software
  17     Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
  18     02111-1307, USA.
  19
  20     The GNU General Public License is contained in the file COPYING.
  21 */
  22
  23
  24
  25 #include <stdio.h>
  26 #include <stdlib.h>
  27 #include <strings.h>
  28
  29 #define POS_NORMAL   0x4000
  30 #define NEG_NORMAL   0x8000
  31 #define POS_INF      0x5000
  32 #define NEG_INF      0x9000
  33 #define POS_ZERO     0x2000
  34 #define NEG_ZERO     0x12000
  35 #define POS_DENORMAL 0x14000
  36 #define NEG_DENORMAL 0x18000
  37 #define NAN          0x11000
  38 #define FPRF_MASK    0x1F000
  39
  40
  41 int main(int argc, char *argv[])
  42 {
  43
  44    double inf, neg0, nan;
  45    union {
  46       double d;
  47       struct {
  48          unsigned int dummy, dummy2:15, fprf:17;
  49       };
  50    } fpscr;
  51
  52    inf = strtod("inf", NULL);
  53    neg0 = strtod("-0", NULL);
  54    nan = strtod("nan", NULL);
  55
  56
  57    /* This set is disabled until fprf is implemented. */
  58    if (0) {
  59       double set[] = { inf, 1.5, 0, neg0, -1.5, -inf, nan };
  60       int i, j, fprf;
  61       for (i = 0; i < 7; ++i) {
  62          for (j = 0; j < 7; ++j) {
  63           asm("fcmpu 1, %1, %2\n\t" "mffs %0\n":"=f"(fpscr.d)
  64           :    "f"(set[i]), "f"(set[j])
  65                 );
  66
  67             if (i == 6 || j == 6) {
  68                fprf = 0x1000;   // Unordered
  69             } else if (i == j || (i == 2 && j == 3) || (i == 3 && j == 2)) {
  70                fprf = 0x2000;   // Equal
  71             } else if (i < j) {
  72                fprf = 0x4000;   // Greater Than
  73             } else if (i > j) {
  74                fprf = 0x8000;   // Less Than
  75             }
  76
  77             printf("fcmpu\t%.1f\t%.1f\t%x\t%s\n", set[i], set[j],
  78                    fpscr.fprf, fpscr.fprf == fprf ? "PASS" : "FAIL");
  79          }
  80       }
  81    }
  82
  83    {
  84       double set[] = { inf, 1.9, 1.1, 0, neg0, -1.1, -1.9, -inf, nan };
  85       double frin[] = { inf, 2.0, 1.0, 0, neg0, -1.0, -2.0, -inf, nan };
  86       double friz[] = { inf, 1.0, 1.0, 0, neg0, -1.0, -1.0, -inf, nan };
  87       double frip[] = { inf, 2.0, 2.0, 0, neg0, -1.0, -1.0, -inf, nan };
  88       double frim[] = { inf, 1.0, 1.0, 0, neg0, -2.0, -2.0, -inf, nan };
  89       int fprf[] = { POS_INF, POS_NORMAL, POS_NORMAL, POS_ZERO, NEG_ZERO,
  90          NEG_NORMAL, NEG_NORMAL, NEG_INF, NAN
  91       };
  92       double set2[] = { 0.9, 0.1, -0.1, -0.9, 1e-40, -1e-40 };
  93       double frin2[] = { 1.0, 0.0, -0.0, -1.0, 0.0, -0.0 };
  94       int frin2rf[] =
  95           { POS_NORMAL, POS_ZERO, NEG_ZERO, NEG_NORMAL, POS_ZERO,
  96  NEG_ZERO };
  97       double friz2[] = { 0.0, 0.0, -0.0, -0.0, 0.0, -0.0 };
  98       int friz2rf[] =
  99           { POS_ZERO, POS_ZERO, NEG_ZERO, NEG_ZERO, POS_ZERO, NEG_ZERO };
 100       double frip2[] = { 1.0, 1.0, -0.0, -0.0, 1.0, -0.0 };
 101       int frip2rf[] =
 102           { POS_NORMAL, POS_NORMAL, NEG_ZERO, NEG_ZERO, POS_NORMAL,
 103  NEG_ZERO };
 104       double frim2[] = { 0.0, 0.0, -1.0, -1.0, 0.0, -1.0 };
 105       int frim2rf[] =
 106           { POS_ZERO, POS_ZERO, NEG_NORMAL, NEG_NORMAL, POS_ZERO,
 107  NEG_NORMAL };
 108       double ret;
 109       int i;
 110
 111 #define DO_TEST(op,in,out,rf)  for (i=0; i<sizeof(in)/sizeof(double); ++i) { \
 112     asm (#op" %0, %2\n\t" \
 113        "mffs %1\n" \
 114     : "=f" (ret), "=f" (fpscr.d) \
 115     : "f" (in[i]) \
 116     ); \
 117     printf(#op"\t%g\t%g\t%x\t%s\n", in[i], ret, fpscr.fprf, \
 118            (!bcmp(&ret, &out[i], sizeof(double))) /*&& (rf[i] == fpscr.fprf)*/ \
 119            ? "PASS" : "FAIL"); \
 120   }
 121       /* Note: fprf check above is disabled until fprf is implemented. */
 122
 123
 124       DO_TEST(frin, set, frin, fprf);
 125       DO_TEST(frin, set2, frin2, frin2rf);
 126       DO_TEST(friz, set, friz, fprf);
 127       DO_TEST(friz, set2, friz2, friz2rf);
 128       DO_TEST(frip, set, frip, fprf);
 129       DO_TEST(frip, set2, frip2, frip2rf);
 130       DO_TEST(frim, set, frim, fprf);
 131       DO_TEST(frim, set2, frim2, frim2rf);
 132    }
 133
 134    /* This set is disabled until fprf is implemented. */
 135    if (0) {
 136       double set1[] = { inf, 0.9, 0.1, 0, neg0, -0.1, -0.9, -inf, nan };
 137       double frsp1[] =
 138           { inf, 0.9f, 0.1f, 0, neg0, -0.1f, -0.9f, -inf, nan };
 139       int fprf1[] =
 140           { POS_INF, POS_NORMAL, POS_NORMAL, POS_ZERO, NEG_ZERO,
 141    NEG_NORMAL,
 142          NEG_NORMAL, NEG_INF, NAN
 143       };
 144       double set2[] =
 145           { 1.2e-38, 1.1e-38, 1e-40, 8e-44, 9e-44, 8e-46, 7e-46 };
 146       double frsp2[] =
 147           { 1.2e-38f, 1.1e-38f, 1e-40f, 8e-44f, 9e-44f, 8e-46f, 0.0 };
 148       int fprf2[] = { POS_NORMAL, POS_DENORMAL, POS_DENORMAL, POS_DENORMAL,
 149          POS_DENORMAL, POS_DENORMAL, POS_ZERO
 150       };
 151       double set3[] =
 152           { -1.2e-38, -1.1e-38, -1e-40, -8e-44, -9e-44, -8e-46, -7e-46 };
 153       double frsp3[] =
 154           { -1.2e-38f, -1.1e-38f, -1e-40f, -8e-44f, -9e-44f, -8e-46f,
 155 -0.0 };
 156       int fprf3[] = { NEG_NORMAL, NEG_DENORMAL, NEG_DENORMAL, NEG_DENORMAL,
 157          NEG_DENORMAL, NEG_DENORMAL, NEG_ZERO
 158       };
 159       double ret;
 160       int i;
 161       DO_TEST(frsp, set1, frsp1, fprf1);
 162       DO_TEST(frsp, set2, frsp2, fprf2);
 163       DO_TEST(frsp, set3, frsp3, fprf3);
 164    }
 165
 166
 167    return 0;
 168 }