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)