28 #include "libavutil/log.h"
29 #include "libavutil/avassert.h"
63 static const double inv[100]={
64 1.0/( 1* 1), 1.0/( 2* 2), 1.0/( 3* 3), 1.0/( 4* 4), 1.0/( 5* 5), 1.0/( 6* 6), 1.0/( 7* 7), 1.0/( 8* 8), 1.0/( 9* 9), 1.0/(10*10),
65 1.0/(11*11), 1.0/(12*12), 1.0/(13*13), 1.0/(14*14), 1.0/(15*15), 1.0/(16*16), 1.0/(17*17), 1.0/(18*18), 1.0/(19*19), 1.0/(20*20),
66 1.0/(21*21), 1.0/(22*22), 1.0/(23*23), 1.0/(24*24), 1.0/(25*25), 1.0/(26*26), 1.0/(27*27), 1.0/(28*28), 1.0/(29*29), 1.0/(30*30),
67 1.0/(31*31), 1.0/(32*32), 1.0/(33*33), 1.0/(34*34), 1.0/(35*35), 1.0/(36*36), 1.0/(37*37), 1.0/(38*38), 1.0/(39*39), 1.0/(40*40),
68 1.0/(41*41), 1.0/(42*42), 1.0/(43*43), 1.0/(44*44), 1.0/(45*45), 1.0/(46*46), 1.0/(47*47), 1.0/(48*48), 1.0/(49*49), 1.0/(50*50),
69 1.0/(51*51), 1.0/(52*52), 1.0/(53*53), 1.0/(54*54), 1.0/(55*55), 1.0/(56*56), 1.0/(57*57), 1.0/(58*58), 1.0/(59*59), 1.0/(60*60),
70 1.0/(61*61), 1.0/(62*62), 1.0/(63*63), 1.0/(64*64), 1.0/(65*65), 1.0/(66*66), 1.0/(67*67), 1.0/(68*68), 1.0/(69*69), 1.0/(70*70),
71 1.0/(71*71), 1.0/(72*72), 1.0/(73*73), 1.0/(74*74), 1.0/(75*75), 1.0/(76*76), 1.0/(77*77), 1.0/(78*78), 1.0/(79*79), 1.0/(80*80),
72 1.0/(81*81), 1.0/(82*82), 1.0/(83*83), 1.0/(84*84), 1.0/(85*85), 1.0/(86*86), 1.0/(87*87), 1.0/(88*88), 1.0/(89*89), 1.0/(90*90),
73 1.0/(91*91), 1.0/(92*92), 1.0/(93*93), 1.0/(94*94), 1.0/(95*95), 1.0/(96*96), 1.0/(97*97), 1.0/(98*98), 1.0/(99*99), 1.0/(10000)
77 for(i=0; v != lastv; i++){
95 int filter_type,
int kaiser_beta){
99 const int center= (tap_count-1)/2;
108 for(ph=0;ph<phase_count;ph++) {
110 for(i=0;i<tap_count;i++) {
111 x =
M_PI * ((double)(i - center) - (double)ph / phase_count) *
factor;
117 x = fabs(((
double)(i - center) - (
double)ph / phase_count) * factor);
118 if(x<1.0) y= 1 - 3*x*x + 2*x*x*x + d*( -x*x + x*x*x);
119 else y= d*(-4 + 8*x - 5*x*x + x*x*x);
122 w = 2.0*x / (factor*tap_count) +
M_PI;
123 y *= 0.3635819 - 0.4891775 * cos(w) + 0.1365995 * cos(2*w) - 0.0106411 * cos(3*w);
126 w = 2.0*x / (factor*tap_count*
M_PI);
140 for(i=0;i<tap_count;i++)
141 ((int16_t*)filter)[ph * alloc + i] = av_clip(
lrintf(tab[i] * scale / norm), INT16_MIN, INT16_MAX);
144 for(i=0;i<tap_count;i++)
145 ((
int32_t*)filter)[ph * alloc + i] = av_clipl_int32(
llrint(tab[i] * scale / norm));
148 for(i=0;i<tap_count;i++)
149 ((
float*)filter)[ph * alloc + i] = tab[i] * scale / norm;
152 for(i=0;i<tap_count;i++)
153 ((
double*)filter)[ph * alloc + i] = tab[i] * scale / norm;
161 double sine[
LEN + tap_count];
162 double filtered[
LEN];
163 double maxff=-2, minff=2, maxsf=-2, minsf=2;
164 for(i=0; i<
LEN; i++){
165 double ss=0, sf=0, ff=0;
166 for(j=0; j<LEN+tap_count; j++)
167 sine[j]= cos(i*j*
M_PI/LEN);
168 for(j=0; j<
LEN; j++){
171 for(k=0; k<tap_count; k++)
172 sum += filter[ph * tap_count + k] * sine[k+j];
173 filtered[j]= sum / (1<<FILTER_SHIFT);
174 ss+= sine[j + center] * sine[j + center];
175 ff+= filtered[j] * filtered[j];
176 sf+= sine[j + center] * filtered[j];
181 maxff=
FFMAX(maxff, ff);
182 minff=
FFMIN(minff, ff);
183 maxsf=
FFMAX(maxsf, sf);
184 minsf=
FFMIN(minsf, sf);
186 av_log(NULL,
AV_LOG_ERROR,
"i:%4d ss:%f ff:%13.6e-%13.6e sf:%13.6e-%13.6e\n", i, ss, maxff, minff, maxsf, minsf);
200 double precision,
int cheby){
201 double cutoff = cutoff0? cutoff0 : 0.97;
202 double factor=
FFMIN(out_rate * cutoff / in_rate, 1.0);
203 int phase_count= 1<<phase_shift;
273 if (compensation_distance)
280 #define TEMPLATE_RESAMPLE_S16
282 #undef TEMPLATE_RESAMPLE_S16
284 #define TEMPLATE_RESAMPLE_S32
286 #undef TEMPLATE_RESAMPLE_S32
288 #define TEMPLATE_RESAMPLE_FLT
290 #undef TEMPLATE_RESAMPLE_FLT
292 #define TEMPLATE_RESAMPLE_DBL
294 #undef TEMPLATE_RESAMPLE_DBL
297 #if HAVE_MMXEXT_INLINE
301 #define TEMPLATE_RESAMPLE_S16_MMX2
303 #undef TEMPLATE_RESAMPLE_S16_MMX2
305 #if HAVE_SSSE3_INLINE
306 #define TEMPLATE_RESAMPLE_S16_SSSE3
308 #undef TEMPLATE_RESAMPLE_S16_SSSE3
311 #endif // HAVE_MMXEXT_INLINE
319 #if HAVE_MMXEXT_INLINE
320 #if HAVE_SSSE3_INLINE
static int multiple_resample(ResampleContext *c, AudioData *dst, int dst_size, AudioData *src, int src_size, int *consumed)
void * av_calloc(size_t nmemb, size_t size) av_malloc_attrib
Allocate a block of nmemb * size bytes with alignment suitable for all memory accesses (including vec...
Audio buffer used for intermediate storage between conversion phases.
int ch_count
number of channels
int swri_resample_float(struct ResampleContext *c, float *dst, const float *src, int *consumed, int src_size, int dst_size, int update_ctx)
SwrFilterType
Resampling Filter Types.
void av_log(void *avcl, int level, const char *fmt,...) av_printf_format(3
Send the specified message to the log if the level is less than or equal to the current av_log_level...
int in_buffer_index
cached buffer position
AudioData in_buffer
cached audio data (convert and resample purpose)
struct ResampleContext * resample
resampling context
void av_freep(void *ptr)
Free a memory block which has been allocated with av_malloc(z)() or av_realloc() and set the pointer ...
int swri_resample_int16_mmx2(struct ResampleContext *c, int16_t *dst, const int16_t *src, int *consumed, int src_size, int dst_size, int update_ctx)
enum AVSampleFormat format
static void resample_free(ResampleContext **c)
static double bessel(double x)
0th order modified bessel function of the first kind.
int swri_realloc_audio(AudioData *a, int count)
int compensation_distance
int swri_resample_double(struct ResampleContext *c, double *dst, const double *src, int *consumed, int src_size, int dst_size, int update_ctx)
static int set_compensation(ResampleContext *c, int sample_delta, int compensation_distance)
enum AVResampleFilterType filter_type
struct Resampler const swri_resampler
int av_get_cpu_flags(void)
Return the flags which specify extensions supported by the CPU.
int in_buffer_count
cached buffer length
#define AV_LOG_ERROR
Something went wrong and cannot losslessly be recovered.
void av_free(void *ptr)
Free a memory block which has been allocated with av_malloc(z)() or av_realloc(). ...
Blackman Nuttall Windowed Sinc.
static int resample_flush(struct SwrContext *s)
#define av_unused
Disable warnings about deprecated features This is useful for sections of code kept for backward comp...
int av_reduce(int *dst_num, int *dst_den, int64_t num, int64_t den, int64_t max)
Reduce a fraction.
static void filter(MpegAudioContext *s, int ch, const short *samples, int incr)
int64_t av_rescale(int64_t a, int64_t b, int64_t c) av_const
Rescale a 64-bit integer with rounding to nearest.
static int64_t get_delay(struct SwrContext *s, int64_t base)
void * av_malloc(size_t size) av_malloc_attrib 1(1)
Allocate a block of size bytes with alignment suitable for all memory accesses (including vectors if ...
int in_sample_rate
input sample rate
static int build_filter(ResampleContext *c, void *filter, double factor, int tap_count, int alloc, int phase_count, int scale, int filter_type, int kaiser_beta)
builds a polyphase filterbank.
static ResampleContext * resample_init(ResampleContext *c, int out_rate, int in_rate, int filter_size, int phase_shift, int linear, double cutoff0, enum AVSampleFormat format, enum SwrFilterType filter_type, int kaiser_beta, double precision, int cheby)
AVSampleFormat
Audio Sample Formats.
int swri_resample_int32(struct ResampleContext *c, int32_t *dst, const int32_t *src, int *consumed, int src_size, int dst_size, int update_ctx)
Describe the class of an AVClass context structure.
static const int factor[16]
#define AV_CPU_FLAG_MMX2
SSE integer functions or AMD MMX ext.
static const uint16_t scale[4]
static struct twinvq_data tab
int planar
1 if planar audio, 0 otherwise
#define av_assert0(cond)
assert() equivalent, that is always enabled.
int swri_resample_int16_ssse3(struct ResampleContext *c, int16_t *dst, const int16_t *src, int *consumed, int src_size, int dst_size, int update_ctx)
int av_get_bytes_per_sample(enum AVSampleFormat sample_fmt)
Return number of bytes per sample.
uint8_t * ch[SWR_CH_MAX]
samples buffer per channel
#define AV_CPU_FLAG_SSSE3
Conroe SSSE3 functions.
int swri_resample_int16(struct ResampleContext *c, int16_t *dst, const int16_t *src, int *consumed, int src_size, int dst_size, int update_ctx)
void * av_mallocz(size_t size) av_malloc_attrib 1(1)
Allocate a block of size bytes with alignment suitable for all memory accesses (including vectors if ...
#define av_assert2(cond)
assert() equivalent, that does lie in speed critical code.