计划接收信号SIGSEGV，分割过错，./sysdeps/ieee754/dbl-64/e_log.c：没有这样的文件或目录

我想用C来写一个简单的超几何测试。此代码适用于小的数字，但无法为大量。计划接收信号SIGSEGV，分割过错，./sysdeps/ieee754/dbl-64/e_log.c：没有这样的文件或目录

例如：hyperG 35 50 90 3400 - 工程

但hyperG ./hyperG 307107 486302 9073845 12147105 --fails

使用我得到一个在线调试器：

计划接收信号SIGSEGV ，分段故障。在../sysdeps/ieee754/dbl-64/e_log.c:76
76 ../sysdeps/ieee754/dbl-64/e_log.c
__ieee754_log_avx（X = 2898563）：没有这样的文件或目录。

任何建议，将不胜感激。

代码如下：

#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 
CC = gcc 
CFLAGS = -Wall -g 
LDFLAGS = -lm 
FILES = hyperG.c 

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

来源

2017-05-04 mplace

我貌似找不到你的数学共享库'libm'。 – jiveturkey

由于其较小的值工作，但失败了较大的值，它可能不多了递归logfact功能 – Pras

堆栈如何让编译器找到库？的lib位于：/usr/lib/x86_64-linux-gnu/libm.a – mplace

由于其较小的值工作，但失败了较大的值，它可能用完堆栈的，如果你改变递归while循环，虽然需要很多时间才能完成执行，但307107 486302 9073845 12147105没有出现段错误

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

    return(fac); 
}

来源

2017-05-04 16:11:13 Pras

计划接收信号SIGSEGV，分割过错，./sysdeps/ieee754/dbl-64/e_log.c：没有这样的文件或目录

回答

相关问题