supercomp11/data/script.r

   1 ## Usage: launch R, and then: source("script.r")
   2
   3 ## load dataset (variable is "ds")
   4 load("data.rda")
   5
   6 ## extract data for plain algorithms, in real mode
   7 #ds.base <- subset(ds, Mode == "R" & grepl("[lt]_plain", Algo))
   8
   9 ## extract data for one experiment
  10 #xx <- subset(ds,
  11 #             Mode == "R" & Distrib == 1 & Ratio == "1:1" &
  12 #             Platform == "cluster" & Size == 16 & Topo == "line")
  13
  14
  15
  16 draw <- function(dset = ds,
  17                  draw_mode = "R",
  18                  draw_distrib = 1,
  19                  draw_ratio = "1:1",
  20                  draw_platform = "cluster",
  21                  draw_topo = "hcube") {
  22
  23   ## extract desired data
  24   dset <- subset(dset,
  25                  Mode == draw_mode & Distrib == draw_distrib &
  26                  Ratio == draw_ratio &
  27                  Platform == draw_platform & Topo == draw_topo)
  28
  29   ## reorder lines (first "plain", then "bookkeeping")
  30   dset <- rbind(subset(dset, grepl("plain", Algo)),
  31                 subset(dset, grepl("bookkeeping", Algo)))
  32   dset <- dset[order(dset$Size), ]
  33
  34   ## finally, add an "Index" attribute(, for further processing
  35   dset <- cbind(dset, Index=seq(1, length(dset$Size))) # FIXME
  36
  37   #print(dset)
  38
  39   ## save the graphical parameters
  40   p <- par(no.readonly=TRUE)
  41
  42 ### TODO
  43 ### - add some spacing between each size
  44 ### - fill in empty rows (e.g. with large size)
  45 ### - add a label telling the size
  46
  47   par(mar=c(13, 4, 4, 4) + 0.1,         # change margins
  48       mgp=c(11, 1, 0))
  49
  50   par(las=3)                           # always draw labels vertically
  51
  52   ## compute the time range: 0..max+20%
  53   time_range <- c(0, range(dset$Conv_max)[2] * 1.2)
  54
  55   ## first draw the bars
  56   barplot(dset$Conv_max, axes=FALSE, col="red",
  57           names=dset$Algo,
  58           ylim=time_range)
  59   barplot(dset$Conv_avg, axes=FALSE, col="green",
  60           add=TRUE)
  61   barplot(dset$Idle_avg, axes=FALSE, col="blue",
  62           add=TRUE)
  63
  64   ## draw the time axis on the left
  65   axis(2, labels=TRUE)
  66   mtext("Time (s)", side=2, line=2.5)
  67
  68   title(xlab="Algorithms")
  69
  70   ## compute the data range
  71   x_data_range <- c(0.6, length(dset$Algo) + 0.4)
  72   data_range <- c(0, range(dset$Data_amnt)[2] * 1.04)
  73
  74   ## draw the data amnt line graphs
  75   for (alg in unique(dset$Algo)) {
  76     print(alg);
  77     dset.da <- subset(dset, Algo == alg)
  78     ## start a new graph, but drawn on the same device
  79     par(new=TRUE)
  80     plot(dset.da$Index, dset.da$Data_amnt, axes=FALSE, xlab="", type="b",
  81          xlim=x_data_range,
  82          yaxs="i", ylim=data_range)
  83   }
  84
  85   ## draw the data amnt axis on the right
  86   axis(4, labels=TRUE)
  87   mtext("Data amount (relative)", side=4, line=2.5)
  88
  89   ## finally, set title
  90   t <- paste0(dset$Mode[1], dset$Distrib[1], " / ", dset$Ratio[1], " / ",
  91               dset$Platform[1], " / ", dset$Topo[1])
  92   title(main=t)
  93
  94   ## restore the graphical parameters
  95   par(p)
  96
  97 }
  98
  99 msg <- function(text,
 100                 wait = TRUE) {
 101   if (wait)
 102     readline("Press <Enter> to continue\n")
 103   message(text)
 104 }
 105
 106 msg("First test, with algorithms \"plain\"...", wait = FALSE);
 107 draw(subset(ds, grepl("[lt]_plain", Algo)))
 108
 109 msg("Second test, with algorithms \"bookkeeping\"...");
 110 draw(subset(ds, grepl("[lt]_bookkeeping", Algo)))
 111
 112 msg("Third test, with all algorithms...");
 113 draw()