Changeset 386

Show
Ignore:
Timestamp:
08/19/08 19:40:12 (3 months ago)
Author:
powell
Message:

One more refining step, better loop_max considering refining, better T

Files:
1 modified

Legend:

Unmodified
Added
Removed
  • trunk/matml/src/ternary/ternary.c

    r383 r386  
    2828int main (int argc, char *argv[]) 
    2929{ 
    30   int i,j, index, loop_max=50, numpoints, *verts, hullnumverts,*hullverts=NULL; 
     30  int i,j, index, loop_max=20, numpoints, *verts, hullnumverts,*hullverts=NULL; 
    3131  char gv_version[100]; 
    3232  FILE *pfd = NULL; 
    33   double T=0.75, P=1.; 
     33  double T=0.7, P=1.; 
    3434  ternary_point *points; 
    3535  /* eparams: R,T0, G1@T0,G2,G3, C1,C2,C3, M1,M2,M3, O12,O13,O23,O123 */ 
     
    8787  if (i=hullCalculate (3, points, numpoints)) 
    8888    { printf ("main: Error %d in hullCalculate\n", i); exit (i); } 
     89  if (i=hullRefine (&points, &numpoints, allparams, 2, T, P, 
     90                    1e-10, 1e-7, 0)) 
     91    { printf ("main: error %d in hullRefine\n", i); exit (i); } 
    8992  if (i=hullRefine (&points, &numpoints, allparams, 2, T, P, 
    9093                    1e-10, 1e-7, 0))