X-Git-Url: https://hackdaworld.org/gitweb/?a=blobdiff_plain;f=pair_correlation_calc.c;h=b3386c737f99b8cc687143f4aa9ac4c7967a1367;hb=32377cc471f36b75710e620289e62284c28b4c76;hp=f6a256fc32eff3fdfae609c5c559a52ee3d50a8f;hpb=e1080fc0dd66b0cf5b7715c5e99e7a34ac04a8cf;p=physik%2Fposic.git

diff --git a/pair_correlation_calc.c b/pair_correlation_calc.c
index f6a256f..b3386c7 100644
--- a/pair_correlation_calc.c
+++ b/pair_correlation_calc.c
@@ -65,27 +65,14 @@ int main(int argc,char **argv) {
 
 	calculate_pair_correlation(&moldyn,dr,stat);
 
-	fd=open("pair_corr_func_ab.txt",
+	fd=open("pair_corr_func.txt",
 	        O_WRONLY|O_CREAT|O_TRUNC,S_IRUSR|S_IWUSR);
-	dprintf(fd,"# type a - type b bonds\n");
+	dprintf(fd,"#r #ab #aa #bb\n");
 	for(i=0;i<slots;i++)
-		dprintf(fd,"%f %f\n",i*dr,stat[i]);
+		dprintf(fd,"%f %f %f %f\n",
+		        i*dr,stat[i],stat[slots+i],stat[2*slots+i]);
 	close(fd);
-		
-	fd=open("pair_corr_func_aa.txt",
-	        O_WRONLY|O_CREAT|O_TRUNC,S_IRUSR|S_IWUSR);
-	dprintf(fd,"# type a - type a bonds\n");
-	for(i=0;i<slots;i++)
-		dprintf(fd,"%f %f\n",i*dr,stat[slots+i]);
-	close(fd);
-		
-	fd=open("pair_corr_func_bb.txt",
-	        O_WRONLY|O_CREAT|O_TRUNC,S_IRUSR|S_IWUSR);
-	dprintf(fd,"# type a - type b bonds\n");
-	for(i=0;i<slots;i++)
-		dprintf(fd,"%f %f\n",i*dr,stat[2*slots+i]);
-	close(fd);
-		
+
 	free(stat);
 
 	return 0;