X-Git-Url: https://hackdaworld.org/gitweb/?a=blobdiff_plain;f=mdrun.c;h=c30ce5ffffca7738a6e851faa56defd3d27ee6da;hb=232fadda0a2f21a081b0a420bd46f767a8627575;hp=34807e76bd76f646f5d587705d09ce0e05840614;hpb=c8a82a10a403a5e8efc9e83ee36f5d6d35ab639e;p=physik%2Fposic.git diff --git a/mdrun.c b/mdrun.c index 34807e7..c30ce5f 100644 --- a/mdrun.c +++ b/mdrun.c @@ -98,6 +98,9 @@ int add_stage(t_mdrun *mdrun,u8 type,void *params) { case STAGE_DISPLACE_ATOM: psize=sizeof(t_displace_atom_params); break; + case STAGE_DEL_ATOMS: + psize=sizeof(t_del_atoms_params); + break; case STAGE_INSERT_ATOMS: psize=sizeof(t_insert_atoms_params); break; @@ -175,6 +178,7 @@ int mdrun_parse_config(t_mdrun *mdrun) { t_set_temp_params stp; t_set_timestep_params stsp; t_fill_params fp; + t_del_atoms_params delp; /* open config file */ fd=open(mdrun->cfile,O_RDONLY); @@ -207,6 +211,7 @@ int mdrun_parse_config(t_mdrun *mdrun) { memset(&stp,0,sizeof(t_set_temp_params)); memset(&stsp,0,sizeof(t_set_timestep_params)); memset(&fp,0,sizeof(t_fill_params)); + memset(&delp,0,sizeof(t_del_atoms_params)); // get command + args wcnt=0; @@ -303,10 +308,12 @@ int mdrun_parse_config(t_mdrun *mdrun) { fp.lz=atoi(word[++i]); fp.lc=atof(word[++i]); mdrun->lc=fp.lc; + continue; } if(!strncmp(word[i],"eb",2)) { fp.fill_element=atoi(word[++i]); fp.fill_brand=atoi(word[++i]); + continue; } if(word[i][0]=='p') { i+=1; @@ -342,6 +349,7 @@ int mdrun_parse_config(t_mdrun *mdrun) { fp.p_params.d.y=atof(word[++i]); fp.p_params.d.z=atof(word[++i]); } + continue; } if(word[i][0]=='d') { switch(word[++i][0]) { @@ -382,6 +390,7 @@ int mdrun_parse_config(t_mdrun *mdrun) { default: break; } + continue; } i+=1; @@ -509,6 +518,13 @@ int mdrun_parse_config(t_mdrun *mdrun) { dap.dz=atof(word[5]); add_stage(mdrun,STAGE_DISPLACE_ATOM,&dap); } + else if(!strncmp(word[1],"del_atoms",9)) { + delp.o.x=atof(word[2]); + delp.o.y=atof(word[3]); + delp.o.z=atof(word[4]); + delp.r=atof(word[5]); + add_stage(mdrun,STAGE_DEL_ATOMS,&delp); + } else if(!strncmp(word[1],"ins_atoms",9)) { iap.ins_steps=atoi(word[2]); iap.ins_atoms=atoi(word[3]); @@ -732,6 +748,31 @@ int displace_atom(t_moldyn *moldyn,t_mdrun *mdrun) { return 0; } +int del_atoms(t_moldyn *moldyn,t_mdrun *mdrun) { + + t_stage *stage; + t_del_atoms_params *delp; + int i; + t_3dvec dist; + + stage=mdrun->stage.current->data; + delp=stage->params; + + for(i=0;icount;i++) { + v3_sub(&dist,&(delp->o),&(moldyn->atom[i].r)); +//printf("%d ----> %f %f %f = %f | %f\n",i,dist.x,dist.y,dist.z,v3_absolute_square(&dist),delp->r*delp->r); + if(v3_absolute_square(&dist)<=(delp->r*delp->r)) { + del_atom(moldyn,moldyn->atom[i].tag); + printf("%s atom deleted: %d %d %d\n",ME, + moldyn->atom[i].tag,moldyn->atom[i].element, + moldyn->atom[i].brand); + } + } + + return 0; + +} + int insert_atoms(t_moldyn *moldyn,t_mdrun *mdrun) { t_insert_atoms_params *iap; @@ -1108,6 +1149,11 @@ int mdrun_hook(void *ptr1,void *ptr2) { displace_atom(moldyn,mdrun); change_stage=TRUE; break; + case STAGE_DEL_ATOMS: + stage_print(" -> del atoms\n\n"); + del_atoms(moldyn,mdrun); + change_stage=TRUE; + break; case STAGE_INSERT_ATOMS: stage_print(" -> insert atoms\n\n"); iap=stage->params; @@ -1213,8 +1259,9 @@ int mdrun_hook(void *ptr1,void *ptr2) { create_lattice(moldyn, fp->lattice,fp->lc, - mdrun->element1, - DEFAULT_ATOM_ATTR,0, + fp->fill_element, + DEFAULT_ATOM_ATTR, + fp->fill_brand, fp->lx,fp->ly,fp->lz, NULL, &(fp->p_params),