00001 { 00002 for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++) 00003 { 00004 fvScalarMatrix tEqn 00005 ( 00006 -fvm::laplacian(K, T) 00007 ); 00008 tEqn.relax(); 00009 eqnResidual = tEqn.solve().initialResidual(); 00010 maxResidual = max(eqnResidual, maxResidual); 00011 00012 } 00013 00014 Info<< "Min/max T:" << min(T).value() << ' ' 00015 << max(T).value() << endl; 00016 }