iccsa-17-gpulab

common.R (5121B)

---

 load_log_file_as_data_frame <- function (prefix, pattern) {
 	runtime <- data.frame()
 	dirs <- sort(list.files(prefix))
 	rownum <- 1
 	run <- 1
 	for (d in dirs) {
 		files <- sort(list.files(file.path(prefix, d), pattern=pattern))
 		for (f in files) {
 			lines <- readLines(file.path(prefix, d, f))
 			t <- lines[grepl("^prfl.*=.*$", lines, perl=TRUE)]
 			routine <- gsub(
 				"^prfl\\s+([a-zA-Z0-9_]+).*=\\s*([0-9]+)us$",
 				"\\1",
 				t,
 				perl=TRUE
 			)
 			t <- gsub(
 				"^prfl\\s+([a-zA-Z0-9_]+).*=\\s*([0-9]+)us$",
 				"\\2",
 				t,
 				perl=TRUE
 			)
 			s <- lines[grepl("^Output grid size", lines, perl=TRUE)]
 			s <- gsub(
 				"Output grid size\\s*=\\s*\\([0-9]+,([0-9]+),[0-9]+\\)",
 				"\\1",
 				s,
 				perl=TRUE
 			)
 			s <- as.numeric(s)
 			# add new column if needed
 			if (!("run" %in% colnames(runtime))) {
 				runtime[, "run"] <- rep(NA, nrow(runtime))
 			}
 			if (!("t" %in% colnames(runtime))) {
 				runtime[, "t"] <- rep(NA, nrow(runtime))
 			}
 			if (!("size" %in% colnames(runtime))) {
 				runtime[, "size"] <- rep(NA, nrow(runtime))
 			}
 			if (!("routine" %in% colnames(runtime))) {
 				runtime[, "routine"] <- rep(NA, nrow(runtime))
 			}
 			# add new row if needed
 			for (i in seq(1, length(t))) {
 				if (!(rownum %in% rownames(runtime))) {
 					runtime[rownum, ] <- rep(NA, ncol(runtime))
 				}
 				runtime[rownum, "run"] <- run
 				runtime[rownum, "t"] <- as.numeric(t[[i]])
 				runtime[rownum, "size"] <- s
 				runtime[rownum, "routine"] <- routine[[i]]
 #				runtime[rownum, ] <- c(as.numeric(t[[i]]), s, routine[[i]])
 				rownum <- rownum + 1
 			}
 			#t <- sum(as.numeric(t))
 			#runtime[r, ] <- c(t, s, routine)
 			#rownum <- rownum + 1
 		}
 		run <- run + 1
 	}
 	runtime
 }
 
 arma.plot_cpu_gpu <- function () {
 	data1 <- load_log_file_as_data_frame("data", "^linear.*\\.log$")
 	data2 <- load_log_file_as_data_frame("data", "^high_amplitude_realtime.*\\.log$")
 	data1 <- aggregate(
 		data1$t,
 		by=list(run=data1$run, size=data1$size),
 		FUN=sum
 	)
 	names(data1) <- c("run", "size", "t") 
 	data2 <- aggregate(
 		data2$t,
 		by=list(run=data2$run, size=data2$size),
 		FUN=sum
 	)
 	names(data2) <- c("run", "size", "t") 
 	#print(data2)
 	data1 <- aggregate(
 		data1$t*1e-6,
 		by=list(size=data1$size),
 		FUN=mean
 	)
 	data2 <- aggregate(
 		data2$t*1e-6,
 		by=list(size=data2$size),
 		FUN=mean
 	)
 	plot.new()
 	plot.window(
 		xlim=range(data1$size),
 		ylim=range(c(data1$x, data2$x))
 	)
 	lines(data1$size, data1$x, lty="solid", type="b")
 	lines(data2$size, data2$x, lty="dashed", type="b")
 	axis(1, at=data1$size, labels=data1$size)
 	axis(2)
 	title(
 		xlab="Wavy surface size",
 		ylab="Time, s"
 	)
 	text(586, 13, "OpenMP", pos=4)
 	text(622, 2, "OpenCL", pos=4)
 	#legend("top", c("OpenMP", "OpenCL"), lty=c("solid", "dashed"))
 	#print(data1)
 	#print(data2)
 }
 
 arma.aggregate_dev <- function (data2) {
 	dev_to_host_opencl <- data2[data2$routine=="dev_to_host_copy", ] 
 	window_function <- data2[data2$routine=="window_function", ] 
 	second_function <- data2[data2$routine=="second_function", ] 
 	fft <- data2[data2$routine=="fft", ] 
 	all_opencl <- data.frame(
 		size=dev_to_host_opencl$size,
 		t_wnd=window_function$x,
 		t_sec=second_function$x,
 		t_fft=fft$x,
 		t_copy=dev_to_host_opencl$x,
 		row.names="size"
 	)
 	all_opencl
 }
 
 arma.plot_aggregate_dev <- function (all_opencl, ...) {
 	args <- list(...)
 	if (args$palette == 'gray') {
 		bar_colours <- gray.colors(ncol(all_opencl))
 	} else if (args$palette == 'spbu') {
 		palette <- colorRampPalette(c("#9F2D20", "#A8ADB4"))
 		bar_colours <- palette(ncol(all_opencl))
 	} else {
 		stop(paste("unknown palette: ", args$palette, sep=" "))
 	}
 	barplot(
 		t(as.matrix(all_opencl)),
 		ylim=args$ylim,
 		axes=FALSE,
 		col=bar_colours
 	)
 	axis(2, at=seq(0, args$ylim[[2]], length.out=6))
 	title(
 		main=args$title,
 		xlab="Wavy surface size",
 		ylab="Time, s"
 	)
 	legend(
 		"top",
 		args$legend,
 		fill=bar_colours
 	)
 }
 
 arma.plot_performance_breakdown <- function (palette="gray") {
 	data1 <- load_log_file_as_data_frame("data", "^linear.*\\.log$")
 	data2 <- load_log_file_as_data_frame("data", "^high_amplitude_realtime.*\\.log$")
 	data1 <- aggregate(
 		data1$t,
 		by=list(run=data1$run, size=data1$size, routine=data1$routine),
 		FUN=sum
 	)
 	names(data1) <- c("run", "size", "routine", "t") 
 	data2 <- aggregate(
 		data2$t,
 		by=list(run=data2$run, size=data2$size, routine=data2$routine),
 		FUN=sum
 	)
 	names(data2) <- c("run", "size", "routine", "t") 
 	#print(data2)
 	data1 <- aggregate(
 		data1$t*1e-6,
 		by=list(size=data1$size, routine=data1$routine),
 		FUN=mean
 	)
 	data2 <- aggregate(
 		data2$t*1e-6,
 		by=list(size=data2$size, routine=data2$routine),
 		FUN=mean
 	)
 	all_openmp <- arma.aggregate_dev(data1)
 	all_openmp$t_copy <- NULL
 	all_opencl <- arma.aggregate_dev(data2)
 	print(all_openmp)
 	print(all_opencl)
 	par(mfrow=c(1,2))
 	arma.plot_aggregate_dev(
 		all_opencl,
 		ylim=c(0,1),
 		title="OpenCL",
 		legend=c(expression("g"[1]), expression("g"[2]), "FFT", "Copy"),
 		palette=palette
 	)
 	arma.plot_aggregate_dev(
 		all_openmp,
 		ylim=c(0,20),
 		title="OpenMP",
 		legend=c(expression("g"[1]), expression("g"[2]), "FFT"),
 		palette=palette
 	)
 }