kcat/openal-soft
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Fix the UHJ decoder parameters in uhjdecoder

Browse Source
This commit is contained in:
Chris Robinson 2022-03-02 21:19:53 -08:00
parent 63ac9c0b5d
commit 240171a53d
1 changed files with 5 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
utils/uhjdecoder.cpp
View File

@ -244,19 +244,19 @@ void UhjDecoder::decode(const float *RESTRICT InSamples, const size_t InChannels
mQ[sFilterDelay+i] = InSamples[i*InChannels + 3]; mQ[sFilterDelay+i] = InSamples[i*InChannels + 3];
} }
/* Precompute j(0.828347*D + 0.767835*T) and store in xoutput. */ /* Precompute j(0.828331*D + 0.767820*T) and store in xoutput. */
auto tmpiter = std::copy(mDTHistory.cbegin(), mDTHistory.cend(), mTemp.begin()); auto tmpiter = std::copy(mDTHistory.cbegin(), mDTHistory.cend(), mTemp.begin());
std::transform(mD.cbegin(), mD.cbegin()+SamplesToDo+sFilterDelay, mT.cbegin(), tmpiter, std::transform(mD.cbegin(), mD.cbegin()+SamplesToDo+sFilterDelay, mT.cbegin(), tmpiter,
[](const float d, const float t) noexcept { return 0.828347f*d + 0.767835f*t; }); [](const float d, const float t) noexcept { return 0.828331f*d + 0.767820f*t; });
std::copy_n(mTemp.cbegin()+SamplesToDo, mDTHistory.size(), mDTHistory.begin()); std::copy_n(mTemp.cbegin()+SamplesToDo, mDTHistory.size(), mDTHistory.begin());
PShift.process({xoutput, SamplesToDo}, mTemp.data()); PShift.process({xoutput, SamplesToDo}, mTemp.data());
for(size_t i{0};i < SamplesToDo;++i) for(size_t i{0};i < SamplesToDo;++i)
{ {
/* W = 0.981532*S + 0.197484*j(0.828347*D + 0.767835*T) */ /* W = 0.981532*S + 0.197484*j(0.828331*D + 0.767820*T) */
woutput[i] = 0.981532f*mS[i] + 0.197484f*xoutput[i]; woutput[i] = 0.981532f*mS[i] + 0.197484f*xoutput[i];
/* X = 0.418504*S - j(0.828347*D + 0.767835*T) */ /* X = 0.418496*S - j(0.828331*D + 0.767820*T) */
xoutput[i] = 0.418504f*mS[i] - xoutput[i]; xoutput[i] = 0.418496f*mS[i] - xoutput[i];
} }
/* Precompute j*S and store in youtput. */ /* Precompute j*S and store in youtput. */