Updated script.r

diff --git a/supercomp11/data/script.r b/supercomp11/data/script.r
index ffac78310bf7bb3d7b7d7ccafd5333d94e281603..f15d382fe81c87571b7ab321855b5276fd0a48a0 100644 (file)
--- a/supercomp11/data/script.r
+++ b/supercomp11/data/script.r
@@ -1,13 +1,92 @@
+# Usage: launch R, and then: source("script.r")
+

 # load dataset (variable is "ds")
 load("data.rda")
 
 # extract data for plain algorithms, in real mode
 # load dataset (variable is "ds")
 load("data.rda")
 
 # extract data for plain algorithms, in real mode

-ds.base <- subset(ds, Mode=="R" & grepl("[lt]_plain", Algo))
+ds.base <- subset(ds, Mode == "R" & grepl("[lt]_plain", Algo))
+
+# extract data for one experiment
+#xx <- subset(ds,
+#             Mode == "R" & Distrib == 1 & Ratio == "1:1" &
+#             Platform == "cluster" & Size == 16 & Topo == "line")
+
+xx <- subset(ds,
+             Mode == "R" & Distrib == "N" & Ratio == "1:1" &
+             Platform == "cluster" & Topo == "hcube")
+
+# only keep plain algorithms with default value for k...
+xx <- subset(xx, grepl("[lt]_plain", Algo))
+
+# reorder lines (first "plain", then "bookkeeping")
+xx <- rbind(subset(xx, grepl("plain", Algo)),
+            subset(xx, grepl("bookkeeping", Algo)))
+xx <- xx[order(xx$Size), ]
+
+# finally, add an "Index" attribute(, for further processing
+xx <- cbind(xx, Index=c(1:length(xx$Size)))
+
+# save the graphical parameters
+p <- par(no.readonly=TRUE)
+
+
+### TODO
+### - add some spacing between each size
+### - fill in empty rows (e.g. with large size)
+### - add a label telling the size
+
+par(mar=c(13, 4, 4, 4) + 0.1,           # change margins
+    mgp=c(11, 1, 0))

 
 

-xx <- subset(ds.base, Ratio == "1:1" & Topo == "line" & Distrib == 1)
-plot(xx$Size, xx$Conv_max, xlab="Size", ylab="Time")
+par(las=3)                              # always draw labels vertically

 
 

-for(pl in c("cluster", "grid")) {
-   yy <- subset(xx, Platform == pl)
-   points(yy$Size, yy$Conv_max, xlab="Size", ylab="Time")
+# compute the time range: 0..max+20%
+time_range <- c(0, range(xx$Conv_max)[2] * 1.2)
+
+# first draw the bars
+barplot(xx$Conv_max, axes=FALSE, col="red",
+        names=xx$Algo,
+        ylim=time_range)
+barplot(xx$Conv_avg, axes=FALSE, col="green",
+        add=TRUE)
+barplot(xx$Idle_avg, axes=FALSE, col="blue",
+        add=TRUE)
+
+# draw the time axis on the left
+axis(2, labels=TRUE)
+mtext("Time (s)", side=2, line=2.5)
+
+title(xlab="Algorithms")
+
+# start a new graph, but drawn on the same device
+
+# compute the data range
+x_data_range <- c(0.6, length(xx$Algo) + 0.4)
+data_range <- c(0, range(xx$Data_amnt)[2] * 1.04)
+
+# draw the data amnt line graph
+#par(new=TRUE)
+#plot(xx$Data_amnt, axes=FALSE, xlab="", type="b",
+#     xlim=x_data_range,
+#     yaxs="i", ylim=data_range)
+
+for (alg in unique(xx$Algo)) {
+  print(alg);
+  yy <- subset(xx, Algo == alg)
+  par(new=TRUE)
+  plot(yy$Index, yy$Data_amnt, axes=FALSE, xlab="", type="b",
+       xlim=x_data_range,
+       yaxs="i", ylim=data_range)

 }
 }

+
+# draw the data amnt axis on the right
+axis(4, labels=TRUE)
+mtext("Data amount (relative)", side=4, line=2.5)
+
+# finally, set title
+t <- paste0(xx$Mode[1], xx$Distrib[1], " / ", xx$Ratio[1], " / ",
+            xx$Platform[1], " / ", xx$Topo[1])
+title(main=t)
+
+# restore the graphical parameters
+par(p)