More fixes for objective function and constraint Jacobian calculation in point equilibrium calculation. Still need to do more work on the sampling scheme to ensure consistent convergence.
Fix for barycentric coordinate computation to estimate the phase fractions.