c-此分段故障背后的原因

Question

这是我计算一个数字的反正弦的函数。它会导致值在0.51-0.8之间的分段故障：

double my_asin(double x)
{
    double sum = x;
    if(x < -1.0 || x > 1.0)
    {
        /* error handling */
    }
    else if(x < -0.5)
    {
        sum = -0.5*PI + my_asin(my_sqrt(1-my_pow(x,2))); // SIG_SEGV
        return sum;
    }
    else if(x > 0.5)
    {
        sum = 0.5*PI - my_asin(my_sqrt(1-my_pow(x,2))); // SIG_SEGV
        return sum;
    }
    /* variable initialization */
    while(my_abs(b - a) > EPSILON2)
    {
        /*code*/
    }
    /* return result */
}

GDB和valgrind都告诉我，函数my_pow中发生了错误，它的作用正是你认为的，所以没有必要在这里发布它。你能看一眼，给我指正确的方向吗？非常感谢。

Answer 1

假设x是sqrt(2)/2（约0.707）。您的函数使用参数sqrt(1-x*x)递归调用自己，该参数等于x。由于无限递归，这将导致堆栈溢出。

对于该值附近的x的其他值也会发生同样的情况。