source: trunk/examples/nardi/shallow_performanceMeasurement/scriptPerformance.sh @ 415

#!/bin/bash
############################## Script which provides the results for the performances on the Shallow water example ######################

#svn checkout svn+ssh://lnalod@forge.ipsl.jussieu.fr/ipsl/forge/projets/yao/svn/trunk/examples/nardi/shallow_performanceMeasurement /tmp/shallow_perfomanceMeasurement
#cd /tmp/shallow_perfomanceMeasurement


##################################### Result for space dimension 50x50 ############################################################
YaoI +O3 shalw_50x50
# Copy the file which allows to make the average of the elapsed times of forward and backward.
cp -f 2shalw_50x50.h Yworkdir/Y2shalw_50x50.h 
YaoI +O3 shalw_50x50
# Serial execution
echo num threads: serial...
time ./shalw_50x50 > resultNotParallelFile_50x50.txt

YaoI -p +O3 shalw_50x50
# 1 thread execution
echo num threads: 1
OMP_NUM_THREADS=1 time ./shalw_50x50 > resultParallelFile_1thread_50x50.txt

# 2 thread execution
echo num threads: 2
OMP_NUM_THREADS=2 time ./shalw_50x50 > resultParallelFile_2thread_50x50.txt

# 4 thread execution
echo num threads: 4
OMP_NUM_THREADS=4 time ./shalw_50x50 > resultParallelFile_4thread_50x50.txt

# 8 thread execution
echo num threads: 8
OMP_NUM_THREADS=8 time ./shalw_50x50 > resultParallelFile_8thread_50x50.txt


##################################### Result for space dimension 100x100 ############################################################
YaoI +O3 shalw_100x100
# Copy the file which allows to make the average of the elapsed times of forward and backward.
cp -f 2shalw_100x100.h Yworkdir/Y2shalw_100x100.h 
YaoI +O3 shalw_100x100
# Serial execution
echo num threads: serial...
time ./shalw_100x100 > resultNotParallelFile_100x100.txt

YaoI -p +O3 shalw_100x100
# 1 thread execution
echo num threads: 1
OMP_NUM_THREADS=1 time ./shalw_100x100 > resultParallelFile_1thread_100x100.txt

# 2 thread execution
echo num threads: 2
OMP_NUM_THREADS=2 time ./shalw_100x100 > resultParallelFile_2thread_100x100.txt

# 4 thread execution
echo num threads: 4
OMP_NUM_THREADS=4 time ./shalw_100x100 > resultParallelFile_4thread_100x100.txt

# 8 thread execution
echo num threads: 8
OMP_NUM_THREADS=8 time ./shalw_100x100 > resultParallelFile_8thread_100x100.txt


##################################### Result for space dimension 200x200 ############################################################
YaoI +O3 shalw_200x200
# Copy the file which allows to make the average of the elapsed times of forward and backward.
cp -f 2shalw_200x200.h Yworkdir/Y2shalw_200x200.h 
YaoI +O3 shalw_200x200
# Serial execution
echo num threads: serial...
time ./shalw_200x200 > resultNotParallelFile_200x200.txt

YaoI -p +O3 shalw_200x200
# 1 thread execution
echo num threads: 1
OMP_NUM_THREADS=1 time ./shalw_200x200 > resultParallelFile_1thread_200x200.txt

# 2 thread execution
echo num threads: 2
OMP_NUM_THREADS=2 time ./shalw_200x200 > resultParallelFile_2thread_200x200.txt

# 4 thread execution
echo num threads: 4
OMP_NUM_THREADS=4 time ./shalw_200x200 > resultParallelFile_4thread_200x200.txt

# 8 thread execution
echo num threads: 8
OMP_NUM_THREADS=8 time ./shalw_200x200 > resultParallelFile_8thread_200x200.txt



###################################################################################################################################
###################################################################################################################################

##### Speedups
# In the Shallow all the code is parallel (we are measuring only the nested loops). We do not need to compute de Amdahl speedup.

# P is the number of cores. P is serial. P4 is 4 cores.
P=1
P1=1
P2=2
P4=4
P8=8

