#!/bin/bash
rm *.txt
rm *.o
rm *.ps
rm ./flutter_optimizer
make


#UNIFORM
./flutter_optimizer 1 1 52 5 code.txt TF_code.txt rendement.txt mean_sdt.txt


gnuplot <<EOF
# set xrange[-1.5:1.5]
set yrange[0.4:1.2]
set lmargin 5
set rmargin 1
set xlabel 'Velocity v'
set ylabel 'Gain factor'
set terminal postscript eps color
set key on  below box title 'Legend'
set key left Left reverse
set output 'rendement.ps
plot "rendement.txt"  using 1:2 title 'Gain factor G(v)' with lines linewidth 5,\
"rendement.txt"  using 1:3 title 'Average gain factor' with lines linewidth 5,\
"rendement.txt"  using 1:4 title 'Velocity Distribution' with dots linewidth 5,\
"rendement.txt"  every 2  using 1:5 title '"Gain=1 line". Above this line: the flutter beats the snapshot in terms of RMSE.' with points linewidth 5 lt 0
EOF


gnuplot <<EOF
set lmargin 5
set rmargin 1
set xlabel 'k'
set ylabel 'Gain'
set sample 1000000
set terminal postscript eps color

set output 'code.ps'
set key on  below box title 'Legend'
set key left Left reverse
plot "code.txt"  with steps title 'Flutter shutter code' linewidth 5
EOF



gnuplot <<EOF
set lmargin 5
set rmargin 1
set xlabel 'xi'
set ylabel 'Fourier tranforms (modulus)'
set sample 1000000
set terminal postscript eps color

set output 'TF_code.ps'
set key on  below box title 'Legend'
set key left Left reverse
plot "TF_code.txt" using 1:3  with lines title "Fourier tranform (modulus) of the optimized flutter shutter function",\
"TF_code.txt" using 1:2  with lines title "Ideal Fourier tranform (modulus)"  linewidth 5
EOF

mv code.txt code_u5.txt
mv TF_code.txt TF_code_u5.txt
mv rendement.txt rendement_u5.txt 
mv code.ps code_u5.ps
mv rendement.ps rendement_u5.ps
mv TF_code.ps TF_code_u5.ps
mv mean_sdt.txt mean_sdt_u5.txt



./flutter_optimizer 1 1 52 10 code.txt TF_code.txt rendement.txt mean_sdt.txt


gnuplot <<EOF
# set xrange[-1.5:1.5]
set yrange[0.4:1.3]
set lmargin 5
set rmargin 1
set xlabel 'Velocity v'
set ylabel 'Gain factor'
set terminal postscript eps color
set key on  below box title 'Legend'
set key left Left reverse
set output 'rendement.ps
plot "rendement.txt"  using 1:2 title 'Gain factor G(v)' with lines linewidth 5,\
"rendement.txt"  using 1:3 title 'Average gain factor' with lines linewidth 5,\
"rendement.txt"  using 1:4 title 'Velocity Distribution' with dots linewidth 5,\
"rendement.txt"  every 2  using 1:5 title '"Gain=1 line". Above this line: the flutter beats the snapshot in terms of RMSE.' with points linewidth 5 lt 0
EOF


gnuplot <<EOF
set lmargin 5
set rmargin 1
set lmargin 5
set rmargin 1
set xlabel 'k'
set ylabel 'Gain'
set sample 1000000
set terminal postscript eps color

set output 'code.ps'
set key on  below box title 'Legend'
set key left Left reverse
plot "code.txt"  with steps title 'Flutter shutter code' linewidth 5
EOF



gnuplot <<EOF
set lmargin 5
set rmargin 1
set xlabel 'xi'
set ylabel 'Fourier tranforms (modulus)'
set sample 1000000
set terminal postscript eps color

set output 'TF_code.ps'
set key on  below box title 'Legend'
set key left Left reverse
plot "TF_code.txt" using 1:3  with lines title "Fourier tranform (modulus) of the optimized flutter shutter function",\
"TF_code.txt" using 1:2  with lines title "Ideal Fourier tranform (modulus)"  linewidth 5
EOF

mv code.txt code_u10.txt
mv TF_code.txt TF_code_u10.txt
mv rendement.txt rendement_u10.txt 
mv code.ps code_u10.ps
mv rendement.ps rendement_u10.ps
mv TF_code.ps TF_code_u10.ps
mv mean_sdt.txt mean_sdtu_10.txt






ps2eps -f -B -l *.ps