99    static void solve(FluidState &fluidState,
 
  100                      ParameterCache ¶mCache,
 
  103        ComponentVector fugVec;
 
  107        paramCache.updatePhase(fluidState, refPhaseIdx);
 
  108        fluidState.setDensity(refPhaseIdx,
 
  109                              FluidSystem::density(fluidState,
 
  112        fluidState.setMolarDensity(refPhaseIdx,
 
  113                                   FluidSystem::molarDensity(fluidState,
 
  118        for (
int compIdx = 0; compIdx < numComponents; ++compIdx)
 
  120            fluidState.setFugacityCoefficient(refPhaseIdx,
 
  122                                              FluidSystem::fugacityCoefficient(fluidState,
 
  126            fugVec[compIdx] = fluidState.fugacity(refPhaseIdx, compIdx);
 
  130        for (
int phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx)
 
  132            if (phaseIdx == refPhaseIdx)