zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

STYLE: avoid floating point comparison with uint8_t in DSMCCloud

Browse Source
This commit is contained in:
Mark Olesen
2019-01-23 22:52:35 +01:00
parent 5bdb6ddf9f
commit c5fd858c14
1 changed files with 19 additions and 22 deletions
Download Patch File Download Diff File Expand all files Collapse all files

41
src/lagrangian/DSMC/clouds/Templates/DSMCCloud/DSMCCloud.C
View File

@ -997,35 +997,32 @@ Foam::scalar Foam::DSMCCloud<ParcelType>::equipartitionInternalEnergy
direction iDof direction iDof
) )
{ {
scalar Ei = 0.0; if (iDof == 0)
if (iDof < SMALL)
{ {
return Ei; return 0;
} }
else if (iDof < 2.0 + SMALL && iDof > 2.0 - SMALL) else if (iDof == 2)
{ {
// Special case for iDof = 2, i.e. diatomics; // Special case for iDof = 2, i.e. diatomics;
Ei = return
-log(rndGen_.sample01<scalar>()) (
*physicoChemical::k.value()*temperature; -log(rndGen_.sample01<scalar>())
*physicoChemical::k.value()*temperature
);
} }
else
const scalar a = 0.5*iDof - 1;
scalar energyRatio = 0;
scalar P = -1;
do
{ {
scalar a = 0.5*iDof - 1; energyRatio = 10*rndGen_.sample01<scalar>();
scalar energyRatio; P = pow((energyRatio/a), a)*exp(a - energyRatio);
scalar P = -1; } while (P < rndGen_.sample01<scalar>());
do return energyRatio*physicoChemical::k.value()*temperature;
{
energyRatio = 10*rndGen_.sample01<scalar>();
P = pow((energyRatio/a), a)*exp(a - energyRatio);
} while (P < rndGen_.sample01<scalar>());
Ei = energyRatio*physicoChemical::k.value()*temperature;
}
return Ei;
} }