tmp<volScalarField> CoEulerDdtScheme<Type>::CorDeltaT() const { surfaceScalarField cofrDeltaT = CofrDeltaT(); tmp<volScalarField> tcorDeltaT ( new volScalarField ( IOobject ( "CorDeltaT", cofrDeltaT.instance(), mesh() ), mesh(), dimensionedScalar("CorDeltaT", cofrDeltaT.dimensions(), 0.0), zeroGradientFvPatchScalarField::typeName ) ); volScalarField& corDeltaT = tcorDeltaT(); const unallocLabelList& owner = mesh().owner(); const unallocLabelList& neighbour = mesh().neighbour(); forAll(owner, faceI) { corDeltaT[owner[faceI]] = max(corDeltaT[owner[faceI]], cofrDeltaT[faceI]); corDeltaT[neighbour[faceI]] = max(corDeltaT[neighbour[faceI]], cofrDeltaT[faceI]); }
tmp<volScalarField> CoEulerDdtScheme<Type>::CorDeltaT() const { const surfaceScalarField cofrDeltaT(CofrDeltaT()); tmp<volScalarField> tcorDeltaT ( new volScalarField ( IOobject ( "CorDeltaT", cofrDeltaT.instance(), mesh() ), mesh(), dimensionedScalar("CorDeltaT", cofrDeltaT.dimensions(), 0.0), extrapolatedCalculatedFvPatchScalarField::typeName ) ); volScalarField& corDeltaT = tcorDeltaT.ref(); const labelUList& owner = mesh().owner(); const labelUList& neighbour = mesh().neighbour(); forAll(owner, faceI) { corDeltaT[owner[faceI]] = max(corDeltaT[owner[faceI]], cofrDeltaT[faceI]); corDeltaT[neighbour[faceI]] = max(corDeltaT[neighbour[faceI]], cofrDeltaT[faceI]); }