[cig-commits] commit: Add plotting utils for benchmarks

[Mercurial]

changeset:   472:38aa5c76ef3f
user:        Walter Landry <wlandry at caltech.edu>
date:        Wed Aug 22 15:03:10 2012 -0700
files:       doc/plotting_utils/sinusoid/extract_heights doc/plotting_utils/sinusoid/gplot doc/plotting_utils/thermal_diffusion/gplot
description:
Add plotting utils for benchmarks


diff -r ba199cee412e -r 38aa5c76ef3f doc/plotting_utils/sinusoid/extract_heights
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/doc/plotting_utils/sinusoid/extract_heights	Wed Aug 22 15:03:10 2012 -0700
@@ -0,0 +1,7 @@
+touch $1
+rm $1
+for i in `seq 1 10`; do
+  h=`h5dump -d /max Mesh.v-mesh.0000$i.h5 | tail -n 4 | head -n 1 | cut -d ' ' -f 6`
+  t=`h5dump -d /currentTime timeInfo.0000$i.h5 | tail -n 4 | head -n 1 | cut -d ' ' -f 5`
+  echo $t $h >> $1
+done
diff -r ba199cee412e -r 38aa5c76ef3f doc/plotting_utils/sinusoid/gplot
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/doc/plotting_utils/sinusoid/gplot	Wed Aug 22 15:03:10 2012 -0700
@@ -0,0 +1,10 @@
+set term postscript eps color
+set output "sinusoid_error.eps"
+set xlabel "Time"
+set ylabel "Error"
+set key top left
+plot [:100] "64/height_64" using ($1):(abs(($2-1)-0.01*exp(-$1/(4*pi)))) w l lw 4 title "64x128", "128/height_128" using ($1):(abs(($2-1)-0.01*exp(-$1/(4*pi)))) w l lw 4 title "128x256", "256/height_256" using ($1):(abs(($2-1)-0.01*exp(-$1/(4*pi)))) w l lw 4 title "256x512"
+
+set output "sinusoid_scaled_error.eps"
+set ylabel "Scaled Error"
+plot [:100] "64/height_64" using ($1):(1*abs(($2-1)-0.01*exp(-$1/(4*pi)))) w l lw 4 title "64x128", "128/height_128" using ($1):(2*abs(($2-1)-0.01*exp(-$1/(4*pi)))) w l lw 4 title "128x256", "256/height_256" using ($1):(4*abs(($2-1)-0.01*exp(-$1/(4*pi)))) w l lw 4 title "256x512"
diff -r ba199cee412e -r 38aa5c76ef3f doc/plotting_utils/thermal_diffusion/gplot
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/doc/plotting_utils/thermal_diffusion/gplot	Wed Aug 22 15:03:10 2012 -0700
@@ -0,0 +1,6 @@
+set term postscript eps color
+set output "thermal_diffusion_error.eps"
+set ylabel "Error"
+set xlabel "h"
+set key off
+plot "data" using (1/$1):($1*abs($2-0.908116048146)) w lp lw 4