research/OpenFOAM-4.x
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

patchInjection: Check for 0-patches

Browse source
This commit is contained in:
Henry Weller 2015-09-30 18:17:32 +01:00
parent df85b892a5
commit fb3eddd2cb
1 changed files with 13 additions and 1 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

14
src/regionModels/surfaceFilmModels/submodels/kinematic/injectionModel/patchInjection/patchInjection.C
View file

@ -82,6 +82,13 @@ patchInjection::patchInjection
patchInjectedMasses_.setSize(pbm.size(), 0); patchInjectedMasses_.setSize(pbm.size(), 0);
} }
if (!patchIDs_.size())
{
FatalErrorIn("patchInjection::patchInjection")
<< "No patches selected"
<< exit(FatalError);
}
} }
@ -100,6 +107,9 @@ void patchInjection::correct
scalarField& diameterToInject scalarField& diameterToInject
) )
{ {
// Do not correct if no patches selected
if (!patchIDs_.size()) return;
const scalarField& delta = owner().delta(); const scalarField& delta = owner().delta();
const scalarField& rho = owner().rho(); const scalarField& rho = owner().rho();
const scalarField& magSf = owner().magSf(); const scalarField& magSf = owner().magSf();
@ -160,6 +170,9 @@ void patchInjection::correct
void patchInjection::patchInjectedMassTotals(scalarField& patchMasses) const void patchInjection::patchInjectedMassTotals(scalarField& patchMasses) const
{ {
// Do not correct if no patches selected
if (!patchIDs_.size()) return;
scalarField patchInjectedMasses scalarField patchInjectedMasses
( (
getModelProperty<scalarField> getModelProperty<scalarField>
@ -170,7 +183,6 @@ void patchInjection::patchInjectedMassTotals(scalarField& patchMasses) const
); );
scalarField patchInjectedMassTotals(patchInjectedMasses_); scalarField patchInjectedMassTotals(patchInjectedMasses_);
//combineReduce(patchInjectedMassTotals, plusEqOp<scalarField>());
Pstream::listCombineGather(patchInjectedMassTotals, plusEqOp<scalar>()); Pstream::listCombineGather(patchInjectedMassTotals, plusEqOp<scalar>());
forAll(patchIDs_, pidi) forAll(patchIDs_, pidi)