With the recent changes in #158 the algorithm for the determination of phase euqilibria in pores needed to be rewritten somewhat. The new implementation uses bulk densities rather than Gibbs energies as iteration variable, which is appealing because it avoids additional iterations of the bulk states. However, it appears to have regressed in terms of robustness.