Changeset 463 for trunk/matml/src
- Timestamp:
- 06/18/2009 02:43:36 PM (3 years ago)
- Location:
- trunk/matml/src/ternary
- Files:
-
- 6 modified
Legend:
- Unmodified
- Added
- Removed
-
trunk/matml/src/ternary/gibbs.h
r462 r463 48 48 etc. +*/ 49 49 SPINODAL /*+ Spinodal boundary +*/ 50 } phase_boundary_ flags;50 } phase_boundary_type; 51 51 52 52 typedef struct { … … 58 58 "current" phase +*/ 59 59 int n_phases; /*+ Number of phases involved +*/ 60 phase_boundary_ flags flags; /*+ See phase_boundary_flags typedef enum +*/60 phase_boundary_type type; /*+ See phase_boundary_flags typedef enum +*/ 61 61 } phase_boundary; /*+ Phase boundary structure which can store either the 62 62 boundary between a single phase and a multi-phase region, -
trunk/matml/src/ternary/Makefile.am
r460 r463 12 12 13 13 libfreenergy_la_SOURCES = freenergy.c 14 libfreenergy_la_LIBADD = -lm 14 15 libfreenergy_la_CFLAGS = -std=c99 #-DDEBUG 15 16 libfreenergy_la_LDFLAGS = -version-info 0:0:0 16 17 17 18 libgibbs_la_SOURCES = geomview.c qhull.c book.c spinodal.c 18 libgibbs_la_LIBADD = -l m -lqhull libfreenergy.la19 libgibbs_la_LIBADD = -lqhull libfreenergy.la 19 20 libgibbs_la_CFLAGS = -std=c99 -DGEOMVIEW=\"@GEOMVIEW@\" #-DDEBUG 20 21 libgibbs_la_LDFLAGS = -version-info 0:0:0 -
trunk/matml/src/ternary/qhull.c
r460 r463 776 776 if ((*boundaries) [b].n_phases == 3) 777 777 if ((*boundaries) [b].compos == 3 && 778 (*boundaries) [b]. flags== TIE_SIMPLICES &&778 (*boundaries) [b].type == TIE_SIMPLICES && 779 779 (*boundaries) [b].phases[0] == phase0 && 780 780 (*boundaries) [b].phases[1] == phase1 && … … 800 800 (*boundaries) [bound].phases [1] = phase1; 801 801 (*boundaries) [bound].phases [2] = phase2; 802 (*boundaries) [bound]. flags= TIE_SIMPLICES;802 (*boundaries) [bound].type = TIE_SIMPLICES; 803 803 (*boundaries) [bound].n_edges = 1; 804 804 if (!((*boundaries) [bound].edges = malloc (3*sizeof (int)))) … … 927 927 for (b=0; b<*n_bounds; b++) 928 928 if ((*boundaries) [b].n_phases == nphases && 929 (*boundaries) [b]. flags== TIE_SIMPLICES &&929 (*boundaries) [b].type == TIE_SIMPLICES && 930 930 (*boundaries) [b].phases[0] == phase0 && 931 931 (*boundaries) [b].phases[nphases-1] == phase1) … … 952 952 (*boundaries) [bound].n_edges = 0; 953 953 (*boundaries) [bound].edges = NULL; 954 (*boundaries) [bound]. flags= TIE_SIMPLICES;954 (*boundaries) [bound].type = TIE_SIMPLICES; 955 955 (*n_bounds) ++; 956 956 e=0; … … 981 981 for (b=0, bound=-1; b<*n_bounds; b++) 982 982 if ((*boundaries) [b].n_phases == nphases && 983 (*boundaries) [b]. flags== ONE_PHASE_EDGE &&983 (*boundaries) [b].type == ONE_PHASE_EDGE && 984 984 (*boundaries) [b].phases[0] == points[p0].efunc && 985 985 (*boundaries) [b].phases[nphases-1] == points[p2].efunc) … … 1006 1006 (*boundaries) [bound].n_edges = 0; 1007 1007 (*boundaries) [bound].edges = NULL; 1008 (*boundaries) [bound]. flags= ONE_PHASE_EDGE;1008 (*boundaries) [bound].type = ONE_PHASE_EDGE; 1009 1009 (*n_bounds) ++; 1010 1010 e=0; … … 1077 1077 for (b=0; b<*n_bounds; b++) 1078 1078 if ((*boundaries) [b].n_phases == 1 && 1079 (*boundaries) [b]. flags== TIE_SIMPLICES &&1079 (*boundaries) [b].type == TIE_SIMPLICES && 1080 1080 (*boundaries) [b].phases[0] == phase) 1081 1081 bound = b; … … 1099 1099 (*boundaries) [bound].n_edges = 0; 1100 1100 (*boundaries) [bound].edges = NULL; 1101 (*boundaries) [bound]. flags= TIE_SIMPLICES;1101 (*boundaries) [bound].type = TIE_SIMPLICES; 1102 1102 (*n_bounds) ++; 1103 1103 e=0; … … 1124 1124 for (b=0, bound=-1; b<*n_bounds; b++) 1125 1125 if ((*boundaries) [b].n_phases == 1 && 1126 (*boundaries) [b]. flags== ONE_PHASE_EDGE &&1126 (*boundaries) [b].type == ONE_PHASE_EDGE && 1127 1127 (*boundaries) [b].phases[0] == phase) 1128 1128 bound = b; … … 1146 1146 (*boundaries) [bound].n_edges = 0; 1147 1147 (*boundaries) [bound].edges = NULL; 1148 (*boundaries) [bound]. flags= ONE_PHASE_EDGE;1148 (*boundaries) [bound].type = ONE_PHASE_EDGE; 1149 1149 (*n_bounds) ++; 1150 1150 e=0; -
trunk/matml/src/ternary/spinodal.c
r460 r463 337 337 *(spinreturn->phases) = efunc; 338 338 spinreturn->n_phases = 1; 339 spinreturn-> flags= SPINODAL;339 spinreturn->type = SPINODAL; 340 340 } 341 341 -
trunk/matml/src/ternary/square.c
r460 r463 135 135 136 136 printf ("Boundary %d: %d-phase %s, phases:", i, (allbounds+i)->n_phases, 137 ((allbounds+i)-> flags==TIE_SIMPLICES) ? "tie-simplices":"edge");137 ((allbounds+i)->type==TIE_SIMPLICES) ? "tie-simplices":"edge"); 138 138 for (ret=0; ret<(allbounds+i)->n_phases; ret++) 139 139 printf (" %d", (allbounds+i)->phases[ret]); -
trunk/matml/src/ternary/ternary.c
r460 r463 134 134 135 135 printf ("Boundary %d: %d-phase %s, phases:", i, (allbounds+i)->n_phases, 136 ((allbounds+i)-> flags==TIE_SIMPLICES) ? "tie-simplices":"edge");136 ((allbounds+i)->type==TIE_SIMPLICES) ? "tie-simplices":"edge"); 137 137 for (ret=0; ret<(allbounds+i)->n_phases; ret++) 138 138 printf (" %d", (allbounds+i)->phases[ret]);