From 0ad947f800bcbe2c488d2d5cbcdb16c46e6d3857 Mon Sep 17 00:00:00 2001 From: Jaron Date: Tue, 8 Nov 2016 08:17:26 -0500 Subject: various changes, including adding central moments and changing the fuse thresholds to long doubles --- src/big_anal_process.c | 11 ++++++++++- 1 file changed, 10 insertions(+), 1 deletion(-) (limited to 'src/big_anal_process.c') diff --git a/src/big_anal_process.c b/src/big_anal_process.c index 8345785..0f7106f 100644 --- a/src/big_anal_process.c +++ b/src/big_anal_process.c @@ -53,6 +53,8 @@ int main(int argc, char *argv[]) { double *errors_c2 = (double *)malloc(nc * sizeof(double)); double *vals_c3 = (double *)malloc(nc * sizeof(double)); double *errors_c3 = (double *)malloc(nc * sizeof(double)); + double *vals_c4 = (double *)malloc(nc * sizeof(double)); + double *errors_c4 = (double *)malloc(nc * sizeof(double)); char *out_filename = (char *)malloc(100 * sizeof(char)); uint_t out_filename_len = 0; @@ -114,6 +116,11 @@ int main(int argc, char *argv[]) { get_mom(dist_len, 3, dist, mom1, mom3); vals_c3[i] = mom3[0]; errors_c3[i] = mom3[1]; + + double mom4[2]; + get_mom(dist_len, 4, dist, mom1, mom4); + vals_c4[i] = mom4[0]; + errors_c4[i] = mom4[1]; } free(dist); } @@ -127,7 +134,7 @@ int main(int argc, char *argv[]) { FILE *out_file = fopen(out_filename, "w"); for (uint_t i = 0; i < nc; i++) { - fprintf(out_file, "%u %g %g %g %g %g %g %g\n", Ls_c[i], betas_c[i], vals_c1[i], errors_c1[i], vals_c2[i], errors_c2[i], vals_c3[i], errors_c3[i]); + fprintf(out_file, "%u %g %g %g %g %g %g %g %g %g\n", Ls_c[i], betas_c[i], vals_c1[i], errors_c1[i], vals_c2[i], errors_c2[i], vals_c3[i], errors_c3[i], vals_c4[i], errors_c4[i]); } fclose(out_file); @@ -141,6 +148,8 @@ int main(int argc, char *argv[]) { free(errors_c2); free(vals_c3); free(errors_c3); + free(vals_c4); + free(errors_c4); return 0; } -- cgit v1.2.3-70-g09d2