aboutsummaryrefslogtreecommitdiffstats
path: root/src/atrac/atrac1_scale.cpp
blob: 461da3a701ef7cde1bd7ad1cb4a13b08e6dbe626 (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
49
50
51
52
53
54
55
#include "atrac1_scale.h"
#include <cmath>
#include <iostream>
#include <algorithm>
namespace NAtrac1 {
using std::vector;
using std::map;

using namespace std;
map<uint32_t, uint8_t> TScaler::ScaleIndex;
static const uint32_t MAX_SCALE = 65536;

static bool absComp(double a, double b) {
    return abs(a) < abs(b);
}

TScaler::TScaler() {
    if (ScaleIndex.empty()) {
        for (int i = 0; i < 64; i++) {
            ScaleIndex[ScaleTable[i] * 256] = i;
        }
    }
}

vector<TScaledBlock> TScaler::Scale(const vector<double>& specs) {
    vector<TScaledBlock> scaledBlocks;
    for (uint8_t bandNum = 0; bandNum < QMF_BANDS; ++bandNum) {
        for (uint8_t blockNum = BlocksPerBand[bandNum]; blockNum < BlocksPerBand[bandNum + 1]; ++blockNum) {
            const uint16_t specNumStart = SpecsStartLong[blockNum]; 
            const uint16_t specNumEnd = SpecsStartLong[blockNum] + 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<uint32_t, uint8_t>::const_iterator scaleIter = ScaleIndex.lower_bound(maxAbsSpec * 256);
            const double scaleFactor = scaleIter->first / 256.0;
            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;
                scaledBlocks.back().Values.push_back(scaledValue);
			}
		}
	}
    return scaledBlocks;
}
}