c - 程序收到信号 SIGSEGV，段错误，./sysdeps/ieee754/dbl-64/e_log.c : No such file or directory

Question

我正在尝试用 c 语言编写一个简单的超几何测试。此代码适用于小数字，但不适用于大数字。

例如:hyperG 35 50 90 3400 -- 有效

但是 hyperG ./hyperG 307107 486302 9073845 12147105 --失败

使用在线调试器我得到:

程序收到信号SIGSEGV，段错误。
__ieee754_log_avx (x=2898563) 位于 ../sysdeps/ieee754/dbl-64/e_log.c:76
76 ../sysdeps/ieee754/dbl-64/e_log.c:没有这样的文件或目录。

如有任何建议，我们将不胜感激。

代码如下:

   #include <stdio.h>
   #include <stdlib.h>
   #include <math.h>
  /*
   * program: hyperG
   * purpose: Calculate a probability based on a hypergeometric distribution
   * input  : k, n, M, N
   * output : Return a probability
   * 
   */
  // Variables
  double k = 0.0;
  double n = 0.0;
  double M = 0.0;
  double N = 0.0;

  double logchoos( double n, double k);
  double logfact( double n);

//Probability of at least k out of n tries of having something
//that occurs M of of N times
double hypergeometric( double k, double n, double M, double N)
{
    double p = 0.0;
    int i = 0;

    for(i = k; i <= n; i++){
        p += exp(logchoos((N-M),(n-i)) + logchoos(M,i) - logchoos(N,n));
    }

    return p;
}

// compute log "n choose k" or n = n!/(n-k)!k!
double logchoos( double n, double k) 
{
    double result = logfact(n) - logfact(n-k) - logfact(k);

    return result;
}

// Calculate the log factorial 
double logfact( double n )
{
    double fac;

    if( n < 1)
        return 0.0;
    else
       fac =  log(n) + logfact(n-1);

    return(fac);
}


int main(int argc, char *argv[])
{
    // Get command line arguments
    k = atof(argv[1]);
    n = atof(argv[2]);
    M = atof(argv[3]);
    N = atof(argv[4]);

    double prob = hypergeometric(k,n,M,N);

    printf("%e\n",prob);

    return 0;
}

我的MakeFile:

Makefile 
CC = gcc
CFLAGS = -Wall -g
LDFLAGS = -lm
FILES = hyperG.c

build:  $(FILES) 
    $(CC) $(FILES) -o hyperG $(CFLAGS) $(LDFLAGS) 
clean:
    rm -f hyperG

Answer 1

由于它适用于较小的值，但适用于较大的值则失败，它可能会耗尽堆栈，如果将递归更改为 while 循环，则不会出现 307107 486302 9073845 12147105 作为输入的段错误，尽管需要很长时间才能完成执行

double logfact( double n )
{
    double fac = 0;
    while(n > 1)
    {
        fac += log(n);
        n = n-1;
    }

    return(fac);
}