Changeset 314
- Timestamp:
- 03/25/08 18:24:18 (8 months ago)
- Files:
-
- 1 modified
-
trunk/matml/src/ternary/ternary.c (modified) (4 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/matml/src/ternary/ternary.c
r313 r314 14 14 char gv_version[100]; 15 15 FILE *pfd = NULL; 16 int i, j, loop_max=50 , vertcolor=1;16 int i, j, loop_max=50; 17 17 coordT Gmin=1.e+100, Gmax=-1.e+100, centerG; 18 18 … … 50 50 fprintf (pfd, "(geometry \"Ternary Phase Diagram\" { : tpd })"); 51 51 fprintf (pfd, "(read geometry { define tpd \n"); 52 fprintf (pfd, "appearance {}\n %sOFF\n", vertcolor ? "C" : "");52 fprintf (pfd, "appearance {}\nCOFF\n"); 53 53 fprintf (pfd, "%d %d 0\n", (loop_max+1)*(loop_max+2)/2, loop_max*loop_max); 54 54 … … 80 80 z = G; 81 81 82 fprintf (pfd, "%g %g %g", x, z, 1.-y); 83 if (vertcolor) 84 fprintf (pfd, " %g %g %g 1", 85 RED(Grel(G)), GREEN(Grel(G)), BLUE(Grel(G))); 86 fprintf (pfd, "\n"); 82 fprintf (pfd, "%g %g %g %g %g %g 1\n", x, z, 1.-y, 83 RED(Grel(G)), GREEN(Grel(G)), BLUE(Grel(G))); 87 84 } 88 85 … … 93 90 for (j=0; j<loop_max-1-i; j++) 94 91 { 95 centerG = free_energy (((coordT)i+0.333333)/loop_max, 96 ((coordT)j+0.333333)/loop_max, 0.); 97 fprintf (pfd, "3 %d %d %d", 92 fprintf (pfd, "3 %d %d %d\n", 98 93 ROWSTART(i)+j, ROWSTART(i)+j+1, ROWSTART(i+1)+j); 99 if (!vertcolor) 100 fprintf (pfd, " %g %g %g 1", 101 RED(Grel(centerG)), GREEN(Grel(centerG)), 102 BLUE(Grel(centerG))); 103 fprintf (pfd, "\n"); 104 105 centerG = free_energy (((coordT)i+0.666667)/loop_max, 106 ((coordT)j+0.666667)/loop_max, 0.); 107 fprintf (pfd, "3 %d %d %d", 94 fprintf (pfd, "3 %d %d %d\n", 108 95 ROWSTART(i)+j+1, ROWSTART(i+1)+j, ROWSTART(i+1)+j+1); 109 if (!vertcolor)110 fprintf (pfd, " %g %g %g 1",111 RED(Grel(centerG)), GREEN(Grel(centerG)),112 BLUE(Grel(centerG)));113 fprintf (pfd, "\n");114 96 } 115 centerG = free_energy (((coordT)i+0.333333)/loop_max, 116 ((coordT)j+0.333333)/loop_max, 0.); 117 fprintf (pfd, "3 %d %d %d", 97 fprintf (pfd, "3 %d %d %d\n", 118 98 ROWSTART(i)+j, ROWSTART(i)+j+1, ROWSTART(i+1)+j); 119 if (!vertcolor)120 fprintf (pfd, " %g %g %g 1",121 RED(Grel(centerG)), GREEN(Grel(centerG)),122 BLUE(Grel(centerG)));123 fprintf (pfd, "\n");124 99 } 125 centerG = free_energy (((coordT)i+0.333333)/loop_max, 126 0.333333/loop_max, 0.); 127 fprintf (pfd, "3 %d %d %d", 100 fprintf (pfd, "3 %d %d %d\n", 128 101 ROWSTART(i), ROWSTART(i)+1, ROWSTART(i+1)); 129 if (!vertcolor)130 fprintf (pfd, " %g %g %g 1",131 RED(Grel(centerG)), GREEN(Grel(centerG)),132 BLUE(Grel(centerG)));133 fprintf (pfd, "\n");134 102 fprintf (pfd, "})\n"); 135 103 fflush (pfd);