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 - report missing octave bands

Browse Source
This commit is contained in:
Andrew Heather
2021-06-10 13:25:28 +01:00
parent 98abafccdb
commit ee4f7feb30
1 changed files with 14 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

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

@ -78,6 +78,7 @@ void Foam::noiseModel::setOctaveBands
// Centre frequencies // Centre frequencies
DynamicList<scalar> fc; DynamicList<scalar> fc;
DynamicList<scalar> missedBins;
// Convert to lower band limit // Convert to lower band limit
fTest /= fRatioL2C; fTest /= fRatioL2C;
@ -91,9 +92,12 @@ void Foam::noiseModel::setOctaveBands
if (f[i] >= fTest) if (f[i] >= fTest)
{ {
// Advance band if appropriate // Advance band if appropriate
label stepi = 0;
while (f[i] > fTest) while (f[i] > fTest)
{ {
if (stepi) missedBins.append(fTest/fRatio*fRatioL2C);
fTest *= fRatio; fTest *= fRatio;
++stepi;
} }
fTest /= fRatio; fTest /= fRatio;
@ -113,6 +117,16 @@ void Foam::noiseModel::setOctaveBands
fBandIDs = bandIDs.sortedToc(); fBandIDs = bandIDs.sortedToc();
if (missedBins.size())
{
label nMiss = missedBins.size();
label nTotal = nMiss + fc.size() - 1;
WarningInFunction
<< "Empty bands found: " << nMiss << " of " << nTotal
<< " with centre frequencies " << flatOutput(missedBins)
<< endl;
}
if (fc.size()) if (fc.size())
{ {
// Remove the last centre frequency (beyond upper frequency limit) // Remove the last centre frequency (beyond upper frequency limit)