University of Iowa

College of Liberal Arts and Sciences

• STAT:7400
• Syllabus
• Outline
• Notes
• Homework
• Resources

HW7 Notes

AddOne OMP Timings

local({
    tmpdir <- tempdir()
    year <- 2023
    tarball <- "AddOne_4.0-2.tar.gz"
    path <- sprintf("classes/STAT7400-%s/examples/%s", year, tarball)
    url <- sprintf("http://www.stat.uiowa.edu/~luke/%s", path)
    download.file(url, file.path(tmpdir, tarball))
    dir.create(file.path(tmpdir, "lib"))
    .libPaths(c(file.path(tmpdir, "lib"), .libPaths()))
    install.packages(file.path(tmpdir, tarball), repos = NULL)
})

library(AddOne)
    
R <- 100000
K <- 3
threads <- c(1, 2, 4)
if (isTRUE(parallel::detectCores() > 4))
    threads <- c(1, 2, 4, 8)
lens <- c(8, 32, 64, 128, 256, 512, 1024, 2048)

n <- rep(rep(lens, length(threads)), K)
P <- rep(rep(threads, each = length(lens)), K)

tm <- function(n, P) {
    x <- rnorm(n)
    system.time(for (i in 1 : R) p.AddOne(x, P = P))[3]
}

system.time(time <- sapply(seq_along(n), function(i) tm(n[i], P[i])))

d <- data.frame(time, n, P = as.factor(P))

library(ggplot2)
ggplot(d, aes(n, time, color = P)) + geom_point() + geom_smooth()

This work is licensed under the CC BY-NC 3.0 Creative Commons License.