22 #include "libavutil/avassert.h"
23 #include "libavutil/channel_layout.h"
25 #define TEMPLATE_REMATRIX_FLT
27 #undef TEMPLATE_REMATRIX_FLT
29 #define TEMPLATE_REMATRIX_DBL
31 #undef TEMPLATE_REMATRIX_DBL
33 #define TEMPLATE_REMATRIX_S16
35 #undef TEMPLATE_REMATRIX_S16
39 #define FRONT_CENTER 2
40 #define LOW_FREQUENCY 3
43 #define FRONT_LEFT_OF_CENTER 6
44 #define FRONT_RIGHT_OF_CENTER 7
49 #define TOP_FRONT_LEFT 12
50 #define TOP_FRONT_CENTER 13
51 #define TOP_FRONT_RIGHT 14
52 #define TOP_BACK_LEFT 15
53 #define TOP_BACK_CENTER 16
54 #define TOP_BACK_RIGHT 17
58 int nb_in, nb_out,
in,
out;
65 for (out = 0; out < nb_out; out++) {
66 for (in = 0; in < nb_in; in++)
67 s->
matrix[out][in] = matrix[in];
76 if(layout&(layout-1))
return 1;
111 double matrix[64][64]={{0}};
112 int64_t unaccounted, in_ch_layout, out_ch_layout;
139 if(in_ch_layout & out_ch_layout & (1ULL<<i))
143 unaccounted= in_ch_layout & ~out_ch_layout;
151 if(in_ch_layout & AV_CH_LAYOUT_STEREO) {
162 if(out_ch_layout & AV_CH_FRONT_CENTER){
165 if(in_ch_layout & AV_CH_FRONT_CENTER)
181 if (unaccounted & (AV_CH_BACK_LEFT | AV_CH_SIDE_LEFT)) {
192 }
else if(out_ch_layout & AV_CH_FRONT_CENTER){
198 if(out_ch_layout & AV_CH_BACK_CENTER){
202 if(in_ch_layout & AV_CH_SIDE_LEFT){
224 }
else if(out_ch_layout & AV_CH_FRONT_CENTER){
232 if(out_ch_layout & AV_CH_BACK_LEFT){
235 if (in_ch_layout & AV_CH_BACK_LEFT) {
242 }
else if(out_ch_layout & AV_CH_BACK_CENTER){
260 }
else if(out_ch_layout & AV_CH_FRONT_CENTER){
271 }
else if(out_ch_layout & AV_CH_FRONT_CENTER){
279 if (out_ch_layout & AV_CH_FRONT_CENTER) {
288 for(out_i=i=0; i<64; i++){
292 s->
matrix[out_i][in_i]= matrix[i][j];
294 sum += fabs(matrix[i][j]);
296 if(in_ch_layout & (1ULL<<j))
299 maxcoef=
FFMAX(maxcoef, sum);
300 if(out_ch_layout & (1ULL<<i))
310 s->
matrix[i][j] /= maxcoef;
345 for (i = 0; i < nb_out; i++)
346 for (j = 0; j < nb_in; j++)
355 for (i = 0; i < nb_out; i++)
356 for (j = 0; j < nb_in; j++)
365 for (i = 0; i < nb_out; i++)
366 for (j = 0; j < nb_in; j++)
397 int out_i, in_i, i, j;
408 off = len1 * out->
bps;
414 for(out_i=0; out_i<out->
ch_count; out_i++){
422 if(s->
matrix[out_i][in_i]!=1.0){
428 memcpy(out->
ch[out_i], in->
ch[in_i], len*out->
bps);
430 out->
ch[out_i]= in->
ch[in_i];
445 for(i=0; i<
len; i++){
449 v+= ((
float*)in->
ch[in_i])[i] * s->
matrix[out_i][in_i];
451 ((
float*)out->
ch[out_i])[i]=
v;
454 for(i=0; i<
len; i++){
458 v+= ((
double*)in->
ch[in_i])[i] * s->
matrix[out_i][in_i];
460 ((
double*)out->
ch[out_i])[i]=
v;
463 for(i=0; i<
len; i++){
467 v+= ((int16_t*)in->
ch[in_i])[i] * s->
matrix32[out_i][in_i];
469 ((int16_t*)out->
ch[out_i])[i]= (v + 16384)>>15;