blob: 1e24cb5d3fe56ff667fe27ce7af3b199ed51e7c9 (
plain) (
blame)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
|
#include "atrac_scale.h"
#include "atrac1.h"
#include <cmath>
#include <iostream>
#include <algorithm>
namespace NAtracDEnc {
using std::vector;
using std::map;
using namespace std;
static const uint32_t MAX_SCALE = 65536;
template<class TBaseData>
vector<TScaledBlock> TScaler<TBaseData>::Scale(const vector<double>& specs, const TBlockSize& blockSize) {
vector<TScaledBlock> scaledBlocks;
for (uint8_t bandNum = 0; bandNum < this->NumQMF; ++bandNum) {
const bool shortWinMode = !!blockSize.LogCount[bandNum];
for (uint8_t blockNum = this->BlocksPerBand[bandNum]; blockNum < this->BlocksPerBand[bandNum + 1]; ++blockNum) {
const uint16_t specNumStart = shortWinMode ? this->SpecsStartShort[blockNum] : this->SpecsStartLong[blockNum];
const uint16_t specNumEnd = specNumStart + this->SpecsPerBlock[blockNum];
double maxAbsSpec = 0;
for (uint16_t curSpec = specNumStart; curSpec < specNumEnd; ++curSpec) {
const double absSpec = abs(specs[curSpec]);
if (absSpec > maxAbsSpec) {
if (absSpec > MAX_SCALE) {
cerr << "got " << absSpec << " value - overflow" << endl;
maxAbsSpec = MAX_SCALE;
} else {
maxAbsSpec = absSpec;
}
}
}
const map<double, uint8_t>::const_iterator scaleIter = ScaleIndex.lower_bound(maxAbsSpec);
const double scaleFactor = scaleIter->first;
const uint8_t scaleFactorIndex = scaleIter->second;
scaledBlocks.push_back(TScaledBlock(scaleFactorIndex));
for (uint16_t specNum = specNumStart; specNum < specNumEnd; ++specNum) {
const double scaledValue = specs[specNum] / scaleFactor;
if (scaledValue > 1.0)
cerr << "got "<< scaledValue << " value - wrong scalling" << endl;
scaledBlocks.back().Values.push_back(scaledValue);
}
}
}
return scaledBlocks;
}
template class TScaler<TAtrac1Data>;
}
|