summaryrefslogtreecommitdiff
path: root/src/libac3/rematrix.c
blob: 7098f0526c20f49357b1d9111b0e54e81f1d3f69 (plain)
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
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
/* 
 *    rematrix.c
 *
 *	Copyright (C) Aaron Holtzman - July 1999
 *
 *  This file is part of ac3dec, a free Dolby AC-3 stream decoder.
 *	
 *  ac3dec is free software; you can redistribute it and/or modify
 *  it under the terms of the GNU General Public License as published by
 *  the Free Software Foundation; either version 2, or (at your option)
 *  any later version.
 *   
 *  ac3dec is distributed in the hope that it will be useful,
 *  but WITHOUT ANY WARRANTY; without even the implied warranty of
 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 *  GNU General Public License for more details.
 *   
 *  You should have received a copy of the GNU General Public License
 *  along with GNU Make; see the file COPYING.  If not, write to
 *  the Free Software Foundation, 675 Mass Ave, Cambridge, MA 02139, USA. 
 *
 *
 */

#include <stdlib.h>
#include <stdio.h>
#include "ac3.h"
#include "ac3_internal.h"


#include "decode.h"
#include "rematrix.h"

struct rematrix_band_s
{
    uint32_t start;
    uint32_t end;
};

struct rematrix_band_s rematrix_band[] = { {13,24}, {25,36}, {37 ,60}, {61,252}};

static inline uint32_t min(uint32_t a,uint32_t b);

static inline uint32_t
min(uint32_t a,uint32_t b)
{
    return (a < b ? a : b);
}

/* This routine simply does stereo rematixing for the 2 channel 
 * stereo mode */
void rematrix(audblk_t *audblk, stream_samples_t samples)
{
    uint32_t num_bands;
    uint32_t start;
    uint32_t end;
    uint32_t i,j;
    float left,right;

    if(!audblk->cplinu || audblk->cplbegf > 2)
	num_bands = 4;
    else if (audblk->cplbegf > 0)
	num_bands = 3;
    else
	num_bands = 2;

    for(i=0;i < num_bands; i++) {
	if(!audblk->rematflg[i])
	    continue;

	start = rematrix_band[i].start;
	end = min(rematrix_band[i].end ,12 * audblk->cplbegf + 36);
	
	for(j=start;j < end; j++) {
	    left  = samples[0][j] + samples[1][j];
	    right = samples[0][j] - samples[1][j];
	    samples[0][j] = left;
	    samples[1][j] = right;
	}
    }
}