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

Bugfix: Avoid zero in log function in dsmcInitialise Authors: Robert Keser & Vanja Skuric

Browse source
This commit is contained in:
Robert Keser 2017-04-19 14:21:43 +02:00
parent 8607f3d224
commit d69eccdc69
1 changed files with 1 additions and 1 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

2
src/lagrangian/dsmc/clouds/Templates/DsmcCloud/DsmcCloudTemplate.C
View file

@ -1144,7 +1144,7 @@ Foam::scalar Foam::DsmcCloud<ParcelType>::equipartitionInternalEnergy
else if (iDof < 2.0 + SMALL && iDof > 2.0 - SMALL)
{
// Special case for iDof = 2, i.e. diatomics;
Ei = -log(rndGen_.scalar01())*kb*temperature;
Ei = -log(max(rndGen_.scalar01(), VSMALL))*kb*temperature;
}
else
{