# E elapsed time of the serial version. E1 for the 1 thread version. E2 for the 2 threads version.
##################################### space 50x50 ############################################################
E_50x50=`awk '/counter 2300/ {print $4}' resultNotParallelFile_50x50.txt`
E1_50x50=`awk '/counter 2300/ {print $4}' resultParallelFile_1thread_50x50.txt`
E2_50x50=`awk '/counter 2300/ {print $4}' resultParallelFile_2thread_50x50.txt`
E4_50x50=`awk '/counter 2300/ {print $4}' resultParallelFile_4thread_50x50.txt`
E8_50x50=`awk '/counter 2300/ {print $4}' resultParallelFile_8thread_50x50.txt`
echo Elapsed time 50x50 serial version
echo $E_50x50
echo Elapsed time 1 thread
echo $E1_50x50
echo Elapsed time 2 threads
echo $E2_50x50
echo Elapsed time 4 threads
echo $E4_50x50
echo Elapsed time 8 threads
echo $E8_50x50
echo 

##################################### space 100x100 ############################################################
E_100x100=`awk '/counter 2100/ {print $4}' resultNotParallelFile_100x100.txt`
E1_100x100=`awk '/counter 2100/ {print $4}' resultParallelFile_1thread_100x100.txt`
E2_100x100=`awk '/counter 2100/ {print $4}' resultParallelFile_2thread_100x100.txt`
E4_100x100=`awk '/counter 2100/ {print $4}' resultParallelFile_4thread_100x100.txt`
E8_100x100=`awk '/counter 2100/ {print $4}' resultParallelFile_8thread_100x100.txt`
echo Elapsed time 100x100 serial version
echo $E_100x100
echo Elapsed time 1 thread
echo $E1_100x100
echo Elapsed time 2 threads
echo $E2_100x100
echo Elapsed time 4 threads
echo $E4_100x100
echo Elapsed time 8 threads
echo $E8_100x100
echo 

##################################### space 200x200 ############################################################
E_200x200=`awk '/counter 2100/ {print $4}' resultNotParallelFile_200x200.txt`
E1_200x200=`awk '/counter 2100/ {print $4}' resultParallelFile_1thread_200x200.txt`
E2_200x200=`awk '/counter 2100/ {print $4}' resultParallelFile_2thread_200x200.txt`
E4_200x200=`awk '/counter 2100/ {print $4}' resultParallelFile_4thread_200x200.txt`
E8_200x200=`awk '/counter 2100/ {print $4}' resultParallelFile_8thread_200x200.txt`
echo Elapsed time 200x200 serial version
echo $E_200x200
echo Elapsed time 1 thread
echo $E1_200x200
echo Elapsed time 2 threads
echo $E2_200x200
echo Elapsed time 4 threads
echo $E4_200x200
echo Elapsed time 8 threads
echo $E8_200x200
echo 


###################################################################################################################################
###################################################################################################################################

# Computation of the real speedup. R is the real speedup of the serial version (is always 1). R2 for the 2 threads version
##################################### space 50x50 ############################################################
R_50x50=1
R1_50x50=`echo "scale=6;  $E_50x50/$E1_50x50" | bc`
R2_50x50=`echo "scale=6;  $E_50x50/$E2_50x50" | bc`
R4_50x50=`echo "scale=6;  $E_50x50/$E4_50x50" | bc`
R8_50x50=`echo "scale=6;  $E_50x50/$E8_50x50" | bc`
echo Real speedup 50x50 serial version
echo $R_50x50
echo Real speedup 1 thread
echo $R1_50x50
echo Real speedup 2 threads
echo $R2_50x50
echo Real speedup 4 threads
echo $R4_50x50
echo Real speedup 8 threads
echo $R8_50x50
echo 


##################################### space 100x100 ############################################################
R_100x100=1
R1_100x100=`echo "scale=6;  $E_100x100/$E1_100x100" | bc`
R2_100x100=`echo "scale=6;  $E_100x100/$E2_100x100" | bc`
R4_100x100=`echo "scale=6;  $E_100x100/$E4_100x100" | bc`
R8_100x100=`echo "scale=6;  $E_100x100/$E8_100x100" | bc`
echo Real speedup 100x100 serial version
echo $R_100x100
echo Real speedup 1 thread
echo $R1_100x100
echo Real speedup 2 threads
echo $R2_100x100
echo Real speedup 4 threads
echo $R4_100x100
echo Real speedup 8 threads
echo $R8_100x100
echo 


