summaryrefslogtreecommitdiffstats
path: root/src/ffts_internal.h
diff options
context:
space:
mode:
authorHaruki Hasegawa <h6a.h4i.0@gmail.com>2016-05-05 13:24:55 +0900
committerHaruki Hasegawa <h6a.h4i.0@gmail.com>2016-05-05 13:24:55 +0900
commitda3213cf045e0c7c4971d8b44272d1d86d689ceb (patch)
tree33e17e0a166f03307ebf11e8ab2891ae1ab90f61 /src/ffts_internal.h
parentfa1780c68593762b1e4bdbc46d83912db3eba27a (diff)
parent944d14c9151f6b20145de0cdae38e366e73c9432 (diff)
downloadffts-da3213cf045e0c7c4971d8b44272d1d86d689ceb.zip
ffts-da3213cf045e0c7c4971d8b44272d1d86d689ceb.tar.gz
Merge remote-tracking branch 'linkotec/master'
Diffstat (limited to 'src/ffts_internal.h')
-rw-r--r--src/ffts_internal.h211
1 files changed, 211 insertions, 0 deletions
diff --git a/src/ffts_internal.h b/src/ffts_internal.h
new file mode 100644
index 0000000..7ae8789
--- /dev/null
+++ b/src/ffts_internal.h
@@ -0,0 +1,211 @@
+/*
+
+This file is part of FFTS -- The Fastest Fourier Transform in the South
+
+Copyright (c) 2012, Anthony M. Blake <amb@anthonix.com>
+Copyright (c) 2012, The University of Waikato
+
+All rights reserved.
+
+Redistribution and use in source and binary forms, with or without
+modification, are permitted provided that the following conditions are met:
+* Redistributions of source code must retain the above copyright
+notice, this list of conditions and the following disclaimer.
+* 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.
+* Neither the name of the organization nor the
+names of its contributors may be used to endorse or promote products
+derived from this software without specific prior written permission.
+
+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 ANTHONY M. BLAKE 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.
+
+*/
+
+#ifndef FFTS_INTERNAL_H
+#define FFTS_INTERNAL_H
+
+//#include "config.h"
+#include "ffts_attributes.h"
+#include "types.h"
+
+#ifdef HAVE_MALLOC_H
+#include <malloc.h>
+#endif
+
+#include <stddef.h>
+
+#ifdef HAVE_STDINT_H
+#include <stdint.h>
+#endif
+
+#ifdef HAVE_STDLIB_H
+#include <stdlib.h>
+#endif
+
+#include <stdio.h>
+
+#ifdef __ANDROID__
+#include <android/log.h>
+#define LOG(s) __android_log_print(ANDROID_LOG_ERROR, "FFTS", s)
+#else
+#define LOG(s) fprintf(stderr, s)
+#endif
+
+struct _ffts_plan_t;
+typedef void (*transform_func_t)(struct _ffts_plan_t *p, const void *in, void *out);
+
+/**
+ * Contains all the Information need to perform FFT
+ *
+ *
+ * DO NOT CHANGE THE ORDER OF MEMBERS
+ * ASSEMBLY CODE USES HARD CODED OFFSETS TO REFERENCE
+ * SOME OF THESE VARIABES!!
+ */
+struct _ffts_plan_t {
+
+ /**
+ *
+ */
+ ptrdiff_t *offsets;
+#ifdef DYNAMIC_DISABLED
+ /**
+ * Twiddle factors
+ */
+ void *ws;
+
+ /**
+ * ee - 2 size x size8
+ * oo - 2 x size4 in parallel
+ * oe -
+ */
+ void *oe_ws, *eo_ws, *ee_ws;
+#else
+ void FFTS_ALIGN(32) *ws;
+ void FFTS_ALIGN(32) *oe_ws, *eo_ws, *ee_ws;
+#endif
+
+ /**
+ * Pointer into an array of precomputed indexes for the input data array
+ */
+ ptrdiff_t *is;
+
+ /**
+ * Twiddle Factor Indexes
+ */
+ size_t *ws_is;
+
+ /**
+ * Size of the loops for the base cases
+ */
+ size_t i0, i1, n_luts;
+
+ /**
+ * Size fo the Transform
+ */
+ size_t N;
+ void *lastlut;
+
+#ifdef __arm__
+ size_t *temporary_fix_as_dynamic_code_assumes_fixed_offset;
+#endif
+
+ /**
+ * Pointer to the dynamically generated function
+ * that will execute the FFT
+ */
+ transform_func_t transform;
+
+ /**
+ * Pointer to the base memory address of
+ * of the transform function
+ */
+ void *transform_base;
+
+ /**
+ * Size of the memory block contain the
+ * generated code
+ */
+ size_t transform_size;
+
+ /**
+ * Points to the cosnant variables used by
+ * the Assembly Code
+ */
+ void *constants;
+
+ // multi-dimensional stuff:
+ struct _ffts_plan_t **plans;
+ int rank;
+ size_t *Ns, *Ms;
+ void *buf;
+
+ void *transpose_buf;
+
+ /**
+ * Pointer to the destroy function
+ * to clean up the plan after use
+ * (differs for real and multi dimension transforms
+ */
+ void (*destroy)(struct _ffts_plan_t *);
+
+ /**
+ * Coefficiants for the real valued transforms
+ */
+ float *A, *B;
+
+ size_t i2;
+};
+
+static FFTS_INLINE void *ffts_aligned_malloc(size_t size)
+{
+#if defined(_WIN32)
+ return _aligned_malloc(size, 32);
+#else
+ return valloc(size);
+#endif
+}
+
+static FFTS_INLINE void ffts_aligned_free(void *p)
+{
+#if defined(_WIN32)
+ _aligned_free(p);
+#else
+ free(p);
+#endif
+}
+
+#if GCC_VERSION_AT_LEAST(3,3)
+#define ffts_ctzl __builtin_ctzl
+#elif defined(_MSC_VER)
+#include <intrin.h>
+#ifdef _M_X64
+#pragma intrinsic(_BitScanForward64)
+static __inline unsigned long ffts_ctzl(size_t N)
+{
+ unsigned long count;
+ _BitScanForward64((unsigned long*) &count, N);
+ return count;
+}
+#else
+#pragma intrinsic(_BitScanForward)
+static __inline unsigned long ffts_ctzl(size_t N)
+{
+ unsigned long count;
+ _BitScanForward((unsigned long*) &count, N);
+ return count;
+}
+#endif /* _WIN64 */
+#endif /* _MSC_VER */
+
+#endif /* FFTS_INTERNAL_H */
OpenPOWER on IntegriCloud