shithub: sox

ref: d6b3a9ef2c8997c28ce52b5482bdbdd804a915d5
dir: /src/mixer.c/

View raw version
/*
 * July 5, 1991
 * Copyright 1991 Lance Norskog And Sundry Contributors
 * This source code is freely redistributable and may be used for
 * any purpose.  This copyright notice must be maintained. 
 * Lance Norskog And Sundry Contributors are not responsible for 
 * the consequences of using this software.
 *
 * Channel duplication code by Graeme W. Gill - 93/5/18
 * General-purpose panning by Geoffrey H. Kuenning -- 2000/11/28
 */

/*
 * libSoX stereo/quad -> mono mixdown effect file.
 * and mono/stereo -> stereo/quad channel duplication.
 */

#include "sox_i.h"
#include <ctype.h>
#include <string.h>
#include <stdlib.h>

typedef struct mixerstuff {
        /* How to generate each output channel.  sources[i][j] */
        /* represents the fraction of channel i that should be passed */
        /* through to channel j on output, and so forth.  Channel 0 is */
        /* left front, channel 1 is right front, and 2 and 3 are left */
        /* and right rear, respectively. (GHK) */
        double  sources[4][4];
        int     num_pans;
        int     mix;                    /* How are we mixing it? */
} *mixer_t;

/* MIX_CENTER is shorthand to mix channels together at 50% each */
#define MIX_CENTER      0
#define MIX_LEFT        1
#define MIX_RIGHT       2
#define MIX_FRONT       3
#define MIX_BACK        4
#define MIX_SPECIFIED   5
#define MIX_LEFT_FRONT  6
#define MIX_RIGHT_FRONT 7
#define MIX_LEFT_BACK   8
#define MIX_RIGHT_BACK  9


/*
 * Process options
 */
static int getopts(sox_effect_t * effp, int n, char **argv) 
{
    mixer_t mixer = (mixer_t) effp->priv;
    double* pans = &mixer->sources[0][0];
    int i;

    for (i = 0;  i < 16;  i++)
        pans[i] = 0.0;
    mixer->mix = MIX_CENTER;
    mixer->num_pans = 0;

    /* Parse parameters.  Since we don't yet know the number of */
    /* input and output channels, we'll record the information for */
    /* later. */
    if (n == 1) {
        if(!strcmp(argv[0], "-l"))
            mixer->mix = MIX_LEFT;
        else if (!strcmp(argv[0], "-r"))
            mixer->mix = MIX_RIGHT;
        else if (!strcmp(argv[0], "-f"))
            mixer->mix = MIX_FRONT;
        else if (!strcmp(argv[0], "-b"))
            mixer->mix = MIX_BACK;
        else if (!strcmp(argv[0], "-1"))
            mixer->mix = MIX_LEFT_FRONT;
        else if (!strcmp(argv[0], "-2"))
            mixer->mix = MIX_RIGHT_FRONT;
        else if (!strcmp(argv[0], "-3"))
            mixer->mix = MIX_LEFT_BACK;
        else if (!strcmp(argv[0], "-4"))
            mixer->mix = MIX_RIGHT_BACK;
        else if (argv[0][0] == '-' && !isdigit((int)argv[0][1])
                && argv[0][1] != '.')
          return sox_usage(effp);
        else {
            int commas;
            char *s;
            mixer->mix = MIX_SPECIFIED;
            pans[0] = atof(argv[0]);
            for (s = argv[0], commas = 0; *s; ++s) {
                if (*s == ',') {
                    ++commas;
                    if (commas >= 16) {
                        sox_fail("mixer can only take up to 16 pan values");
                        return (SOX_EOF);
                    }
                    pans[commas] = atof(s+1);
                }
            }
            mixer->num_pans = commas + 1;
        }
    }
    else if (n == 0) {
        mixer->mix = MIX_CENTER;
    }
    else
      return sox_usage(effp);

    return (SOX_SUCCESS);
}