##################################### space 200x200 ############################################################
R_200x200=1
R1_200x200=`echo "scale=6;  $E_200x200/$E1_200x200" | bc`
R2_200x200=`echo "scale=6;  $E_200x200/$E2_200x200" | bc`
R4_200x200=`echo "scale=6;  $E_200x200/$E4_200x200" | bc`
R8_200x200=`echo "scale=6;  $E_200x200/$E8_200x200" | bc`
echo Real speedup 200x200 serial version
echo $R_200x200
echo Real speedup 1 thread
echo $R1_200x200
echo Real speedup 2 threads
echo $R2_200x200
echo Real speedup 4 threads
echo $R4_200x200
echo Real speedup 8 threads
echo $R8_200x200
echo 

###################################################################################################################################
###################################################################################################################################


# Computation of the efficency. EFF is the efficency of the serial version (is always 100). EFF2 for the 2 threads version
# (Formule Speedup / number of cores) * 100
##################################### space 50x50 ############################################################
EFF_50x50=100
EFF1_50x50=`echo "scale=6;  ($R1_50x50/$P1)*100" | bc`
EFF2_50x50=`echo "scale=6;  ($R2_50x50/$P2)*100" | bc`
EFF4_50x50=`echo "scale=6;  ($R4_50x50/$P4)*100" | bc`
EFF8_50x50=`echo "scale=6;  ($R8_50x50/$P8)*100" | bc`
echo Efficency 50x50 serial version
echo $EFF_50x50
echo Efficency 1 thread
echo $EFF1_50x50
echo Efficency 2 threads
echo $EFF2_50x50
echo Efficency 4 threads
echo $EFF4_50x50
echo Efficency 8 threads
echo $EFF8_50x50
echo 


##################################### space 100x100 ############################################################
EFF_100x100=100
EFF1_100x100=`echo "scale=6;  ($R1_100x100/$P1)*100" | bc`
EFF2_100x100=`echo "scale=6;  ($R2_100x100/$P2)*100" | bc`
EFF4_100x100=`echo "scale=6;  ($R4_100x100/$P4)*100" | bc`
EFF8_100x100=`echo "scale=6;  ($R8_100x100/$P8)*100" | bc`
echo Efficency 100x100 serial version
echo $EFF_100x100
echo Efficency 1 thread
echo $EFF1_100x100
echo Efficency 2 threads
echo $EFF2_100x100
echo Efficency 4 threads
echo $EFF4_100x100
echo Efficency 8 threads
echo $EFF8_100x100
echo 


##################################### space 200x200 ############################################################
EFF_200x200=100
EFF1_200x200=`echo "scale=6;  ($R1_200x200/$P1)*100" | bc`
EFF2_200x200=`echo "scale=6;  ($R2_200x200/$P2)*100" | bc`
EFF4_200x200=`echo "scale=6;  ($R4_200x200/$P4)*100" | bc`
EFF8_200x200=`echo "scale=6;  ($R8_200x200/$P8)*100" | bc`
echo Efficency 200x200 serial version
echo $EFF_200x200
echo Efficency 1 thread
echo $EFF1_200x200
echo Efficency 2 threads
echo $EFF2_200x200
echo Efficency 4 threads
echo $EFF4_200x200
echo Efficency 8 threads
echo $EFF8_200x200
echo 




# Generate the matlab file for the elapsed time
echo "A=[1, 1, 1; 2, 2, 2; 4, 4, 4; 8, 8, 8];" > shalw_matlabFile.m
#echo "A=[1 1 1 1 1 1; 2 2 2 2 2 2; 4 4 4 4 4 4; 8 8 8 8 8 8];" >> shalw_matlabFile.m
echo "B=[$E_50x50, $E_100x100, $E_200x200; $E2_50x50, $E2_100x100, $E2_200x200; $E4_50x50, $E4_100x100, $E4_200x200; $E8_50x50, $E8_100x100, $E8_200x200]; " >> shalw_matlabFile.m

