fft_flt.cc

用systemC制作的快速付立叶变换

// fft_flt.cc
// Pierre Bomel - LESTER/UBS - janvier 2005
//    Formation chez France Telecom R&D
#include <systemc.h>
#include "gen.h"
#include "sense.h"
#include "pack.h"
#include "fft_flt.h"

int const points = 64;
int const log2points = 6;
int const sine_waves = 8;
int const sampling_period_ns = 10;
int const nwaves = 1;
int const sampling_length = (nwaves * points + 1) * sampling_period_ns;

int sc_main(int argc, char * argv[])
{
	cout << "Points = " << points << endl;
	cout << "Waves  = " << sine_waves << endl;
	gen<double, points, sine_waves, sampling_period_ns> GEN("gen");
	sense<double, 0> SENSE("sense");
	pack<double, points> PACK("pack");
	fft_flt<points, log2points, false, true> FFT("fft");
	sc_signal<double> sreal, simag;
	sc_signal<double> sreals[points], simags[points];
	sc_signal<double> freals[points], fimags[points];

	GEN.real(sreal);
	GEN.imag(simag);
	SENSE.real(sreal);
	SENSE.imag(simag);

	PACK.real(sreal);
	PACK.imag(simag);

	for(int i = 0; i < points; i++){
		PACK.reals[i](sreals[i]);
		PACK.imags[i](simags[i]);

		FFT.sreals[i](sreals[i]);
		FFT.simags[i](simags[i]);

		FFT.freals[i](freals[i]);
		FFT.fimags[i](fimags[i]);
	}

	sc_start(sampling_length, SC_NS);

	return 0;
}