#ifndef __CP_SSE_H__
#define __CP_SSE_H__

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <complex.h>
#include <stddef.h>
#include <stdint.h>
#include <stdalign.h>

#include "codegen.h"

typedef complex float cdata_t;
typedef alignas(16) float data_t;

#define W(N,k) (cexp(-2.0f * M_PI * I * (float)(k) / (float)(N)))

typedef size_t transform_index_t;

typedef void (*transform_func_t)(float *data, size_t N, float *LUT);

struct _ffts_plan_t {
	ptrdiff_t *is;
	ptrdiff_t *offsets;
	void __attribute__ ((aligned(32))) *ws;
	size_t *ws_is;
	void (*firstpass)(const float * restrict, float * restrict, struct _ffts_plan_t * restrict);
	size_t i0, i1, n_luts;
	size_t N;
	void *lastlut;
	transform_index_t *transforms; 
	transform_func_t transform;
};

typedef struct _ffts_plan_t ffts_plan_t;

#endif