various changes, including adding central moments and changing the fuse thresholds to long doubles

1 files changed, 10 insertions, 1 deletions

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;
 }