Changeset 436

Show
Ignore:
Timestamp:
03/16/2009 03:38:26 PM (3 years ago)
Author:
powell
Message:

Important fixes to tie line/phase edge test.

Files:
1 modified

Legend:

Unmodified
Added
Removed

  • trunk/matml/src/ternary/qhull.c

    r434 r436  
    855855                    nphases=1; 
    856856 
    857                     if (free_energy ((C02+C12)/2,(C03,C13)/2,T,P,eparams+phase0) < 
    858                         G0+G1/2) 
     857                    if (free_energy ((C02+C12)/2,(C03+C13)/2,T,P,eparams+phase0) < 
     858                        (G0+G1)/2) 
    859859                      { 
    860860                        printf ("v2 alone, Gav=%g, Gmid=%g\n", 
    861                                 G0+G1/2, free_energy ((C02+C12)/2,(C03,C13)/2,T,P,eparams+phase0)); 
     861                                (G0+G1)/2, free_energy ((C02+C12)/2,(C03+C13)/2,T,P,eparams+phase0)); 
    862862                        p0 = MIN (v0, v1); 
    863863                        p1 = MAX (v0, v1); 
    864864                        p2 = v2; 
    865865                      } 
    866                     else if (free_energy ((C02+C22)/2,(C03,C23)/2,T,P,eparams+phase0) < 
    867                         G0+G2/2) 
     866                    else if (free_energy ((C02+C22)/2,(C03+C23)/2,T,P,eparams+phase0) < 
     867                             (G0+G2)/2) 
    868868                      { 
    869869                        printf ("v1 alone, Gav=%g, Gmid=%g\n", 
    870                                 G0+G2/2, free_energy ((C02+C22)/2,(C03,C23)/2,T,P,eparams+phase0)); 
     870                                (G0+G2)/2, free_energy ((C02+C22)/2,(C03+C23)/2,T,P,eparams+phase0)); 
    871871                        p0 = MIN (v0, v2); 
    872872                        p1 = MAX (v0, v2); 
    … …  
    876876                      { 
    877877                        printf ("v0 alone, Gav=%g, Gmid=%g\n", 
    878                                 G1+G2/2, free_energy ((C22+C12)/2,(C23,C13)/2,T,P,eparams+phase0)); 
     878                                (G1+G2)/2, free_energy ((C22+C12)/2,(C23+C13)/2,T,P,eparams+phase0)); 
    879879                        p0 = MIN (v1, v2); 
    880880                        p1 = MAX (v1, v2);