X-Git-Url: https://www.hackdaworld.org/gitweb/?a=blobdiff_plain;f=sic.c;h=82b371018cfc11fa99bef055030b541b079b03d2;hb=a33fe7f5d4d7e09b5b52f5f4dadd462cf868ec2f;hp=2cd32eb04853a4c14518561999f9d6a429d9e5a1;hpb=b9696b3202fea9d936da68c322517ecbae994597;p=physik%2Fposic.git diff --git a/sic.c b/sic.c index 2cd32eb..82b3710 100644 --- a/sic.c +++ b/sic.c @@ -86,6 +86,7 @@ printf("%d\n",sizeof(t_atom)); /* cutoff radius */ printf("[sic] setting cutoff radius\n"); set_cutoff(&md,TM_S_SI); + //set_cutoff(&md,1.0*LC_SI); /* set (initial) dimensions of simulation volume */ printf("[sic] setting dimensions\n"); @@ -99,6 +100,7 @@ printf("%d\n",sizeof(t_atom)); printf("[sic] creating atoms\n"); create_lattice(&md,DIAMOND,LC_SI,SI,M_SI, ATOM_ATTR_1BP|ATOM_ATTR_2BP|ATOM_ATTR_3BP, + //ATOM_ATTR_2BP, 0,4,4,4); /* setting a nearest neighbour distance for the moldyn checks */ @@ -106,7 +108,7 @@ printf("%d\n",sizeof(t_atom)); /* set temperature */ printf("[sic] setting temperature\n"); - set_temperature(&md,273.0+450.0); + set_temperature(&md,10.0); /* initial thermal fluctuations of particles */ printf("[sic] thermal init\n"); @@ -114,12 +116,12 @@ printf("%d\n",sizeof(t_atom)); /* create the simulation schedule */ printf("[sic] adding schedule\n"); - moldyn_add_schedule(&md,10000,1.0e-15); + moldyn_add_schedule(&md,1000,1.0e-15); /* activate logging */ printf("[sic] activate logging\n"); - moldyn_set_log(&md,LOG_TOTAL_ENERGY,"saves/test-energy",100); - moldyn_set_log(&md,VISUAL_STEP,"saves/test-visual",100); + moldyn_set_log(&md,LOG_TOTAL_ENERGY,"saves/test-energy",10); + moldyn_set_log(&md,VISUAL_STEP,"saves/test-visual",10); /* * let's do the actual md algorithm now