echo " " >> shalw_matlabFile.m
echo "p1 =  plot(A(:,1), B(:,1), '-k'); % shallow 50x50 " >> shalw_matlabFile.m
echo "  hold on" >> shalw_matlabFile.m
echo "p2 =  plot(A(:,2), B(:,2), ':k'); % shallow 100x100" >> shalw_matlabFile.m
echo "  hold on" >> shalw_matlabFile.m
echo "p3 =  plot(A(:,3), B(:,3), '--k'); % shallow 200x200" >> shalw_matlabFile.m
echo "  hold on" >> shalw_matlabFile.m
#echo "p4 =  plot(A(:,4), B(:,4), '--ko'); % shallow 100x100 Amdahl" >> shalw_matlabFile.m
#echo "  hold on" >> shalw_matlabFile.m
#echo "p5 =  plot(A(:,5), B(:,5), ':kx'); % shallow 200x200 Real" >> shalw_matlabFile.m
#echo "  hold on" >> shalw_matlabFile.m
#echo "p6 =  plot(A(:,6), B(:,6), ':ko'); % shallow 200x200 Amdahl" >> shalw_matlabFile.m
#echo "  hold on" >> shalw_matlabFile.m
echo " " >> shalw_matlabFile.m
echo "xlabel('Cores' , 'FontSize' ,22);" >> shalw_matlabFile.m
echo "ylabel('Time (seconds)' , 'FontSize', 22);" >> shalw_matlabFile.m
echo " " >> shalw_matlabFile.m
echo "legend([p1, p2, p3],'50x50','100x100','200x200',1);" >> shalw_matlabFile.m
#echo "set(h1, 'Location', 'NorthWest');" >> shalw_matlabFile.m
#echo "legend([p1,p2,p3,p4,p5,p6],'50x50 real','50x50 Amdahl','100x100 real','100x100 Amdahl','200x200 real','200x200 Amdahl',1);" >> shalw_50x50_matlabFile.m
echo " " >> shalw_matlabFile.m
echo "saveas(gcf, 'shalw_matlabFile_print', 'fig');" >> shalw_matlabFile.m
echo "saveas(gcf, 'shalw_matlabFile_print', 'eps');" >> shalw_matlabFile.m
echo "quit;" >> shalw_matlabFile.m


# Generate the matlab file for the speedup time
echo "A=[1, 1, 1; 2, 2, 2; 4, 4, 4; 8, 8, 8];" > shalw_matlabFile_speedup.m
#echo "A=[1 1 1 1 1 1; 2 2 2 2 2 2; 4 4 4 4 4 4; 8 8 8 8 8 8];" >> shalw_matlabFile_speedup.m
echo "B=[$R_50x50, $R_100x100, $R_200x200; $R2_50x50, $R2_100x100, $R2_200x200; $R4_50x50, $R4_100x100, $R4_200x200; $R8_50x50, $R8_100x100, $R8_200x200]; " >> shalw_matlabFile_speedup.m

echo " " >> shalw_matlabFile_speedup.m
echo "p1 =  plot(A(:,1), B(:,1), '-k'); % shalw 50x50 real speedup" >> shalw_matlabFile_speedup.m
echo "  hold on" >> shalw_matlabFile_speedup.m
echo "p2 =  plot(A(:,2), B(:,2), ':k'); % shalw 100x100 real speedup" >> shalw_matlabFile_speedup.m
echo "  hold on" >> shalw_matlabFile_speedup.m
echo "p3 =  plot(A(:,3), B(:,3), '--k'); % shallow 200x200 real speedup" >> shalw_matlabFile_speedup.m
echo "  hold on" >> shalw_matlabFile_speedup.m
#echo "p4 =  plot(A(:,4), B(:,4), '--ko'); % shallow 100x100 Amdahl" >> shalw_50x50_matlabFile_speedup.m
#echo "  hold on" >> shalw_50x50_matlabFile_speedup.m
#echo "p5 =  plot(A(:,5), B(:,5), ':kx'); % shallow 200x200 Real" >> shalw_50x50_matlabFile_speedup.m
#echo "  hold on" >> shalw_50x50_matlabFile_speedup.m
#echo "p6 =  plot(A(:,6), B(:,6), ':ko'); % shallow 200x200 Amdahl" >> shalw_50x50_matlabFile_speedup.m
#echo "  hold on" >> shalw_50x50_matlabFile_speedup.m
echo " " >> shalw_matlabFile_speedup.m
echo "xlabel('Cores' , 'FontSize' ,22);" >> shalw_matlabFile_speedup.m
echo "ylabel('Speedup' , 'FontSize', 22);" >> shalw_matlabFile_speedup.m
echo " " >> shalw_matlabFile_speedup.m
echo "h1 = legend([p1, p2, p3],'50x50','100x100','200x200',1);" >> shalw_matlabFile_speedup.m
echo "set(h1, 'Location', 'NorthWest');" >> shalw_matlabFile_speedup.m
#echo "legend([p1,p2,p3,p4,p5,p6],'50x50 real','50x50 Amdahl','100x100 real','100x100 Amdahl','200x200 real','200x200 Amdahl',1);" >> shalw_50x50_matlabFile_speedup.m
echo " " >> shalw_matlabFile_speedup.m
echo "saveas(gcf, 'shalw_matlabFile_speedup_print', 'fig');" >> shalw_matlabFile_speedup.m
echo "saveas(gcf, 'shalw_matlabFile_speedup_print', 'eps');" >> shalw_matlabFile_speedup.m
echo "quit;" >> shalw_matlabFile_speedup.m






