s2xyz.c

   1 /*
   2  * calcultae pair correlation function
   3  *
   4  * author: frank.zirkelbach@physik.uni-augsburg.de
   5  *
   6  */
   7
   8 #define _GNU_SOURCE
   9 #include <stdio.h>
  10 #include <stdlib.h>
  11 #include <unistd.h>
  12 #include <string.h>
  13 #include <sys/types.h>
  14 #include <sys/stat.h>
  15 #include <fcntl.h>
  16
  17 #include "moldyn.h"
  18
  19 #define ME "[s2xyz]"
  20
  21 int main(int argc,char **argv) {
  22
  23         t_moldyn moldyn;
  24         int i,ret;
  25         char color[2][64];
  26         char type[2][8];
  27
  28         memset(&moldyn,0,sizeof(t_moldyn));
  29
  30         strcpy(color[0],"Yellow");
  31         strcpy(type[0],"Si");
  32         strcpy(color[1],"Gray");
  33         strcpy(type[1],"C");
  34
  35         ret=moldyn_read_save_file(&moldyn,argv[1]);
  36         if(ret) {
  37                 printf("%s exit!\n",ME);
  38                 return ret;
  39         }
  40
  41
  42         printf("# [P] %d 0 <0,0,0>\n",moldyn.count);
  43
  44         for(i=0;i<moldyn.count;i++)
  45                 printf("%s %f %f %f %s %f\n",
  46                        type[moldyn.atom[i].brand],
  47                        moldyn.atom[i].r.x,
  48                        moldyn.atom[i].r.y,
  49                        moldyn.atom[i].r.z,
  50                        color[moldyn.atom[i].brand],
  51                        moldyn.atom[i].ekin/EV);
  52
  53         return 0;
  54 }
  55