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