# Generate the matlab file for the efficency
echo "A=[1, 1, 1; 2, 2, 2; 4, 4, 4; 8, 8, 8];" > shalw_matlabFile_efficency.m
#echo "A=[1 1 1 1 1 1; 2 2 2 2 2 2; 4 4 4 4 4 4; 8 8 8 8 8 8];" >> shalw_50x50_matlabFile_efficency.m

echo "B=[$EFF_50x50, $EFF_100x100, $EFF_200x200; $EFF2_50x50, $EFF2_100x100, $EFF2_200x200; $EFF4_50x50, $EFF4_100x100, $EFF4_200x200; $EFF8_50x50, $EFF8_100x100, $EFF8_200x200]; " >> shalw_matlabFile_efficency.m

echo " " >> shalw_matlabFile_efficency.m
echo "p1 =  plot(A(:,1), B(:,1), '-k'); % shalw 50x50 efficency" >> shalw_matlabFile_efficency.m
echo "  hold on" >> shalw_matlabFile_efficency.m
echo "p2 =  plot(A(:,2), B(:,2), ':k'); % shalw 100x100 efficency" >> shalw_matlabFile_efficency.m
echo "  hold on" >> shalw_matlabFile_efficency.m
echo "p3 =  plot(A(:,3), B(:,3), '--k'); % shallow 200x200 efficency" >> shalw_matlabFile_efficency.m
echo "  hold on" >> shalw_matlabFile_efficency.m
#echo "p4 =  plot(A(:,4), B(:,4), '--ko'); % shallow 100x100 Amdahl" >> shalw_50x50_matlabFile_efficency.m
#echo "  hold on" >> shalw_50x50_matlabFile_efficency.m
#echo "p5 =  plot(A(:,5), B(:,5), ':kx'); % shallow 200x200 Real" >> shalw_50x50_matlabFile_efficency.m
#echo "  hold on" >> shalw_50x50_matlabFile_efficency.m
#echo "p6 =  plot(A(:,6), B(:,6), ':ko'); % shallow 200x200 Amdahl" >> shalw_50x50_matlabFile_efficency.m
#echo "  hold on" >> shalw_50x50_matlabFile_efficency.m
echo " " >> shalw_matlabFile_efficency.m
echo "xlabel('Cores' , 'FontSize' ,22);" >> shalw_matlabFile_efficency.m
echo "ylabel('Efficency' , 'FontSize', 22);" >> shalw_matlabFile_efficency.m
echo " " >> shalw_matlabFile_efficency.m
echo "legend([p1, p2, p3],'50x50','100x100','200x200',1);" >> shalw_matlabFile_efficency.m
#echo "set(h1, 'Location', 'NorthWest');" >> shalw_matlabFile_efficency.m
#echo "legend([p1,p2,p3,p4,p5,p6],'50x50 real','50x50 Amdahl','100x100 real','100x100 Amdahl','200x200 real','200x200 Amdahl',1);" >> shalw_50x50_matlabFile_efficency.m
echo " " >> shalw_matlabFile_efficency.m
echo "saveas(gcf, 'shalw_matlabFile_efficency_print', 'fig');" >> shalw_matlabFile_efficency.m
echo "saveas(gcf, 'shalw_matlabFile_efficency_print', 'eps');" >> shalw_matlabFile_efficency.m
echo "quit;" >> shalw_matlabFile_efficency.m





# Lanciare matlab
#matlab
matlab -r "shalw_matlabFile"
matlab -r "shalw_matlabFile_speedup"
matlab -r "shalw_matlabFile_efficency"

showfoto shalw_matlabFile_print.eps &
showfoto shalw_matlabFile_speedup_print.eps &
showfoto shalw_matlabFile_efficency_print.eps &