1 files changed, 167 insertions, 0 deletions
diff --git a/private/fp32/tran/mips/powt.c b/private/fp32/tran/mips/powt.c
new file mode 100644
index 000000000..9609a7ea8
--- /dev/null
+++ b/private/fp32/tran/mips/powt.c
@@ -0,0 +1,167 @@
+#include <stdio.h>
+#include <math.h>
+#include <errno.h>
+
+int test1(void);
+int test2(void);
+
+int main(int argc, char **argv)
+{
+	int k = 0;
+
+	k += test1();
+
+	k += test2();
+
+	if (k) {
+		printf("\n\tFailed %d tests...\n", k);
+	} else {
+		printf("\n\tPassed all tests...\n", k);
+	}
+
+	return k;
+}
+
+int test1()
+{
+    double x, y, result, answer;
+    char str[80];
+    int i;
+    int k = 0;
+
+    x = y = 0.0;
+    answer = 1.0;
+    result = pow(x,y);
+
+    if (result != answer) {
+        printf("pow(%g,%g) = %g, should be %g\n", x, y, result, answer);
+	k++;
+    }
+
+    x = 0.0;
+    y = -1.0;
+    result = pow(x,y);
+
+    sprintf(str, "%g", result);
+    if (strcmp(str, "1.#INF")) {
+        printf("pow(%g,%g) = %g, should be %s\n", x, y, result, "1.#INF");
+	k++;
+    }
+
+    x = -1.1e300;
+    y = 21.0;
+    result = pow(x,y);
+
+    sprintf(str, "%le", result);
+    if (strcmp(str, "-1.#INF00e+000")) {
+        printf("pow(%g,%g) = %g, should be %s\n", x, y, result, "-1.#INF00e+000");
+	k++;
+    }
+
+    for (i = 1, x = 0.0; i < 1000; i++) {
+        y = i;
+        answer = 0.0;
+        result = pow(x,y);
+
+        if (result != answer) {
+            printf("pow(%g,%g) = %g, should be %g\n", x, y, result, answer);
+		k++;
+        }
+    }
+
+    return (k);
+}
+
+
+typedef union   {
+    long lng[2];
+    double dbl;
+    } dbl;
+
+dbl d_inf = { 0x0, 0x7ff00000 };
+
+#define D_INF  (d_inf.dbl)
+
+typedef struct {
+	double x;
+	double y;
+	double answer;
+} _test;
+
+int test2(void)
+{
+	_test tests[] = {
+		{ 21.0,		-1.1e300,	0.0	},	// -D_INF???
+		{ 21.0,		1.1e-300,	1.0	},	// +D_INF???
+		{ -1.1e300,	21.0,		-D_INF	},
+		{ 21.0,		1.1e300,	D_INF	},
+		{ 1.0e100,	21.0,		D_INF	},
+		{ 21.0,		1.0e100,	D_INF	},
+		{ 1.0e100,	0.0,		1.0	},
+		{ 1.0e100,	1.0,		1.0e100	},
+		{ 1.0e100,	2.0,		1.0e200	},
+		{ 1.0e300,	1.0,		1.0e300	},
+		{ 1.0e300,	2.0,		D_INF	},
+	};
+	double result;
+	int i;
+	int k = 0;
+	char buf[BUFSIZ];
+
+	dbl foo = {0x78b58c40, 0x4415af1d};
+
+	for (i = 0; i < sizeof(tests) / sizeof(tests[0]); i++) {
+		result = pow(tests[i].x,tests[i].y);
+		if (result != tests[i].answer) {
+			// sprintf(buf, "%e", result);
+			// result = atof(buf);
+			// if (result != tests[i].answer) {
+				k++;
+				printf("pow(%e,%e) = %e, should be %e\n",
+					tests[i].x,
+					tests[i].y,
+					result,
+					tests[i].answer	
+					);
+			// }
+		}
+	}
+
+	return(k);
+}
+
+
+int test3()
+{
+	_test tests[] = {
+		{ 0.0,		-2.0,	D_INF	},
+		{ 0.0,		-1.0,	D_INF	},
+		{ 0.0,		0.0,	1.0	},
+		{ 0.0,		1.0,	0.0	},
+		{ 0.0,		2.0,	0.0	},
+		{ 0.0,		3.0,	0.0	},
+		{ 0.0,		4.0,	0.0	},
+	};
+	double result;
+	int i;
+	int k = 0;
+	char buf[BUFSIZ];
+
+	for (i = 0; i < sizeof(tests) / sizeof(tests[0]); i++) {
+		result = pow(tests[i].x,tests[i].y);
+		if (result != tests[i].answer) {
+			if (result != tests[i].answer) {
+				k++;
+				printf("pow(%e,%e) = %e, should be %e\n",
+					tests[i].x,
+					tests[i].y,
+					result,
+					tests[i].answer	
+					);
+			}
+		}
+	}
+
+	return(k);
+}
+