[tests][material] update material tests with molar density

test/material/fluidsystems/checkfluidsystem.hh

@@ -86,6 +86,7 @@ public:
             BaseFluidState::setSaturation(phaseIdx, 1.0 / numPhases);
             BaseFluidState::setPressure(phaseIdx, 1e5);
             BaseFluidState::setDensity(phaseIdx, 1.0);
+            BaseFluidState::setMolarDensity(phaseIdx, 1.0);
 
             for (int compIdx = 0; compIdx < numComponents; ++compIdx)
             {

test/material/immiscibleflash/test_immiscibleflash.cc

@@ -142,7 +142,9 @@ void completeReferenceFluidState(FluidState &fs,
     paramCache.updateAll(fs);
     for (int phaseIdx = 0; phaseIdx < numPhases; ++ phaseIdx) {
         Scalar rho = FluidSystem::density(fs, paramCache, phaseIdx);
+        Scalar rhoMolar = FluidSystem::molarDensity(fs, paramCache, phaseIdx);
         fs.setDensity(phaseIdx, rho);
+        fs.setMolarDensity(phaseIdx, rhoMolar);
     }
 }
 

test/material/pengrobinson/test_pengrobinson.cc

@@ -71,6 +71,7 @@ Scalar bringOilToSurface(FluidState &surfaceFluidState, Scalar alpha, const Flui
                                                   reservoirFluidState.moleFraction(phaseIdx, compIdx));
             }
             surfaceFluidState.setDensity(phaseIdx, reservoirFluidState.density(phaseIdx));
+            surfaceFluidState.setMolarDensity(phaseIdx, reservoirFluidState.molarDensity(phaseIdx));
             surfaceFluidState.setPressure(phaseIdx, reservoirFluidState.pressure(phaseIdx));
             surfaceFluidState.setSaturation(phaseIdx, 0.0);
         }
@@ -189,7 +190,9 @@ int main(int argc, char** argv)
     // density
     paramCache.updatePhase(fluidState, oPhaseIdx);
     Scalar rho = FluidSystem::density(fluidState, paramCache, oPhaseIdx);
+    Scalar rhoMolar = FluidSystem::molarDensity(fluidState, paramCache, oPhaseIdx);
     fluidState.setDensity(oPhaseIdx, rho);
+    fluidState.setMolarDensity(oPhaseIdx, rhoMolar);
 
     ////////////
     // Calculate the total molarities of the components