diff options
Diffstat (limited to 'src/test.c')
-rw-r--r-- | src/test.c | 19 |
1 files changed, 10 insertions, 9 deletions
@@ -1257,21 +1257,22 @@ int test_norm_cdf(char *err) static double gsl_time_pdf(double x, void *params) { - double mu_noise, mu_indirect, mu_direct, mu_shower, ts, ts_shower, tmean, tsigma, ts_sigma; + double mu_noise, mu_indirect, tmean; + double mu[2], ts[2], ts_sigma[2]; double *pars = (double *) params; mu_noise = pars[0]; mu_indirect = pars[1]; - mu_direct = pars[2]; - mu_shower = pars[3]; - ts = pars[4]; - ts_shower = pars[5]; + mu[0] = pars[2]; + mu[1] = pars[3]; + ts[0] = pars[4]; + ts[1] = pars[5]; tmean = pars[6]; - tsigma = pars[7]; - ts_sigma = pars[8]; + ts_sigma[0] = pars[7]; + ts_sigma[1] = pars[8]; - return time_pdf(x,mu_noise,mu_indirect,&mu_direct,&mu_shower,1,&ts,&ts_shower,tmean,tsigma,&ts_sigma); + return time_pdf(x,mu_noise,mu_indirect,mu,2,ts,tmean,ts_sigma); } int test_time_pdf_norm(char *err) @@ -1378,7 +1379,7 @@ int test_time_cdf(char *err) goto err; } - expected = time_cdf(t,params[0],params[1],¶ms[2],¶ms[3],1,¶ms[4],¶ms[5],params[6],params[7],¶ms[8]); + expected = time_cdf(t,params[0],params[1],¶ms[2],2,¶ms[4],params[6],¶ms[7]); if (!isclose(result, expected, 1e-2, 0)) { sprintf(err, "integral of time_pdf = %.5f, but expected %.5f", result, expected); |