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

1 files changed, 8 insertions, 1 deletions

diff --git a/src/anal_process.c b/src/anal_process.c
index c66aea0..0ebec84 100644
--- a/src/anal_process.c
+++ b/src/anal_process.c
@@ -37,6 +37,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;
@@ -94,6 +96,11 @@ int main(int argc, char *argv[]) {
 			mom(dist_len, 3, dist, mom3);
 			vals_c3[i] = mom3[0];
 			errors_c3[i] = mom3[1];
+
+			double mom4[2];
+			mom(dist_len, 4, dist, mom4);
+			vals_c4[i] = mom4[0];
+			errors_c4[i] = mom4[1];
 		}
 		free(dist);
 	}
@@ -107,7 +114,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);