/* GStreamer
 * (c) 2011 Stefan Kost <ensonic@users.sf.net>
 *
 * This library is free software; you can redistribute it and/or
 * modify it under the terms of the GNU Library General Public
 * License as published by the Free Software Foundation; either
 * version 2 of the License, or (at your option) any later version.
 *
 * This library 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
 * Library General Public License for more details.
 *
 * You should have received a copy of the GNU Library General Public
 * License along with this library; if not, write to the
 * Free Software Foundation, Inc., 59 Temple Place - Suite 330,
 * Boston, MA 02111-1307, USA.
 */

#include <stdio.h>
#include <gst/fft/gstffts16.h>
#include <gst/fft/gstffts32.h>
#include <gst/fft/gstfftf32.h>
#include <gst/fft/gstfftf64.h>

/* effectively max range seems to be 1/4 of what it should be */

#define MAKE_I_TEST(_g_,_G_,_t_,_T_,_f_)                        \
static void                                                     \
test_##_t_ (const gchar *test_name, gint num_freq, gint window) \
{                                                               \
  GstFFT ##_T_ *ctx;                                            \
  GstFFT ##_T_ ##Complex *fdata;                                \
  _g_ *adata;                                                   \
  _g_ maxfr = 0, maxfi = 0;                                     \
  gint num_samples = num_freq * 2 - 2;                          \
  gint s, f;                                                    \
                                                                \
  ctx = gst_fft_ ##_t_ ##_new (num_samples, FALSE);             \
  fdata = g_new (GstFFT ##_T_ ##Complex, num_freq);             \
  adata = g_new (_g_, num_samples);                             \
                                                                \
  for (s = 0; s < num_samples;) {                               \
    adata[s++]=G_MIN##_G_;                                      \
    adata[s++]=G_MAX##_G_;                                      \
  }                                                             \
                                                                \
  gst_fft_ ##_t_ ##_window (ctx, adata, window);                \
  gst_fft_ ##_t_ ##_fft (ctx, adata, fdata);                    \
                                                                \
  for (f = 0; f < num_freq; f++) {                              \
    if (fdata[1+f].r > maxfr)                                   \
      maxfr = fdata[1+f].r;                                     \
    if (fdata[1+f].i > maxfi)                                   \
      maxfi = fdata[1+f].i;                                     \
  }                                                             \
                                                                \
  printf (#_t_" %-15s: maxfr: %"_f_" %10.5f maxfi: %"_f_" %10.5f\n",\
    test_name,                                                  \
    maxfr, (gfloat)G_MAX##_G_/maxfr,                            \
    maxfi, (gfloat)G_MAX##_G_/maxfi);                           \
                                                                \
  gst_fft_ ##_t_ ##_free (ctx);                                 \
  g_free (fdata);                                               \
  g_free (adata);                                               \
}

MAKE_I_TEST (gint16, INT16, s16, S16, "6d");
MAKE_I_TEST (gint32, INT32, s32, S32, "9d");

#define MAKE_F_TEST(_g_,_G_,_t_,_T_,_f_)                        \
static void                                                     \
test_##_t_ (const gchar *test_name, gint num_freq, gint window) \
{                                                               \
  GstFFT ##_T_ *ctx;                                            \
  GstFFT ##_T_ ##Complex *fdata;                                \
  _g_ *adata;                                                   \
  _g_ maxfr = 0, maxfi = 0;                                     \
  gint num_samples = num_freq * 2 - 2;                          \
  gint s, f;                                                    \
                                                                \
  ctx = gst_fft_ ##_t_ ##_new (num_samples, FALSE);             \
  fdata = g_new (GstFFT ##_T_ ##Complex, num_freq);             \
  adata = g_new (_g_, num_samples);                             \
                                                                \
  for (s = 0; s < num_samples;) {                               \
    adata[s++]=-1.0;                                            \
    adata[s++]=+1.0;                                            \
  }                                                             \
                                                                \
  gst_fft_ ##_t_ ##_window (ctx, adata, window);                \
  gst_fft_ ##_t_ ##_fft (ctx, adata, fdata);                    \
                                                                \
  for (f = 0; f < num_freq; f++) {                              \
    if (fdata[1+f].r > maxfr)                                   \
      maxfr = fdata[1+f].r;                                     \
    if (fdata[1+f].i > maxfi)                                   \
      maxfi = fdata[1+f].i;                                     \
  }                                                             \
                                                                \
  printf (#_t_" %-15s: maxfr: %"_f_" %10.5f maxfi: %"_f_" %10.5f\n",\
    test_name,                                                  \
    maxfr, (gfloat)1.0/maxfr,                                   \
    maxfi, (gfloat)1.0/maxfi);                                  \
                                                                \
  gst_fft_ ##_t_ ##_free (ctx);                                 \
  g_free (fdata);                                               \
  g_free (adata);                                               \
}

MAKE_F_TEST (gfloat, FLOAT, f32, F32, "10.5f");
MAKE_F_TEST (gdouble, DOUBLE, f64, F64, "10.5f");

gint
main (gint argc, gchar * argv[])
{
  gint num_bands;

  gst_init (&argc, &argv);

  num_bands = 200;
  test_s16 ("200, none", num_bands, GST_FFT_WINDOW_RECTANGULAR);
  test_s16 ("200, hamming", num_bands, GST_FFT_WINDOW_HAMMING);
  test_s16 ("200, hann", num_bands, GST_FFT_WINDOW_HANN);
  test_s16 ("200, bartlett", num_bands, GST_FFT_WINDOW_BARTLETT);
  test_s16 ("200, blackman", num_bands, GST_FFT_WINDOW_BLACKMAN);
  puts ("");

  num_bands = 300;
  test_s16 ("300, none", num_bands, GST_FFT_WINDOW_RECTANGULAR);
  test_s16 ("300, hamming", num_bands, GST_FFT_WINDOW_HAMMING);
  test_s16 ("300, hann", num_bands, GST_FFT_WINDOW_HANN);
  test_s16 ("300, bartlett", num_bands, GST_FFT_WINDOW_BARTLETT);
  test_s16 ("300, blackman", num_bands, GST_FFT_WINDOW_BLACKMAN);
  puts ("\n");

  num_bands = 200;
  test_s32 ("200, none", num_bands, GST_FFT_WINDOW_RECTANGULAR);
  test_s32 ("200, hamming", num_bands, GST_FFT_WINDOW_HAMMING);
  test_s32 ("200, hann", num_bands, GST_FFT_WINDOW_HANN);
  test_s32 ("200, bartlett", num_bands, GST_FFT_WINDOW_BARTLETT);
  test_s32 ("200, blackman", num_bands, GST_FFT_WINDOW_BLACKMAN);
  puts ("");

  num_bands = 300;
  test_s32 ("300, none", num_bands, GST_FFT_WINDOW_RECTANGULAR);
  test_s32 ("300, hamming", num_bands, GST_FFT_WINDOW_HAMMING);
  test_s32 ("300, hann", num_bands, GST_FFT_WINDOW_HANN);
  test_s32 ("300, bartlett", num_bands, GST_FFT_WINDOW_BARTLETT);
  test_s32 ("300, blackman", num_bands, GST_FFT_WINDOW_BLACKMAN);
  puts ("\n");

  num_bands = 200;
  test_f32 ("200, none", num_bands, GST_FFT_WINDOW_RECTANGULAR);
  test_f32 ("200, hamming", num_bands, GST_FFT_WINDOW_HAMMING);
  test_f32 ("200, hann", num_bands, GST_FFT_WINDOW_HANN);
  test_f32 ("200, bartlett", num_bands, GST_FFT_WINDOW_BARTLETT);
  test_f32 ("200, blackman", num_bands, GST_FFT_WINDOW_BLACKMAN);
  puts ("");

  num_bands = 300;
  test_f32 ("300, none", num_bands, GST_FFT_WINDOW_RECTANGULAR);
  test_f32 ("300, hamming", num_bands, GST_FFT_WINDOW_HAMMING);
  test_f32 ("300, hann", num_bands, GST_FFT_WINDOW_HANN);
  test_f32 ("300, bartlett", num_bands, GST_FFT_WINDOW_BARTLETT);
  test_f32 ("300, blackman", num_bands, GST_FFT_WINDOW_BLACKMAN);
  puts ("\n");

  num_bands = 200;
  test_f64 ("200, none", num_bands, GST_FFT_WINDOW_RECTANGULAR);
  test_f64 ("200, hamming", num_bands, GST_FFT_WINDOW_HAMMING);
  test_f64 ("200, hann", num_bands, GST_FFT_WINDOW_HANN);
  test_f64 ("200, bartlett", num_bands, GST_FFT_WINDOW_BARTLETT);
  test_f64 ("200, blackman", num_bands, GST_FFT_WINDOW_BLACKMAN);
  puts ("");

  num_bands = 300;
  test_f64 ("300, none", num_bands, GST_FFT_WINDOW_RECTANGULAR);
  test_f64 ("300, hamming", num_bands, GST_FFT_WINDOW_HAMMING);
  test_f64 ("300, hann", num_bands, GST_FFT_WINDOW_HANN);
  test_f64 ("300, bartlett", num_bands, GST_FFT_WINDOW_BARTLETT);
  test_f64 ("300, blackman", num_bands, GST_FFT_WINDOW_BLACKMAN);
  puts ("\n");
}
