diff options
| author | John Glover <john@john-glovers-macbook.local> | 2010-12-24 19:08:42 +0000 |
|---|---|---|
| committer | John Glover <john@john-glovers-macbook.local> | 2010-12-24 19:08:42 +0000 |
| commit | 5daa78a0630c865519e557fde5aed57dcdee4aea (patch) | |
| tree | a1726e2f87d7eb10573350fa3aeebea4eb024536 | |
| parent | 5b884d40c155fb47ed5ebf7dbbea2cff46f0a1c4 (diff) | |
| download | simpl-5daa78a0630c865519e557fde5aed57dcdee4aea.tar.gz simpl-5daa78a0630c865519e557fde5aed57dcdee4aea.tar.bz2 simpl-5daa78a0630c865519e557fde5aed57dcdee4aea.zip | |
Fixed a bug in converting SMS peak amplitudes from DB to linear
| -rw-r--r-- | sms/analysis.c | 7 | ||||
| -rw-r--r-- | tests/sms.py | 2 |
2 files changed, 5 insertions, 4 deletions
diff --git a/sms/analysis.c b/sms/analysis.c index 345b95f..efd98f4 100644 --- a/sms/analysis.c +++ b/sms/analysis.c @@ -197,7 +197,7 @@ int sms_findPeaks(int sizeWaveform, sfloat *pWaveform, SMS_AnalParams *pAnalPara if(i < pSpectralPeaks->nPeaksFound) { pSpectralPeaks->pSpectralPeaks[i].fMag = - pow(10.0, 0.05*(pAnalParams->ppFrames[iCurrentFrame]->pSpectralPeaks[i].fMag)); + sms_dBToMag(pAnalParams->ppFrames[iCurrentFrame]->pSpectralPeaks[i].fMag); pSpectralPeaks->pSpectralPeaks[i].fFreq = pAnalParams->ppFrames[iCurrentFrame]->pSpectralPeaks[i].fFreq; pSpectralPeaks->pSpectralPeaks[i].fPhase = @@ -210,6 +210,7 @@ int sms_findPeaks(int sizeWaveform, sfloat *pWaveform, SMS_AnalParams *pAnalPara pSpectralPeaks->pSpectralPeaks[i].fPhase = 0.0; } } + /*printf("\n");*/ return pSpectralPeaks->nPeaks; } else @@ -243,7 +244,7 @@ void sms_setPeaks(SMS_AnalParams *pAnalParams, int numamps, sfloat* amps, for(i = 0; i < numamps; i++) { /* copy current peaks data */ - frame->pSpectralPeaks[i].fMag = 20.0 * log10(amps[i]); + frame->pSpectralPeaks[i].fMag = sms_magToDB(amps[i]); frame->pSpectralPeaks[i].fFreq = freqs[i]; frame->pSpectralPeaks[i].fPhase = phases[i]; } @@ -515,7 +516,7 @@ int sms_analyze(int sizeWaveform, sfloat *pWaveform, SMS_Data *pSmsData, SMS_Ana for(i = 0; i < numTracks; i++) { pSmsData->pFSinFreq[i] = pAnalParams->ppFrames[iCurrentFrame]->pSpectralPeaks[i].fFreq; - pSmsData->pFSinAmp[i] = pAnalParams->ppFrames[iCurrentFrame]->pSpectralPeaks[i].fMag; + pSmsData->pFSinAmp[i] = sms_dBToMag(pAnalParams->ppFrames[iCurrentFrame]->pSpectralPeaks[i].fMag); if(pAnalParams->iFormat == SMS_FORMAT_HP || pAnalParams->iFormat == SMS_FORMAT_IHP) { diff --git a/tests/sms.py b/tests/sms.py index a1b8fa9..a5fb702 100644 --- a/tests/sms.py +++ b/tests/sms.py @@ -502,7 +502,7 @@ class TestSimplSMS(object): if simplsms_amps[i]: p = simpl.Peak() # convert amplitude back to linear - p.amplitude = 10**(simplsms_amps[i]/20.0) + p.amplitude = simplsms_amps[i] p.frequency = simplsms_freqs[i] p.phase = simplsms_phases[i] frame_peaks.append(p) |