4 main(int argc, char *argv[])
9 float mu1=0. , mu2=6., density, afac, s;
11 FILE *fp = (FILE *)NULL;
13 fprintf(stderr, "Usage: %s iview nview phmfilename\n", argv[0]);
18 nview = atoi(argv[2]);
19 sprintf(filename, "%s", argv[3]);
21 s = (float)view/((float)(nview-1));
25 else if ( s > 9./16. )
28 afac = ( (s - 7./16.) / 2./16.);
29 density = log(1/((1-afac)*exp(-mu1) + afac * exp(-mu2)));
33 density = mu1 + (mu2-mu1)*s;
39 if ( (fp = fopen(filename, "w")) == (FILE *)NULL){
40 fprintf(stderr,"Error, can not open file \"tmpphmfile\"\n");
43 fprintf(fp, "rectangle 0 0 11.5 11.5 0 0\n");
44 fprintf(fp, "ellipse 0 0 11.4 11.4 0 1\n");
45 fprintf(fp, "ellipse 0 0 1.25 1.25 0 %f\n", density);