#!/bin/csh

set levels=35

/bin/rm -f dmean
touch dmean
foreach flavor ( vanilla chocolate )
  /bin/rm -f $flavor.out
  echo -n " Running rk_scalar_tend_driver_$flavor for " 

  echo 7 | ../rk_scalar_tend_driver_$flavor   >& $flavor.out

  grep 'allocating big scalar arrays for' $flavor.out | sed 's/allocating big scalar arrays for *//'

  echo -n " $flavor microseconds: "
  grep 'timings:' $flavor.out | tail -1 | sed 's/timings://'

  grep dmean $flavor.out | awk '{print $4}' | paste - dmean > tmp$$
  /bin/mv tmp$$ dmean
end

skippy:
set len_a = `wc dmean`
awk '{ x = ( $1 - $2 ) / $1 ; print x >= 0 ? x : -x  } ' dmean | awk -v N=$len_a[1] '{ sum += $1} END { print 100 * sum / N " percent"  }'

/bin/rm -f tmp$$

exit

set len_a = `wc dmean`

while ( $len_a[1] > 0 ) 
   head -$levels dmean > tmp$$
   set len_a = `wc tmp$$`

   if ( $len_a[1] == $levels ) then
     awk '{ x = ( $1 - $2 ) / $1 ; print x >= 0 ? x : -x  } ' tmp$$ | awk -v N=$len_a[1] '{ sum += $1} END { print sum / N  }'
   endif

   sed "1,${levels}d" dmean > tmp$$
   /bin/mv tmp$$ dmean


end