/*
 * Start processing
 */
static int start(sox_effect_t * effp)
{
    /*
       Hmmm, this is tricky.  Lemme think:
       channel orders are [0][0],[0][1], etc.
       i.e., 0->0, 0->1, 0->2, 0->3, 1->0, 1->1, ...
       trailing zeros are omitted
       L/R balance is x= -1 for left only, 1 for right only
       1->1 channel effects:
       changing volume by x is x,0,0,0
       1->2 channel effects:
       duplicating everywhere is 1,1,0,0
       1->4 channel effects:
       duplicating everywhere is 1,1,1,1
       2->1 channel effects:
       left only is 1,0,0,0 0,0,0,0
       right only is 0,0,0,0 1,0,0,0
       left+right is 0.5,0,0,0 0.5,0,0,0
       left-right is 1,0,0,0 -1,0,0,0
       2->2 channel effects:
       L/R balance can be done several ways.  The standard stereo
       way is both the easiest and the most sensible:
       min(1-x,1),0,0,0 0,min(1+x,1),0,0
       left to both is 1,1,0,0
       right to both is 0,0,0,0 1,1,0,0
       left+right to both is 0.5,0.5,0,0 0.5,0.5,0,0
       left-right to both is 1,1,0,0 -1,-1,0,0
       left-right to left, right-left to right is 1,-1,0,0 -1,1,0,0
       2->4 channel effects:
       front duplicated into rear is 1,0,1,0 0,1,0,1
       front swapped into rear (why?) is 1,0,0,1 0,1,1,0
       front put into rear as mono (why?) is 1,0,0.5,0.5 0,1,0.5,0.5
       4->1 channel effects:
       left front only is 1,0,0,0
       left only is 0.5,0,0,0 0,0,0,0 0.5,0,0,0
       etc.
       4->2 channel effects:
       merge front/back is 0.5,0,0,0 0,0.5,0,0 0.5,0,0,0 0,0.5,0,0
       selections similar to above
       4->4 channel effects:
       left front to all is 1,1,1,1 0,0,0,0
       right front to all is 0,0,0,0 1,1,1,1
       left f/r to all f/r is 1,1,0,0 0,0,0,0 0,0,1,1 0,0,0,0
       etc.

       The interesting cases from above (deserving of abbreviations of
       less than 16 numbers) are:

       0) n->n volume change (1 number)
       1) 1->n duplication (0 numbers)
       2) 2->1 mixdown (0 or 2 numbers)
       3) 2->2 balance (1 number)
       4) 2->2 fully general mix (4 numbers)
       5) 2->4 duplication (0 numbers)
       6) 4->1 mixdown (0 or 4 numbers)
       7) 4->2 mixdown (0, or 2 numbers)
       8) 4->4 balance (1 or 2 numbers)

       The above has one ambiguity: n->n volume change conflicts with
       n->n balance for n != 1.  In such a case, we'll prefer
       balance, since there is already a volume effect in vol.c.

       GHK 2000/11/28
     */
     mixer_t mixer = (mixer_t) effp->priv;
     double pans[16];
     int i, j;
     int ichan, ochan;

     for (i = 0;  i < 16;  i++)
         pans[i] = ((double*)&mixer->sources[0][0])[i];

     ichan = effp->ininfo.channels;
     ochan = effp->outinfo.channels;
     if (ochan == -1) {
         sox_fail("Output must have known number of channels to use mixer effect");
         return(SOX_EOF);
     }

     if ((ichan != 1 && ichan != 2 && ichan != 4 &&
          mixer->mix != MIX_CENTER && ochan != 1)
             ||  (ochan != 1 && ochan != 2 && ochan != 4)) {
         sox_fail("Can't average %d channels into %d channels",
                 ichan, ochan);
         return (SOX_EOF);
     }

     /* Handle the special-case flags */
     switch (mixer->mix) {
         case MIX_CENTER:
             if (ichan == ochan)
               return SOX_EFF_NULL;
             break;             /* Code below will handle this case */
         case MIX_LEFT:
             if (ichan == 2 && ochan == 1)
             {
                 pans[0] = 1.0;
                 pans[1] = 0.0;
                 mixer->num_pans = 2;
             }
             else if (ichan == 4 && ochan == 1)
             {
                 pans[0] = 0.5;
                 pans[1] = 0.0;
                 pans[2] = 0.5;
                 pans[3] = 0.0;
                 mixer->num_pans = 4;
             }
             else
             {
                 sox_fail("Can't average %d channels into %d channels",
                         ichan, ochan);
                 return SOX_EOF;
             }
             break;
         case MIX_RIGHT:
             if (ichan == 2 && ochan == 1)
             {
                 pans[0] = 0.0;
                 pans[1] = 1.0;
                 mixer->num_pans = 2;
             }
             else if (ichan == 4 && ochan == 1)
             {
                 pans[0] = 0.0;
                 pans[1] = 0.5;
                 pans[2] = 0.0;
                 pans[3] = 0.5;
                 mixer->num_pans = 4;
             }
             else
             {
                 sox_fail("Can't average %d channels into %d channels",
                         ichan, ochan);
                 return SOX_EOF;
             }
             break;
         case MIX_FRONT:
             if (ichan == 4 && ochan == 2)
             {
                 pans[0] = 1.0;
                 pans[1] = 0.0;
                 mixer->num_pans = 2;
             }
             else
             {
                 sox_fail("-f option requires 4 channels input and 2 channel output");
                 return SOX_EOF;
             }
             break;
         case MIX_BACK:
             if (ichan == 4 && ochan == 2)
             {
                 pans[0] = 0.0;
                 pans[1] = 1.0;
                 mixer->num_pans = 2;
             }
             else
             {
                 sox_fail("-b option requires 4 channels input and 2 channel output");
                 return SOX_EOF;
             }
             break;
         case MIX_LEFT_FRONT:
             if (ichan == 2 && ochan == 1)
             {
                 pans[0] = 1.0;
                 pans[1] = 0.0;
                 mixer->num_pans = 2;
             }
             else if (ichan == 4 && ochan == 1)
             {
                 pans[0] = 1.0;
                 pans[1] = 0.0;
                 pans[2] = 0.0;
                 pans[3] = 0.0;
                 mixer->num_pans = 4;
             }
             else
             {
                 sox_fail("-1 option requires 4 channels input and 1 channel output");
                 return SOX_EOF;
             }
             break;
         case MIX_RIGHT_FRONT:
             if (ichan == 2 && ochan == 1)
             {
                 pans[0] = 0.0;
                 pans[1] = 1.0;
                 mixer->num_pans = 2;
             }
             else if (ichan == 4 && ochan == 1)
             {
                 pans[0] = 0.0;
                 pans[1] = 1.0;
                 pans[2] = 0.0;
                 pans[3] = 0.0;
                 mixer->num_pans = 4;
             }
             else
             {
                 sox_fail("-2 option requires 4 channels input and 1 channel output");
                 return SOX_EOF;
             }
             break;
         case MIX_LEFT_BACK:
             if (ichan == 4 && ochan == 1)
             {
                 pans[0] = 0.0;
                 pans[1] = 0.0;
                 pans[2] = 1.0;
                 pans[3] = 0.0;
                 mixer->num_pans = 4;
             }
             else
             {
                 sox_fail("-3 option requires 4 channels input and 1 channel output");
                 return SOX_EOF;
             }
         case MIX_RIGHT_BACK:
             if (ichan == 4 && ochan == 1)
             {
                 pans[0] = 0.0;
                 pans[1] = 0.0;
                 pans[2] = 0.0;
                 pans[3] = 1.0;
                 mixer->num_pans = 4;
             }
             else
             {
                 sox_fail("-4 option requires 4 channels input and 1 channel output");
                 return SOX_EOF;
             }

         case MIX_SPECIFIED:
             break;
         default:
             sox_fail("Unknown mix option in average effect");
             return SOX_EOF;
     }

     /* If number of pans is 4 or less then its a shorthand
      * representation.  If user specified it, then we have
      * garbage in our sources[][] array.  Need to clear that
      * now that all data is stored in pans[] array.
      */
     if (mixer->num_pans <= 4)
     {
         for (i = 0; i < ichan; i++)
         {
             for (j = 0; j < ochan; j++) 
             {
                 mixer->sources[i][j] = 0;
             }
         }
     }

     /* If the number of pans given is 4 or fewer, handle the special */
     /* cases listed in the comments above.  The code is lengthy but */
     /* straightforward. */
     if (mixer->num_pans == 0) {
         /* CASE 1 */
         if (ichan == 1 && ochan > ichan) {
             mixer->sources[0][0] = 1.0;
             mixer->sources[0][1] = 1.0;
             mixer->sources[0][2] = 1.0;
             mixer->sources[0][3] = 1.0;
         }
         /* CASE 2, 6 */
         else if (ochan == 1) {
             mixer->sources[0][0] = 1.0 / ichan;
         }
         /* CASE 5 */
         else if (ichan == 2 && ochan == 4) {
             mixer->sources[0][0] = 1.0;
             mixer->sources[0][2] = 1.0;
             mixer->sources[1][1] = 1.0;
             mixer->sources[1][3] = 1.0;
         }
         /* CASE 7 */
         else if (ichan == 4 && ochan == 2) {
             mixer->sources[0][0] = 0.5;
             mixer->sources[1][1] = 0.5;
             mixer->sources[2][0] = 0.5;
             mixer->sources[3][1] = 0.5;
         }
         else {
             sox_fail("You must specify at least one mix level when using mixer with an unusual number of channels.");
             return(SOX_EOF);
         }
     }
     else if (mixer->num_pans == 1) {
         /* Might be volume change or balance change */
         /* CASE 3 and CASE 8 */
         if ((ichan == 2 || ichan == 4) &&  ichan == ochan) {
             /* -1 is left only, 1 is right only */
             if (pans[0] <= 0.0) {
                 mixer->sources[1][1] = pans[0] + 1.0;
                 if (mixer->sources[1][1] < 0.0)
                     mixer->sources[1][1] = 0.0;
                 mixer->sources[0][0] = 1.0;
             }
             else {
                 mixer->sources[0][0] = 1.0 - pans[0];
                 if (mixer->sources[0][0] < 0.0)
                     mixer->sources[0][0] = 0.0;
                 mixer->sources[1][1] = 1.0;
             }
             if (ichan == 4) {
                 mixer->sources[2][2] = mixer->sources[0][0];
                 mixer->sources[3][3] = mixer->sources[1][1];
             }
         }
         else
         {
             sox_fail("Invalid options specified to mixer while not mixing");
             return SOX_EOF;
         }
     }
     else if (mixer->num_pans == 2) {
         /* CASE 2 */
         if (ichan == 2 && ochan == 1) {
             mixer->sources[0][0] = pans[0];
             mixer->sources[1][0] = pans[1];
         }
         /* CASE 7 */
         else if (ichan == 4 && ochan == 2) {
             mixer->sources[0][0] = pans[0];
             mixer->sources[1][1] = pans[0];
             mixer->sources[2][0] = pans[1];
             mixer->sources[3][1] = pans[1];
         }
         /* CASE 8 */
         else if (ichan == 4 && ochan == 4) {
             /* pans[0] is front -> front, pans[1] is for back */
             mixer->sources[0][0] = pans[0];
             mixer->sources[1][1] = pans[0];
             mixer->sources[2][2] = pans[1];
             mixer->sources[3][3] = pans[1];
         }
         else
         {
             sox_fail("Invalid options specified to mixer for this channel combination");
             return SOX_EOF;
         }
     }
     else if (mixer->num_pans == 4) {
         /* CASE 4 */
         if (ichan == 2 && ochan == 2) {
             /* Shorthand for 2-channel case */
             mixer->sources[0][0] = pans[0];
             mixer->sources[0][1] = pans[1];
             mixer->sources[1][0] = pans[2];
             mixer->sources[1][1] = pans[3];
         }
         /* CASE 6 */
         else if (ichan == 4 && ochan == 1) {
             mixer->sources[0][0] = pans[0];
             mixer->sources[1][0] = pans[1];
             mixer->sources[2][0] = pans[2];
             mixer->sources[3][0] = pans[3];
         }
         else
         {
             sox_fail("Invalid options specified to mixer for this channel combination");
             return SOX_EOF;
         }
     }
     else
     {
         sox_fail("Invalid options specified to mixer while not mixing");
         return SOX_EOF;
     }

#if 0  /* TODO: test the following: */
     if (effp->ininfo.channels != effp->outinfo.channels)
       return SOX_SUCCESS;

     for (i = 0; i < (int)effp->ininfo.channels; ++i)
       for (j = 0; j < (int)effp->outinfo.channels; ++j)
         if (avg->sources[i][j] != (i == j))
           return SOX_SUCCESS;

     return SOX_EFF_NULL;
#else
     return SOX_SUCCESS;
#endif
}

