gamma replaced with cwt_gamma

This commit is contained in:
Wei-Hsuan Yang 2019-04-15 13:33:55 +08:00
parent 1e597aaffa
commit d76bd20573

View File

@ -111,7 +111,7 @@ void gauss(int N,int p,double lb,double ub,double *psi,double *t) {
t[i] = lb + delta * i;
}
den = sqrt(gamma(p+0.5));
den = sqrt(cwt_gamma(p+0.5));
if ((p+1)%2 == 0) {
num = 1.0;