Compare commits
1 Commits
feature-me
...
issue-3348
| Author | SHA1 | Date | |
|---|---|---|---|
| bd41a2abdb |
@ -36,11 +36,12 @@ Description
|
||||
if (adjustTimeStep)
|
||||
{
|
||||
scalar maxDeltaTFact = maxCo/(CoNum + StCoNum + SMALL);
|
||||
scalar deltaTFact = min(min(maxDeltaTFact, 1.0 + 0.1*maxDeltaTFact), 1.2);
|
||||
scalar deltaTFact =
|
||||
Foam::min(Foam::min(maxDeltaTFact, 1.0 + 0.1*maxDeltaTFact), 1.2);
|
||||
|
||||
runTime.setDeltaT
|
||||
(
|
||||
min
|
||||
Foam::min
|
||||
(
|
||||
deltaTFact*runTime.deltaTValue(),
|
||||
maxDeltaT
|
||||
|
||||
@ -1,5 +1,5 @@
|
||||
if (adjustTimeStep)
|
||||
{
|
||||
runTime.setDeltaT(min(dtChem, maxDeltaT));
|
||||
runTime.setDeltaT(Foam::min(dtChem, maxDeltaT));
|
||||
Info<< "deltaT = " << runTime.deltaTValue() << endl;
|
||||
}
|
||||
|
||||
@ -54,9 +54,18 @@ if (adjustTimeStep)
|
||||
|
||||
runTime.setDeltaT
|
||||
(
|
||||
min
|
||||
Foam::min
|
||||
(
|
||||
dt0*min(min(TFactorFluid, min(TFactorFilm, TFactorSolid)), 1.2),
|
||||
dt0
|
||||
* Foam::min
|
||||
(
|
||||
Foam::min
|
||||
(
|
||||
TFactorFluid,
|
||||
Foam::min(TFactorFilm, TFactorSolid)
|
||||
),
|
||||
1.2
|
||||
),
|
||||
maxDeltaT
|
||||
)
|
||||
);
|
||||
|
||||
@ -75,8 +75,10 @@
|
||||
)
|
||||
{
|
||||
rDeltaT =
|
||||
(
|
||||
rDeltaT0
|
||||
*max(rDeltaT/rDeltaT0, scalar(1) - rDeltaTDampingCoeff);
|
||||
* Foam::max(rDeltaT/rDeltaT0, scalar(1) - rDeltaTDampingCoeff)
|
||||
);
|
||||
|
||||
Info<< "Damped flow time scale min/max = "
|
||||
<< gMin(1/rDeltaT.primitiveField())
|
||||
|
||||
@ -1,7 +1,7 @@
|
||||
scalar CoNum = -GREAT;
|
||||
forAll(fluidRegions, regionI)
|
||||
{
|
||||
CoNum = max
|
||||
CoNum = Foam::max
|
||||
(
|
||||
compressibleCourantNo
|
||||
(
|
||||
|
||||
@ -31,7 +31,7 @@
|
||||
);
|
||||
|
||||
|
||||
CoNum =
|
||||
scalar maxPhiCoNum =
|
||||
0.5*gMax
|
||||
(
|
||||
sumPhi/fluidRegions[regioni].V().field()
|
||||
@ -41,9 +41,9 @@
|
||||
(
|
||||
fvc::surfaceSum(mag(phi1 - phi2))().primitiveField()
|
||||
/ fluidRegions[regioni].V().field()
|
||||
)*runTime.deltaTValue(),
|
||||
)*runTime.deltaTValue();
|
||||
|
||||
CoNum = max(UrCoNum, CoNum);
|
||||
CoNum = Foam::max(CoNum, Foam::max(maxPhiCoNum, UrCoNum));
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
@ -2,7 +2,7 @@
|
||||
|
||||
forAll(fluidRegions, regioni)
|
||||
{
|
||||
CoNum = max
|
||||
CoNum = Foam::max
|
||||
(
|
||||
compressibleCourantNo
|
||||
(
|
||||
@ -17,7 +17,7 @@
|
||||
/*
|
||||
forAll(porousFluidRegions, porousi)
|
||||
{
|
||||
CoNum = max
|
||||
CoNum = Foam::max
|
||||
(
|
||||
compressibleCourantNo
|
||||
(
|
||||
|
||||
@ -47,10 +47,10 @@ if (adjustTimeStep)
|
||||
|
||||
runTime.setDeltaT
|
||||
(
|
||||
min
|
||||
Foam::min
|
||||
(
|
||||
min(maxCo/CoNum, maxDi/DiNum)*runTime.deltaTValue(),
|
||||
min(runTime.deltaTValue(), maxDeltaT)
|
||||
Foam::min(maxCo/CoNum, maxDi/DiNum)*runTime.deltaTValue(),
|
||||
Foam::min(runTime.deltaTValue(), maxDeltaT)
|
||||
)
|
||||
);
|
||||
Info<< "deltaT = " << runTime.deltaTValue() << endl;
|
||||
|
||||
@ -49,17 +49,17 @@ if (adjustTimeStep)
|
||||
scalar maxDeltaTSolid = maxDi/(DiNum + SMALL);
|
||||
|
||||
scalar deltaTFluid =
|
||||
min
|
||||
Foam::min
|
||||
(
|
||||
min(maxDeltaTFluid, 1.0 + 0.1*maxDeltaTFluid),
|
||||
Foam::min(maxDeltaTFluid, 1.0 + 0.1*maxDeltaTFluid),
|
||||
1.2
|
||||
);
|
||||
|
||||
runTime.setDeltaT
|
||||
(
|
||||
min
|
||||
Foam::min
|
||||
(
|
||||
min(deltaTFluid, maxDeltaTSolid)*runTime.deltaTValue(),
|
||||
Foam::min(deltaTFluid, maxDeltaTSolid)*runTime.deltaTValue(),
|
||||
maxDeltaT
|
||||
)
|
||||
);
|
||||
|
||||
@ -22,7 +22,7 @@ forAll(solidRegions, i)
|
||||
tmp<volScalarField> trho = thermo.rho();
|
||||
const volScalarField& rho = trho();
|
||||
|
||||
DiNum = max
|
||||
DiNum = Foam::max
|
||||
(
|
||||
solidRegionDiffNo
|
||||
(
|
||||
|
||||
@ -17,7 +17,7 @@ scalar DiNum = -GREAT;
|
||||
tmp<volScalarField> trho = thermo.rho();
|
||||
const volScalarField& rho = trho();
|
||||
|
||||
DiNum = max
|
||||
DiNum = Foam::max
|
||||
(
|
||||
solidRegionDiffNo
|
||||
(
|
||||
|
||||
@ -36,13 +36,16 @@ Description
|
||||
if (adjustTimeStep)
|
||||
{
|
||||
const scalar maxDeltaTFact =
|
||||
min(maxCo/(CoNum + SMALL), maxCo/(surfaceFilm.CourantNumber() + SMALL));
|
||||
Foam::min
|
||||
(
|
||||
maxCo/(CoNum + SMALL), maxCo/(surfaceFilm.CourantNumber() + SMALL)
|
||||
);
|
||||
const scalar deltaTFact =
|
||||
min(min(maxDeltaTFact, 1.0 + 0.1*maxDeltaTFact), 1.2);
|
||||
Foam::min(Foam::min(maxDeltaTFact, 1.0 + 0.1*maxDeltaTFact), 1.2);
|
||||
|
||||
runTime.setDeltaT
|
||||
(
|
||||
min
|
||||
Foam::min
|
||||
(
|
||||
deltaTFact*runTime.deltaTValue(),
|
||||
maxDeltaT
|
||||
|
||||
@ -36,13 +36,14 @@ Description
|
||||
if (adjustTimeStep)
|
||||
{
|
||||
scalar maxDeltaTFact =
|
||||
min(maxCo/(CoNum + SMALL), maxAlphaCo/(alphaCoNum + SMALL));
|
||||
Foam::min(maxCo/(CoNum + SMALL), maxAlphaCo/(alphaCoNum + SMALL));
|
||||
|
||||
scalar deltaTFact = min(min(maxDeltaTFact, 1.0 + 0.1*maxDeltaTFact), 1.2);
|
||||
scalar deltaTFact =
|
||||
Foam::min(Foam::min(maxDeltaTFact, 1.0 + 0.1*maxDeltaTFact), 1.2);
|
||||
|
||||
runTime.setDeltaT
|
||||
(
|
||||
min
|
||||
Foam::min
|
||||
(
|
||||
deltaTFact*runTime.deltaTValue(),
|
||||
maxDeltaT
|
||||
|
||||
@ -36,13 +36,14 @@ Description
|
||||
if (adjustTimeStep)
|
||||
{
|
||||
scalar maxDeltaTFact =
|
||||
min(maxCo/(CoNum + SMALL), maxAcousticCo/(acousticCoNum + SMALL));
|
||||
Foam::min(maxCo/(CoNum + SMALL), maxAcousticCo/(acousticCoNum + SMALL));
|
||||
|
||||
scalar deltaTFact = min(min(maxDeltaTFact, 1.0 + 0.1*maxDeltaTFact), 1.2);
|
||||
scalar deltaTFact =
|
||||
Foam::min(Foam::min(maxDeltaTFact, 1.0 + 0.1*maxDeltaTFact), 1.2);
|
||||
|
||||
runTime.setDeltaT
|
||||
(
|
||||
min
|
||||
Foam::min
|
||||
(
|
||||
deltaTFact*runTime.deltaTValue(),
|
||||
maxDeltaT
|
||||
|
||||
@ -37,11 +37,15 @@ if (adjustTimeStep)
|
||||
if (CoNum > SMALL)
|
||||
{
|
||||
scalar maxDeltaTFact =
|
||||
min(maxCo/(CoNum + SMALL), maxAcousticCo/(acousticCoNum + SMALL));
|
||||
Foam::min
|
||||
(
|
||||
maxCo/(CoNum + SMALL),
|
||||
maxAcousticCo/(acousticCoNum + SMALL)
|
||||
);
|
||||
|
||||
runTime.setDeltaT
|
||||
(
|
||||
min
|
||||
Foam::min
|
||||
(
|
||||
maxDeltaTFact*runTime.deltaTValue(),
|
||||
maxDeltaT
|
||||
|
||||
@ -26,12 +26,12 @@ forAll(dgdt, celli)
|
||||
{
|
||||
if (dgdt[celli] > 0.0)
|
||||
{
|
||||
Sp[celli] -= dgdt[celli]/max(1.0 - alpha1[celli], 1e-4);
|
||||
Su[celli] += dgdt[celli]/max(1.0 - alpha1[celli], 1e-4);
|
||||
Sp[celli] -= dgdt[celli]/Foam::max(1.0 - alpha1[celli], 1e-4);
|
||||
Su[celli] += dgdt[celli]/Foam::max(1.0 - alpha1[celli], 1e-4);
|
||||
}
|
||||
else if (dgdt[celli] < 0.0)
|
||||
{
|
||||
Sp[celli] += dgdt[celli]/max(alpha1[celli], 1e-4);
|
||||
Sp[celli] += dgdt[celli]/Foam::max(alpha1[celli], 1e-4);
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
@ -26,12 +26,12 @@ forAll(dgdt, celli)
|
||||
{
|
||||
if (dgdt[celli] > 0.0)
|
||||
{
|
||||
Sp[celli] -= dgdt[celli]/max(1.0 - alpha1[celli], 1e-4);
|
||||
Su[celli] += dgdt[celli]/max(1.0 - alpha1[celli], 1e-4);
|
||||
Sp[celli] -= dgdt[celli]/Foam::max(1.0 - alpha1[celli], 1e-4);
|
||||
Su[celli] += dgdt[celli]/Foam::max(1.0 - alpha1[celli], 1e-4);
|
||||
}
|
||||
else if (dgdt[celli] < 0.0)
|
||||
{
|
||||
Sp[celli] += dgdt[celli]/max(alpha1[celli], 1e-4);
|
||||
Sp[celli] += dgdt[celli]/Foam::max(alpha1[celli], 1e-4);
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
@ -26,12 +26,12 @@ forAll(dgdt, celli)
|
||||
{
|
||||
if (dgdt[celli] > 0.0)
|
||||
{
|
||||
Sp[celli] -= dgdt[celli]/max(1.0 - alpha1[celli], 1e-4);
|
||||
Su[celli] += dgdt[celli]/max(1.0 - alpha1[celli], 1e-4);
|
||||
Sp[celli] -= dgdt[celli]/Foam::max(1.0 - alpha1[celli], 1e-4);
|
||||
Su[celli] += dgdt[celli]/Foam::max(1.0 - alpha1[celli], 1e-4);
|
||||
}
|
||||
else if (dgdt[celli] < 0.0)
|
||||
{
|
||||
Sp[celli] += dgdt[celli]/max(alpha1[celli], 1e-4);
|
||||
Sp[celli] += dgdt[celli]/Foam::max(alpha1[celli], 1e-4);
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
@ -36,13 +36,13 @@ Description
|
||||
if (adjustTimeStep)
|
||||
{
|
||||
scalar maxDeltaTFact =
|
||||
min
|
||||
Foam::min
|
||||
(
|
||||
maxCo/(CoNum + SMALL),
|
||||
min
|
||||
Foam::min
|
||||
(
|
||||
maxAlphaCo/(alphaCoNum + SMALL),
|
||||
min
|
||||
Foam::min
|
||||
(
|
||||
maxAlphaDdt/(ddtAlphaNum + SMALL),
|
||||
maxDi/(DiNum + SMALL)
|
||||
@ -50,11 +50,12 @@ if (adjustTimeStep)
|
||||
)
|
||||
);
|
||||
|
||||
scalar deltaTFact = min(min(maxDeltaTFact, 1.0 + 0.1*maxDeltaTFact), 1.2);
|
||||
scalar deltaTFact =
|
||||
Foam::min(Foam::min(maxDeltaTFact, 1.0 + 0.1*maxDeltaTFact), 1.2);
|
||||
|
||||
runTime.setDeltaT
|
||||
(
|
||||
min
|
||||
Foam::min
|
||||
(
|
||||
deltaTFact*runTime.deltaTValue(),
|
||||
maxDeltaT
|
||||
|
||||
@ -8,5 +8,5 @@
|
||||
|
||||
Info<< "Max Ur Courant Number = " << UrCoNum << endl;
|
||||
|
||||
CoNum = max(CoNum, UrCoNum);
|
||||
CoNum = Foam::max(CoNum, UrCoNum);
|
||||
}
|
||||
|
||||
@ -8,5 +8,5 @@
|
||||
|
||||
Info<< "Max Ur Courant Number = " << UrCoNum << endl;
|
||||
|
||||
CoNum = max(CoNum, UrCoNum);
|
||||
CoNum = Foam::max(CoNum, UrCoNum);
|
||||
}
|
||||
|
||||
@ -108,13 +108,13 @@ Input
|
||||
if (numNbrs == 1)
|
||||
{
|
||||
//if (pointSetPtr) pointSetPtr->insert(mesh.edges()[meshEdgei]);
|
||||
labelTyp = max(labelTyp, TopoType::OPEN);
|
||||
labelTyp = Foam::max(labelTyp, TopoType::OPEN);
|
||||
}
|
||||
else if (numNbrs == 0 || numNbrs > 2)
|
||||
{
|
||||
if (pointSetPtr) pointSetPtr->insert(mesh.edges()[meshEdgei]);
|
||||
if (badEdgesPtr) badEdgesPtr->insert(edgei);
|
||||
labelTyp = max(labelTyp, TopoType::ILLEGAL);
|
||||
labelTyp = Foam::max(labelTyp, TopoType::ILLEGAL);
|
||||
}
|
||||
}
|
||||
reduce(labelTyp, maxOp<label>());
|
||||
|
||||
@ -368,7 +368,8 @@ void subsetTopoSets
|
||||
|
||||
Info<< "Subsetting " << set.type() << ' ' << set.name() << endl;
|
||||
|
||||
labelHashSet subset(2*min(set.size(), map.size()));
|
||||
labelHashSet subset;
|
||||
subset.reserve(Foam::min(set.size(), map.size()));
|
||||
|
||||
// Map the data
|
||||
forAll(map, i)
|
||||
|
||||
@ -142,26 +142,26 @@ scalar getEdgeStats(const primitiveMesh& mesh, const direction excludeCmpt)
|
||||
|
||||
if (mag(eVec & x) > 1-edgeTol)
|
||||
{
|
||||
minX = min(minX, eMag);
|
||||
maxX = max(maxX, eMag);
|
||||
minX = Foam::min(minX, eMag);
|
||||
maxX = Foam::max(maxX, eMag);
|
||||
nX++;
|
||||
}
|
||||
else if (mag(eVec & y) > 1-edgeTol)
|
||||
{
|
||||
minY = min(minY, eMag);
|
||||
maxY = max(maxY, eMag);
|
||||
minY = Foam::min(minY, eMag);
|
||||
maxY = Foam::max(maxY, eMag);
|
||||
nY++;
|
||||
}
|
||||
else if (mag(eVec & z) > 1-edgeTol)
|
||||
{
|
||||
minZ = min(minZ, eMag);
|
||||
maxZ = max(maxZ, eMag);
|
||||
minZ = Foam::min(minZ, eMag);
|
||||
maxZ = Foam::max(maxZ, eMag);
|
||||
nZ++;
|
||||
}
|
||||
else
|
||||
{
|
||||
minOther = min(minOther, eMag);
|
||||
maxOther = max(maxOther, eMag);
|
||||
minOther = Foam::min(minOther, eMag);
|
||||
maxOther = Foam::max(maxOther, eMag);
|
||||
}
|
||||
}
|
||||
|
||||
@ -179,19 +179,19 @@ scalar getEdgeStats(const primitiveMesh& mesh, const direction excludeCmpt)
|
||||
|
||||
if (excludeCmpt == 0)
|
||||
{
|
||||
return min(minY, min(minZ, minOther));
|
||||
return Foam::min(minY, Foam::min(minZ, minOther));
|
||||
}
|
||||
else if (excludeCmpt == 1)
|
||||
{
|
||||
return min(minX, min(minZ, minOther));
|
||||
return Foam::min(minX, Foam::min(minZ, minOther));
|
||||
}
|
||||
else if (excludeCmpt == 2)
|
||||
{
|
||||
return min(minX, min(minY, minOther));
|
||||
return Foam::min(minX, Foam::min(minY, minOther));
|
||||
}
|
||||
else
|
||||
{
|
||||
return min(minX, min(minY, min(minZ, minOther)));
|
||||
return Foam::min(minX, Foam::min(minY, Foam::min(minZ, minOther)));
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
@ -271,7 +271,7 @@ int main(int argc, char *argv[])
|
||||
if (blockPFacePointi != blockPFacePointi2)
|
||||
{
|
||||
sqrMergeTol =
|
||||
min
|
||||
Foam::min
|
||||
(
|
||||
sqrMergeTol,
|
||||
magSqr
|
||||
@ -338,16 +338,16 @@ int main(int argc, char *argv[])
|
||||
blockNFacePoints[blockNFacePointi]
|
||||
+ blockOffsets[blockNlabel];
|
||||
|
||||
label minPN = min(PpointLabel, NpointLabel);
|
||||
label minPN = Foam::min(PpointLabel, NpointLabel);
|
||||
|
||||
if (pointMergeList[PpointLabel] != -1)
|
||||
{
|
||||
minPN = min(minPN, pointMergeList[PpointLabel]);
|
||||
minPN = Foam::min(minPN, pointMergeList[PpointLabel]);
|
||||
}
|
||||
|
||||
if (pointMergeList[NpointLabel] != -1)
|
||||
{
|
||||
minPN = min(minPN, pointMergeList[NpointLabel]);
|
||||
minPN = Foam::min(minPN, pointMergeList[NpointLabel]);
|
||||
}
|
||||
|
||||
pointMergeList[PpointLabel]
|
||||
@ -411,7 +411,7 @@ int main(int argc, char *argv[])
|
||||
|
||||
pointMergeList[PpointLabel]
|
||||
= pointMergeList[NpointLabel]
|
||||
= min
|
||||
= Foam::min
|
||||
(
|
||||
pointMergeList[PpointLabel],
|
||||
pointMergeList[NpointLabel]
|
||||
|
||||
@ -350,7 +350,7 @@ mtype {space}"MTYPE:"{space}
|
||||
// Find out how many labels are expected. If less or equal to
|
||||
// seven, read them all and finish with it. If there is more,
|
||||
// set read of the next line
|
||||
label labelsToRead = min(8, nVertices);
|
||||
label labelsToRead = Foam::min(8, nVertices);
|
||||
label labelI = 0;
|
||||
for (; labelI < labelsToRead; labelI++)
|
||||
{
|
||||
@ -387,7 +387,7 @@ mtype {space}"MTYPE:"{space}
|
||||
|
||||
labelList& curLabels = cellLabels[curNumberOfCells];
|
||||
|
||||
label labelsToRead = min
|
||||
label labelsToRead = Foam::min
|
||||
(
|
||||
(nCellContinuationLines + 1)*7,
|
||||
curLabels.size()
|
||||
|
||||
@ -269,7 +269,7 @@ void readCells
|
||||
label maxUnvPoint = 0;
|
||||
forAll(unvPointID, pointi)
|
||||
{
|
||||
maxUnvPoint = max(maxUnvPoint, unvPointID[pointi]);
|
||||
maxUnvPoint = Foam::max(maxUnvPoint, unvPointID[pointi]);
|
||||
}
|
||||
labelList unvToFoam(invert(maxUnvPoint+1, unvPointID));
|
||||
|
||||
@ -784,7 +784,7 @@ int main(int argc, char *argv[])
|
||||
label maxUnvPoint = 0;
|
||||
forAll(unvPointID, pointi)
|
||||
{
|
||||
maxUnvPoint = max(maxUnvPoint, unvPointID[pointi]);
|
||||
maxUnvPoint = Foam::max(maxUnvPoint, unvPointID[pointi]);
|
||||
}
|
||||
labelList unvToFoam(invert(maxUnvPoint+1, unvPointID));
|
||||
|
||||
|
||||
@ -8,10 +8,10 @@
|
||||
{
|
||||
if (kivaVersion == kiva3v)
|
||||
{
|
||||
regionIndex = max
|
||||
regionIndex = Foam::max
|
||||
(
|
||||
max(idface[quadFace[0]], idface[quadFace[1]]),
|
||||
max(idface[quadFace[2]], idface[quadFace[3]])
|
||||
Foam::max(idface[quadFace[0]], idface[quadFace[1]]),
|
||||
Foam::max(idface[quadFace[2]], idface[quadFace[3]])
|
||||
);
|
||||
|
||||
if (regionIndex > 0)
|
||||
|
||||
@ -148,7 +148,7 @@ for (label i=0; i<nPoints; i++)
|
||||
end = pointMap[end];
|
||||
}
|
||||
|
||||
label minLabel = min(start, end);
|
||||
label minLabel = Foam::min(start, end);
|
||||
|
||||
pointMap[start] = pointMap[end] = minLabel;
|
||||
}
|
||||
@ -331,7 +331,7 @@ if
|
||||
{
|
||||
forAll(pf, pfi)
|
||||
{
|
||||
minz = min(minz, points[pf[pfi]].z());
|
||||
minz = Foam::min(minz, points[pf[pfi]].z());
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
@ -188,7 +188,7 @@ int main(int argc, char *argv[])
|
||||
}
|
||||
|
||||
|
||||
maxPatch = max(maxPatch, patchi);
|
||||
maxPatch = Foam::max(maxPatch, patchi);
|
||||
|
||||
triFace tri(readLabel(str)-1, readLabel(str)-1, readLabel(str)-1);
|
||||
|
||||
|
||||
@ -556,8 +556,8 @@ void calcEdgeMinMaxZone
|
||||
forAll(eFaces, i)
|
||||
{
|
||||
label zoneI = mappedZoneID[eFaces[i]];
|
||||
minZoneID[edgeI] = min(minZoneID[edgeI], zoneI);
|
||||
maxZoneID[edgeI] = max(maxZoneID[edgeI], zoneI);
|
||||
minZoneID[edgeI] = Foam::min(minZoneID[edgeI], zoneI);
|
||||
maxZoneID[edgeI] = Foam::max(maxZoneID[edgeI], zoneI);
|
||||
}
|
||||
}
|
||||
}
|
||||
@ -813,8 +813,8 @@ void addCoupledPatches
|
||||
forAll(eFaces, i)
|
||||
{
|
||||
label proci = procID[eFaces[i]];
|
||||
minProcID[edgeI] = min(minProcID[edgeI], proci);
|
||||
maxProcID[edgeI] = max(maxProcID[edgeI], proci);
|
||||
minProcID[edgeI] = Foam::min(minProcID[edgeI], proci);
|
||||
maxProcID[edgeI] = Foam::max(maxProcID[edgeI], proci);
|
||||
}
|
||||
}
|
||||
}
|
||||
@ -1291,7 +1291,7 @@ void extrudeGeometricProperties
|
||||
label celli = regionMesh.faceOwner()[facei];
|
||||
if (regionMesh.isInternalFace(facei))
|
||||
{
|
||||
celli = max(celli, regionMesh.faceNeighbour()[facei]);
|
||||
celli = Foam::max(celli, regionMesh.faceNeighbour()[facei]);
|
||||
}
|
||||
|
||||
// Calculate layer from cell numbering (see createShellMesh)
|
||||
@ -2210,8 +2210,8 @@ int main(int argc, char *argv[])
|
||||
|
||||
if (addSidePatches && (zone0 != zone1)) // || (cos(angle) > blabla))
|
||||
{
|
||||
label minZone = min(zone0,zone1);
|
||||
label maxZone = max(zone0,zone1);
|
||||
label minZone = Foam::min(zone0,zone1);
|
||||
label maxZone = Foam::max(zone0,zone1);
|
||||
label index = minZone*zoneNames.size()+maxZone;
|
||||
|
||||
ePatches.setSize(eFaces.size());
|
||||
|
||||
@ -26,7 +26,7 @@ void maxFaceToCell
|
||||
const cell& cFaces = cells[cellI];
|
||||
forAll(cFaces, i)
|
||||
{
|
||||
cellFld[cellI] = max(cellFld[cellI], faceData[cFaces[i]]);
|
||||
cellFld[cellI] = Foam::max(cellFld[cellI], faceData[cFaces[i]]);
|
||||
}
|
||||
}
|
||||
|
||||
@ -56,7 +56,7 @@ void minFaceToCell
|
||||
const cell& cFaces = cells[cellI];
|
||||
forAll(cFaces, i)
|
||||
{
|
||||
cellFld[cellI] = min(cellFld[cellI], faceData[cFaces[i]]);
|
||||
cellFld[cellI] = Foam::min(cellFld[cellI], faceData[cFaces[i]]);
|
||||
}
|
||||
}
|
||||
|
||||
@ -87,8 +87,8 @@ void minFaceToCell
|
||||
// Internal faces
|
||||
forAll(own, facei)
|
||||
{
|
||||
cellFld[own[facei]] = min(cellFld[own[facei]], faceData[facei]);
|
||||
cellFld[nei[facei]] = min(cellFld[nei[facei]], faceData[facei]);
|
||||
cellFld[own[facei]] = Foam::min(cellFld[own[facei]], faceData[facei]);
|
||||
cellFld[nei[facei]] = Foam::min(cellFld[nei[facei]], faceData[facei]);
|
||||
}
|
||||
|
||||
// Patch faces
|
||||
@ -99,7 +99,7 @@ void minFaceToCell
|
||||
|
||||
forAll(fc, i)
|
||||
{
|
||||
cellFld[fc[i]] = min(cellFld[fc[i]], fvp[i]);
|
||||
cellFld[fc[i]] = Foam::min(cellFld[fc[i]], fvp[i]);
|
||||
}
|
||||
}
|
||||
|
||||
@ -178,7 +178,14 @@ void Foam::writeFields
|
||||
(
|
||||
radToDeg
|
||||
(
|
||||
Foam::acos(min(scalar(1), max(scalar(-1), faceOrthogonality)))
|
||||
Foam::acos
|
||||
(
|
||||
Foam::min
|
||||
(
|
||||
scalar(1),
|
||||
Foam::max(scalar(-1), faceOrthogonality)
|
||||
)
|
||||
)
|
||||
)
|
||||
);
|
||||
|
||||
@ -534,7 +541,7 @@ void Foam::writeFields
|
||||
ownCc,
|
||||
fc
|
||||
).quality();
|
||||
ownVol = min(ownVol, tetQual);
|
||||
ownVol = Foam::min(ownVol, tetQual);
|
||||
}
|
||||
}
|
||||
if (mesh.isInternalFace(facei))
|
||||
@ -550,7 +557,7 @@ void Foam::writeFields
|
||||
fc,
|
||||
neiCc
|
||||
).quality();
|
||||
neiVol = min(neiVol, tetQual);
|
||||
neiVol = Foam::min(neiVol, tetQual);
|
||||
}
|
||||
}
|
||||
}
|
||||
@ -602,8 +609,10 @@ void Foam::writeFields
|
||||
// Internal faces
|
||||
forAll(own, facei)
|
||||
{
|
||||
cellFld[own[facei]] = min(cellFld[own[facei]], ownPyrVol[facei]);
|
||||
cellFld[nei[facei]] = min(cellFld[nei[facei]], neiPyrVol[facei]);
|
||||
cellFld[own[facei]] =
|
||||
Foam::min(cellFld[own[facei]], ownPyrVol[facei]);
|
||||
cellFld[nei[facei]] =
|
||||
Foam::min(cellFld[nei[facei]], neiPyrVol[facei]);
|
||||
}
|
||||
|
||||
// Patch faces
|
||||
@ -614,7 +623,8 @@ void Foam::writeFields
|
||||
forAll(fc, i)
|
||||
{
|
||||
const label meshFacei = fvp.patch().start();
|
||||
cellFld[fc[i]] = min(cellFld[fc[i]], ownPyrVol[meshFacei]);
|
||||
cellFld[fc[i]] =
|
||||
Foam::min(cellFld[fc[i]], ownPyrVol[meshFacei]);
|
||||
}
|
||||
}
|
||||
|
||||
@ -625,7 +635,7 @@ void Foam::writeFields
|
||||
if (writeFaceFields)
|
||||
{
|
||||
scalarField minFacePyrVol(neiPyrVol);
|
||||
minFacePyrVol = min
|
||||
minFacePyrVol = Foam::min
|
||||
(
|
||||
minFacePyrVol,
|
||||
SubField<scalar>(ownPyrVol, mesh.nInternalFaces())
|
||||
|
||||
@ -100,26 +100,26 @@ void printEdgeStats(const polyMesh& mesh)
|
||||
|
||||
if (mag(eVec & x) > 1-edgeTol)
|
||||
{
|
||||
minX = min(minX, eMag);
|
||||
maxX = max(maxX, eMag);
|
||||
minX = Foam::min(minX, eMag);
|
||||
maxX = Foam::max(maxX, eMag);
|
||||
nX++;
|
||||
}
|
||||
else if (mag(eVec & y) > 1-edgeTol)
|
||||
{
|
||||
minY = min(minY, eMag);
|
||||
maxY = max(maxY, eMag);
|
||||
minY = Foam::min(minY, eMag);
|
||||
maxY = Foam::max(maxY, eMag);
|
||||
nY++;
|
||||
}
|
||||
else if (mag(eVec & z) > 1-edgeTol)
|
||||
{
|
||||
minZ = min(minZ, eMag);
|
||||
maxZ = max(maxZ, eMag);
|
||||
minZ = Foam::min(minZ, eMag);
|
||||
maxZ = Foam::max(maxZ, eMag);
|
||||
nZ++;
|
||||
}
|
||||
else
|
||||
{
|
||||
minOther = min(minOther, eMag);
|
||||
maxOther = max(maxOther, eMag);
|
||||
minOther = Foam::min(minOther, eMag);
|
||||
maxOther = Foam::max(maxOther, eMag);
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
@ -465,7 +465,7 @@ labelList getRegionFaceOrder
|
||||
|
||||
// Do region interfaces
|
||||
{
|
||||
const label nRegions = max(cellToRegion)+1;
|
||||
const label nRegions = Foam::max(cellToRegion)+1;
|
||||
|
||||
// Sort in increasing region
|
||||
SortableList<label> sortKey(mesh.nInternalFaces(), labelMax);
|
||||
@ -478,8 +478,10 @@ labelList getRegionFaceOrder
|
||||
if (ownRegion != neiRegion)
|
||||
{
|
||||
sortKey[facei] =
|
||||
min(ownRegion, neiRegion)*nRegions
|
||||
+max(ownRegion, neiRegion);
|
||||
(
|
||||
Foam::min(ownRegion, neiRegion)*nRegions
|
||||
+ Foam::max(ownRegion, neiRegion)
|
||||
);
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
@ -320,10 +320,10 @@ bool doCommand
|
||||
const globalMeshData& parData = mesh.globalData();
|
||||
|
||||
label typSize =
|
||||
max
|
||||
Foam::max
|
||||
(
|
||||
parData.nTotalCells(),
|
||||
max
|
||||
Foam::max
|
||||
(
|
||||
parData.nTotalFaces(),
|
||||
parData.nTotalPoints()
|
||||
@ -375,7 +375,7 @@ bool doCommand
|
||||
|
||||
topoSet& currentSet = currentSetPtr();
|
||||
// Presize it according to current mesh data.
|
||||
currentSet.reserve(max(currentSet.size(), typSize));
|
||||
currentSet.reserve(Foam::max(currentSet.size(), typSize));
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
@ -269,11 +269,7 @@ void addToInterface
|
||||
EdgeMap<Map<label>>& regionsToSize
|
||||
)
|
||||
{
|
||||
edge interface
|
||||
(
|
||||
min(ownRegion, neiRegion),
|
||||
max(ownRegion, neiRegion)
|
||||
);
|
||||
const edge interface(ownRegion, neiRegion, true); // sort=true
|
||||
|
||||
auto iter = regionsToSize.find(interface);
|
||||
|
||||
@ -509,11 +505,7 @@ void getInterfaceSizes
|
||||
zoneID = mesh.faceZones().whichZone(facei);
|
||||
}
|
||||
|
||||
edge interface
|
||||
(
|
||||
min(ownRegion, neiRegion),
|
||||
max(ownRegion, neiRegion)
|
||||
);
|
||||
const edge interface(ownRegion, neiRegion, true); // sort=true
|
||||
|
||||
faceToInterface[facei] = regionsToInterface[interface][zoneID];
|
||||
}
|
||||
@ -532,11 +524,7 @@ void getInterfaceSizes
|
||||
zoneID = mesh.faceZones().whichZone(facei);
|
||||
}
|
||||
|
||||
edge interface
|
||||
(
|
||||
min(ownRegion, neiRegion),
|
||||
max(ownRegion, neiRegion)
|
||||
);
|
||||
const edge interface(ownRegion, neiRegion, true); // sort=true
|
||||
|
||||
faceToInterface[facei] = regionsToInterface[interface][zoneID];
|
||||
}
|
||||
|
||||
@ -295,7 +295,8 @@ void subsetTopoSets
|
||||
|
||||
Info<< "Subsetting " << set.type() << " " << set.name() << endl;
|
||||
|
||||
labelHashSet subset(2*min(set.size(), map.size()));
|
||||
labelHashSet subset;
|
||||
subset.reserve(Foam::min(set.size(), map.size()));
|
||||
|
||||
// Map the data
|
||||
forAll(map, i)
|
||||
|
||||
@ -1131,7 +1131,7 @@ int main(int argc, char *argv[])
|
||||
for
|
||||
(
|
||||
label addedI=next;
|
||||
addedI<min(proci+step, nProcs);
|
||||
addedI < Foam::min(proci+step, nProcs);
|
||||
addedI++
|
||||
)
|
||||
{
|
||||
|
||||
@ -403,11 +403,11 @@ void printMeshData(const polyMesh& mesh)
|
||||
<< nBndFaces-nProcFaces << endl;
|
||||
}
|
||||
|
||||
maxProcCells = max(maxProcCells, nLocalCells);
|
||||
maxProcCells = Foam::max(maxProcCells, nLocalCells);
|
||||
totProcFaces += nProcFaces;
|
||||
totProcPatches += nei.size();
|
||||
maxProcFaces = max(maxProcFaces, nProcFaces);
|
||||
maxProcPatches = max(maxProcPatches, nei.size());
|
||||
maxProcFaces = Foam::max(maxProcFaces, nProcFaces);
|
||||
maxProcPatches = Foam::max(maxProcPatches, nei.size());
|
||||
}
|
||||
|
||||
// Summary stats
|
||||
|
||||
@ -131,7 +131,7 @@ int main(int argc, char *argv[])
|
||||
args.readIfPresent("format", setFormat);
|
||||
|
||||
args.readIfPresent("stride", sampleFrequency);
|
||||
sampleFrequency = max(1, sampleFrequency); // sanity
|
||||
sampleFrequency = Foam::max(1, sampleFrequency); // sanity
|
||||
|
||||
// Setup the writer
|
||||
auto writerPtr =
|
||||
@ -179,7 +179,7 @@ int main(int argc, char *argv[])
|
||||
maxIds.resize(origProc+1, -1);
|
||||
}
|
||||
|
||||
maxIds[origProc] = max(maxIds[origProc], origId);
|
||||
maxIds[origProc] = Foam::max(maxIds[origProc], origId);
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
@ -1056,7 +1056,7 @@ void calc_drag_etc
|
||||
const scalar expon =
|
||||
(
|
||||
br > 0.0
|
||||
? min(max((surr_br / br - 0.25) * 4.0 / 3.0, scalar(0)), scalar(1))
|
||||
? Foam::clamp((surr_br / br - 0.25) * 4.0 / 3.0, Foam::zero_one{})
|
||||
: 0.0
|
||||
);
|
||||
|
||||
@ -1114,16 +1114,16 @@ void Foam::PDRarrays::blockageSummary() const
|
||||
totVolBlock += v_block(ijk) * pdrBlock.V(ijk);
|
||||
totArea += surf(ijk);
|
||||
|
||||
totCount += max(0, obs_count(ijk));
|
||||
totCount += Foam::max(0, obs_count(ijk));
|
||||
|
||||
totDrag.x() += max(0, drag_s(ijk).xx());
|
||||
totDrag.y() += max(0, drag_s(ijk).yy());
|
||||
totDrag.z() += max(0, drag_s(ijk).zz());
|
||||
totDrag.x() += Foam::max(0, drag_s(ijk).xx());
|
||||
totDrag.y() += Foam::max(0, drag_s(ijk).yy());
|
||||
totDrag.z() += Foam::max(0, drag_s(ijk).zz());
|
||||
|
||||
for (direction cmpt=0; cmpt < vector::nComponents; ++cmpt)
|
||||
{
|
||||
totBlock[cmpt] += max(0, area_block_s(ijk)[cmpt]);
|
||||
totBlock[cmpt] += max(0, area_block_r(ijk)[cmpt]);
|
||||
totBlock[cmpt] += Foam::max(0, area_block_s(ijk)[cmpt]);
|
||||
totBlock[cmpt] += Foam::max(0, area_block_r(ijk)[cmpt]);
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
@ -302,7 +302,7 @@ void Foam::PDRutils::circle_overlap
|
||||
scalar da = ac - 0.5 * (a1 + a2);
|
||||
scalar db = bc - 0.5 * (b1 + b2);
|
||||
scalar dc = std::hypot(da, db);
|
||||
scalar rat1 = min(max((dc / sqrt(area) - 0.3) * 1.4, 0), 1);
|
||||
scalar rat1 = Foam::min(Foam::max((dc / sqrt(area) - 0.3) * 1.4, 0), 1);
|
||||
scalar drg0 = c_drag(ia,ib).xx();
|
||||
scalar drg1 = c_drag(ia,ib).yy();
|
||||
scalar drg = std::hypot(drg0, drg1);
|
||||
@ -449,8 +449,8 @@ scalar block_overlap
|
||||
{
|
||||
PDRobstacle over;
|
||||
|
||||
over.pt = max(blk1.pt, blk2.pt);
|
||||
over.span = min(max1, max2) - over.pt;
|
||||
over.pt = Foam::max(blk1.pt, blk2.pt);
|
||||
over.span = Foam::min(max1, max2) - over.pt;
|
||||
|
||||
assert(cmptProduct(over.span) > 0.0);
|
||||
|
||||
@ -603,11 +603,11 @@ scalar block_cylinder_overlap
|
||||
|
||||
over.x() = a_centre - 0.5 * a_lblk;
|
||||
over.y() = b_centre - 0.5 * b_lblk;
|
||||
over.z() = max(blk1.z(), cyl2.z());
|
||||
over.z() = Foam::max(blk1.z(), cyl2.z());
|
||||
|
||||
over.span.x() = a_lblk;
|
||||
over.span.y() = b_lblk;
|
||||
over.span.z() = min(max1.z(), cyl2.z() + cyl2.len()) - over.z();
|
||||
over.span.z() = Foam::min(max1.z(), cyl2.z() + cyl2.len()) - over.z();
|
||||
assert(over.x() > -200.0);
|
||||
assert(over.x() < 2000.0);
|
||||
}
|
||||
@ -668,11 +668,11 @@ scalar block_cylinder_overlap
|
||||
|
||||
over.z() = a_centre - a_lblk * 0.5;
|
||||
over.x() = b_centre - b_lblk * 0.5;
|
||||
over.y() = max(blk1.y(), cyl2.y());
|
||||
over.y() = Foam::max(blk1.y(), cyl2.y());
|
||||
|
||||
over.span.z() = a_lblk;
|
||||
over.span.x() = b_lblk;
|
||||
over.span.y() = min(max1.y(), cyl2.y() + cyl2.len()) - over.y();
|
||||
over.span.y() = Foam::min(max1.y(), cyl2.y() + cyl2.len()) - over.y();
|
||||
}
|
||||
break;
|
||||
|
||||
@ -734,11 +734,11 @@ scalar block_cylinder_overlap
|
||||
|
||||
over.y() = a_centre - a_lblk * 0.5;
|
||||
over.z() = b_centre - b_lblk * 0.5;
|
||||
over.x() = max(blk1.x(), cyl2.x());
|
||||
over.x() = Foam::max(blk1.x(), cyl2.x());
|
||||
|
||||
over.span.y() = a_lblk;
|
||||
over.span.z() = b_lblk;
|
||||
over.span.x() = min(max1.x(), cyl2.x() + cyl2.len()) - over.x();
|
||||
over.span.x() = Foam::min(max1.x(), cyl2.x() + cyl2.len()) - over.x();
|
||||
}
|
||||
break;
|
||||
}
|
||||
|
||||
@ -71,9 +71,9 @@ int main(int argc, char *argv[])
|
||||
{
|
||||
scalar b = j1(swirlProfile*r/cylinderRadius).value();
|
||||
scalar vEff = omega*b;
|
||||
r = max(r, SMALL);
|
||||
r = Foam::max(r, SMALL);
|
||||
U[celli] = ((vEff/r)*(c & yT))*xT + (-(vEff/r)*(c & xT))*yT;
|
||||
Umax = max(Umax, mag(U[celli]));
|
||||
Umax = Foam::max(Umax, mag(U[celli]));
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
@ -887,7 +887,7 @@ int main(int argc, char *argv[])
|
||||
writeParts
|
||||
(
|
||||
surf,
|
||||
min(outputThreshold, numZones),
|
||||
Foam::min(outputThreshold, numZones),
|
||||
faceZone,
|
||||
surfFilePath,
|
||||
surfFileStem
|
||||
@ -953,7 +953,7 @@ int main(int argc, char *argv[])
|
||||
writeParts
|
||||
(
|
||||
surf,
|
||||
min(outputThreshold, numNormalZones),
|
||||
Foam::min(outputThreshold, numNormalZones),
|
||||
normalZone,
|
||||
surfFilePath,
|
||||
surfFileStem + "_normal"
|
||||
|
||||
@ -280,7 +280,7 @@ int main(int argc, char *argv[])
|
||||
const IOdictionary dict(dictIO);
|
||||
|
||||
const scalar dist(args.get<scalar>(1));
|
||||
const scalar matchTolerance(max(1e-6*dist, SMALL));
|
||||
const scalar matchTolerance(Foam::max(1e-6*dist, SMALL));
|
||||
const label maxIters = 100;
|
||||
|
||||
Info<< "Hooking distance = " << dist << endl;
|
||||
|
||||
@ -276,7 +276,7 @@ label detectIntersectionPoints
|
||||
|
||||
// 1. Extrusion offset vectors intersecting new surface location
|
||||
{
|
||||
scalar tol = max(tolerance, 10*s.tolerance());
|
||||
scalar tol = Foam::max(tolerance, 10*s.tolerance());
|
||||
|
||||
// Collect all the edge vectors. Slightly shorten the edges to prevent
|
||||
// finding lots of intersections. The fast triangle intersection routine
|
||||
@ -296,7 +296,7 @@ label detectIntersectionPoints
|
||||
&& !localFaces[hits[pointI].index()].found(pointI)
|
||||
)
|
||||
{
|
||||
scale[pointI] = max(0.0, scale[pointI]-0.2);
|
||||
scale[pointI] = Foam::max(0.0, scale[pointI]-0.2);
|
||||
|
||||
isPointOnHitEdge.set(pointI);
|
||||
nHits++;
|
||||
@ -330,7 +330,7 @@ label detectIntersectionPoints
|
||||
<< pt
|
||||
<< endl;
|
||||
|
||||
scale[e[0]] = max(0.0, scale[e[0]]-0.2);
|
||||
scale[e[0]] = Foam::max(0.0, scale[e[0]]-0.2);
|
||||
nHits++;
|
||||
}
|
||||
if (isPointOnHitEdge.set(e[1]))
|
||||
@ -342,7 +342,7 @@ label detectIntersectionPoints
|
||||
<< pt
|
||||
<< endl;
|
||||
|
||||
scale[e[1]] = max(0.0, scale[e[1]]-0.2);
|
||||
scale[e[1]] = Foam::max(0.0, scale[e[1]]-0.2);
|
||||
nHits++;
|
||||
}
|
||||
}
|
||||
@ -418,7 +418,7 @@ void minSmooth
|
||||
const edge& e = edges[edgeI];
|
||||
scalar w = mag(points[mp[e[0]]]-points[mp[e[1]]]);
|
||||
|
||||
edgeWeights[edgeI] = 1.0/(max(w, SMALL));
|
||||
edgeWeights[edgeI] = 1.0/(Foam::max(w, SMALL));
|
||||
}
|
||||
|
||||
tmp<scalarField> tavgFld = avg(s, fld, edgeWeights);
|
||||
@ -429,7 +429,7 @@ void minSmooth
|
||||
{
|
||||
if (isAffectedPoint.test(pointI))
|
||||
{
|
||||
newFld[pointI] = min
|
||||
newFld[pointI] = Foam::min
|
||||
(
|
||||
fld[pointI],
|
||||
0.5*fld[pointI] + 0.5*avgFld[pointI]
|
||||
|
||||
@ -109,7 +109,7 @@ void writeOBJ
|
||||
const auto& constraints = ppp.constraints();
|
||||
forAll(constraints, i)
|
||||
{
|
||||
maxConstraint = max(maxConstraint, constraints[i].first());
|
||||
maxConstraint = Foam::max(maxConstraint, constraints[i].first());
|
||||
}
|
||||
reduce(maxConstraint, maxOp<label>());
|
||||
|
||||
|
||||
@ -180,8 +180,8 @@ int main(int argc, char *argv[])
|
||||
|
||||
scalar o2 = (1.0/equiv)*stoicO2;
|
||||
scalar n2 = (0.79/0.21)*o2;
|
||||
scalar fres = max(1.0 - 1.0/equiv, 0.0);
|
||||
scalar ores = max(1.0/equiv - 1.0, 0.0);
|
||||
scalar fres = Foam::max(1.0 - 1.0/equiv, 0.0);
|
||||
scalar ores = Foam::max(1.0/equiv - 1.0, 0.0);
|
||||
scalar fburnt = 1.0 - fres;
|
||||
|
||||
thermo reactants
|
||||
|
||||
@ -196,12 +196,12 @@ int main(int argc, char *argv[])
|
||||
// Number of moles of species for one mole of fuel
|
||||
scalar o2 = (1.0/equiv)*stoicO2;
|
||||
scalar n2 = (0.79/0.21)*o2;
|
||||
scalar fres = max(1.0 - 1.0/equiv, 0.0);
|
||||
scalar fres = Foam::max(1.0 - 1.0/equiv, 0.0);
|
||||
scalar fburnt = 1.0 - fres;
|
||||
|
||||
// Initial guess for number of moles of product species
|
||||
// ignoring product dissociation
|
||||
scalar oresInit = max(1.0/equiv - 1.0, 0.0)*stoicO2;
|
||||
scalar oresInit = Foam::max(1.0/equiv - 1.0, 0.0)*stoicO2;
|
||||
scalar co2Init = fburnt*stoicCO2;
|
||||
scalar h2oInit = fburnt*stoicH2O;
|
||||
|
||||
@ -231,18 +231,18 @@ int main(int argc, char *argv[])
|
||||
if (j > 0)
|
||||
{
|
||||
co = co2*
|
||||
min
|
||||
Foam::min
|
||||
(
|
||||
CO2BreakUp.Kn(P, equilibriumFlameTemperature, N)
|
||||
/::sqrt(max(ores, 0.001)),
|
||||
/::sqrt(Foam::max(ores, 0.001)),
|
||||
1.0
|
||||
);
|
||||
|
||||
h2 = h2o*
|
||||
min
|
||||
Foam::min
|
||||
(
|
||||
H2OBreakUp.Kn(P, equilibriumFlameTemperature, N)
|
||||
/::sqrt(max(ores, 0.001)),
|
||||
/::sqrt(Foam::max(ores, 0.001)),
|
||||
1.0
|
||||
);
|
||||
|
||||
|
||||
@ -36,11 +36,12 @@ Description
|
||||
if (adjustTimeStep)
|
||||
{
|
||||
scalar maxDeltaTFact = maxCo/(CoNum + SMALL);
|
||||
scalar deltaTFact = min(min(maxDeltaTFact, 1.0 + 0.1*maxDeltaTFact), 1.2);
|
||||
scalar deltaTFact =
|
||||
Foam::min(Foam::min(maxDeltaTFact, 1.0 + 0.1*maxDeltaTFact), 1.2);
|
||||
|
||||
runTime.setDeltaT
|
||||
(
|
||||
min
|
||||
Foam::min
|
||||
(
|
||||
deltaTFact*runTime.deltaTValue(),
|
||||
maxDeltaT
|
||||
|
||||
@ -38,10 +38,10 @@ if (adjustTimeStep)
|
||||
{
|
||||
runTime.setDeltaT
|
||||
(
|
||||
min
|
||||
Foam::min
|
||||
(
|
||||
maxCo*runTime.deltaTValue()/CoNum,
|
||||
min(runTime.deltaTValue(), maxDeltaT)
|
||||
Foam::min(runTime.deltaTValue(), maxDeltaT)
|
||||
)
|
||||
);
|
||||
}
|
||||
|
||||
@ -163,7 +163,8 @@ scalar pyrolysisModelCollection::maxDiff() const
|
||||
scalar maxDiff = 0.0;
|
||||
forAll(*this, i)
|
||||
{
|
||||
maxDiff = max(maxDiff, this->operator[](i).maxDiff());
|
||||
maxDiff =
|
||||
Foam::max(maxDiff, this->operator[](i).maxDiff());
|
||||
}
|
||||
|
||||
return maxDiff;
|
||||
@ -175,7 +176,8 @@ scalar pyrolysisModelCollection::solidRegionDiffNo() const
|
||||
scalar totalDiNum = GREAT;
|
||||
forAll(*this, i)
|
||||
{
|
||||
totalDiNum = min(totalDiNum, this->operator[](i).solidRegionDiffNo());
|
||||
totalDiNum =
|
||||
Foam::min(totalDiNum, this->operator[](i).solidRegionDiffNo());
|
||||
}
|
||||
|
||||
return totalDiNum;
|
||||
|
||||
Reference in New Issue
Block a user