X-Git-Url: https://hackdaworld.org/gitweb/?p=physik%2Fnlsop.git;a=blobdiff_plain;f=nlsop.c;h=40803d5fe5e75b2d4ff23c13e69a475cf2ea4bab;hp=d51d24cfdc274199c5e18ba856877bb7f5310e89;hb=9a0a1b10e8e13f1cc1117eaa4c25a02112a36117;hpb=c80606ee95b42185168c1ae7ce9c6b24ea8085bc diff --git a/nlsop.c b/nlsop.c index d51d24c..40803d5 100644 --- a/nlsop.c +++ b/nlsop.c @@ -61,9 +61,9 @@ int usage(void) printf("-s \t steps (default %d)\n",STEPS); printf("-d \t refresh display (default %d)\n",REFRESH); printf("-r \t amorphous influence range (default %d)\n",RANGE); - printf("-f \t pressure = * 1/distance^2 (default %f)\n",A_AP); - printf("-p \t pressure offset (default %f)\n",B_AP); - printf("-F \t proportionality constant between c conc and ability to get amorphous (default %f)\n",A_CP); + printf("-f \t stress induced amorphization influence (default %f)\n",S_D); + printf("-p \t ballistic amorphization influence (default %f)\n",B_D); + printf("-F \t carbon induced amorphization influence (default %f)\n",C_D); printf("-D \t diffusion rate from cryst to amorph cells (default %f)\n",DR_AC); printf("-c \t diffusion rate in cryst cells (default %f)\n",DR_CC); printf("-e \t do diffusion every steps (default %d)\n",DIFF_RATE); @@ -91,7 +91,7 @@ int process_cell(d3_lattice *d3_l,u32 x,u32 y,u32 z,info *my_info,u32 nel_z) thiz=d3_l->status+x+y*d3_l->max_x+z*d3_l->max_x*d3_l->max_y; conc=d3_l->extra+x+y*d3_l->max_x+z*d3_l->max_x*d3_l->max_y; - p=my_info->b_ap*nel_z; + p=my_info->b*nel_z; for(i=-(my_info->range);i<=my_info->range;i++) { for(j=-(my_info->range);j<=my_info->range;j++) @@ -99,11 +99,11 @@ int process_cell(d3_lattice *d3_l,u32 x,u32 y,u32 z,info *my_info,u32 nel_z) if(!(i==0 && j==0)) { off=((x+d3_l->max_x+i)%d3_l->max_x)+((y+d3_l->max_y+j)%d3_l->max_x)*d3_l->max_x+z*d3_l->max_x*d3_l->max_y; - if(*(d3_l->status+off)&AMORPH) p+=my_info->a_ap*(*(d3_l->extra+off))*URAND_MAX/(i*i+j*j); + if(*(d3_l->status+off)&AMORPH) p+=my_info->s*(*(d3_l->extra+off))*URAND_MAX/(i*i+j*j); } } } - p+=*conc*my_info->a_cp*URAND_MAX; + p+=*conc*my_info->c*URAND_MAX; if(!(*thiz&AMORPH)) { if(get_rand(URAND_MAX)<=p) MAKE_AMORPH(thiz); @@ -341,11 +341,14 @@ int write_ac_distr(d3_lattice *d3_l,int ac_distr) { int fd,x,y,z; int count,offset; - char file[16]; - - if(ac_distr==1) strcpy(file,"a.plot"); - if(ac_distr==2) strcpy(file,"c.plot"); - if(ac_distr==3) strcpy(file,"b.plot"); + char file[32]; + int si_count; + + si_count=d3_l->max_x*d3_l->max_y*SI_PER_VOLUME; + if(ac_distr==1) strcpy(file,"carbon_in_av.plot"); + if(ac_distr==2) strcpy(file,"carbon_in_cv.plot"); + if(ac_distr==3) strcpy(file,"carbon.plot"); + if(ac_distr==4) strcpy(file,"amorphous_volumes.plot"); if((fd=open(file,O_WRONLY|O_CREAT))<0) { @@ -366,9 +369,12 @@ int write_ac_distr(d3_lattice *d3_l,int ac_distr) if(ac_distr==2) if(!(*(d3_l->status+offset)&AMORPH)) count+=*(d3_l->extra+offset); if(ac_distr==3) count+=*(d3_l->extra+offset); + if(ac_distr==4) + if(*(d3_l->status+offset)&AMORPH) count+=1; } } - dprintf(fd,"%d %d\n",z,count); + if(ac_distr==4) dprintf(fd,"%d %d\n",z*CELL_LENGTH,count); + else dprintf(fd,"%d %f\n",z*CELL_LENGTH,100.0*count/si_count); } close(fd); @@ -430,7 +436,54 @@ int write_bmp(d3_lattice *d3_l,int window,u32 x,u32 y,u32 z,int max) width=d3_l->max_x; height=d3_l->max_y; } - + if(window==7) + { + sprintf(bmpfile,"x-z_a_cdistr_%d.bmp",y); + foo=3*d3_l->max_x; + size=(foo+(4-foo%4))*d3_l->max_z; + height=d3_l->max_z; + width=d3_l->max_x; + } + if(window==8) + { + sprintf(bmpfile,"y-z_a_cdistr_%d.bmp",x); + foo=3*d3_l->max_y; + size=(foo+(4-foo%4))*d3_l->max_z; + height=d3_l->max_z; + width=d3_l->max_y; + } + if(window==9) + { + sprintf(bmpfile,"x-y_a_cdistr_%d.bmp",z); + foo=3*d3_l->max_x; + size=(foo+(4-foo%4))*d3_l->max_y; + height=d3_l->max_y; + width=d3_l->max_x; + } + if(window==10) + { + sprintf(bmpfile,"x-z_c_cdistr_%d.bmp",y); + foo=3*d3_l->max_x; + size=(foo+(4-foo%4))*d3_l->max_z; + height=d3_l->max_z; + width=d3_l->max_x; + } + if(window==11) + { + sprintf(bmpfile,"y-z_c_cdistr_%d.bmp",x); + foo=3*d3_l->max_y; + size=(foo+(4-foo%4))*d3_l->max_z; + height=d3_l->max_z; + width=d3_l->max_y; + } + if(window==12) + { + sprintf(bmpfile,"x-y_c_cdistr_%d.bmp",z); + foo=3*d3_l->max_x; + size=(foo+(4-foo%4))*d3_l->max_y; + height=d3_l->max_y; + width=d3_l->max_x; + } if((fd=open(bmpfile,O_WRONLY|O_CREAT))<0) { puts("cannot open bmp file"); @@ -609,7 +662,7 @@ int write_bmp(d3_lattice *d3_l,int window,u32 x,u32 y,u32 z,int max) } if(window==6) { - for(j=0;jmax_z;j++) + for(j=0;jmax_y;j++) { for(i=0;imax_x;i++) { @@ -633,6 +686,168 @@ int write_bmp(d3_lattice *d3_l,int window,u32 x,u32 y,u32 z,int max) } } } + if(window==7) + { + for(j=0;jmax_z;j++) + { + for(i=0;imax_x;i++) + { + if(*(d3_l->status+i+y*d3_l->max_x+(d3_l->max_z-j-1)*d3_l->max_x*d3_l->max_y)&RED) sum=*(d3_l->extra+i+y*d3_l->max_x+(d3_l->max_z-j-1)*d3_l->max_x*d3_l->max_y); + else sum=0; + sum=sum*255/max; + memset(buf,(unsigned char)sum,3); + if(write(fd,buf,3)<3) + { + puts("failed writing rgb values to bmp file"); + return-1; + } + } + if(foo%4) + { + memset(buf,0,4-foo%4); + if(write(fd,buf,4-foo%4)<4-foo%4) + { + puts("failed writing 4 byte ending"); + return -1; + } + } + } + } + if(window==8) + { + for(j=0;jmax_z;j++) + { + for(i=0;imax_x;i++) + { + if(*(d3_l->status+x+i*d3_l->max_x+(d3_l->max_z-j-1)*d3_l->max_x*d3_l->max_y)&RED) sum=*(d3_l->extra+x+i*d3_l->max_x+(d3_l->max_z-j-1)*d3_l->max_x*d3_l->max_y); + else sum=0; + sum=sum*255/max; + memset(buf,(unsigned char)sum,3); + if(write(fd,buf,3)<3) + { + puts("failed writing rgb values to bmp file"); + return-1; + } + } + if(foo%4) + { + memset(buf,0,4-foo%4); + if(write(fd,buf,4-foo%4)<4-foo%4) + { + puts("failed writing 4 byte ending"); + return -1; + } + } + } + } + if(window==9) + { + for(j=0;jmax_y;j++) + { + for(i=0;imax_x;i++) + { + if(*(d3_l->status+i+(d3_l->max_y-j-1)*d3_l->max_x+z*d3_l->max_x*d3_l->max_y)&RED) sum=*(d3_l->extra+i+(d3_l->max_y-j-1)*d3_l->max_x+z*d3_l->max_x*d3_l->max_y); + else sum=0; + sum=sum*255/max; + memset(buf,(unsigned char)sum,3); + if(write(fd,buf,3)<3) + { + puts("failed writing rgb values to bmp file"); + return-1; + } + } + if(foo%4) + { + memset(buf,0,4-foo%4); + if(write(fd,buf,4-foo%4)<4-foo%4) + { + puts("failed writing 4 byte ending"); + return -1; + } + } + } + } + if(window==10) + { + for(j=0;jmax_z;j++) + { + for(i=0;imax_x;i++) + { + if(!(*(d3_l->status+i+y*d3_l->max_x+(d3_l->max_z-j-1)*d3_l->max_x*d3_l->max_y)&RED)) sum=*(d3_l->extra+i+y*d3_l->max_x+(d3_l->max_z-j-1)*d3_l->max_x*d3_l->max_y); + else sum=0; + sum=sum*255/max; + memset(buf,(unsigned char)sum,3); + if(write(fd,buf,3)<3) + { + puts("failed writing rgb values to bmp file"); + return-1; + } + } + if(foo%4) + { + memset(buf,0,4-foo%4); + if(write(fd,buf,4-foo%4)<4-foo%4) + { + puts("failed writing 4 byte ending"); + return -1; + } + } + } + } + if(window==11) + { + for(j=0;jmax_z;j++) + { + for(i=0;imax_x;i++) + { + if(!(*(d3_l->status+x+i*d3_l->max_x+(d3_l->max_z-j-1)*d3_l->max_x*d3_l->max_y)&RED)) sum=*(d3_l->extra+x+i*d3_l->max_x+(d3_l->max_z-j-1)*d3_l->max_x*d3_l->max_y); + else sum=0; + sum=sum*255/max; + memset(buf,(unsigned char)sum,3); + if(write(fd,buf,3)<3) + { + puts("failed writing rgb values to bmp file"); + return-1; + } + } + if(foo%4) + { + memset(buf,0,4-foo%4); + if(write(fd,buf,4-foo%4)<4-foo%4) + { + puts("failed writing 4 byte ending"); + return -1; + } + } + } + } + if(window==12) + { + for(j=0;jmax_y;j++) + { + for(i=0;imax_x;i++) + { + if(!(*(d3_l->status+i+(d3_l->max_y-j-1)*d3_l->max_x+z*d3_l->max_x*d3_l->max_y)&RED)) sum=*(d3_l->extra+i+(d3_l->max_y-j-1)*d3_l->max_x+z*d3_l->max_x*d3_l->max_y); + else sum=0; + sum=sum*255/max; + memset(buf,(unsigned char)sum,3); + if(write(fd,buf,3)<3) + { + puts("failed writing rgb values to bmp file"); + return-1; + } + } + if(foo%4) + { + memset(buf,0,4-foo%4); + if(write(fd,buf,4-foo%4)<4-foo%4) + { + puts("failed writing 4 byte ending"); + return -1; + } + } + } + } close(fd); return 1; @@ -866,7 +1081,7 @@ u32 get_reject_graph(info *my_info,d3_lattice *d3_l,char *file,u32 *graph) { free(flag); #ifdef DEBUG_INTERPOL_PROFILE - printf("debug: (interpolated profile)\n"); + printf("debug: %s (interpolated profile)\n",file); for(i=0;imax_z;i++) printf("%d %d\n",i,graph[i]); #endif @@ -930,9 +1145,9 @@ int main(int argc,char **argv) resave=RESAVE; my_info.z_diff=0; my_info.c_diff=1; - my_info.a_ap=A_AP; - my_info.b_ap=B_AP; - my_info.a_cp=A_CP; + my_info.s=S_D; + my_info.b=B_D; + my_info.c=C_D; my_info.cc=CC; my_info.dr_ac=DR_AC; my_info.dr_cc=DR_CC; @@ -993,13 +1208,13 @@ int main(int argc,char **argv) my_info.range=atoi(argv[++i]); break; case 'f': - my_info.a_ap=atof(argv[++i]); + my_info.s=atof(argv[++i]); break; case 'p': - my_info.b_ap=atof(argv[++i]); + my_info.b=atof(argv[++i]); break; case 'F': - my_info.a_cp=atof(argv[++i]); + my_info.c=atof(argv[++i]); break; case 'W': resave=atoi(argv[++i]); @@ -1108,9 +1323,9 @@ int main(int argc,char **argv) sprintf(s_txt,"steps: %d",my_info.steps); sprintf(dose_txt,"dose: %.2fe+17 C/cm²",my_info.steps*1.0/(d3_l.max_x*d3_l.max_y*CELL_LENGTH*CELL_LENGTH*1000)); sprintf(r_txt,"pressure range: %d",my_info.range); - sprintf(stress_txt,"stress term: %f",my_info.a_ap); - sprintf(ballistic_txt,"ballistic term: %f",my_info.b_ap); - sprintf(carbon_txt,"carbon term: %f",my_info.a_cp); + sprintf(stress_txt,"stress term: %f",my_info.s); + sprintf(ballistic_txt,"ballistic term: %f",my_info.b); + sprintf(carbon_txt,"carbon term: %f",my_info.c); sprintf(dr_ac_txt,"a/c diffusion rate: %f",my_info.dr_ac); if(my_info.c_diff!=0) sprintf(dr_cc_txt,"c/c diffusion rate: %f",my_info.dr_cc); else sprintf(dr_cc_txt,"c/c diffusion rate: none"); @@ -1146,6 +1361,20 @@ int main(int argc,char **argv) arg_v[25]=NULL; #endif + /* compute graphs for random number rejection method */ + if((c_profile=(u32 *)malloc(d3_l.max_z*sizeof(unsigned int)))==NULL) + { + puts("failed allocating memory for carbon profile graph"); + return -1; + } + if((n_e_loss=(u32 *)malloc(d3_l.max_z*sizeof(unsigned int)))==NULL) + { + puts("failed allocating memory for nuclear energy loss graph"); + return -1; + } + ip_max=get_reject_graph(&my_info,&d3_l,p_file,c_profile); + ne_max=get_reject_graph(&my_info,&d3_l,n_e_file,n_e_loss); + if((!strcmp(l_file,""))||(c_step)) { /* calculate ratio of c_simwindow / c_total */ @@ -1154,21 +1383,9 @@ int main(int argc,char **argv) puts("failed calculating ratio"); return -1; } - /* compute graphs for random number rejection method */ - if((c_profile=(u32 *)malloc(d3_l.max_z*sizeof(unsigned int)))==NULL) - { - puts("failed allocating memory for carbon profile graph"); - return -1; - } - if((n_e_loss=(u32 *)malloc(d3_l.max_z*sizeof(unsigned int)))==NULL) - { - puts("failed allocating memory for nuclear energy loss graph"); - return -1; - } - ip_max=get_reject_graph(&my_info,&d3_l,p_file,c_profile); - ne_max=get_reject_graph(&my_info,&d3_l,n_e_file,n_e_loss); #ifdef DEBUG_RAND + i=0; while(1) { #ifdef DEBUG_CP @@ -1180,6 +1397,8 @@ int main(int argc,char **argv) #ifdef DEBUG_NORM printf("%d\n",get_rand(d3_l.max_z)); #endif + if(i==10000000) return 1; + i++; } #endif @@ -1204,7 +1423,7 @@ int main(int argc,char **argv) sprintf(conc_txt,"conc: %d",*(d3_l.extra+x+y*d3_l.max_x+z*d3_l.max_x*d3_l.max_y)); sprintf(steps_txt,"step: %d",i); sprintf(cc_txt,"total c: %d",my_info.cc); - d3_lattice_draw(&d3_l,x,y,z,25,arg_v,mode,0); + d3_lattice_draw(&d3_l,x,y,z,25,arg_v,mode,0,NULL,0,NULL,0); } #endif if(i%resave==0 && strcmp(s_file,"") && resave!=0 && i!=0) @@ -1242,6 +1461,7 @@ int main(int argc,char **argv) if(switchmode==0) mode=0; if(switchmode==1) mode=1; if(switchmode==2) mode=2; + if(switchmode==3) mode=3; /* end of bahh! */ sprintf(xyz_txt,"x: %d y: %d z: %d",x+1,y+1,z+1); sprintf(status_txt,"status: %c",(*(d3_l.status+x+y*d3_l.max_x+z*d3_l.max_x*d3_l.max_y)&AMORPH)?'a':'c'); @@ -1251,7 +1471,8 @@ int main(int argc,char **argv) if(switchmode==0) strcpy(mode_txt,"view: a/c mode"); if(switchmode==1) strcpy(mode_txt,"view: c conc mode"); if(switchmode==2) strcpy(mode_txt,"view: a pressure mode"); - d3_lattice_draw(&d3_l,x,y,z,25,arg_v,mode,max_extra); + if(switchmode==3) strcpy(mode_txt,"view: a/c + profiles mode"); + d3_lattice_draw(&d3_l,x,y,z,25,arg_v,mode,max_extra,c_profile,ip_max,n_e_loss,ne_max); bmp=0; ac_distr=0; scan_event(&d3_l,&x,&y,&z,&quit,&escape,&switchmode,&bmp,&ac_distr);