Homec4science
Diffusion SCITAS application benchmarks (master) Browse

R195/comparison.R
8b2b21cc8ff2master

Last Change
Raw File

comparison.R
View Options

library(ggplot2)
library(reshape2)
rougeflash="#E3051A"
rougesc="#AE0010"
ltgray="#8E938F"
dkgray="#49423A"
#mycol = "#552683"
mycol = rougeflash
col4=rougesc
col3=ltgray
col2=dkgray
#bgcol="#E2E2E3"
bgcol="#ffffff"
clusters = c('cluster1', 'cluster2', 'cluster3')
bench=c("cpmd", "gear", "qe", "miniFE")
for (b in bench){
assign(b,NULL)
}
for (fn in clusters){
path=paste(fn,"user.Rdata",sep="/")
myenv=new.env(hash = TRUE, parent = globalenv(), size = 29L)
load(path,envir=myenv)
for (b in bench){
assign(b,rbind(get(b),get(b,envir=myenv)))
}
rm(myenv)
}
cpmdplot <- ggplot(data=cpmd, aes(x=cores,y=value,group=machine,shape=machine),color=mycol) + geom_point(color=mycol) +
geom_line(aes(linetype=machine),color=mycol) + scale_y_log10(breaks=2^(0:5)) + scale_x_log10(breaks=2^(0:4)) + ggtitle("CPMD average time per iteration")
gearplot <- ggplot(gear, aes(x=cores,y=value,group=machine,shape=machine)) +
geom_point(aes(shape=machine),color=mycol) +
scale_y_log10(breaks=0.1*2^(-3:4)) + scale_x_log10(breaks=2^(0:4)) +
geom_line(aes(linetype=machine),color=mycol) + xlab('np')+ylab('Steps per time unit') + ggtitle("GEAR Steps per time unit")
feplot <- ggplot(miniFE, aes(x=cores,y=Total_Program_Time,group=machine)) +
geom_point(aes(shape=machine),color=mycol) + facet_grid(. ~ size , scales="free_y") +
scale_y_log10(breaks=10*2^(-2:10)) + scale_x_log10(breaks=2^(0:4)) +
geom_line(aes(linetype=machine),color=mycol) + xlab('nodes') + ggtitle("MiniFE total execution time")
qeplot <- ggplot(qe, aes(x=cores,y=total,group=machine,color=machine)) +
geom_point(aes(shape=machine),color=mycol) + facet_grid(. ~ bench, scales = "free_y") +
scale_y_log10(breaks=10*2^(-4:10)) + scale_x_log10(breaks=2^(0:4)) +
geom_line(aes(linetype=machine),color=mycol) + xlab('nodes')+ylab('Total Execution Time') + ggtitle("Quantum Espresso total execution time")
##
library(extrafont)
font_import() # Import all fonts
fonts() # Print list of all fonts
# Configure Theme
kobe_theme <- function() {
theme(
plot.background = element_rect(fill = bgcol, colour = bgcol),
panel.background = element_rect(fill = bgcol),
panel.background = element_rect(fill = "white"),
axis.text = element_text(colour = col2, family = "Impact"),
plot.title = element_text(colour = mycol, face = "bold", size = 18, vjust = 1, family = "Impact"),
axis.title = element_text(colour = mycol, face = "bold", size = 13, family = "Impact"),
panel.grid.major.x = element_line(colour = col2),
panel.grid.minor.x = element_blank(),
panel.grid.major.y = element_blank(),
panel.grid.minor.y = element_blank(),
strip.text = element_text(family = "Impact", colour = "white"),
strip.background = element_rect(fill = col2),
axis.ticks = element_line(colour = col2)
)
}
x_id <- rep(12:1, 3) # use this index for reordering the x ticks
kobe_theme2 <- function() {
theme(
legend.position = "bottom", legend.title = element_text(family = "Impact", colour = mycol, size = 10),
legend.background = element_rect(fill = bgcol),
legend.key = element_rect(fill = bgcol, colour = bgcol),
legend.text = element_text(family = "Impact", colour = col2, size = 10),
plot.background = element_rect(fill = bgcol, colour = bgcol),
panel.background = element_rect(fill = bgcol),
panel.background = element_rect(fill = "white"),
axis.text = element_text(colour = col2, family = "Impact"),
plot.title = element_text(colour = mycol, face = "bold", size = 18, vjust = 1, family = "Impact"),
axis.title = element_text(colour = mycol, face = "bold", size = 13, family = "Impact"),
panel.grid.major.y = element_line(colour = col2),
panel.grid.minor.y = element_blank(),
panel.grid.major.x = element_blank(),
panel.grid.minor.x = element_blank(),
strip.text = element_text(family = "Impact", colour = "white"),
strip.background = element_rect(fill = col2),
axis.ticks = element_line(colour = col2)
)
}
p0 <- cpmdplot + kobe_theme2()
p1 <- gearplot + kobe_theme2()
p2 <- feplot + kobe_theme2()
p3 <- qeplot + kobe_theme2()
require(useful)
# Generate Infographic in PDF format
library(grid)
pdf("user.pdf", width = 10, height = 20)
grid.newpage()
pushViewport(viewport(layout = grid.layout(5, 3)))
grid.rect(gp = gpar(fill = bgcol, col = bgcol))
# grid.text: y=unit(%height,"npc")
grid.text("BENCHMARKS", y = unit(1, "npc"), x = unit(0.5, "npc"), vjust = 1, hjust = .5, gp = gpar(fontfamily = "Impact", col = col3, cex = 11.2, alpha = 0.3))
grid.text("SCITAS", y = unit(0.97, "npc"), gp = gpar(fontfamily = "Impact", col = col4, cex = 6.4))
grid.text("scitas.epfl.ch", vjust = 0, y = unit(0.936, "npc"), gp = gpar(fontfamily = "Impact", col = mycol, cex = 0.8))
grid.text("Fidis mini-tender", vjust = 0, hjust = 0,
x = unit(0.01, "npc"), y = unit(0.92, "npc"), gp = gpar(fontfamily = "Impact", col = mycol, cex = 1.2))
# version intelMPI et icc.
# orange rect for "APPLICATIONS"
grid.rect(gp = gpar(fill = col2, col = col2), x = unit(0.5, "npc"), y = unit(0.85, "npc"), width = unit(1, "npc"), height = unit(0.11, "npc"))
grid.text("APPLICATIONS", y = unit(0.85, "npc"), x = unit(0.5, "npc"), vjust = .5, hjust = .5,
gp = gpar(fontfamily = "Impact", col = col3, cex = 10.9, alpha = 0.3))
# grid.text("MACHINE INFO", vjust = 0, hjust = 0, x = unit(0.01, "npc"), y = unit(0.89, "npc"), gp = gpar(fontfamily = "Impact", col = "white", cex = 1.2))
# grid.text(paste(
# "Machine ID",
# "Number of nodes",
# "Sockets per node",
# "CPU model",
# "Cores per socket",
# "Frequency",
# "Compilers and MPI",
# "Date", sep = "\n"), vjust = 0, hjust = 0, x = unit(0.01, "npc"), y = unit(0.81, "npc"), gp = gpar(fontfamily = "Impact", col = mycol, cex = 0.8))
# grid.text(paste(
# "cluster1",
# 1,
# 2,
# "Intel(R) Xeon(R) CPU E5-2680 v3",
# 12,
# "2.5 GHz",
# "icc (ICC) 15.0.2 , Intel(R) MPI Library 5.1.3",
# "April 2016", sep = "\n"), vjust = 0, hjust = 0, x = unit(0.15, "npc"), y = unit(0.81, "npc"), gp = gpar(fontfamily = "Impact", col = mycol, cex = 0.8))
# grid.text(paste(
# "eltanin",
# 1,
# 2,
# "Intel(R) Xeon(R) CPU E5-2680 v3",
# 12,
# "2.5 GHz",
# " ",
# "April 2016", sep = "\n"), vjust = 0, hjust = 0, x = unit(0.34, "npc"), y = unit(0.81, "npc"), gp = gpar(fontfamily = "Impact", col = mycol, cex = 0.8))
print(p0, vp = vplayout(5, 1:3))
print(p3, vp = vplayout(4, 1:3))
print(p1, vp = vplayout(3, 1:3))
print(p2, vp = vplayout(2, 1:3))
dev.off()
c4science · Help