【发布时间】：2018-12-06 06:34:51
【问题描述】：
所以我正在尝试实现我自己的 pow 函数。为了让事情变得更简单，我可以使用 int 类型作为指数。所以代码正在工作......主要是。我有一个名为 xmath.c 的文件，我在其中实现了我的 pow 函数和另一个 test.c 来测试限制。它不工作的唯一情况是pow(-111, INT_MAX)。我怎样才能改进我的代码？还告诉我是否也应该发布 test.c 代码。 xfabs 是另一个我必须自己实现的函数，而不是使用 fabs。
忘记在我的帖子中包含 xfabs：
#define xfabs(x) (((x) > 0) ? (x) : (-(x)))
#include <stdio.h>
#include <errno.h>
#include <math.h>
#include <float.h>
#include "xmath.h"
double xpow(double x, int y)
{
errno = 0;
double product = 1;
double prod = 1;
int i;
if (y == 0) {
return 1.0;
} else if (x == 0 && y < 0) {
errno = EDOM;
return 0.0;
}
if (xfabs(x) >= 1) {
if (y > 0) {
for (i = 0; i < y; i++) {
if (xfabs(product) > DBL_MAX / xfabs(x)) {
errno = ERANGE;
return HUGE_VAL;
}
product *= x;
}
return product;
} else if (y < 0) {
for (i = 0; i > y; i--) {
if (xfabs(product) < DBL_MIN / xfabs(x)) {
errno = ERANGE;
return -0.0;
}
product *= 1 / x;
}
return product;
}
} else if (xfabs(x) < 1) {
if (y > 0) {
for (i = 0; i < y; i++) {
if (xfabs(prod) < DBL_MIN / xfabs(x)) {
errno = ERANGE;
return +0.0;
}
prod *= x;
}
return prod;
} else if (y < 0) {
for (i = 0; i > y; i--) {
if (xfabs(prod) < DBL_MAX / xfabs(x)) {
errno = ERANGE;
return -HUGE_VAL;
}
prod *= 1 / x;
}
return prod;
}
}
}
这里是测试功能。我还在为我在 test.c 中使用的 xmath.c 创建一个目标文件：
#include <stdio.h>
#include <float.h>
#include <math.h>
#include <errno.h>
#include <limits.h>
#include "xmath.h"
void xpow_tests(void)
{
printf("XPOW TESTS (SUCCESS == 1 / FAILURE == 0)\n");
printf("=================================================\n");
printf("Test (xpow(2, 5)): \t\t%i\t", xpow(2, 5) == pow(2, 5));
xpow(2, 5);
printf("errno: %i\n", errno);
printf("Test (xpow(-2, -4)): \t\t%i\t", xpow(-2, -4) == pow(-2, -4));
xpow(-2, -4);
printf("errno: %i\n", errno);
printf("Test (xpow(1, 15)): \t\t%i\t", xpow(1, 15) == pow(1, 15));
xpow(1, 15);
printf("errno: %i\n", errno);
printf("Test (xpow(-4, 414)): \t\t%i\t", xpow(-4, 414) == pow(-4, 414));
xpow(-4, 414);
printf("errno: %i\n", errno);
printf("Test (xpow(-5, 303)): \t\t%i\t", xpow(-5, 303) == pow(-5, 303));
xpow(-5, 303);
printf("errno: %i\n", errno);
printf("Test (xpow(0, -3)): \t\t%i\t", xpow(0, -3) == 0);
xpow(0, -3);
printf("errno: %i\n", errno);
printf("Test (xpow(1.0e-10, 100)): \t%i\t", xpow(1.0e-10, 100) == pow(1.0e-10, 100));
xpow(1.0e-10, 100);
printf("errno: %i\n", errno);
printf("Test (xpow(-1.0e-10, 101)): \t%i\t", xpow(-1.0e-10, 101) == pow(-1.0e-10, 101));
xpow(-1.0e-10, 101);
printf("errno: %i\n", errno);
printf("Test (xpow(-111, INT_MAX)): \t%i\t", xpow(-111, INT_MAX) == pow(-111, INT_MAX));
xpow(-111, INT_MAX);
printf("errno: %i\n", errno);
printf("Test (xpow(-111, INT_MAX-1)): \t%i\t", xpow(-111, INT_MAX-1) == pow(-111, INT_MAX-1));
xpow(-111, INT_MAX-1);
printf("errno: %i\n", errno);
printf("Test (xpow(-111, -INT_MAX)): \t%i\t", xpow(-111, -INT_MAX) == pow(-111, -INT_MAX));
xpow(-111, -INT_MAX);
printf("errno: %i\n", errno);
printf("Test (xpow(-111, -INT_MAX+1)): \t%i\t", xpow(-111, -INT_MAX+1) == pow(-111, -INT_MAX+1));
xpow(-111, -INT_MAX+1);
printf("errno: %i\n", errno);
printf("=================================================\n");
printf("errno (EDOM: %i / ERANGE: %i)\n\n", EDOM, ERANGE);
}
int main(void)
{
xpow_tests();
return 0;
}
【问题讨论】：
-
“不工作”是什么意思？请在这里给我们一些具体的东西。
-
所以函数根本不计算 pow(-111, INT_MAX) 这意味着我的函数 xpow 没有为这些值提供任何结果，或者它们不等于实际 pow 函数的值。 .这有帮助吗？
-
无论如何它都会给出
double结果。您是否有此代码的测试用例，您只是随意敲打它吗？测试驱动开发，您可以针对所有边界条件进行明确定义的测试，这是这里的关键。 -
是的，我将编辑测试功能...我知道该功能有点复杂...我不知道如何改进它，因为我还在学习...
-
你熟悉溢出的概念吗？ double 可以准确表示的最大尺寸是多少？我建议从较低的指数开始并增加它直到你的结果是垃圾，然后在调试器中逐步确定原因。