利用oprofile分析fortran并用gnuplot画图的shell脚本

最经利用oprofile分析fortran代码，一个小插曲就是当你的fortran程序运行时间太短的话，oprofile根本就分析不到。先 贴下我的脚本代码；

!/bin/sh

llog=$1

if [ "$#" -ne 1 ]

then

    echo "Usage:oprofile.sh interger(300,500,800)"

    exit

fi

opcontrol --no-vmlinux

opcontrol --init

opcontrol --reset

opcontrol --setup --event=CPU_CLK_UNHALTED:6000

opcontrol --start

./4_1_1

./4_1_2

./4_1_3

./4_1_4

./4_1_7

./4_1_8

opcontrol --dump

opcontrol --stop

opcontrol --shutdown

opreport -l ./4_1_1 > 1-report-$llog

opreport -l ./4_1_2 > 2-report-$llog

opreport -l ./4_1_3 > 3-report-$llog

opreport -l ./4_1_4 > 4-report-$llog

opreport -l ./4_1_7 > 7-report-$llog

opreport -l ./4_1_8 > 8-report-$llog

opannotate --source ./4_1_1 > 1-annotate-$llog

opannotate --source ./4_1_2 > 2-annotate-$llog

opannotate --source ./4_1_3 > 3-annotate-$llog

opannotate --source ./4_1_4 > 4-annotate-$llog

opannotate --source ./4_1_7 > 7-annotate-$llog

opannotate --source ./4_1_8 > 8-annotate-$llog

--------------------------------------------------------------------------------

上述脚本中的4_1_1至4_1_8是我编译后的fortran可执行文件。100，300，500，800是我里面的数组规模，我这里面的意思是分别在数组是100，300，500，800的时候，测试优化前后的代码的CPU_CLK_UNHALTED,然后把执行结果分别保存在log里面。log的名称唯一的不同就在于后面跟的llog变量的值，他是数组规模。
截取部分3-annotate-300代码：

--------------------------------------------------------------------------------

               :      Program test1 /* MAIN__ total:   9016 99.9335 */
               :        implicit none
               :        integer,parameter::scale1=300
               :       ! real:: res(scale1,scale1,scale1)=2.3
               :       ! real:: x(scale1,scale1,scale1)=3.6
其中的9016是占用的CPU周期数。

--------------------------------------------------------------------------------

当得到上述结果之后，可以画出柱状图来的到直观的比较。画图脚本如下：

--------------------------------------------------------------------------------

#!/bin/sh
lvqq_data=$1
lvqq_title=$2
lvqq_xlabel=$3
lvqq_ylabel=$4
lvqq_before=$5
lvqq_after=$6
if [ $# -ne 6 ]
then
    echo "Usage:`basename $0 ` data-file title xlabel ylabel first_histogram_name second_histogram_name "
    exit
fi
gnuplot <<CMDS
set title "$lvqq_title"
set key left top Left reverse width 0 box 3
set xlabel "$lvqq_xlabel" 0,0
set ylabel "$lvqq_ylabel" 0,0
set boxwidth 0.9 absolute
set style fill pattern 3 border -1
set style histogram clustered gap 1 title offset character 0, 0, 0
set datafile missing '-'
set style data histograms
set xtics   ("100" 0, "300" 1, "500" 2, "800" 3)
set yrange [0:200] noreverse nowriteback
plot "$lvqq_data" using 2 ti col title "$lvqq_before" , '' u 3 ti col title "$lvqq_after"
set output "test"
set terminal png
replot
quit
CMDS

因为我需要画好些图，每次都在命令行里面输入，很麻烦，所以就把他们整合进一个shell脚本，只有有数据文件，直接带上参数运行这个脚本就可以得到柱状图。