2 * Copyright (c) 2012 Mans Rullgard <mans@mansr.com>
4 * This file is part of FFmpeg.
6 * FFmpeg is free software; you can redistribute it and/or
7 * modify it under the terms of the GNU Lesser General Public
8 * License as published by the Free Software Foundation; either
9 * version 2.1 of the License, or (at your option) any later version.
11 * FFmpeg is distributed in the hope that it will be useful,
12 * but WITHOUT ANY WARRANTY; without even the implied warranty of
13 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
14 * Lesser General Public License for more details.
16 * You should have received a copy of the GNU Lesser General Public
17 * License along with FFmpeg; if not, write to the Free Software
18 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
21 #include "libavutil/attributes.h"
22 #include "libavutil/samplefmt.h"
26 #define SAMPLE_SIZE 16
28 #include "flacdsp_template.c"
29 #include "flacdsp_lpc_template.c"
33 #include "flacdsp_template.c"
37 #define SAMPLE_SIZE 32
39 #include "flacdsp_template.c"
40 #include "flacdsp_lpc_template.c"
44 #include "flacdsp_template.c"
46 static void flac_lpc_16_c(int32_t *decoded
, const int coeffs
[32],
47 int pred_order
, int qlevel
, int len
)
51 for (i
= pred_order
; i
< len
- 1; i
+= 2, decoded
+= 2) {
55 for (j
= 1; j
< pred_order
; j
++) {
62 d
= decoded
[j
] += s0
>> qlevel
;
64 decoded
[j
+ 1] += s1
>> qlevel
;
68 for (j
= 0; j
< pred_order
; j
++)
69 sum
+= coeffs
[j
] * decoded
[j
];
70 decoded
[j
] += sum
>> qlevel
;
74 static void flac_lpc_32_c(int32_t *decoded
, const int coeffs
[32],
75 int pred_order
, int qlevel
, int len
)
79 for (i
= pred_order
; i
< len
; i
++, decoded
++) {
81 for (j
= 0; j
< pred_order
; j
++)
82 sum
+= (int64_t)coeffs
[j
] * decoded
[j
];
83 decoded
[j
] += sum
>> qlevel
;
88 av_cold
void ff_flacdsp_init(FLACDSPContext
*c
, enum AVSampleFormat fmt
, int channels
,
92 c
->lpc
= flac_lpc_32_c
;
93 c
->lpc_encode
= flac_lpc_encode_c_32
;
95 c
->lpc
= flac_lpc_16_c
;
96 c
->lpc_encode
= flac_lpc_encode_c_16
;
100 case AV_SAMPLE_FMT_S32
:
101 c
->decorrelate
[0] = flac_decorrelate_indep_c_32
;
102 c
->decorrelate
[1] = flac_decorrelate_ls_c_32
;
103 c
->decorrelate
[2] = flac_decorrelate_rs_c_32
;
104 c
->decorrelate
[3] = flac_decorrelate_ms_c_32
;
107 case AV_SAMPLE_FMT_S32P
:
108 c
->decorrelate
[0] = flac_decorrelate_indep_c_32p
;
109 c
->decorrelate
[1] = flac_decorrelate_ls_c_32p
;
110 c
->decorrelate
[2] = flac_decorrelate_rs_c_32p
;
111 c
->decorrelate
[3] = flac_decorrelate_ms_c_32p
;
114 case AV_SAMPLE_FMT_S16
:
115 c
->decorrelate
[0] = flac_decorrelate_indep_c_16
;
116 c
->decorrelate
[1] = flac_decorrelate_ls_c_16
;
117 c
->decorrelate
[2] = flac_decorrelate_rs_c_16
;
118 c
->decorrelate
[3] = flac_decorrelate_ms_c_16
;
121 case AV_SAMPLE_FMT_S16P
:
122 c
->decorrelate
[0] = flac_decorrelate_indep_c_16p
;
123 c
->decorrelate
[1] = flac_decorrelate_ls_c_16p
;
124 c
->decorrelate
[2] = flac_decorrelate_rs_c_16p
;
125 c
->decorrelate
[3] = flac_decorrelate_ms_c_16p
;
130 ff_flacdsp_init_arm(c
, fmt
, channels
, bps
);
132 ff_flacdsp_init_x86(c
, fmt
, channels
, bps
);