+# Usage: launch R, and then: source("script.r")
+
# 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)
