dist/ntp/scripts/plot_summary.in

   1 #! @PATH_PERL@ -w
   2 # $Id: plot_summary.in,v 1.2 2003/12/04 16:23:38 drochner Exp $
   3 #
   4 # Use Gnuplot to display data in summary files produced by summary.pl.
   5 # This script requires GNUPLOT 3.7!
   6 #
   7 # Copyright (c) 1997, 1999 by Ulrich Windl <Ulrich.Windl@rz.uni-regensburg.de>
   8 #
   9 # This program is free software; you can redistribute it and/or modify
  10 # it under the terms of the GNU General Public License as published by
  11 # the Free Software Foundation; either version 2 of the License, or
  12 # (at your option) any later version.
  13 #
  14 # This program is distributed in the hope that it will be useful, but
  15 # WITHOUT ANY WARRANTY; without even the implied warranty of
  16 # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  17 # General Public License for more details.
  18 #
  19 # You should have received a copy of the GNU General Public License
  20 # along with this program; if not, write to the Free Software
  21 # Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
  22
  23 require 5.003; # "never tested with any other version of Perl"
  24 use strict;
  25
  26 use Time::Local;
  27 use Getopt::Long;
  28
  29 # parse command line
  30 my $summary_dir = "/tmp";
  31 my $identifier = "host " . `hostname`;  # origin of these data
  32 chomp $identifier;                      # remove newline
  33 my $offset_limit = 0.128;               # limit of absolute offset
  34 my $output_file = "";                   # output file defaults to stdout
  35 my $output_file_number = 1;             # numbering of output files
  36 my $gnuplot_terminal = $ENV{DISPLAY} ? "x11" : "dumb";
  37 my $wait_after_plot = 1;
  38 my @peer_list = ();
  39
  40 my %options = ("directory|input-directory=s" => \$summary_dir,
  41                "identifier=s" => \$identifier,
  42                "offset-limit=f" => \$offset_limit,
  43                "output-file=s" => \$output_file,
  44                "peer=s@" => \@peer_list,
  45                "plot-term|gnuplot-term=s" => \$gnuplot_terminal,
  46                "wait-after-plot!" => \$wait_after_plot,
  47                );
  48
  49 if ( !GetOptions(%options) )
  50 {
  51     print STDERR "valid options for $0 are:\n";
  52     my $opt;
  53     foreach $opt (sort(keys %options)) {
  54         print STDERR "\t--$opt\t(default is ";
  55         if ( ref($options{$opt}) eq "ARRAY" ) {
  56             print STDERR join(", ",  map { "'$_'" } @{$options{$opt}});
  57         } else {
  58             print STDERR "'${$options{$opt}}'";
  59         }
  60         print STDERR ")\n";
  61     }
  62     print STDERR "\n";
  63     die;
  64 }
  65
  66 chomp $identifier;
  67 die "illegal offset-limit: $offset_limit" unless $offset_limit > 0.0;
  68 $offset_limit *= 1e6;                   # scale to microseconds
  69
  70 # return the smallest value in the given list
  71 sub min
  72 {
  73     my ($result, @rest) = @_;
  74     map { $result = $_ if ($_ < $result) } @rest;
  75     return($result);
  76 }
  77
  78 # return the largest value in the given list
  79 sub max
  80 {
  81     my ($result, @rest) = @_;
  82     map { $result = $_ if ($_ > $result) } @rest;
  83     return($result);
  84 }
  85
  86 # maybe open alternate output file
  87 sub open_output
  88 {
  89     my $file;
  90     if ($output_file) {
  91         while ( -r ($file = "$output_file$output_file_number") ) {
  92             ++$output_file_number;
  93         }
  94         open TOUCH, ">$file" and close TOUCH or die "$file: $!";
  95         print "set output \"$file\"\n";
  96     }
  97 }
  98
  99 # make Gnuplot wait
 100 sub maybe_add_pause
 101 {
 102     print "pause -1 \"Press key to continue...\"\n" if $wait_after_plot;
 103 }
 104
 105 # plot data from loop summary
 106 sub do_loop
 107 {
 108     my $fname = shift;
 109     my $line;
 110     my $out_file = "/tmp/tempdata$$";
 111     my $cmd_file = "/tmp/tempcmd$$";
 112     my ($first_day, $day_out) = ("", 0);
 113     my ($lower_bound, $upper_bound, $rms);
 114     my ($min_offs, $max_offs) = (1e9, -1e9);
 115     my ($min_rms, $max_rms) = (1e9, -1e9);
 116     open INPUT, "$fname" or die "$fname: $!";
 117     open OUTPUT, ">$out_file" or die "$out_file: $!";
 118     my @Fld;
 119     while (<INPUT>) {
 120         chop;   # strip record separator
 121         @Fld = split;
 122         if ($#Fld == 0) {
 123 # loops.19960405
 124             $_ = $Fld[0]; s/.*([12]\d{3}[01]\d[0-3]\d)$/$1/;
 125             m/(\d{4})(\d{2})(\d{2})/;
 126             $line = timegm(59, 59, 23, $3, $2 - 1, $1 - 1900, 0, 0, 0);
 127             $line = int $line / 86400;  # days relative to 1970
 128             $first_day = "$1-$2-$3 ($line)" unless $day_out;
 129             next;
 130         }
 131         if ($#Fld != 8) {
 132             warn "Illegal number of fields in file $fname, line $.";
 133             next;
 134         }
 135 # loop 216, 856106+/-874041.5, rms 117239.8, freq 67.52+/-10.335, var 4.850
 136         $_ = $Fld[1]; s/,/ /; $line .= " $_";
 137         $_ = $Fld[2]; m:(.+?)\+/-(.+),:;
 138         $lower_bound = $1 - $2;
 139         $upper_bound = $1 + $2;
 140         $line .= "$1 $lower_bound $upper_bound";
 141         $min_offs = min($min_offs, $lower_bound);
 142         $max_offs = max($max_offs, $upper_bound);
 143         $_ = $Fld[4]; s/,/ /; $rms = $_;
 144         $min_rms = min($min_rms, $rms);
 145         $max_rms = max($max_rms, $rms);
 146         $line .= " $rms";
 147         $_ = $Fld[6]; m:(.+?)\+/-(.+),:;
 148         $line .= " $1 " . ($1-$2) . " " . ($1+$2);
 149         $line .= " $Fld[8]";
 150         print OUTPUT "$line\n";
 151         $day_out = 1;
 152 # 9621 216 856106 -17935.5 1730147.5 117239.8  67.52 57.185 77.855 4.850
 153     }
 154     close INPUT;
 155     close OUTPUT or die "close failed on $out_file: $!";
 156     my $ylimit = "[";
 157     if ($min_offs < -$offset_limit) {
 158         $ylimit .= "-$offset_limit";
 159     }
 160     $ylimit .= ":";
 161     if ($max_offs > $offset_limit) {
 162         $ylimit .= "$offset_limit";
 163     }
 164     if ( $ylimit eq "[:" ) {
 165         $ylimit = "";
 166     } else {
 167         $ylimit = "[] $ylimit]";
 168     }
 169 # build command file for GNUplot
 170     open OUTPUT, "> $cmd_file" or die "$cmd_file: $!";
 171     my $oldfh = select OUTPUT;
 172     print "set term $gnuplot_terminal\n";
 173     open_output;
 174     print "set grid\n";
 175     print "set title \"Loop Summary for $identifier: " .
 176         "Daily mean values since $first_day\\n" .
 177         "(Offset limit is $offset_limit microseconds)\"\n";
 178     print "set ylabel \"[us]\"\n";
 179     print "set data style yerrorbars\n";
 180     print "set multiplot\n";
 181     print "set size 1, 0.5\n";
 182     print "set lmargin 8\n";
 183     print "set origin 0, 0.5\n";
 184     print "plot $ylimit \"$out_file\"" .
 185         " using 1:3:4:5 title \"mean offset\", ";
 186     print "\"$out_file\" using 1:(\$3-\$6/2) " .
 187         "title \"(sigma low)\" with lines, ";
 188     print "\"$out_file\" using 1:3 smooth bezier " .
 189         "title \"(Bezier med)\" with lines, ";
 190     print "\"$out_file\" using 1:(\$3+\$6/2) " .
 191         "title \"(sigma high)\" with lines\n";
 192     print "set ylabel \"[ppm]\"\n";
 193     print "set origin 0, 0.0\n";
 194     print "set title\n";
 195     print "set xlabel \"Days relative to 1970\"\n";
 196     print "plot \"$out_file\" using 1:7:8:9 title \"mean frequency\", ";
 197     print "\"$out_file\" using 1:(\$7-\$10/2) " .
 198         "title \"(sigma low)\" with lines, ";
 199     print "\"$out_file\" using 1:7 smooth bezier " .
 200         "title \"(Bezier med)\" with lines, ";
 201     print "\"$out_file\" using 1:(\$7+\$10/2) " .
 202         "title \"(sigma high)\" with lines\n";
 203     print "set nomultiplot\n";
 204     maybe_add_pause;
 205
 206     $ylimit = "[";
 207     if ($min_rms < -$offset_limit) {
 208         $ylimit .= "-$offset_limit";
 209     }
 210     $ylimit .= ":";
 211     if ($max_rms > $offset_limit) {
 212         $ylimit .= "$offset_limit";
 213     }
 214     if ( $ylimit eq "[:" ) {
 215         $ylimit ="";
 216     } else {
 217         $ylimit = "[] $ylimit]";
 218     }
 219
 220     open_output;
 221     print "set title \"Loop Summary for $identifier: " .
 222         "Standard deviation since $first_day\\n" .
 223         "(Offset limit is $offset_limit microseconds)\"\n";
 224     print "set xlabel\n";
 225     print "set ylabel \"[us]\"\n";
 226     print "set origin 0, 0.5\n";
 227     print "set data style linespoints\n";
 228     print "set multiplot\n";
 229     print "plot $ylimit \"$out_file\" using 1:6 title \"Offset\", ";
 230     print "\"$out_file\" using 1:6 smooth bezier " .
 231         "title \"(Bezier)\" with lines\n";
 232     print "set title\n";
 233     print "set origin 0, 0.0\n";
 234     print "set xlabel \"Days relative to 1970\"\n";
 235     print "set ylabel \"[ppm]\"\n";
 236     print "plot \"$out_file\" using 1:10 title \"Frequency\", ";
 237     print "\"$out_file\" using 1:10 smooth bezier " .
 238         "title \"(Bezier)\" with lines\n";
 239     print "set nomultiplot\n";
 240     maybe_add_pause;
 241
 242     close OUTPUT or die "close failed on $cmd_file: $!";
 243     select $oldfh;
 244     print `gnuplot $cmd_file`;
 245     unlink $cmd_file;
 246     unlink $out_file;
 247 }
 248
 249 # plot data form peer summary
 250 sub do_peer
 251 {
 252     my $fname = shift;
 253     my $peer = shift;
 254     my $out_file = "/tmp/tempdata$$";
 255     my $cmd_file = "/tmp/tempcmd$$";
 256     my $line;
 257     my ($first_day, $day_out) = ("", 0);
 258     open INPUT, "$fname" or die "$fname: $!";
 259     open OUTPUT, ">$out_file" or die "$out_file: $!";
 260     my @Fld;
 261     while (<INPUT>) {
 262         chop;   # strip record separator
 263         @Fld = split;
 264         if ($#Fld == 0) {
 265 # peers.19960405
 266             $_ = $Fld[0]; s/.*([12]\d{3}[01]\d[0-3]\d)$/$1/;
 267             m/(\d{4})(\d{2})(\d{2})/ or next;
 268             $line = timegm(59, 59, 23, $3, $2 - 1, $1 - 1900, 0, 0, 0);
 269             $line = int $line / 86400;  # days relative to 1970
 270             $first_day = "$1-$2-$3 ($line)" unless $day_out;
 271             next;
 272         }
 273         if ($#Fld != 7) {
 274             warn "Illegal number of fields in file $fname, line $.";
 275             next;
 276         }
 277         next if ($Fld[0] ne $peer);
 278 #       ident     cnt     mean     rms      max     delay     dist     disp
 279 # 127.127.8.1       38   30.972  189.867 1154.607    0.000  879.760  111.037
 280         $Fld[0] = $line;
 281         print OUTPUT join(' ', @Fld) . "\n";
 282 # 9969 38 30.972 189.867 1154.607 0.000 879.760 111.037
 283         $day_out = 1;
 284     }
 285     close INPUT;
 286     close OUTPUT or die "close failed on $out_file: $!";
 287     die "no data found for peer $peer" if !$day_out;
 288     open OUTPUT, "> $cmd_file" or die "$cmd_file: $!";
 289     my $oldfh = select OUTPUT;
 290     print "set term $gnuplot_terminal\n";
 291     open_output;
 292     print "set grid\n";
 293     print "set multiplot\n";
 294     print "set lmargin 8\n";
 295     print "set size 1, 0.34\n";
 296     print "set origin 0, 0.66\n";
 297     print "set title " .
 298         "\"Peer Summary for $peer on $identifier since $first_day\"\n";
 299     print "set data style linespoints\n";
 300     print "set ylabel \"[us]\"\n";
 301     print "plot \"$out_file\" using 1:3 title \"mean offset\", ";
 302     print "\"$out_file\" using 1:3 smooth bezier " .
 303         "title \"(Bezier)\" with lines, ";
 304     print "\"$out_file\" using 1:(\$3-\$7/2) " .
 305         "title \"(sigma low)\" with lines, ";
 306     print "\"$out_file\" using 1:(\$3+\$7/2) " .
 307         "title \"(sigma high)\" with lines\n";
 308     print "set title\n";
 309     print "set origin 0, 0.34\n";
 310     print "set size 1, 0.32\n";
 311     print "set ylabel\n";
 312     print "plot \"$out_file\" using 1:7 title \"dist\", ";
 313     print "\"$out_file\" using 1:7 smooth bezier " .
 314         "title \"(Bezier)\" with lines\n";
 315     print "set origin 0, 0.00\n";
 316     print "set size 1, 0.35\n";
 317     print "set xlabel \"Days relative to 1970\"\n";
 318     print "plot \"$out_file\" using 1:8 title \"disp\", ";
 319     print "\"$out_file\" using 1:8 smooth bezier " .
 320         "title \"(Bezier)\" with lines\n";
 321     print "set nomultiplot\n";
 322     maybe_add_pause;
 323
 324     select $oldfh;
 325     close OUTPUT or die "close failed on $cmd_file: $!";
 326     print `gnuplot $cmd_file`;
 327     unlink $cmd_file;
 328     unlink $out_file;
 329 }
 330
 331
 332 my $loop_summary ="$summary_dir/loop_summary";
 333 my $peer_summary ="$summary_dir/peer_summary";
 334 my $clock_summary="$summary_dir/clock_summary";
 335
 336 do_loop $loop_summary;
 337 map { do_peer $peer_summary, $_ } @peer_list;