introduce a 2body bond function with callback, modified pair corr calc

[physik/posic.git] / bond_analyze.c

diff --git a/bond_analyze.c b/bond_analyze.c
new file mode 100644 (file)
index 0000000..0feffef
--- /dev/null
+++ b/bond_analyze.c
@@ -0,0 +1,54 @@
+/*
+ * bonding analyzation code
+ *
+ * author: frank.zirkelbach@physik.uni-augsburg.de
+ *
+ */
+
+#define _GNU_SOURCE
+#include <stdio.h>
+//#include <stdlib.h>
+//#include <unistd.h>
+//#include <string.h>
+//#include <sys/types.h>
+//#include <sys/stat.h>
+//#include <fcntl.h>
+
+#include "moldyn.h"
+
+int usage(char *prog) {
+
+       printf("\nusage:\n");
+       printf("  %s <save file>\n\n",prog);
+
+       return -1;
+}
+
+int main(int argc,char **argv) {
+
+       t_moldyn moldyn;
+       int ret;
+       double quality;
+
+       if(argc!=2) {
+               usage(argv[0]);
+               return -1;
+       }
+
+       memset(&moldyn,0,sizeof(t_moldyn));
+
+       printf("[bond analyze] reading save file ...\n");
+       ret=moldyn_read_save_file(&moldyn,argv[1]);
+       if(ret) {
+               printf("[bond analyze] exit!\n");
+               return ret;
+       }
+
+       bond_analyze(&moldyn,&quality);
+
+       printf("[bond analyze] quality = %f\n",quality);
+
+       moldyn_free_save_file(&moldyn);
+
+       return 0;
+}