summaryrefslogtreecommitdiffstats
path: root/src
diff options
context:
space:
mode:
Diffstat (limited to 'src')
-rw-r--r--src/DabMod.cpp12
-rw-r--r--src/OfdmGenerator.cpp104
-rw-r--r--src/OfdmGenerator.h14
-rw-r--r--src/Resampler.cpp164
-rw-r--r--src/Resampler.h19
-rw-r--r--src/Utils.cpp4
-rw-r--r--src/kiss_fftsimd.c44
-rw-r--r--src/kiss_fftsimd.h60
8 files changed, 13 insertions, 408 deletions
diff --git a/src/DabMod.cpp b/src/DabMod.cpp
index 43550c5..a60bdd6 100644
--- a/src/DabMod.cpp
+++ b/src/DabMod.cpp
@@ -291,18 +291,6 @@ int launch_modulator(int argc, char* argv[])
#endif
<< std::endl;
- std::cerr << "Using FFT library " <<
-#if defined(USE_FFTW)
- "FFTW" <<
-#endif
-#if defined(USE_KISS_FFT)
- "Kiss FFT" <<
-#endif
-#if defined(USE_SIMD)
- " (with fft_simd)" <<
-#endif
- "\n";
-
std::cerr << "Compiled with features: " <<
#if defined(HAVE_ZEROMQ)
"zeromq " <<
diff --git a/src/OfdmGenerator.cpp b/src/OfdmGenerator.cpp
index 3dc85a9..55dba0c 100644
--- a/src/OfdmGenerator.cpp
+++ b/src/OfdmGenerator.cpp
@@ -2,7 +2,7 @@
Copyright (C) 2005, 2006, 2007, 2008, 2009, 2010, 2011 Her Majesty
the Queen in Right of Canada (Communications Research Center Canada)
- Copyright (C) 2014
+ Copyright (C) 2016
Matthias P. Braendli, matthias.braendli@mpb.li
http://opendigitalradio.org
@@ -26,12 +26,8 @@
#include "OfdmGenerator.h"
#include "PcDebug.h"
-#if USE_FFTW
-# include "fftw3.h"
-# define FFT_TYPE fftwf_complex
-#else
-# include "kiss_fftsimd.h"
-#endif
+#include "fftw3.h"
+#define FFT_TYPE fftwf_complex
#include <stdio.h>
#include <string.h>
@@ -48,11 +44,7 @@ OfdmGenerator::OfdmGenerator(size_t nbSymbols,
ModCodec(ModFormat(nbSymbols * nbCarriers * sizeof(FFT_TYPE)),
ModFormat(nbSymbols * spacing * sizeof(FFT_TYPE))),
myFftPlan(NULL),
-#if USE_FFTW
myFftIn(NULL), myFftOut(NULL),
-#else
- myFftBuffer(NULL),
-#endif
myNbSymbols(nbSymbols),
myNbCarriers(nbCarriers),
mySpacing(spacing)
@@ -93,7 +85,6 @@ OfdmGenerator::OfdmGenerator(size_t nbSymbols,
PDEBUG(" myZeroDst: %u\n", myZeroDst);
PDEBUG(" myZeroSize: %u\n", myZeroSize);
-#if USE_FFTW
const int N = mySpacing; // The size of the FFT
myFftIn = (FFT_TYPE*)fftwf_malloc(sizeof(FFT_TYPE) * N);
myFftOut = (FFT_TYPE*)fftwf_malloc(sizeof(FFT_TYPE) * N);
@@ -107,11 +98,6 @@ OfdmGenerator::OfdmGenerator(size_t nbSymbols,
throw std::runtime_error(
"OfdmGenerator::process complexf size is not FFT_TYPE size!");
}
-#else
- myFftPlan = kiss_fft_alloc(mySpacing, 1, NULL, NULL);
- myFftBuffer = (FFT_TYPE*)memalign(16, mySpacing * sizeof(FFT_TYPE));
-#endif
-
}
@@ -119,7 +105,6 @@ OfdmGenerator::~OfdmGenerator()
{
PDEBUG("OfdmGenerator::~OfdmGenerator() @ %p\n", this);
-#if USE_FFTW
if (myFftIn) {
fftwf_free(myFftIn);
}
@@ -131,18 +116,6 @@ OfdmGenerator::~OfdmGenerator()
if (myFftPlan) {
fftwf_destroy_plan(myFftPlan);
}
-
-#else
- if (myFftPlan != NULL) {
- kiss_fft_free(myFftPlan);
- }
-
- if (myFftBuffer != NULL) {
- free(myFftBuffer);
- }
-
- kiss_fft_cleanup();
-#endif
}
int OfdmGenerator::process(Buffer* const dataIn, Buffer* dataOut)
@@ -175,8 +148,6 @@ int OfdmGenerator::process(Buffer* const dataIn, Buffer* dataOut)
"OfdmGenerator::process output size not valid!");
}
-#if USE_FFTW
- // No SIMD/no-SIMD distinction, it's too early to optimize anything
for (size_t i = 0; i < myNbSymbols; ++i) {
myFftIn[0][0] = 0;
myFftIn[0][1] = 0;
@@ -194,75 +165,6 @@ int OfdmGenerator::process(Buffer* const dataIn, Buffer* dataOut)
in += myNbCarriers;
out += mySpacing;
}
-#else
-# ifdef USE_SIMD
- for (size_t i = 0, j = 0; i < sizeIn; ) {
- // Pack 4 fft operations
- typedef struct {
- float r[4];
- float i[4];
- } fft_data;
- assert(sizeof(FFT_TYPE) == sizeof(fft_data));
- complexf *cplxIn = (complexf*)in;
- complexf *cplxOut = (complexf*)out;
- fft_data *dataBuffer = (fft_data*)myFftBuffer;
-
- FFT_REAL(myFftBuffer[0]) = _mm_setzero_ps();
- FFT_IMAG(myFftBuffer[0]) = _mm_setzero_ps();
- for (size_t k = 0; k < myZeroSize; ++k) {
- FFT_REAL(myFftBuffer[myZeroDst + k]) = _mm_setzero_ps();
- FFT_IMAG(myFftBuffer[myZeroDst + k]) = _mm_setzero_ps();
- }
- for (int k = 0; k < 4; ++k) {
- if (i < sizeIn) {
- for (size_t l = 0; l < myPosSize; ++l) {
- dataBuffer[myPosDst + l].r[k] = cplxIn[i + myPosSrc + l].real();
- dataBuffer[myPosDst + l].i[k] = cplxIn[i + myPosSrc + l].imag();
- }
- for (size_t l = 0; l < myNegSize; ++l) {
- dataBuffer[myNegDst + l].r[k] = cplxIn[i + myNegSrc + l].real();
- dataBuffer[myNegDst + l].i[k] = cplxIn[i + myNegSrc + l].imag();
- }
- i += myNbCarriers;
- }
- else {
- for (size_t l = 0; l < myNbCarriers; ++l) {
- dataBuffer[l].r[k] = 0.0f;
- dataBuffer[l].i[k] = 0.0f;
- }
- }
- }
-
- kiss_fft(myFftPlan, myFftBuffer, myFftBuffer);
-
- for (int k = 0; k < 4; ++k) {
- if (j < sizeOut) {
- for (size_t l = 0; l < mySpacing; ++l) {
- cplxOut[j + l] = complexf(dataBuffer[l].r[k], dataBuffer[l].i[k]);
- }
- j += mySpacing;
- }
- }
- }
-# else
- for (size_t i = 0; i < myNbSymbols; ++i) {
- FFT_REAL(myFftBuffer[0]) = 0;
- FFT_IMAG(myFftBuffer[0]) = 0;
- bzero(&myFftBuffer[myZeroDst], myZeroSize * sizeof(FFT_TYPE));
- memcpy(&myFftBuffer[myPosDst], &in[myPosSrc],
- myPosSize * sizeof(FFT_TYPE));
- memcpy(&myFftBuffer[myNegDst], &in[myNegSrc],
- myNegSize * sizeof(FFT_TYPE));
-
- kiss_fft(myFftPlan, myFftBuffer, out);
-
- in += myNbCarriers;
- out += mySpacing;
-
- }
-# endif
-#endif
-
return sizeOut;
}
diff --git a/src/OfdmGenerator.h b/src/OfdmGenerator.h
index ec9f14a..061ad91 100644
--- a/src/OfdmGenerator.h
+++ b/src/OfdmGenerator.h
@@ -2,7 +2,7 @@
Copyright (C) 2005, 2006, 2007, 2008, 2009, 2010, 2011 Her Majesty
the Queen in Right of Canada (Communications Research Center Canada)
- Copyright (C) 2014
+ Copyright (C) 2016
Matthias P. Braendli, matthias.braendli@mpb.li
http://opendigitalradio.org
@@ -34,12 +34,7 @@
#include "porting.h"
#include "ModCodec.h"
-#if USE_FFTW
-# include "fftw3.h"
-#else
-# include "kiss_fftsimd.h"
-# include <kiss_fft.h>
-#endif
+#include "fftw3.h"
#include <sys/types.h>
@@ -57,13 +52,8 @@ public:
const char* name() { return "OfdmGenerator"; }
protected:
-#if USE_FFTW
fftwf_plan myFftPlan;
fftwf_complex *myFftIn, *myFftOut;
-#else
- FFT_PLAN myFftPlan;
- FFT_TYPE *myFftBuffer;
-#endif
size_t myNbSymbols;
size_t myNbCarriers;
size_t mySpacing;
diff --git a/src/Resampler.cpp b/src/Resampler.cpp
index f09a58d..38fea33 100644
--- a/src/Resampler.cpp
+++ b/src/Resampler.cpp
@@ -33,10 +33,8 @@
#include <stdexcept>
#include <assert.h>
-#if USE_FFTW
-# define FFT_REAL(x) x[0]
-# define FFT_IMAG(x) x[1]
-#endif
+#define FFT_REAL(x) x[0]
+#define FFT_IMAG(x) x[1]
template<class T>
T gcd(T a, T b)
@@ -64,44 +62,6 @@ Resampler::Resampler(size_t inputRate, size_t outputRate, size_t resolution) :
{
PDEBUG("Resampler::Resampler(%zu, %zu) @ %p\n", inputRate, outputRate, this);
-#if USE_FFTW
- fprintf(stderr, "This software uses the FFTW library.\n\n");
-#else
- fprintf(stderr, "This software uses KISS FFT.\n\n");
- fprintf(stderr, "Copyright (c) 2003-2004 Mark Borgerding\n"
- "\n"
- "All rights reserved.\n"
- "\n"
- "Redistribution and use in source and binary forms, with or "
- "without modification, are permitted provided that the following "
- "conditions are met:\n"
- "\n"
- " * Redistributions of source code must retain the above "
- "copyright notice, this list of conditions and the following "
- "disclaimer.\n"
- " * Redistributions in binary form must reproduce the above "
- "copyright notice, this list of conditions and the following "
- "disclaimer in the documentation and/or other materials provided "
- "with the distribution.\n"
- " * Neither the author nor the names of any contributors may be "
- "used to endorse or promote products derived from this software "
- "without specific prior written permission.\n"
- "\n"
- "THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND "
- "CONTRIBUTORS \"AS IS\" AND ANY EXPRESS OR IMPLIED WARRANTIES, "
- "INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF "
- "MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE "
- "DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS "
- "BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, "
- "EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED "
- "TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, "
- "DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON "
- "ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, "
- "OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY "
- "OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE "
- "POSSIBILITY OF SUCH DAMAGE.\n");
-#endif
-
size_t divisor = gcd(inputRate, outputRate);
L = outputRate / divisor;
M = inputRate / divisor;
@@ -128,7 +88,6 @@ Resampler::Resampler(size_t inputRate, size_t outputRate, size_t resolution) :
PDEBUG("Window[%zu] = %f\n", i, myWindow[i]);
}
-#if USE_FFTW
myFftIn = (FFT_TYPE*)fftwf_malloc(sizeof(FFT_TYPE) * myFftSizeIn);
myFront = (FFT_TYPE*)fftwf_malloc(sizeof(FFT_TYPE) * myFftSizeIn);
myFftPlan1 = fftwf_plan_dft_1d(myFftSizeIn,
@@ -143,16 +102,6 @@ Resampler::Resampler(size_t inputRate, size_t outputRate, size_t resolution) :
myBufferIn = (complexf*)fftwf_malloc(sizeof(FFT_TYPE) * myFftSizeIn / 2);
myBufferOut = (complexf*)fftwf_malloc(sizeof(FFT_TYPE) * myFftSizeOut / 2);
-#else
- myFftIn = (FFT_TYPE*)memalign(16, myFftSizeIn * sizeof(FFT_TYPE));
- myFftOut = (FFT_TYPE*)memalign(16, myFftSizeOut * sizeof(FFT_TYPE));
- myBufferIn = (complexf*)memalign(16, myFftSizeIn / 2 * sizeof(FFT_TYPE));
- myBufferOut = (complexf*)memalign(16, myFftSizeOut / 2 * sizeof(FFT_TYPE));
- myFront = (FFT_TYPE*)memalign(16, myFftSizeIn * sizeof(FFT_TYPE));
- myBack = (FFT_TYPE*)memalign(16, myFftSizeOut * sizeof(FFT_TYPE));
- myFftPlan1 = kiss_fft_alloc(myFftSizeIn, 0, NULL, NULL);
- myFftPlan2 = kiss_fft_alloc(myFftSizeOut, 1, NULL, NULL);
-#endif
memset(myBufferIn, 0, myFftSizeIn / 2 * sizeof(FFT_TYPE));
memset(myBufferOut, 0, myFftSizeOut / 2 * sizeof(FFT_TYPE));
@@ -163,7 +112,6 @@ Resampler::~Resampler()
{
PDEBUG("Resampler::~Resampler() @ %p\n", this);
-#if USE_FFTW
if (myFftPlan1 != NULL) { fftwf_free(myFftPlan1); }
if (myFftPlan2 != NULL) { fftwf_free(myFftPlan2); }
if (myFftIn != NULL) { fftwf_free(myFftIn); }
@@ -175,18 +123,6 @@ Resampler::~Resampler()
if (myWindow != NULL) { fftwf_free(myWindow); }
fftwf_destroy_plan(myFftPlan1);
fftwf_destroy_plan(myFftPlan2);
-#else
- if (myFftPlan1 != NULL) { free(myFftPlan1); }
- if (myFftPlan2 != NULL) { free(myFftPlan2); }
- if (myFftIn != NULL) { free(myFftIn); }
- if (myFftOut != NULL) { free(myFftOut); }
- if (myBufferIn != NULL) { free(myBufferIn); }
- if (myBufferOut != NULL) { free(myBufferOut); }
- if (myFront != NULL) { free(myFront); }
- if (myBack != NULL) { free(myBack); }
- if (myWindow != NULL) { free(myWindow); }
- kiss_fft_cleanup();
-#endif
}
@@ -201,93 +137,6 @@ int Resampler::process(Buffer* const dataIn, Buffer* dataOut)
FFT_TYPE* out = reinterpret_cast<FFT_TYPE*>(dataOut->getData());
size_t sizeIn = dataIn->getLength() / sizeof(complexf);
-#if defined(USE_SIMD) && !USE_FFTW
- size_t sizeOut = dataOut->getLength() / sizeof(complexf);
-
- typedef struct {
- float r[4];
- float i[4];
- } fft_data;
- assert(sizeof(FFT_TYPE) == sizeof(fft_data));
- fft_data *fftDataIn = (fft_data*)myFftIn;
- fft_data *fftDataOut = (fft_data*)myFftOut;
- complexf *cplxIn = (complexf*)in;
- complexf *cplxOut = (complexf*)out;
- for (size_t i = 0, j = 0; i < sizeIn; ) {
- for (int k = 0; k < 4; ++k) {
- if (i < sizeIn) {
- for (size_t l = 0; l < myFftSizeIn / 2; ++l) {
- fftDataIn[l].r[k] = myBufferIn[l].real();
- fftDataIn[l].i[k] = myBufferIn[l].imag();
- fftDataIn[myFftSizeIn / 2 + l].r[k] = cplxIn[i + l].real();
- fftDataIn[myFftSizeIn / 2 + l].i[k] = cplxIn[i + l].imag();
- }
- memcpy(myBufferIn, cplxIn + i, myFftSizeIn / 2 * sizeof(complexf));
- i += myFftSizeIn / 2;
- } else {
- for (size_t l = 0; l < myFftSizeIn; ++l) {
- fftDataIn[l].r[k] = 0.0f;
- fftDataIn[l].i[k] = 0.0f;
- }
- }
- }
- for (size_t k = 0; k < myFftSizeIn; ++ k) {
- FFT_REAL(myFftIn[k]) = _mm_mul_ps(FFT_REAL(myFftIn[k]), _mm_set_ps1(myWindow[k]));
- FFT_IMAG(myFftIn[k]) = _mm_mul_ps(FFT_IMAG(myFftIn[k]), _mm_set_ps1(myWindow[k]));
- }
-
- kiss_fft(myFftPlan1, myFftIn, myFront);
-
- if (myFftSizeOut > myFftSizeIn) {
- memset(myBack, 0, myFftSizeOut * sizeof(FFT_TYPE));
- memcpy(myBack, myFront, myFftSizeIn / 2 * sizeof(FFT_TYPE));
- memcpy(&myBack[myFftSizeOut - (myFftSizeIn / 2)],
- &myFront[myFftSizeIn / 2],
- myFftSizeIn / 2 * sizeof(FFT_TYPE));
- // Copy input Fs
- FFT_REAL(myBack[myFftSizeIn / 2]) =
- FFT_REAL(myFront[myFftSizeIn / 2]);
- FFT_IMAG(myBack[myFftSizeIn / 2]) =
- FFT_IMAG(myFront[myFftSizeIn / 2]);
- } else {
- memcpy(myBack, myFront, myFftSizeOut / 2 * sizeof(FFT_TYPE));
- memcpy(&myBack[myFftSizeOut / 2],
- &myFront[myFftSizeIn - (myFftSizeOut / 2)],
- myFftSizeOut / 2 * sizeof(FFT_TYPE));
- // Average output Fs from input
- FFT_REAL(myBack[myFftSizeOut / 2]) =
- _mm_add_ps(FFT_REAL(myBack[myFftSizeOut / 2]),
- FFT_REAL(myFront[myFftSizeOut / 2]));
- FFT_IMAG(myBack[myFftSizeOut / 2]) =
- _mm_add_ps(FFT_IMAG(myBack[myFftSizeOut / 2]),
- FFT_IMAG(myFront[myFftSizeOut / 2]));
- FFT_REAL(myBack[myFftSizeOut / 2]) =
- _mm_mul_ps(FFT_REAL(myBack[myFftSizeOut / 2]), _mm_set_ps1(0.5f));
- FFT_IMAG(myBack[myFftSizeOut / 2]) =
- _mm_mul_ps(FFT_IMAG(myBack[myFftSizeOut / 2]), _mm_set_ps1(0.5f));
- }
- for (size_t k = 0; k < myFftSizeOut; ++k) {
- FFT_REAL(myBack[k]) = _mm_mul_ps(FFT_REAL(myBack[k]), _mm_set_ps1(myFactor));
- FFT_IMAG(myBack[k]) = _mm_mul_ps(FFT_IMAG(myBack[k]), _mm_set_ps1(myFactor));
- }
-
- kiss_fft(myFftPlan2, myBack, myFftOut);
-
- for (size_t k = 0; k < 4; ++k) {
- if (j < sizeOut) {
- for (size_t l = 0; l < myFftSizeOut / 2; ++l) {
- cplxOut[j + l] = complexf(myBufferOut[l].real() + fftDataOut[l].r[k],
- myBufferOut[l].imag() + fftDataOut[l].i[k]);
- myBufferOut[l] = complexf(fftDataOut[myFftSizeOut / 2 + l].r[k],
- fftDataOut[myFftSizeOut / 2 + l].i[k]);
- }
- }
- j += myFftSizeOut / 2;
- }
- }
-#endif
-
-#if USE_FFTW || (!defined(USE_SIMD))
for (size_t i = 0, j = 0; i < sizeIn; i += myFftSizeIn / 2, j += myFftSizeOut / 2) {
memcpy(myFftIn, myBufferIn, myFftSizeIn / 2 * sizeof(FFT_TYPE));
memcpy(myFftIn + (myFftSizeIn / 2), in + i, myFftSizeIn / 2 * sizeof(FFT_TYPE));
@@ -297,11 +146,7 @@ int Resampler::process(Buffer* const dataIn, Buffer* dataOut)
FFT_IMAG(myFftIn[k]) *= myWindow[k];
}
-#if USE_FFTW
fftwf_execute(myFftPlan1);
-#else
- kiss_fft(myFftPlan1, myFftIn, myFront);
-#endif
if (myFftSizeOut > myFftSizeIn) {
memset(myBack, 0, myFftSizeOut * sizeof(FFT_TYPE));
@@ -332,11 +177,7 @@ int Resampler::process(Buffer* const dataIn, Buffer* dataOut)
FFT_IMAG(myBack[k]) *= myFactor;
}
-#if USE_FFTW
fftwf_execute(myFftPlan2);
-#else
- kiss_fft(myFftPlan2, myBack, myFftOut);
-#endif
for (size_t k = 0; k < myFftSizeOut / 2; ++k) {
FFT_REAL(out[j + k]) = myBufferOut[k].real() + FFT_REAL(myFftOut[k]);
@@ -344,7 +185,6 @@ int Resampler::process(Buffer* const dataIn, Buffer* dataOut)
}
memcpy(myBufferOut, myFftOut + (myFftSizeOut / 2), (myFftSizeOut / 2) * sizeof(FFT_TYPE));
}
-#endif
return 1;
}
diff --git a/src/Resampler.h b/src/Resampler.h
index 392d0a6..f8388ee 100644
--- a/src/Resampler.h
+++ b/src/Resampler.h
@@ -33,19 +33,12 @@
#include "porting.h"
#include "ModCodec.h"
-#if USE_FFTW
-# include <sys/types.h>
-# include <fftw3.h>
-
-# define FFT_TYPE fftwf_complex
-# define FFT_PLAN fftwf_plan
-
-#else
-# include "kiss_fftsimd.h"
-# include <sys/types.h>
-# include <kiss_fft.h>
-# include <tools/kiss_fftr.h>
-#endif
+#include <sys/types.h>
+#include <fftw3.h>
+
+#define FFT_TYPE fftwf_complex
+#define FFT_PLAN fftwf_plan
+
#include <complex>
typedef std::complex<float> complexf;
diff --git a/src/Utils.cpp b/src/Utils.cpp
index 1560e07..4805374 100644
--- a/src/Utils.cpp
+++ b/src/Utils.cpp
@@ -110,10 +110,6 @@ void printVersion(void)
" You should have received a copy of the GNU General Public License along\n"
" with ODR-DabMod. If not, see <http://www.gnu.org/licenses/>.\n"
"\n"
-#if USE_KISS_FFT
- "ODR-DabMod makes use of the following open source packages:\n"
- " Kiss FFT v1.2.9 (Revised BSD) - http://kissfft.sourceforge.net/\n"
-#endif
);
}
diff --git a/src/kiss_fftsimd.c b/src/kiss_fftsimd.c
deleted file mode 100644
index 83f2e4a..0000000
--- a/src/kiss_fftsimd.c
+++ /dev/null
@@ -1,44 +0,0 @@
-/*
- Copyright (C) 2011, 2012
- Her Majesty the Queen in Right of Canada (Communications Research
- Center Canada)
- */
-/*
- This file is part of ODR-DabMod.
-
- ODR-DabMod 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 3 of the
- License, or (at your option) any later version.
-
- ODR-DabMod 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 ODR-DabMod. If not, see <http://www.gnu.org/licenses/>.
- */
-
-#include "kiss_fftsimd.h"
-
-
-#ifdef USE_SIMD
-void kiss_fft_pack(kiss_fft_complex *in, size_t in_size, size_t in_offset,
- kiss_fft_complex *out, size_t out_size, size_t out_offset,
- size_t stride, size_t n)
-{
- size_t i, j;
- complex_float *in_cplx = (complex_float*)in;
-
- for (i = 0; i < 4; ++i) {
- if (in_offset < in_size) {
- for (j = 0; j < n; ++j) {
- out[out_offset + j].r[i] = in_cplx[in_offset + j].r;
- out[out_offset + j].i[i] = in_cplx[in_offset + j].i;
- }
- in_offset += stride;
- }
- }
-}
-#endif
diff --git a/src/kiss_fftsimd.h b/src/kiss_fftsimd.h
deleted file mode 100644
index 90ee435..0000000
--- a/src/kiss_fftsimd.h
+++ /dev/null
@@ -1,60 +0,0 @@
-/*
- Copyright (C) 2011, 2012
- Her Majesty the Queen in Right of Canada (Communications Research
- Center Canada)
- */
-/*
- This file is part of ODR-DabMod.
-
- ODR-DabMod 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 3 of the
- License, or (at your option) any later version.
-
- ODR-DabMod 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 ODR-DabMod. If not, see <http://www.gnu.org/licenses/>.
- */
-
-#ifndef KISS_FFTSIMD_H
-#define KISS_FFTSIMD_H
-
-#ifdef HAVE_CONFIG_H
-# include "config.h"
-#endif
-
-#include <kiss_fft.h>
-
-
-#define FFT_TYPE kiss_fft_cpx
-#define FFT_PLAN kiss_fft_cfg
-#define FFT_REAL(a) (a).r
-#define FFT_IMAG(a) (a).i
-
-
-
-#ifdef USE_SIMD
-
-typedef struct {
- float r[4];
- float i[4];
-} kiss_fft_complex;
-
-typedef struct {
- float r;
- float i;
-} complex_float;
-
-void kiss_fft_pack(kiss_fft_complex *in, size_t in_size, size_t in_offset,
- kiss_fft_complex *out, size_t out_size, size_t out_offset,
- size_t stride, size_t n);
-void kiss_fft_unpack(kiss_fft_complex *dst, kiss_fft_complex *src, size_t n, size_t offset, size_t stride);
-
-#endif
-
-
-#endif // KISS_FFTSIMD_H
generated by cgit v1.2.3 (git 2.39.1) at 2024-11-11 23:16:45 +0000