exaFOAM/foam-extend4.1-coherent-io
Archived
1
1
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

Array access optimisation

Browse source
This commit is contained in:
Hrvoje Jasak 2015-04-07 13:51:12 +01:00
parent 4ade918e1c
commit 52fbcd4029
1 changed files with 4 additions and 7 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

11
src/finiteVolume/finiteVolume/gradSchemes/limitedGradSchemes/cellLimitedGrad/cellLimitedGrads.C
View file

@ -411,7 +411,7 @@ tmp<BlockLduSystem<vector, vector> > cellLimitedGrad<scalar>::fvmGrad
) const ) const
{ {
// Consider doing a calculateLimiter member function since both fvmGrad and // Consider doing a calculateLimiter member function since both fvmGrad and
// grad use almost the same procedure to calculate limiter. VV, 9/June/2014. // grad use almost the same procedure to calculate limiter. VV, 9/June/2014
const fvMesh& mesh = vsf.mesh(); const fvMesh& mesh = vsf.mesh();
tmp<BlockLduSystem<vector, vector> > tbs = basicGradScheme_().fvmGrad(vsf); tmp<BlockLduSystem<vector, vector> > tbs = basicGradScheme_().fvmGrad(vsf);
@ -494,9 +494,6 @@ tmp<BlockLduSystem<vector, vector> > cellLimitedGrad<scalar>::fvmGrad
scalarField maxMinVsf = (1.0/k_ - 1.0)*(maxVsf - minVsf); scalarField maxMinVsf = (1.0/k_ - 1.0)*(maxVsf - minVsf);
maxVsf += maxMinVsf; maxVsf += maxMinVsf;
minVsf -= maxMinVsf; minVsf -= maxMinVsf;
//maxVsf *= 1.0/k_;
//minVsf *= 1.0/k_;
} }
@ -600,6 +597,8 @@ tmp<BlockLduSystem<vector, vector> > cellLimitedGrad<scalar>::fvmGrad
const fvPatchScalarField& pf = vsf.boundaryField()[patchI]; const fvPatchScalarField& pf = vsf.boundaryField()[patchI];
const fvPatch& patch = pf.patch(); const fvPatch& patch = pf.patch();
const labelList& fc = patch.faceCells();
if (patch.coupled()) if (patch.coupled())
{ {
CoeffField<vector>::linearTypeField& pcoupleUpper = CoeffField<vector>::linearTypeField& pcoupleUpper =
@ -612,9 +611,7 @@ tmp<BlockLduSystem<vector, vector> > cellLimitedGrad<scalar>::fvmGrad
forAll(pf, faceI) forAll(pf, faceI)
{ {
label cellI = patch.faceCells()[faceI]; pcoupleUpper[faceI] *= lfIn[fc[faceI]];
pcoupleUpper[faceI] *= lfIn[cellI];
pcoupleLower[faceI] *= lfNei[faceI]; pcoupleLower[faceI] *= lfNei[faceI];
} }
} }