applications/solvers/multiphase/cavitatingFoam/gammaPsi.H

   1 {
   2     gamma = max(min((rho - rholSat)/(rhovSat - rholSat), scalar(1)), scalar(0));
   3
   4     Info<< "max-min gamma: " << max(gamma).value()
   5         << " " << min(gamma).value() << endl;
   6
   7     psiModel->correct();
   8
   9     //Info<< "min a: " << 1.0/sqrt(max(psi)).value() << endl;
  10 }