libs/ardour/mix.cc

   1 /*
   2     Copyright (C) 2000-2005 Paul Davis
   3
   4     This program is free software; you can redistribute it and/or modify
   5     it under the terms of the GNU General Public License as published by
   6     the Free Software Foundation; either version 2 of the License, or
   7     (at your option) any later version.
   8
   9     This program is distributed in the hope that it will be useful,
  10     but WITHOUT ANY WARRANTY; without even the implied warranty of
  11     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  12     GNU General Public License for more details.
  13
  14     You should have received a copy of the GNU General Public License
  15     along with this program; if not, write to the Free Software
  16     Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
  17
  18     $Id$
  19 */
  20
  21 #include <cmath>
  22 #include <ardour/types.h>
  23 #include <ardour/utils.h>
  24 #include <ardour/mix.h>
  25 #include <stdint.h>
  26
  27 #if defined (ARCH_X86) && defined (BUILD_SSE_OPTIMIZATIONS)
  28
  29 // Debug wrappers
  30
  31 float
  32 debug_compute_peak (ARDOUR::Sample *buf, nframes_t nsamples, float current)
  33 {
  34         if ( ((intptr_t)buf % 16) != 0) {
  35                 cerr << "compute_peak(): buffer unaligned!" << endl;
  36         }
  37
  38         return x86_sse_compute_peak(buf, nsamples, current);
  39 }
  40
  41 void
  42 debug_apply_gain_to_buffer (ARDOUR::Sample *buf, nframes_t nframes, float gain)
  43 {
  44         if ( ((intptr_t)buf % 16) != 0) {
  45                 cerr << "apply_gain_to_buffer(): buffer unaligned!" << endl;
  46         }
  47
  48         x86_sse_apply_gain_to_buffer(buf, nframes, gain);
  49 }
  50
  51 void
  52 debug_mix_buffers_with_gain (ARDOUR::Sample *dst, ARDOUR::Sample *src, nframes_t nframes, float gain)
  53 {
  54         if ( ((intptr_t)dst & 15) != 0) {
  55                 cerr << "mix_buffers_with_gain(): dst unaligned!" << endl;
  56         }
  57
  58         if ( ((intptr_t)dst & 15) != ((intptr_t)src & 15) ) {
  59                 cerr << "mix_buffers_with_gain(): dst & src don't have the same alignment!" << endl;
  60                 mix_buffers_with_gain(dst, src, nframes, gain);
  61         } else {
  62                 x86_sse_mix_buffers_with_gain(dst, src, nframes, gain);
  63         }
  64 }
  65
  66 void
  67 debug_mix_buffers_no_gain (ARDOUR::Sample *dst, ARDOUR::Sample *src, nframes_t nframes)
  68 {
  69         if ( ((intptr_t)dst & 15) != 0) {
  70                 cerr << "mix_buffers_no_gain(): dst unaligned!" << endl;
  71         }
  72
  73         if ( ((intptr_t)dst & 15) != ((intptr_t)src & 15) ) {
  74                 cerr << "mix_buffers_no_gain(): dst & src don't have the same alignment!" << endl;
  75                 mix_buffers_no_gain(dst, src, nframes);
  76         } else {
  77                 x86_sse_mix_buffers_no_gain(dst, src, nframes);
  78         }
  79 }
  80
  81 #endif
  82
  83
  84 float
  85 compute_peak (ARDOUR::Sample *buf, nframes_t nsamples, float current)
  86 {
  87         for (nframes_t i = 0; i < nsamples; ++i) {
  88                 current = f_max (current, fabsf (buf[i]));
  89         }
  90
  91         return current;
  92 }
  93
  94 void
  95 apply_gain_to_buffer (ARDOUR::Sample *buf, nframes_t nframes, float gain)
  96 {
  97         for (nframes_t i=0; i<nframes; i++)
  98                 buf[i] *= gain;
  99 }
 100
 101 void
 102 mix_buffers_with_gain (ARDOUR::Sample *dst, ARDOUR::Sample *src, nframes_t nframes, float gain)
 103 {
 104         for (nframes_t i = 0; i < nframes; i++) {
 105                 dst[i] += src[i] * gain;
 106         }
 107 }
 108
 109 void
 110 mix_buffers_no_gain (ARDOUR::Sample *dst, ARDOUR::Sample *src, nframes_t nframes)
 111 {
 112         for (nframes_t i=0; i < nframes; i++) {
 113                 dst[i] += src[i];
 114         }
 115 }
 116
 117 #if defined (__APPLE__) && defined (BUILD_VECLIB_OPTIMIZATIONS)
 118 #include <Accelerate/Accelerate.h>
 119
 120 float
 121 veclib_compute_peak (ARDOUR::Sample *buf, nframes_t nsamples, float current)
 122 {
 123         float tmpmax = 0.0f;
 124         vDSP_maxmgv(buf, 1, &tmpmax, nsamples);
 125         return f_max(current, tmpmax);
 126 }
 127
 128 void
 129 veclib_apply_gain_to_buffer (ARDOUR::Sample *buf, nframes_t nframes, float gain)
 130 {
 131         vDSP_vsmul(buf, 1, &gain, buf, 1, nframes);
 132 }
 133
 134 void
 135 veclib_mix_buffers_with_gain (ARDOUR::Sample *dst, ARDOUR::Sample *src, nframes_t nframes, float gain)
 136 {
 137         vDSP_vsma(src, 1, &gain, dst, 1, dst, 1, nframes);
 138 }
 139
 140 void
 141 veclib_mix_buffers_no_gain (ARDOUR::Sample *dst, ARDOUR::Sample *src, nframes_t nframes)
 142 {
 143         // It seems that a vector mult only operation does not exist...
 144         float gain = 1.0f;
 145         vDSP_vsma(src, 1, &gain, dst, 1, dst, 1, nframes);
 146 }
 147
 148 #endif
 149
 150