]> www.hackdaworld.org Git - physik/posic.git/blob - s2xyz.c
oups, forgot the basis trafo stuff
[physik/posic.git] / 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