00001 { 00002 #include <finiteVolume/rhoEqn.H> 00003 } 00004 { 00005 scalar sumLocalContErr = 00006 sum 00007 ( 00008 mag(rho.internalField() - (psi*p)().internalField()) 00009 )/sum(rho.internalField()); 00010 00011 scalar globalContErr = 00012 sum(rho.internalField() - (psi*p)().internalField()) 00013 /sum(rho.internalField()); 00014 00015 cumulativeContErr += globalContErr; 00016 00017 Info<< "time step continuity errors : sum local = " << sumLocalContErr 00018 << ", global = " << globalContErr 00019 << ", cumulative = " << cumulativeContErr << endl; 00020 }