+
+ atom=moldyn->atom;
+
+ v.x=0; v.y=0; v.z=0;
+
+ for(j=0;j<INS_ATOMS;j++) {
+ run=1;
+ while(run) {
+ // tetrahedral
+ /*
+ r.x=0.0;
+ r.y=0.0;
+ r.z=0.0;
+ */
+ // hexagonal
+ /*
+ r.x=-1.0/8.0*ALBE_LC_SI;
+ r.y=-1.0/8.0*ALBE_LC_SI;
+ r.z=1.0/8.0*ALBE_LC_SI;
+ */
+ // 110 dumbbell
+ /*
+ r.x=(-0.5+0.25+0.125)*ALBE_LC_SI;
+ r.y=(-0.5+0.25+0.125)*ALBE_LC_SI;
+ r.z=(-0.5+0.25)*ALBE_LC_SI;
+ md->atom[4372].r.x=(-0.5+0.125+0.125)*ALBE_LC_SI;
+ md->atom[4372].r.y=(-0.5+0.125+0.125)*ALBE_LC_SI;
+ */
+ // random
+ //
+ r.x=(rand_get_double(&(moldyn->random))-0.5)*INS_LENX;
+ r.y=(rand_get_double(&(moldyn->random))-0.5)*INS_LENY;
+ r.z=(rand_get_double(&(moldyn->random))-0.5)*INS_LENZ;
+ //
+ // offset
+ r.x+=INS_OFFSET;
+ r.y+=INS_OFFSET;
+ r.z+=INS_OFFSET;
+ /* assume valid coordinates */
+ run=0;
+ for(i=0;i<moldyn->count;i++) {
+ atom=&(moldyn->atom[i]);
+ v3_sub(&dist,&(atom->r),&r);
+ check_per_bound(moldyn,&dist);
+ d=v3_absolute_square(&dist);
+ /* reject coordinates */
+ if(d<INS_R_C) {
+ //printf("atom %d - %f\n",i,d);
+ run=1;
+ break;
+ }
+ }
+ }
+ add_atom(moldyn,INS_TYPE,INS_MASS,INS_BRAND,
+ ATOM_ATTR_1BP|ATOM_ATTR_2BP|ATOM_ATTR_3BP|\
+ //ATOM_ATTR_HB|ATOM_ATTR_VB,
+ ATOM_ATTR_HB,
+ &r,&v);
+ }
+
+ return 0;
+}
+
+int sic_hook(void *moldyn,void *hook_params) {
+