nlsop.c

   1 /*
   2  * nlsop.c
   3  *
   4  * this program tries helping to understand the amorphous depuration
   5  * and recrystallization of SiCx while ion implantation at temperatures
   6  * below 400 degree celsius.
   7  * hopefully the program will simulate the stabilization of the
   8  * selforganizing lamella structure in the observed behaviour.
   9  *
  10  * refs:
  11  *  - J. K. N. Lindner. Habilationsschrift, Universitaet Augsburg.
  12  *  - Maik Haeberlen. Diplomarbeit, Universitaet Augsburg.
  13  */
  14
  15 #define _GNU_SOURCE
  16 #include <stdio.h>
  17 #include <stdlib.h>
  18 #include <string.h>
  19 #include <sys/types.h>
  20 #include <sys/stat.h>
  21 #include <fcntl.h>
  22 #include <unistd.h>
  23
  24 #include "nlsop.h"
  25
  26 #include "dfbapi.h"
  27 #include "random.h"
  28
  29 #define MAKE_AMORPH(N) *(N)|=AMORPH
  30 #define MAKE_CRYST(N) *(N)&=~AMORPH
  31
  32 int usage(void)
  33 {
  34  puts("usage:");
  35  puts("-h \t\t help");
  36  puts("-n \t\t no user interaction");
  37  puts("-Z \t\t cryst -> amorph c diffusion in z direction");
  38  printf("-a <value> \t slope of nuclear energy loss (default %f)\n",A_EL);
  39  printf("-b <value> \t nuclear energy loss offset (default %f)\n",B_EL);
  40  printf("-x <value> \t # x cells (default %d)\n",X);
  41  printf("-y <value> \t # y cells (default %d)\n",Y);
  42  printf("-z <value> \t # z cells (default %d)\n",Z);
  43  printf("-s <value> \t steps (default %d)\n",STEPS);
  44  printf("-d <value> \t refresh display (default %d)\n",REFRESH);
  45  printf("-r <value> \t amorphous influence range (default %d)\n",RANGE);
  46  printf("-f <value> \t pressure = <value> * 1/distance^2 (default %f)\n",A_AP);
  47  printf("-p <value> \t pressure offset (default %f)\n",B_AP);
  48  printf("-F <value> \t proportionality constant between c conc and ability to get amorphous (default %f)\n",A_CP);
  49  printf("-A <value> \t slope of linear c distribution (default %f)\n",A_CD);
  50  printf("-B <value> \t linear c distribution offset (default %f)\n",B_CD);
  51  printf("-D <value> \t diffusion rate from cryst to amorph cells (default %f)\n",DR_AC);
  52  printf("-c <value> \t diffusion rate in cryst cells (default %f)\n",DR_CC);
  53  printf("-e <value> \t do diffusion every <value> steps (default %d)\n",DIFF_RATE);
  54  printf("-W <value> \t write every <value> steps to save file (default %d)\n",RESAVE);
  55  puts("-C <file> \t convert file to gnuplot format");
  56  puts("-L <file> \t load from file");
  57  puts("-S <file> \t save to file");
  58  puts("-R <file> \t read from random file");
  59  puts("-P <file> \t specify implantatin profile file");
  60
  61  return 1;
  62 }
  63
  64 int process_cell(d3_lattice *d3_l,u32 x,u32 y,u32 z,info *my_info)
  65 {
  66  unsigned char *thiz;
  67  int *conc;
  68  int i,j;
  69  int off;
  70  double p;
  71
  72  thiz=d3_l->status+x+y*d3_l->max_x+z*d3_l->max_x*d3_l->max_y;
  73  conc=d3_l->extra+x+y*d3_l->max_x+z*d3_l->max_x*d3_l->max_y;
  74  p=my_info->b_ap*URAND_MAX;
  75  for(i=-(my_info->range);i<=my_info->range;i++)
  76  {
  77   for(j=-(my_info->range);j<=my_info->range;j++)
  78   {
  79    if(!(i==0 && j==0))
  80    {
  81     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;
  82     if(*(d3_l->status+off)&AMORPH) p+=my_info->a_ap*(*(d3_l->extra+off))*URAND_MAX/(i*i+j*j);
  83    }
  84   }
  85  }
  86  p+=*conc*my_info->a_cp*URAND_MAX;
  87  printf("debug: %f - %u\n",p,URAND_MAX);
  88  if(!(*thiz&AMORPH))
  89  {
  90   if(get_rand(URAND_MAX)<=p) MAKE_AMORPH(thiz);
  91  } else
  92  {
  93   /* assume 1-p probability */
  94   if(get_rand(URAND_MAX)>p) MAKE_CRYST(thiz);
  95  }
  96
  97  my_info->cc+=1;
  98
  99  return 1;
 100 }
 101
 102 int distrib_c(d3_lattice *d3_l,info *my_info,int step)
 103 {
 104  u32 x,y,z;
 105  int i,j,k,c;
 106  int offset,off;
 107  int carry;
 108
 109  /* put one c ion somewhere in the lattice */
 110  x=get_rand(d3_l->max_x);
 111  y=get_rand(d3_l->max_y);
 112  z=get_rand_lgp(d3_l->max_z,my_info->a_cd,my_info->b_cd);
 113  *(d3_l->extra+x+y*d3_l->max_x+z*d3_l->max_x*d3_l->max_y)+=1;
 114
 115  if(step%my_info->diff_rate==0)
 116  {
 117
 118  for(i=0;i<d3_l->max_x;i++)
 119  {
 120   for(j=0;j<d3_l->max_y;j++)
 121   {
 122    for(k=0;k<d3_l->max_z;k++)
 123    {
 124     offset=i+j*d3_l->max_x+k*d3_l->max_x*d3_l->max_y;
 125     /* case amorph - amorph <-> cryst diffusion */
 126     if(*(d3_l->status+offset)&AMORPH)
 127     {
 128      for(c=-1;c<=1;c++)
 129      {
 130       if(c!=0)
 131       {
 132        off=((i+d3_l->max_x+c)%d3_l->max_x)+j*d3_l->max_x+k*d3_l->max_x*d3_l->max_y;
 133        carry=0;
 134        if(!(*(d3_l->status+off)&AMORPH)) carry=(int)(my_info->dr_ac*(*(d3_l->extra+off)));
 135        if(carry!=0)
 136        {
 137         *(d3_l->extra+offset)+=carry;
 138         *(d3_l->extra+off)-=carry;
 139        }
 140       }
 141      }
 142      for(c=-1;c<=1;c++)
 143      {
 144       if(c!=0)
 145       {
 146        off=i+((j+c+d3_l->max_y)%d3_l->max_y)*d3_l->max_x+k*d3_l->max_x*d3_l->max_y;
 147        carry=0;
 148        if(!(*(d3_l->status+off)&AMORPH)) carry=(int)(my_info->dr_ac*(*(d3_l->extra+off)));
 149        if(carry!=0)
 150        {
 151         *(d3_l->extra+offset)+=carry;
 152         *(d3_l->extra+off)-=carry;
 153        }
 154       }
 155      }
 156      if(my_info->z_diff)
 157      {
 158       for(c=-1;c<=1;c++)
 159       {
 160        if(c!=0)
 161        {
 162         off=i+j*d3_l->max_x+((k+c+d3_l->max_z)%d3_l->max_z)*d3_l->max_x*d3_l->max_y;
 163         carry=0;
 164         if(!*(d3_l->status+off)&AMORPH) carry=(int)(my_info->dr_ac*(*(d3_l->extra+off)));
 165         if(carry!=0)
 166         {
 167          *(d3_l->extra+off)-=carry;
 168          *(d3_l->extra+offset)+=carry;
 169         }
 170        }
 171       }
 172      }
 173     } else
 174     /* case not amorph - cryst <-> cryst diffusion */
 175     {
 176      for(c=-1;c<=1;c++)
 177      {
 178       if(c!=0)
 179       {
 180        off=i+((j+c+d3_l->max_y)%d3_l->max_y)*d3_l->max_x+k*d3_l->max_x*d3_l->max_y;
 181        carry=0;
 182        if(!(*(d3_l->status+off)&AMORPH))
 183        {
 184         carry=(*(d3_l->extra+off)-*(d3_l->extra+offset))/2;
 185         if(carry!=0)
 186         {
 187          *(d3_l->extra+offset)+=carry;
 188          *(d3_l->extra+off)-=carry;
 189         }
 190        }
 191       }
 192      }
 193      for(c=-1;c<=1;c++)
 194      {
 195       if(c!=0)
 196       {
 197        off=((i+c+d3_l->max_x)%d3_l->max_x)+j*d3_l->max_x+k*d3_l->max_x*d3_l->max_y;
 198        carry=0;
 199        if(!(*(d3_l->status+off)&AMORPH))
 200        {
 201         carry=(*(d3_l->extra+off)-*(d3_l->extra+offset))/2;
 202         if(carry!=0)
 203         {
 204          *(d3_l->extra+offset)+=carry;
 205          *(d3_l->extra+off)-=carry;
 206         }
 207        }
 208       }
 209      }
 210     }
 211    } /* for z */
 212   } /* for y */
 213  } /* for x */
 214
 215  } /* if step modulo diff_rate == 0 */
 216
 217  return 1;
 218 }
 219
 220 int calc_pressure(d3_lattice *d3_l,int range)
 221 {
 222  int i,j,off;
 223  double count;
 224  int x,y,z;
 225
 226  for(x=0;x<d3_l->max_x;x++)
 227  {
 228   for(y=0;y<d3_l->max_y;y++)
 229   {
 230    for(z=0;z<d3_l->max_z;z++)
 231    {
 232     count=0;
 233     for(i=-range;i<=range;i++)
 234     {
 235      for(j=-range;j<=range;j++)
 236      {
 237       if(i!=0 && j!=0)
 238       {
 239        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;
 240        if(*(d3_l->status+off)&AMORPH) count+=((double)*(d3_l->extra+off))/(i*i+j*j);
 241       }
 242      }
 243     }
 244     *(unsigned char *)(d3_l->v_ptr+x+y*d3_l->max_x+z*d3_l->max_x*d3_l->max_y)=(unsigned char)(count*255/MAX_VPTR);
 245    }
 246   }
 247  }
 248
 249  return 1;
 250 }
 251
 252 int save_to_file(char *sf,d3_lattice *d3_l,info *my_inf)
 253 {
 254  int sf_fd,c;
 255
 256  if((sf_fd=open(sf,O_WRONLY|O_CREAT))<0)
 257  {
 258   puts("cannot open save file");
 259   return -1;
 260  }
 261  if(write(sf_fd,d3_l,sizeof(d3_lattice))<sizeof(d3_lattice))
 262  {
 263   puts("failed saving d3 lattice struct");
 264   return -1;
 265  }
 266  if(write(sf_fd,my_inf,sizeof(info))<sizeof(info))
 267  {
 268   puts("failed saving info struct");
 269   return-1;
 270  }
 271  c=d3_l->max_x*d3_l->max_y*d3_l->max_z;
 272  if(write(sf_fd,d3_l->status,c*sizeof(unsigned char))<c*sizeof(unsigned char))
 273  {
 274   puts("failed saving status of d3 lattice sites");
 275   return -1;
 276  }
 277  if(write(sf_fd,d3_l->extra,c*sizeof(int))<c*sizeof(int))
 278  {
 279   puts("failed saving sites concentration");
 280   return -1;
 281  }
 282  close(sf_fd);
 283
 284  return 1;
 285 }
 286
 287 int load_from_file(char *lf,d3_lattice *d3_l,info *my_inf)
 288 {
 289  int lf_fd,c;
 290
 291  if((lf_fd=open(lf,O_RDONLY))<0)
 292  {
 293   puts("cannot open load file");
 294   return -1;
 295  }
 296  if(read(lf_fd,d3_l,sizeof(d3_lattice))<sizeof(d3_lattice))
 297  {
 298   puts("failed reading d3 lattice struct");
 299   return -1;
 300  }
 301  if(read(lf_fd,my_inf,sizeof(info))<sizeof(info))
 302  {
 303   puts("failed reading info struct");
 304   return-1;
 305  }
 306  c=d3_l->max_x*d3_l->max_y*d3_l->max_z;
 307  if((d3_l->status=(unsigned char*)malloc(c*sizeof(unsigned char)))==NULL)
 308  {
 309   puts("cannot allocate status buffer");
 310   return -1;
 311  }
 312  if((d3_l->extra=(int *)malloc(c*sizeof(int)))==NULL)
 313  {
 314   puts("cannot allocate concentration buffer");
 315   return -1;
 316  }
 317  if(read(lf_fd,d3_l->status,c*sizeof(unsigned char))<c*sizeof(unsigned char))
 318  {
 319   puts("failed reading status of d3 lattice sites");
 320   return -1;
 321  }
 322  if(read(lf_fd,d3_l->extra,c*sizeof(int))<c*sizeof(int))
 323  {
 324   puts("failed reading sites concentration");
 325   return -1;
 326  }
 327  close(lf_fd);
 328
 329  return 1;
 330 }
 331
 332 int convert_file(char *cf,d3_lattice *d3_l)
 333 {
 334  int x,y,z;
 335  int c_fd;
 336
 337  if((c_fd=open(cf,O_WRONLY|O_CREAT))<0)
 338  {
 339   puts("cannot open convert file");
 340   return -1;
 341  }
 342  dprintf(c_fd,"# created by nlsop (gnuplot format)\n");
 343  for(x=0;x<d3_l->max_x;x++)
 344  {
 345   for(y=0;y<d3_l->max_y;y++)
 346   {
 347    for(z=0;z<d3_l->max_z;z++)
 348    {
 349     if(*(d3_l->status+x+y*d3_l->max_x+z*d3_l->max_x*d3_l->max_y)&AMORPH) dprintf(c_fd,"%d %d %d\n",x,y,z);
 350    }
 351   }
 352  }
 353  close(c_fd);
 354
 355  return 1;
 356 }
 357
 358 int main(int argc,char **argv)
 359 {
 360  u32 x,y,z,x_c,y_c,z_c;
 361  int i,quit,escape,switchmode,nowait;
 362  int refresh,resave;
 363  char s_file[MAX_CHARS];
 364  char s_file_tmp[MAX_CHARS];
 365  char l_file[MAX_CHARS];
 366  char c_file[MAX_CHARS];
 367  char p_file[MAX_CHARS];
 368  char convert;
 369  char r_file[MAX_CHARS];
 370 #ifdef USE_DFB_API
 371  char xyz_txt[MAX_TXT];
 372  char status_txt[MAX_TXT];
 373  char conc_txt[MAX_TXT];
 374  char steps_txt[MAX_TXT];
 375  char cc_txt[MAX_TXT];
 376  char a_txt[MAX_TXT];
 377  char s_txt[MAX_TXT];
 378  char ap_txt[MAX_TXT];
 379  char el_txt[MAX_TXT];
 380  char cd_txt[MAX_TXT];
 381  char r_txt[MAX_TXT];
 382  char cp_txt[MAX_TXT];
 383  char zdiff_txt[MAX_TXT];
 384  char diff_txt[MAX_TXT];
 385  char dr_ac_txt[MAX_TXT];
 386  char dr_cc_txt[MAX_TXT];
 387  char mode_txt[MAX_TXT];
 388  char *arg_v[MAX_ARGV];
 389 #endif
 390  d3_lattice d3_l;
 391  info my_info;
 392  unsigned char mode;
 393
 394  d3_l.max_x=X;
 395  d3_l.max_y=Y;
 396  d3_l.max_z=Z;
 397  my_info.steps=STEPS;
 398  my_info.range=RANGE;
 399  refresh=REFRESH;
 400  resave=RESAVE;
 401  my_info.z_diff=0;
 402  my_info.a_el=A_EL;
 403  my_info.b_el=B_EL;
 404  my_info.a_cd=A_CD;
 405  my_info.b_cd=B_CD;
 406  my_info.a_ap=A_AP;
 407  my_info.b_ap=B_AP;
 408  my_info.a_cp=A_CP;
 409  my_info.cc=CC;
 410  my_info.dr_ac=DR_AC;
 411  my_info.dr_cc=DR_CC;
 412  my_info.diff_rate=DIFF_RATE;
 413  nowait=0;
 414  quit=0;
 415  escape=0;
 416  switchmode=0;
 417  strcpy(s_file,"");
 418  strcpy(l_file,"");
 419  strcpy(c_file,"");
 420  strcpy(p_file,IMP_PROFILE);
 421  convert=0;
 422  strcpy(r_file,"");
 423  mode=0;
 424
 425  for(i=1;i<argc;i++)
 426  {
 427   if(argv[i][0]=='-')
 428   {
 429    switch(argv[i][1])
 430    {
 431     case 'h':
 432      usage();
 433      return -1;
 434     case 'n':
 435      nowait=1;
 436      break;
 437     case 'a':
 438      my_info.a_el=atof(argv[++i]);
 439      break;
 440     case 'b':
 441      my_info.b_el=atof(argv[++i]);
 442      break;
 443     case 'x':
 444      d3_l.max_x=atoi(argv[++i]);
 445      break;
 446     case 'y':
 447      d3_l.max_y=atoi(argv[++i]);
 448      break;
 449     case 'z':
 450      d3_l.max_z=atoi(argv[++i]);
 451      break;
 452     case 'Z':
 453      my_info.z_diff=1;
 454      break;
 455     case 's':
 456      my_info.steps=atoi(argv[++i]);
 457      break;
 458     case 'd':
 459      refresh=atoi(argv[++i]);
 460      break;
 461     case 'r':
 462      my_info.range=atoi(argv[++i]);
 463      break;
 464     case 'f':
 465      my_info.a_ap=atof(argv[++i]);
 466      break;
 467     case 'p':
 468      my_info.b_ap=atof(argv[++i]);
 469      break;
 470     case 'F':
 471      my_info.a_cp=atof(argv[++i]);
 472      break;
 473     case 'A':
 474      my_info.a_cd=atof(argv[++i]);
 475      break;
 476     case 'B':
 477      my_info.b_cd=atof(argv[++i]);
 478      break;
 479     case 'W':
 480      resave=atoi(argv[++i]);
 481      break;
 482     case 'C':
 483      strcpy(l_file,argv[++i]);
 484      if(i<argc-1) if(argv[i+1][0]!='-') strcpy(c_file,argv[++i]);
 485      convert=1;
 486      break;
 487     case 'D':
 488      my_info.dr_ac=atof(argv[++i]);
 489      break;
 490     case 'c':
 491      my_info.dr_cc=atof(argv[++i]);
 492      break;
 493     case 'e':
 494      my_info.diff_rate=atoi(argv[++i]);
 495      break;
 496     case 'L':
 497      strcpy(l_file,argv[++i]);
 498      break;
 499     case 'S':
 500      strcpy(s_file,argv[++i]);
 501      break;
 502     case 'R':
 503      strcpy(r_file,argv[++i]);
 504      break;
 505     case 'P':
 506      strcpy(p_file,argv[++i]);
 507      break;
 508     default:
 509      usage();
 510      return -1;
 511    }
 512   } else usage();
 513  }
 514
 515  x=d3_l.max_x/2-1;
 516  y=d3_l.max_y/2-1;
 517  z=d3_l.max_z/2-1;
 518
 519 #ifdef NODFB
 520  if(!strcmp(s_file,""))
 521  {
 522   puts("NODFB defined, run with -S option");
 523   return -1;
 524  }
 525 #endif
 526
 527  if(!strcmp(r_file,"")) rand_init(NULL);
 528  else rand_init(r_file);
 529
 530  if(!strcmp(l_file,""))
 531  {
 532   i=d3_l.max_x*d3_l.max_y*d3_l.max_z;
 533 #ifdef USE_DFB_API
 534   d3_lattice_init(&argc,argv,&d3_l);
 535 #endif
 536   if((d3_l.status=(unsigned char *)malloc(i*sizeof(unsigned char)))==NULL)
 537   {
 538    puts("failed allocating status buffer");
 539    return -1;
 540   }
 541   memset(d3_l.status,0,i*sizeof(unsigned char));
 542   if((d3_l.extra=(int *)malloc(i*sizeof(int)))==NULL)
 543   {
 544    puts("failed allocating concentration buffer");
 545    return -1;
 546   }
 547   memset(d3_l.extra,0,i*sizeof(int));
 548  } else
 549  {
 550   load_from_file(l_file,&d3_l,&my_info);
 551   if(convert)
 552   {
 553    if(!strcmp(c_file,"")) sprintf(c_file,"%s_gnuplot",l_file);
 554    printf("converting file %s to %s\n",l_file,c_file);
 555    convert_file(c_file,&d3_l);
 556    puts("done");
 557    return 1;
 558   }
 559 #ifdef USE_DFB_API
 560     else d3_lattice_init(&argc,argv,&d3_l);
 561 #endif
 562  }
 563
 564 #ifdef USE_DFB_API
 565  d3_event_init(&d3_l);
 566 #endif
 567
 568 #ifdef USE_DFB_API
 569  strcpy(a_txt,"args:");
 570  sprintf(s_txt,"steps: %d",my_info.steps);
 571  sprintf(r_txt,"pressure range: %d",my_info.range);
 572  sprintf(ap_txt,"a_ap: %.2f  b_ap: %.3f",my_info.a_ap,my_info.b_ap);
 573  sprintf(el_txt,"a_el: %.2f  b_el: %.3f",my_info.a_el,my_info.b_el);
 574  sprintf(cd_txt,"a_cd: %.2f  b_cd: %.3f",my_info.a_cd,my_info.b_cd);
 575  sprintf(cp_txt,"a_cp: %.4f",my_info.a_cp);
 576  sprintf(dr_ac_txt,"a/c diffusion rate: %.4f",my_info.dr_ac);
 577  sprintf(dr_cc_txt,"c/c diffusion rate: %.4f",my_info.dr_cc);
 578  sprintf(zdiff_txt,"diffusion in z direction: %c",my_info.z_diff?'y':'n');
 579  sprintf(diff_txt,"diffusion every %d steps",my_info.diff_rate);
 580  strcpy(mode_txt,"view: a/c mode");
 581  arg_v[1]=xyz_txt;
 582  arg_v[2]=NULL;
 583  arg_v[3]=status_txt;
 584  arg_v[4]=conc_txt;
 585  arg_v[5]=NULL;
 586  arg_v[6]=mode_txt;
 587  arg_v[7]=NULL;
 588  arg_v[8]=steps_txt;
 589  arg_v[9]=cc_txt;
 590  arg_v[10]=NULL;
 591  arg_v[11]=diff_txt;
 592  arg_v[12]=zdiff_txt;
 593  arg_v[13]=NULL;
 594  arg_v[14]=a_txt;
 595  arg_v[15]=s_txt;
 596  arg_v[16]=r_txt;
 597  arg_v[17]=ap_txt;
 598  arg_v[18]=el_txt;
 599  arg_v[19]=cd_txt;
 600  arg_v[20]=cp_txt;
 601  arg_v[21]=dr_ac_txt;
 602  arg_v[22]=dr_cc_txt;
 603  arg_v[23]=NULL;
 604  arg_v[24]=NULL;
 605 #endif
 606
 607  if(!strcmp(l_file,""))
 608  {
 609   i=0;
 610   while((i<my_info.steps) && (quit==0) && (escape==0))
 611   {
 612    x_c=get_rand(d3_l.max_x);
 613    y_c=get_rand(d3_l.max_y);
 614    z_c=get_rand_lgp(d3_l.max_z,my_info.a_el,my_info.b_el);
 615    distrib_c(&d3_l,&my_info);
 616    process_cell(&d3_l,x_c,y_c,z_c,&my_info);
 617 #ifdef USE_DFB_API
 618    if(i%refresh==0)
 619    {
 620     sprintf(xyz_txt,"x: %d  y: %d  z: %d",x+1,y+1,z+1);
 621     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');
 622     sprintf(conc_txt,"conc: %d",*(d3_l.extra+x+y*d3_l.max_x+z*d3_l.max_x*d3_l.max_y));
 623     sprintf(steps_txt,"step: %d",i);
 624     sprintf(cc_txt,"total c: %d",my_info.cc);
 625     d3_lattice_draw(&d3_l,x,y,z,24,arg_v,mode);
 626    }
 627 #endif
 628    if(i%resave==0 && strcmp(s_file,"") && resave!=0 && i!=0)
 629    {
 630     sprintf(s_file_tmp,"%s_%d_of_%d",s_file,i,my_info.steps);
 631     save_to_file(s_file_tmp,&d3_l,&my_info);
 632 #ifdef NODFB
 633     printf("saved %s\n",s_file_tmp);
 634 #endif
 635    }
 636    i++;
 637   }
 638  }
 639
 640  if(strcmp(s_file,""))
 641  {
 642    printf("saved %s\n",s_file);
 643    save_to_file(s_file,&d3_l,&my_info);
 644  }
 645
 646 #ifdef USE_DFB_API
 647  /* allocating buffer for pressure values */
 648  if((d3_l.v_ptr=malloc(d3_l.max_x*d3_l.max_y*d3_l.max_z*sizeof(unsigned char)))==NULL)
 649  {
 650   puts("cannot allocate buffer for pressure values");
 651   return -1;
 652  }
 653  /* calc values */
 654  calc_pressure(&d3_l,my_info.range);
 655
 656  while((quit==0) && (escape==0) && (nowait==0))
 657  {
 658   /* bahh! */
 659   if(switchmode==0) mode=0;
 660   if(switchmode==1) mode=1;
 661   if(switchmode==2) mode=2;
 662   /* end of bahh! */
 663   sprintf(xyz_txt,"x: %d  y: %d  z: %d",x+1,y+1,z+1);
 664   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');
 665   sprintf(conc_txt,"conc: %d",*(d3_l.extra+x+y*d3_l.max_x+z*d3_l.max_x*d3_l.max_y));
 666   strcpy(steps_txt,"step: end!");
 667   sprintf(cc_txt,"total c: %d",my_info.cc);
 668   if(switchmode==0) strcpy(mode_txt,"view: a/c mode");
 669   if(switchmode==1) strcpy(mode_txt,"view: c conc mode");
 670   if(switchmode==2) strcpy(mode_txt,"view: a pressure mode");
 671   d3_lattice_draw(&d3_l,x,y,z,24,arg_v,mode);
 672   scan_event(&d3_l,&x,&y,&z,&quit,&escape,&switchmode);
 673  }
 674
 675  d3_lattice_release(&d3_l);
 676 #endif
 677
 678  return 1;
 679 }