/*
 * Process either isamp or osamp samples, whichever is smaller.
 */

static int flow(sox_effect_t * effp, const sox_ssample_t *ibuf, sox_ssample_t *obuf, 
                sox_size_t *isamp, sox_size_t *osamp)
{
    mixer_t mixer = (mixer_t) effp->priv;
    sox_size_t len, done;
    int ichan, ochan;
    int i, j;
    double samp;

    ichan = effp->ininfo.channels;
    ochan = effp->outinfo.channels;
    len = *isamp / ichan;
    if (len > *osamp / ochan)
        len = *osamp / ochan;
    for (done = 0; done < len; done++, ibuf += ichan, obuf += ochan) {
        for (j = 0; j < ochan; j++) {
            samp = 0.0;
            for (i = 0; i < ichan; i++)
                samp += ibuf[i] * mixer->sources[mixer->mix == MIX_CENTER? 0 : i][j];
            SOX_SAMPLE_CLIP_COUNT(samp, effp->clips);
            obuf[j] = samp;
        }
    }
    *isamp = len * ichan;
    *osamp = len * ochan;
    return (SOX_SUCCESS);
}

sox_effect_handler_t const * sox_mixer_effect_fn(void)
{
  static sox_effect_handler_t handler = {
    "mixer",
    "[ -l | -r | -f | -b | -1 | -2 | -3 | -4 | n,n,n...,n ]",
    SOX_EFF_MCHAN | SOX_EFF_CHAN,
    getopts, start, flow, 0, 0, 0
  };
  return &handler;
}

sox_effect_handler_t const * sox_avg_effect_fn(void)
{
  static sox_effect_handler_t handler;
  handler = *sox_mixer_effect_fn();
  handler.name = "avg";
  handler.flags |= SOX_EFF_DEPRECATED;
  return &handler;
}

sox_effect_handler_t const * sox_pick_effect_fn(void)
{
  static sox_effect_handler_t handler;
  handler = *sox_avg_effect_fn();
  handler.name = "pick";
  return &handler;
}

static int oops_getopts(sox_effect_t * effp, int argc, char * * argv UNUSED) 
{
  char * args[] = {"1,1,-1,-1"};
  return argc? sox_usage(effp) : sox_mixer_effect_fn()->getopts(effp, array_length(args), args);
}

sox_effect_handler_t const * sox_oops_effect_fn(void)
{
  static sox_effect_handler_t handler;
  handler = *sox_mixer_effect_fn();
  handler.name = "oops";
  handler.usage = NULL;
  handler.getopts = oops_getopts;
  return &handler;
}