00001 { 00002 fvScalarMatrix rhoEqn 00003 ( 00004 fvm::ddt(rho) 00005 + fvm::div(phiv, rho) 00006 ); 00007 00008 rhoEqn.solve(); 00009 00010 phi = rhoEqn.flux(); 00011 00012 Info<< "max-min rho: " << max(rho).value() 00013 << " " << min(rho).value() << endl; 00014 00015 rho == max(rho, rhoMin); 00016 }