No-op; fix GPL address and use the explicit-program-name version.

[dcpomatic.git] / src / lib / audio_filter.cc

diff --git a/src/lib/audio_filter.cc b/src/lib/audio_filter.cc
index 2cf1cf6047ac4b3e00adb5d6812ac8bd1ad5061f..44345fc9d11f4e817865a402266d2b2bedff1b41 100644 (file)
--- a/src/lib/audio_filter.cc
+++ b/src/lib/audio_filter.cc
@@ -1,19 +1,20 @@
 /*
     Copyright (C) 2014 Carl Hetherington <cth@carlh.net>
 
-    This program is free software; you can redistribute it and/or modify
+    This file is part of DCP-o-matic.
+
+    DCP-o-matic is free software; you can redistribute it and/or modify
     it under the terms of the GNU General Public License as published by
     the Free Software Foundation; either version 2 of the License, or
     (at your option) any later version.
 
-    This program is distributed in the hope that it will be useful,
+    DCP-o-matic is distributed in the hope that it will be useful,
     but WITHOUT ANY WARRANTY; without even the implied warranty of
     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
     GNU General Public License for more details.
 
     You should have received a copy of the GNU General Public License
-    along with this program; if not, write to the Free Software
-    Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
+    along with DCP-o-matic.  If not, see <http://www.gnu.org/licenses/>.
 
 */
 
@@ -21,14 +22,14 @@
 #include "audio_buffers.h"
 #include <cmath>
 
-using std::vector;
 using std::min;
 using boost::shared_ptr;
 
-vector<float>
+/** @return array of floats which the caller must destroy with delete[] */
+float *
 AudioFilter::sinc_blackman (float cutoff, bool invert) const
 {
-       vector<float> ir (_M + 1);
+       float* ir = new float[_M + 1];
 
        /* Impulse response */
 
@@ -66,8 +67,13 @@ AudioFilter::sinc_blackman (float cutoff, bool invert) const
        return ir;
 }
 
+AudioFilter::~AudioFilter ()
+{
+       delete[] _ir;
+}
+
 shared_ptr<AudioBuffers>
-AudioFilter::run (shared_ptr<AudioBuffers> in)
+AudioFilter::run (shared_ptr<const AudioBuffers> in)
 {
        shared_ptr<AudioBuffers> out (new AudioBuffers (in->channels(), in->frames()));
 
@@ -76,18 +82,24 @@ AudioFilter::run (shared_ptr<AudioBuffers> in)
                _tail->make_silent ();
        }
 
-       for (int i = 0; i < in->channels(); ++i) {
-               for (int j = 0; j < in->frames(); ++j) {
+       int const channels = in->channels ();
+       int const frames = in->frames ();
+
+       for (int i = 0; i < channels; ++i) {
+               float* tail_p = _tail->data (i);
+               float* in_p = in->data (i);
+               float* out_p = out->data (i);
+               for (int j = 0; j < frames; ++j) {
                        float s = 0;
                        for (int k = 0; k <= _M; ++k) {
                                if ((j - k) < 0) {
-                                       s += _tail->data(i)[j - k + _M + 1] * _ir[k];
+                                       s += tail_p[j - k + _M + 1] * _ir[k];
                                } else {
-                                       s += in->data(i)[j - k] * _ir[k];
+                                       s += in_p[j - k] * _ir[k];
                                }
                        }
 
-                       out->data(i)[j] = s;
+                       out_p[j] = s;
                }
        }
 
@@ -122,14 +134,18 @@ HighPassAudioFilter::HighPassAudioFilter (float transition_bandwidth, float cuto
 BandPassAudioFilter::BandPassAudioFilter (float transition_bandwidth, float lower, float higher)
        : AudioFilter (transition_bandwidth)
 {
-       vector<float> lpf = sinc_blackman (lower, false);
-       vector<float> hpf = sinc_blackman (higher, true);
+       float* lpf = sinc_blackman (lower, false);
+       float* hpf = sinc_blackman (higher, true);
 
-       _ir.resize (_M + 1);
+       delete[] _ir;
+       _ir = new float[_M + 1];
        for (int i = 0; i <= _M; ++i) {
                _ir[i] = lpf[i] + hpf[i];
        }
 
+       delete[] lpf;
+       delete[] hpf;
+
        /* We now have a band-stop, so invert for band-pass */
        for (int i = 0; i <= _M; ++i) {
                _ir[i] = -_ir[i];