CFLAGS = -Ofast -march=native -Wall
LDLIBS = -lm

all:Out.noise.svg Out.timing.svg

Out.noise.svg: out.noise.data Makefile
	echo '\
	set term svg size 640,480 background "white" font "Arial,18";\
	set out "$@";\
	set title "Noise filtering using DFT";\
	set xlabel "time";\
	set ylabel "signal";\
	plot \
	 "$<" index 0 with lines title "original signal"\
	,"$<" index 1 with lines title "signal with random noise"\
	,"$<" index 2 with lines title "amplitude filtered"\
	,"$<" index 3 with lines title "frequency filtered"\
	'|gnuplot

Out.timing.svg: out.timing.data Makefile
	echo '\
	set term svg size 640,480 background "white" font "Arial,18";\
	set out "$@";\
	set title "FFT time as function of data set size";\
	set key left;\
	set xlabel "data set size N";\
	set ylabel "time to perform FFT";\
	a=1;\
	f(x)=a*x*log(x);\
	fit f(x) "$<"  via a;\
	plot "$<" with points title "data", f(x) with lines title "N*log(N)";\
	'|gnuplot

#N = 16384 32768 65536 131072 #262144 segmentation fault
#N = $(shell echo '2^14; 2^15; 2^16; 2^17'|bc -l)
N = $(shell echo '3*2^12; 3*2^13; 3*2^14; 3*2^15'|bc -l)
out.timing.data:main Makefile
	>$@
	for n in $(N); do time -f "$$n %U" -ao $@ ./$< $$n; done

out.noise.data:main
	./$< > $@

main: main.o fft2.o

clean:
	find . -type f -executable -delete
	$(RM) *.o [Oo]ut* *.log

TEST=$(shell echo '2^11; 2^12; 2^13'|bc -l)
test:
	echo $(TEST)