gst/audioconvert/gstchannelmix.c

   1 /* GStreamer
   2  * Copyright (C) 2004 Ronald Bultje <rbultje@ronald.bitfreak.net>
   3  *
   4  * gstchannelmix.c: setup of channel conversion matrices
   5  *
   6  * This library is free software; you can redistribute it and/or
   7  * modify it under the terms of the GNU Library General Public
   8  * License as published by the Free Software Foundation; either
   9  * version 2 of the License, or (at your option) any later version.
  10  *
  11  * This library 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  * Library General Public License for more details.
  15  *
  16  * You should have received a copy of the GNU Library General Public
  17  * License along with this library; if not, write to the
  18  * Free Software Foundation, Inc., 59 Temple Place - Suite 330,
  19  * Boston, MA 02111-1307, USA.
  20  */
  21
  22 #ifdef HAVE_CONFIG_H
  23 #include "config.h"
  24 #endif
  25
  26 #include <math.h>
  27 #include <string.h>
  28 #include <gst/audio/multichannel.h>
  29
  30 #include "gstchannelmix.h"
  31
  32 /*
  33  * Channel matrix functions.
  34  */
  35
  36 void
  37 gst_channel_mix_unset_matrix (AudioConvertCtx * this)
  38 {
  39   gint i;
  40
  41   /* don't access if nothing there */
  42   if (!this->matrix)
  43     return;
  44
  45   /* free */
  46   for (i = 0; i < this->in.channels; i++)
  47     g_free (this->matrix[i]);
  48   g_free (this->matrix);
  49
  50   this->matrix = NULL;
  51   g_free (this->tmp);
  52   this->tmp = NULL;
  53 }
  54
  55 /*
  56  * Detect and fill in identical channels. E.g.
  57  * forward the left/right front channels in a
  58  * 5.1 to 2.0 conversion.
  59  */
  60
  61 static void
  62 gst_channel_mix_fill_identical (AudioConvertCtx * this)
  63 {
  64   gint ci, co;
  65
  66   /* Apart from the compatible channel assignments, we can also have
  67    * same channel assignments. This is much simpler, we simply copy
  68    * the value from source to dest! */
  69   for (co = 0; co < this->out.channels; co++) {
  70     /* find a channel in input with same position */
  71     for (ci = 0; ci < this->in.channels; ci++) {
  72       if (this->in.pos[ci] == this->out.pos[co]) {
  73         this->matrix[ci][co] = 1.0;
  74       }
  75     }
  76   }
  77 }
  78
  79 /*
  80  * Detect and fill in compatible channels. E.g.
  81  * forward left/right front to mono (or the other
  82  * way around) when going from 2.0 to 1.0.
  83  */
  84
  85 static void
  86 gst_channel_mix_fill_compatible (AudioConvertCtx * this)
  87 {
  88   /* Conversions from one-channel to compatible two-channel configs */
  89   struct
  90   {
  91     GstAudioChannelPosition pos1[2];
  92     GstAudioChannelPosition pos2[1];
  93   } conv[] = {
  94     /* front: mono <-> stereo */
  95     { {
  96     GST_AUDIO_CHANNEL_POSITION_FRONT_LEFT,
  97             GST_AUDIO_CHANNEL_POSITION_FRONT_RIGHT}, {
  98     GST_AUDIO_CHANNEL_POSITION_FRONT_MONO}},
  99         /* front center: 2 <-> 1 */
 100     { {
 101     GST_AUDIO_CHANNEL_POSITION_FRONT_LEFT_OF_CENTER,
 102             GST_AUDIO_CHANNEL_POSITION_FRONT_RIGHT_OF_CENTER}, {
 103     GST_AUDIO_CHANNEL_POSITION_FRONT_CENTER}},
 104         /* rear: 2 <-> 1 */
 105     { {
 106     GST_AUDIO_CHANNEL_POSITION_REAR_LEFT,
 107             GST_AUDIO_CHANNEL_POSITION_REAR_RIGHT}, {
 108     GST_AUDIO_CHANNEL_POSITION_REAR_CENTER}}, { {
 109     GST_AUDIO_CHANNEL_POSITION_INVALID}}
 110   };
 111   gint c;
 112
 113   /* conversions from compatible (but not the same) channel schemes. This
 114    * goes two ways: if the sink has both pos1[0,1] and src has pos2[0] or
 115    * if the src has both pos1[0,1] and sink has pos2[0], then we do the
 116    * conversion. We hereby assume that the existance of pos1[0,1] and
 117    * pos2[0] are mututally exclusive. There are no checks for that,
 118    * unfortunately. This shouldn't lead to issues (like crashes or so),
 119    * though. */
 120   for (c = 0; conv[c].pos1[0] != GST_AUDIO_CHANNEL_POSITION_INVALID; c++) {
 121     gint pos1_0 = -1, pos1_1 = -1, pos2_0 = -1, n;
 122
 123     /* Try to go from the given 2 channels to the given 1 channel */
 124     for (n = 0; n < this->in.channels; n++) {
 125       if (this->in.pos[n] == conv[c].pos1[0])
 126         pos1_0 = n;
 127       else if (this->in.pos[n] == conv[c].pos1[1])
 128         pos1_1 = n;
 129     }
 130     for (n = 0; n < this->out.channels; n++) {
 131       if (this->out.pos[n] == conv[c].pos2[0])
 132         pos2_0 = n;
 133     }
 134
 135     if (pos1_0 != -1 && pos1_1 != -1 && pos2_0 != -1) {
 136       this->matrix[pos1_0][pos2_0] = 1.0;
 137       this->matrix[pos1_1][pos2_0] = 1.0;
 138     }
 139
 140     /* Try to go from the given 1 channel to the given 2 channels */
 141     pos1_0 = -1;
 142     pos1_1 = -1;
 143     pos2_0 = -1;
 144
 145     for (n = 0; n < this->out.channels; n++) {
 146       if (this->out.pos[n] == conv[c].pos1[0])
 147         pos1_0 = n;
 148       else if (this->out.pos[n] == conv[c].pos1[1])
 149         pos1_1 = n;
 150     }
 151     for (n = 0; n < this->in.channels; n++) {
 152       if (this->in.pos[n] == conv[c].pos2[0])
 153         pos2_0 = n;
 154     }
 155
 156     if (pos1_0 != -1 && pos1_1 != -1 && pos2_0 != -1) {
 157       this->matrix[pos2_0][pos1_0] = 1.0;
 158       this->matrix[pos2_0][pos1_1] = 1.0;
 159     }
 160   }
 161 }
 162
 163 /*
 164  * Detect and fill in channels not handled by the
 165  * above two, e.g. center to left/right front in
 166  * 5.1 to 2.0 (or the other way around).
 167  *
 168  * Unfortunately, limited to static conversions
 169  * for now.
 170  */
 171
 172 static void
 173 gst_channel_mix_detect_pos (AudioConvertFmt * caps,
 174     gint * f, gboolean * has_f,
 175     gint * c, gboolean * has_c, gint * r, gboolean * has_r,
 176     gint * s, gboolean * has_s, gint * b, gboolean * has_b)
 177 {
 178   gint n;
 179
 180   for (n = 0; n < caps->channels; n++) {
 181     switch (caps->pos[n]) {
 182       case GST_AUDIO_CHANNEL_POSITION_FRONT_MONO:
 183       case GST_AUDIO_CHANNEL_POSITION_FRONT_LEFT:
 184       case GST_AUDIO_CHANNEL_POSITION_FRONT_RIGHT:
 185         *has_f = TRUE;
 186         if (f[0] == -1)
 187           f[0] = n;
 188         else
 189           f[1] = n;
 190         break;
 191       case GST_AUDIO_CHANNEL_POSITION_FRONT_CENTER:
 192       case GST_AUDIO_CHANNEL_POSITION_FRONT_LEFT_OF_CENTER:
 193       case GST_AUDIO_CHANNEL_POSITION_FRONT_RIGHT_OF_CENTER:
 194         *has_c = TRUE;
 195         if (c[0] == -1)
 196           c[0] = n;
 197         else
 198           c[1] = n;
 199         break;
 200       case GST_AUDIO_CHANNEL_POSITION_REAR_CENTER:
 201       case GST_AUDIO_CHANNEL_POSITION_REAR_LEFT:
 202       case GST_AUDIO_CHANNEL_POSITION_REAR_RIGHT:
 203         *has_r = TRUE;
 204         if (r[0] == -1)
 205           r[0] = n;
 206         else
 207           r[1] = n;
 208         break;
 209       case GST_AUDIO_CHANNEL_POSITION_SIDE_LEFT:
 210       case GST_AUDIO_CHANNEL_POSITION_SIDE_RIGHT:
 211         *has_s = TRUE;
 212         if (s[0] == -1)
 213           s[0] = n;
 214         else
 215           s[1] = n;
 216         break;
 217       case GST_AUDIO_CHANNEL_POSITION_LFE:
 218         *has_b = TRUE;
 219         b[0] = n;
 220         break;
 221       default:
 222         break;
 223     }
 224   }
 225 }
 226
 227 static void
 228 gst_channel_mix_fill_one_other (gfloat ** matrix,
 229     AudioConvertFmt * from_caps, gint * from_idx,
 230     GstAudioChannelPosition from_pos_l,
 231     GstAudioChannelPosition from_pos_r,
 232     GstAudioChannelPosition from_pos_c,
 233     AudioConvertFmt * to_caps, gint * to_idx,
 234     GstAudioChannelPosition to_pos_l,
 235     GstAudioChannelPosition to_pos_r,
 236     GstAudioChannelPosition to_pos_c, gfloat ratio)
 237 {
 238   gfloat in_r, out_r[2] = { 0.f, 0.f };
 239
 240   /*
 241    * The idea is that we add up from the input (which means that if we
 242    * have stereo input, we divide their sum by two) and put that in
 243    * the matrix for their output ratio (given in $ratio).
 244    * For left channels, we need to invert the signal sign (* -1).
 245    */
 246
 247   if (from_caps->pos[from_idx[0]] == from_pos_c)
 248     in_r = 1.0;
 249   else
 250     in_r = 0.5;
 251
 252   if (to_caps->pos[to_idx[0]] == to_pos_l)
 253     out_r[0] = in_r * -ratio;
 254   else
 255     out_r[0] = in_r * ratio;
 256
 257   if (to_idx[1] != -1) {
 258     if (to_caps->pos[to_idx[1]] == to_pos_l)
 259       out_r[1] = in_r * -ratio;
 260     else
 261       out_r[1] = in_r * ratio;
 262   }
 263
 264   matrix[from_idx[0]][to_idx[0]] = out_r[0];
 265   if (to_idx[1] != -1)
 266     matrix[from_idx[0]][to_idx[1]] = out_r[1];
 267   if (from_idx[1] != -1) {
 268     matrix[from_idx[1]][to_idx[0]] = out_r[0];
 269     if (to_idx[1] != -1)
 270       matrix[from_idx[1]][to_idx[1]] = out_r[1];
 271   }
 272 }
 273
 274 #define RATIO_FRONT_CENTER (1.0 / sqrt (2.0))
 275 #define RATIO_FRONT_REAR (1.0 / sqrt (2.0))
 276 #define RATIO_FRONT_BASS (1.0)
 277 #define RATIO_REAR_BASS (1.0 / sqrt (2.0))
 278 #define RATIO_CENTER_BASS (1.0 / sqrt (2.0))
 279
 280 static void
 281 gst_channel_mix_fill_others (AudioConvertCtx * this)
 282 {
 283   gboolean in_has_front = FALSE, out_has_front = FALSE,
 284       in_has_center = FALSE, out_has_center = FALSE,
 285       in_has_rear = FALSE, out_has_rear = FALSE,
 286       in_has_side = FALSE, out_has_side = FALSE,
 287       in_has_bass = FALSE, out_has_bass = FALSE;
 288   gint in_f[2] = { -1, -1 }, out_f[2] = {
 289   -1, -1}, in_c[2] = {
 290   -1, -1}, out_c[2] = {
 291   -1, -1}, in_r[2] = {
 292   -1, -1}, out_r[2] = {
 293   -1, -1}, in_s[2] = {
 294   -1, -1}, out_s[2] = {
 295   -1, -1}, in_b[2] = {
 296   -1, -1}, out_b[2] = {
 297   -1, -1};
 298
 299   /* First see where (if at all) the various channels from/to
 300    * which we want to convert are located in our matrix/array. */
 301   gst_channel_mix_detect_pos (&this->in,
 302       in_f, &in_has_front,
 303       in_c, &in_has_center, in_r, &in_has_rear,
 304       in_s, &in_has_side, in_b, &in_has_bass);
 305   gst_channel_mix_detect_pos (&this->out,
 306       out_f, &out_has_front,
 307       out_c, &out_has_center, out_r, &out_has_rear,
 308       out_s, &out_has_side, out_b, &out_has_bass);
 309
 310   /* center/front */
 311   if (!in_has_center && in_has_front && out_has_center) {
 312     gst_channel_mix_fill_one_other (this->matrix,
 313         &this->in, in_f,
 314         GST_AUDIO_CHANNEL_POSITION_FRONT_LEFT,
 315         GST_AUDIO_CHANNEL_POSITION_FRONT_RIGHT,
 316         GST_AUDIO_CHANNEL_POSITION_FRONT_MONO,
 317         &this->out, out_c,
 318         GST_AUDIO_CHANNEL_POSITION_FRONT_LEFT_OF_CENTER,
 319         GST_AUDIO_CHANNEL_POSITION_FRONT_RIGHT_OF_CENTER,
 320         GST_AUDIO_CHANNEL_POSITION_FRONT_CENTER, RATIO_FRONT_CENTER);
 321   } else if (in_has_center && !out_has_center && out_has_front) {
 322     gst_channel_mix_fill_one_other (this->matrix,
 323         &this->in, in_c,
 324         GST_AUDIO_CHANNEL_POSITION_FRONT_LEFT_OF_CENTER,
 325         GST_AUDIO_CHANNEL_POSITION_FRONT_RIGHT_OF_CENTER,
 326         GST_AUDIO_CHANNEL_POSITION_FRONT_CENTER,
 327         &this->out, out_f,
 328         GST_AUDIO_CHANNEL_POSITION_FRONT_LEFT,
 329         GST_AUDIO_CHANNEL_POSITION_FRONT_RIGHT,
 330         GST_AUDIO_CHANNEL_POSITION_FRONT_MONO, RATIO_FRONT_CENTER);
 331   }
 332
 333   /* rear/front */
 334   if (!in_has_rear && in_has_front && out_has_rear) {
 335     gst_channel_mix_fill_one_other (this->matrix,
 336         &this->in, in_f,
 337         GST_AUDIO_CHANNEL_POSITION_FRONT_LEFT,
 338         GST_AUDIO_CHANNEL_POSITION_FRONT_RIGHT,
 339         GST_AUDIO_CHANNEL_POSITION_FRONT_MONO,
 340         &this->out, out_r,
 341         GST_AUDIO_CHANNEL_POSITION_REAR_LEFT,
 342         GST_AUDIO_CHANNEL_POSITION_REAR_RIGHT,
 343         GST_AUDIO_CHANNEL_POSITION_REAR_CENTER, RATIO_FRONT_REAR);
 344   } else if (in_has_rear && !out_has_rear && out_has_front) {
 345     gst_channel_mix_fill_one_other (this->matrix,
 346         &this->in, in_r,
 347         GST_AUDIO_CHANNEL_POSITION_REAR_LEFT,
 348         GST_AUDIO_CHANNEL_POSITION_REAR_RIGHT,
 349         GST_AUDIO_CHANNEL_POSITION_REAR_CENTER,
 350         &this->out, out_f,
 351         GST_AUDIO_CHANNEL_POSITION_FRONT_LEFT,
 352         GST_AUDIO_CHANNEL_POSITION_FRONT_RIGHT,
 353         GST_AUDIO_CHANNEL_POSITION_FRONT_MONO, RATIO_FRONT_REAR);
 354   }
 355
 356   /* bass/any */
 357   if (in_has_bass && !out_has_bass) {
 358     if (out_has_front) {
 359       gst_channel_mix_fill_one_other (this->matrix,
 360           &this->in, in_b,
 361           GST_AUDIO_CHANNEL_POSITION_INVALID,
 362           GST_AUDIO_CHANNEL_POSITION_INVALID,
 363           GST_AUDIO_CHANNEL_POSITION_LFE,
 364           &this->out, out_f,
 365           GST_AUDIO_CHANNEL_POSITION_FRONT_LEFT,
 366           GST_AUDIO_CHANNEL_POSITION_FRONT_RIGHT,
 367           GST_AUDIO_CHANNEL_POSITION_FRONT_MONO, RATIO_FRONT_BASS);
 368     }
 369     if (out_has_center) {
 370       gst_channel_mix_fill_one_other (this->matrix,
 371           &this->in, in_b,
 372           GST_AUDIO_CHANNEL_POSITION_INVALID,
 373           GST_AUDIO_CHANNEL_POSITION_INVALID,
 374           GST_AUDIO_CHANNEL_POSITION_LFE,
 375           &this->out, out_c,
 376           GST_AUDIO_CHANNEL_POSITION_FRONT_LEFT_OF_CENTER,
 377           GST_AUDIO_CHANNEL_POSITION_FRONT_RIGHT_OF_CENTER,
 378           GST_AUDIO_CHANNEL_POSITION_FRONT_CENTER, RATIO_CENTER_BASS);
 379     }
 380     if (out_has_rear) {
 381       gst_channel_mix_fill_one_other (this->matrix,
 382           &this->in, in_b,
 383           GST_AUDIO_CHANNEL_POSITION_INVALID,
 384           GST_AUDIO_CHANNEL_POSITION_INVALID,
 385           GST_AUDIO_CHANNEL_POSITION_LFE,
 386           &this->out, out_r,
 387           GST_AUDIO_CHANNEL_POSITION_REAR_LEFT,
 388           GST_AUDIO_CHANNEL_POSITION_REAR_RIGHT,
 389           GST_AUDIO_CHANNEL_POSITION_REAR_CENTER, RATIO_REAR_BASS);
 390     }
 391   } else if (!in_has_bass && out_has_bass) {
 392     if (in_has_front) {
 393       gst_channel_mix_fill_one_other (this->matrix,
 394           &this->in, in_f,
 395           GST_AUDIO_CHANNEL_POSITION_FRONT_LEFT,
 396           GST_AUDIO_CHANNEL_POSITION_FRONT_RIGHT,
 397           GST_AUDIO_CHANNEL_POSITION_FRONT_MONO,
 398           &this->out, out_b,
 399           GST_AUDIO_CHANNEL_POSITION_INVALID,
 400           GST_AUDIO_CHANNEL_POSITION_INVALID,
 401           GST_AUDIO_CHANNEL_POSITION_LFE, RATIO_FRONT_BASS);
 402     }
 403     if (in_has_center) {
 404       gst_channel_mix_fill_one_other (this->matrix,
 405           &this->in, in_c,
 406           GST_AUDIO_CHANNEL_POSITION_FRONT_LEFT_OF_CENTER,
 407           GST_AUDIO_CHANNEL_POSITION_FRONT_RIGHT_OF_CENTER,
 408           GST_AUDIO_CHANNEL_POSITION_FRONT_CENTER,
 409           &this->out, out_b,
 410           GST_AUDIO_CHANNEL_POSITION_INVALID,
 411           GST_AUDIO_CHANNEL_POSITION_INVALID,
 412           GST_AUDIO_CHANNEL_POSITION_LFE, RATIO_CENTER_BASS);
 413     }
 414     if (in_has_rear) {
 415       gst_channel_mix_fill_one_other (this->matrix,
 416           &this->in, in_r,
 417           GST_AUDIO_CHANNEL_POSITION_REAR_LEFT,
 418           GST_AUDIO_CHANNEL_POSITION_REAR_RIGHT,
 419           GST_AUDIO_CHANNEL_POSITION_REAR_CENTER,
 420           &this->out, out_b,
 421           GST_AUDIO_CHANNEL_POSITION_INVALID,
 422           GST_AUDIO_CHANNEL_POSITION_INVALID,
 423           GST_AUDIO_CHANNEL_POSITION_LFE, RATIO_REAR_BASS);
 424     }
 425   }
 426
 427   /* FIXME: side */
 428 }
 429
 430 /*
 431  * Normalize output values.
 432  */
 433
 434 static void
 435 gst_channel_mix_fill_normalize (AudioConvertCtx * this)
 436 {
 437   gfloat sum, top = 0;
 438   gint i, j;
 439
 440   for (j = 0; j < this->out.channels; j++) {
 441     /* calculate sum */
 442     sum = 0.0;
 443     for (i = 0; i < this->in.channels; i++) {
 444       sum += fabs (this->matrix[i][j]);
 445     }
 446     if (sum > top) {
 447       top = sum;
 448     }
 449   }
 450
 451   /* normalize to this */
 452   for (j = 0; j < this->out.channels; j++) {
 453     for (i = 0; i < this->in.channels; i++) {
 454       this->matrix[i][j] /= top;
 455     }
 456   }
 457 }
 458
 459 /*
 460  * Automagically generate conversion matrix.
 461  */
 462
 463 static void
 464 gst_channel_mix_fill_matrix (AudioConvertCtx * this)
 465 {
 466   gst_channel_mix_fill_identical (this);
 467
 468   if (!this->in.unpositioned_layout) {
 469     gst_channel_mix_fill_compatible (this);
 470     gst_channel_mix_fill_others (this);
 471     gst_channel_mix_fill_normalize (this);
 472   }
 473 }
 474
 475 /* only call after this->out and this->in are filled in */
 476 void
 477 gst_channel_mix_setup_matrix (AudioConvertCtx * this)
 478 {
 479   gint i, j;
 480   GString *s;
 481
 482   /* don't lose memory */
 483   gst_channel_mix_unset_matrix (this);
 484
 485   /* temp storage */
 486   if (this->in.is_int || this->out.is_int) {
 487     this->tmp = (gpointer) g_new (gint32, this->out.channels);
 488   } else {
 489     this->tmp = (gpointer) g_new (gdouble, this->out.channels);
 490   }
 491
 492   /* allocate */
 493   this->matrix = g_new0 (gfloat *, this->in.channels);
 494   for (i = 0; i < this->in.channels; i++) {
 495     this->matrix[i] = g_new (gfloat, this->out.channels);
 496     for (j = 0; j < this->out.channels; j++)
 497       this->matrix[i][j] = 0.;
 498   }
 499
 500   /* setup the matrix' internal values */
 501   gst_channel_mix_fill_matrix (this);
 502
 503   /* debug */
 504   s = g_string_new ("Matrix for");
 505   g_string_append_printf (s, " %d -> %d: ",
 506       this->in.channels, this->out.channels);
 507   g_string_append (s, "{");
 508   for (i = 0; i < this->in.channels; i++) {
 509     if (i != 0)
 510       g_string_append (s, ",");
 511     g_string_append (s, " {");
 512     for (j = 0; j < this->out.channels; j++) {
 513       if (j != 0)
 514         g_string_append (s, ",");
 515       g_string_append_printf (s, " %f", this->matrix[i][j]);
 516     }
 517     g_string_append (s, " }");
 518   }
 519   g_string_append (s, " }");
 520   GST_DEBUG (s->str);
 521   g_string_free (s, TRUE);
 522 }
 523
 524 gboolean
 525 gst_channel_mix_passthrough (AudioConvertCtx * this)
 526 {
 527   gint i;
 528
 529   /* only NxN matrices can be identities */
 530   if (this->in.channels != this->out.channels)
 531     return FALSE;
 532
 533   /* this assumes a normalized matrix */
 534   for (i = 0; i < this->in.channels; i++)
 535     if (this->matrix[i][i] != 1.)
 536       return FALSE;
 537
 538   return TRUE;
 539 }
 540
 541 /* IMPORTANT: out_data == in_data is possible, make sure to not overwrite data
 542  * you might need later on! */
 543 void
 544 gst_channel_mix_mix_int (AudioConvertCtx * this,
 545     gint32 * in_data, gint32 * out_data, gint samples)
 546 {
 547   gint in, out, n;
 548   gint64 res;
 549   gboolean backwards;
 550   gint inchannels, outchannels;
 551   gint32 *tmp = (gint32 *) this->tmp;
 552
 553   g_return_if_fail (this->matrix != NULL);
 554   g_return_if_fail (this->tmp != NULL);
 555
 556   inchannels = this->in.channels;
 557   outchannels = this->out.channels;
 558   backwards = outchannels > inchannels;
 559
 560   /* FIXME: use liboil here? */
 561   for (n = (backwards ? samples - 1 : 0); n < samples && n >= 0;
 562       backwards ? n-- : n++) {
 563     for (out = 0; out < outchannels; out++) {
 564       /* convert */
 565       res = 0;
 566       for (in = 0; in < inchannels; in++) {
 567         res += in_data[n * inchannels + in] * this->matrix[in][out];
 568       }
 569
 570       /* clip (shouldn't we use doubles instead as intermediate format?) */
 571       if (res < G_MININT32)
 572         res = G_MININT32;
 573       else if (res > G_MAXINT32)
 574         res = G_MAXINT32;
 575       tmp[out] = res;
 576     }
 577     memcpy (&out_data[n * outchannels], this->tmp,
 578         sizeof (gint32) * outchannels);
 579   }
 580 }
 581
 582 void
 583 gst_channel_mix_mix_float (AudioConvertCtx * this,
 584     gdouble * in_data, gdouble * out_data, gint samples)
 585 {
 586   gint in, out, n;
 587   gdouble res;
 588   gboolean backwards;
 589   gint inchannels, outchannels;
 590   gdouble *tmp = (gdouble *) this->tmp;
 591
 592   g_return_if_fail (this->matrix != NULL);
 593   g_return_if_fail (this->tmp != NULL);
 594
 595   inchannels = this->in.channels;
 596   outchannels = this->out.channels;
 597   backwards = outchannels > inchannels;
 598
 599   /* FIXME: use liboil here? */
 600   for (n = (backwards ? samples - 1 : 0); n < samples && n >= 0;
 601       backwards ? n-- : n++) {
 602     for (out = 0; out < outchannels; out++) {
 603       /* convert */
 604       res = 0.0;
 605       for (in = 0; in < inchannels; in++) {
 606         res += in_data[n * inchannels + in] * this->matrix[in][out];
 607       }
 608
 609       /* clip (shouldn't we use doubles instead as intermediate format?) */
 610       if (res < -1.0)
 611         res = -1.0;
 612       else if (res > 1.0)
 613         res = 1.0;
 614       tmp[out] = res;
 615     }
 616     memcpy (&out_data[n * outchannels], this->tmp,
 617         sizeof (gdouble) * outchannels);
 618   }
 619 }