X-Git-Url: https://www.hackdaworld.org/gitweb/?p=physik%2Fposic.git;a=blobdiff_plain;f=bond_analyze.c;h=39928941fbf3aae13c2588c35e0a1a902edb6d8c;hp=0feffef13399501587f61bb4f7d76f4769295755;hb=HEAD;hpb=2346d076f3d0955fb3a7bfb47f2b17b9eae6dd5b diff --git a/bond_analyze.c b/bond_analyze.c index 0feffef..3992894 100644 --- a/bond_analyze.c +++ b/bond_analyze.c @@ -15,6 +15,7 @@ //#include #include "moldyn.h" +#include "potentials/albe.h" int usage(char *prog) { @@ -28,7 +29,7 @@ int main(int argc,char **argv) { t_moldyn moldyn; int ret; - double quality; + double quality[2]; if(argc!=2) { usage(argv[0]); @@ -44,9 +45,18 @@ int main(int argc,char **argv) { return ret; } - bond_analyze(&moldyn,&quality); + /* potential setting */ + set_potential(&moldyn,MOLDYN_POTENTIAL_AM); + albe_mult_set_params(&moldyn,SI,C); - printf("[bond analyze] quality = %f\n",quality); + /* link cell init */ + link_cell_init(&moldyn,VERBOSE); + + /* analyzing ... */ + bond_analyze(&moldyn,quality); + + printf("[bond analyze] t: %f | quality = %f | %f | %d\n", + moldyn.time,quality[0],quality[1],moldyn.count); moldyn_free_save_file(&moldyn);