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

ENH: noiseModel - protected log10 calls

Browse Source
This commit is contained in:
Andrew Heather
2021-08-05 13:03:03 +01:00
parent d8eecc5e2a
commit b78dccdd51
1 changed files with 23 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

26
src/randomProcesses/noise/noiseModels/noiseModel/noiseModel.C
View File

@ -137,6 +137,26 @@ void Foam::noiseModel::setOctaveBands
}
namespace Foam
{
tmp<scalarField> safeLog10(const scalarField& fld)
{
auto tresult = tmp<scalarField>::New(fld.size(), Zero);
auto& result = tresult.ref();
forAll(result, i)
{
if (fld[i] > 0)
{
result[i] = log10(fld[i]);
}
}
return tresult;
}
}
// * * * * * * * * * * * * * Protected Member Functions * * * * * * * * * * //
void Foam::noiseModel::readWriteOption
@ -688,7 +708,7 @@ Foam::tmp<Foam::scalarField> Foam::noiseModel::PSD
const scalarField& PSDf
) const
{
return 10*log10(PSDf/sqr(2e-5));
return 10*safeLog10(PSDf/sqr(2e-5));
}
@ -698,7 +718,7 @@ Foam::tmp<Foam::scalarField> Foam::noiseModel::SPL
const scalar f
) const
{
tmp<scalarField> tspl(10*log10(Prms2/sqr(2e-5)));
tmp<scalarField> tspl(10*safeLog10(Prms2/sqr(2e-5)));
scalarField& spl = tspl.ref();
switch (SPLweighting_)
@ -745,7 +765,7 @@ Foam::tmp<Foam::scalarField> Foam::noiseModel::SPL
const scalarField& f
) const
{
tmp<scalarField> tspl(10*log10(Prms2/sqr(2e-5)));
tmp<scalarField> tspl(10*safeLog10(Prms2/sqr(2e-5)));
scalarField& spl = tspl.ref();
switch (SPLweighting_)