libaf/af_lavcresample.c

   1 /*
   2  * Copyright (c) 2004 Michael Niedermayer <michaelni@gmx.at>
   3  *
   4  * This file is part of MPlayer.
   5  *
   6  * MPlayer is free software; you can redistribute it and/or modify
   7  * it under the terms of the GNU General Public License as published by
   8  * the Free Software Foundation; either version 2 of the License, or
   9  * (at your option) any later version.
  10  *
  11  * MPlayer 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
  14  * GNU General Public License for more details.
  15  *
  16  * You should have received a copy of the GNU General Public License along
  17  * with MPlayer; if not, write to the Free Software Foundation, Inc.,
  18  * 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
  19  */
  20
  21 #include <stdio.h>
  22 #include <stdlib.h>
  23 #include <string.h>
  24 #include <inttypes.h>
  25
  26 #include "config.h"
  27 #include "af.h"
  28 #include "libavcodec/avcodec.h"
  29 #include "libavutil/rational.h"
  30
  31 // Data for specific instances of this filter
  32 typedef struct af_resample_s{
  33     struct AVResampleContext *avrctx;
  34     int16_t *in[AF_NCH];
  35     int in_alloc;
  36     int index;
  37
  38     int filter_length;
  39     int linear;
  40     int phase_shift;
  41     double cutoff;
  42 }af_resample_t;
  43
  44
  45 // Initialization and runtime control
  46 static int control(struct af_instance_s* af, int cmd, void* arg)
  47 {
  48   af_resample_t* s   = (af_resample_t*)af->setup;
  49   af_data_t *data= (af_data_t*)arg;
  50   int out_rate, test_output_res; // helpers for checking input format
  51
  52   switch(cmd){
  53   case AF_CONTROL_REINIT:
  54     if((af->data->rate == data->rate) || (af->data->rate == 0))
  55         return AF_DETACH;
  56
  57     af->data->nch    = data->nch;
  58     if (af->data->nch > AF_NCH) af->data->nch = AF_NCH;
  59     af->data->format = AF_FORMAT_S16_NE;
  60     af->data->bps    = 2;
  61     af->mul = (double)af->data->rate / data->rate;
  62     af->delay = af->data->nch * s->filter_length / min(af->mul, 1); // *bps*.5
  63
  64     if(s->avrctx) av_resample_close(s->avrctx);
  65     s->avrctx= av_resample_init(af->data->rate, /*in_rate*/data->rate, s->filter_length, s->phase_shift, s->linear, s->cutoff);
  66
  67     // hack to make af_test_output ignore the samplerate change
  68     out_rate = af->data->rate;
  69     af->data->rate = data->rate;
  70     test_output_res = af_test_output(af, (af_data_t*)arg);
  71     af->data->rate = out_rate;
  72     return test_output_res;
  73   case AF_CONTROL_COMMAND_LINE:{
  74     s->cutoff= 0.0;
  75     sscanf((char*)arg,"%d:%d:%d:%d:%lf", &af->data->rate, &s->filter_length, &s->linear, &s->phase_shift, &s->cutoff);
  76     if(s->cutoff <= 0.0) s->cutoff= max(1.0 - 6.5/(s->filter_length+8), 0.80);
  77     return AF_OK;
  78   }
  79   case AF_CONTROL_RESAMPLE_RATE | AF_CONTROL_SET:
  80     af->data->rate = *(int*)arg;
  81     return AF_OK;
  82   }
  83   return AF_UNKNOWN;
  84 }
  85
  86 // Deallocate memory
  87 static void uninit(struct af_instance_s* af)
  88 {
  89     if(af->data)
  90         free(af->data->audio);
  91     free(af->data);
  92     if(af->setup){
  93         int i;
  94         af_resample_t *s = af->setup;
  95         if(s->avrctx) av_resample_close(s->avrctx);
  96         for (i=0; i < AF_NCH; i++)
  97             free(s->in[i]);
  98         free(s);
  99     }
 100 }
 101
 102 // Filter data through filter
 103 static af_data_t* play(struct af_instance_s* af, af_data_t* data)
 104 {
 105   af_resample_t *s = af->setup;
 106   int i, j, consumed, ret;
 107   int16_t *in = (int16_t*)data->audio;
 108   int16_t *out;
 109   int chans   = data->nch;
 110   int in_len  = data->len/(2*chans);
 111   int out_len = in_len * af->mul + 10;
 112   int16_t tmp[AF_NCH][out_len];
 113
 114   if(AF_OK != RESIZE_LOCAL_BUFFER(af,data))
 115       return NULL;
 116
 117   out= (int16_t*)af->data->audio;
 118
 119   out_len= min(out_len, af->data->len/(2*chans));
 120
 121   if(s->in_alloc < in_len + s->index){
 122       s->in_alloc= in_len + s->index;
 123       for(i=0; i<chans; i++){
 124           s->in[i]= realloc(s->in[i], s->in_alloc*sizeof(int16_t));
 125       }
 126   }
 127
 128   if(chans==1){
 129       memcpy(&s->in[0][s->index], in, in_len * sizeof(int16_t));
 130   }else if(chans==2){
 131       for(j=0; j<in_len; j++){
 132           s->in[0][j + s->index]= *(in++);
 133           s->in[1][j + s->index]= *(in++);
 134       }
 135   }else{
 136       for(j=0; j<in_len; j++){
 137           for(i=0; i<chans; i++){
 138               s->in[i][j + s->index]= *(in++);
 139           }
 140       }
 141   }
 142   in_len += s->index;
 143
 144   for(i=0; i<chans; i++){
 145       ret= av_resample(s->avrctx, tmp[i], s->in[i], &consumed, in_len, out_len, i+1 == chans);
 146   }
 147   out_len= ret;
 148
 149   s->index= in_len - consumed;
 150   for(i=0; i<chans; i++){
 151       memmove(s->in[i], s->in[i] + consumed, s->index*sizeof(int16_t));
 152   }
 153
 154   if(chans==1){
 155       memcpy(out, tmp[0], out_len*sizeof(int16_t));
 156   }else if(chans==2){
 157       for(j=0; j<out_len; j++){
 158           *(out++)= tmp[0][j];
 159           *(out++)= tmp[1][j];
 160       }
 161   }else{
 162       for(j=0; j<out_len; j++){
 163           for(i=0; i<chans; i++){
 164               *(out++)= tmp[i][j];
 165           }
 166       }
 167   }
 168
 169   data->audio = af->data->audio;
 170   data->len   = out_len*chans*2;
 171   data->rate  = af->data->rate;
 172   return data;
 173 }
 174
 175 static int af_open(af_instance_t* af){
 176   af_resample_t *s = calloc(1,sizeof(af_resample_t));
 177   af->control=control;
 178   af->uninit=uninit;
 179   af->play=play;
 180   af->mul=1;
 181   af->data=calloc(1,sizeof(af_data_t));
 182   s->filter_length= 16;
 183   s->cutoff= max(1.0 - 6.5/(s->filter_length+8), 0.80);
 184   s->phase_shift= 10;
 185 //  s->setup = RSMP_INT | FREQ_SLOPPY;
 186   af->setup=s;
 187   return AF_OK;
 188 }
 189
 190 af_info_t af_info_lavcresample = {
 191   "Sample frequency conversion using libavcodec",
 192   "lavcresample",
 193   "Michael Niedermayer",
 194   "",
 195   AF_FLAGS_REENTRANT,
 196   af_open
 197 };