gistfile1.txt

# Script to the video: https://youtu.be/erlWsquoHlM

set.seed(1234)

x <- abs(rnorm(1000))


in_base_R <- `^`(x, 2)

# code from https://helloacm.com/exponentiation-by-squaring/
fast_power <- function(base, power){
  result <- 1
  while(power > 0){
    if(power %% 2 == 1){
      result <- result * base
      power <- power - 1
    } else {
      base <- base * base
      power <- power / 2
    }
  }
  result
}


benchmark <- microbenchmark::microbenchmark(
                  `^`(x, 2),
                  fast_power(x, 2)
                )

benchmark

boxplot(benchmark)