计划接收信号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
答
由于其较小的值工作,但失败了较大的值,它可能用完堆栈的,如果你改变递归while循环,虽然需要很多时间才能完成执行,但307107 486302 9073845 12147105没有出现段错误
double logfact(double n)
{
double fac = 0;
while(n > 1)
{
fac += log(n);
n = n-1;
}
return(fac);
}
我貌似找不到你的数学共享库'libm'。 – jiveturkey
由于其较小的值工作,但失败了较大的值,它可能不多了递归logfact功能 – Pras
堆栈如何让编译器找到库?的lib位于:/usr/lib/x86_64-linux-gnu/libm.a – mplace