solution_bd <- function(i, XX, K, XKX) {
xkx <- XKX[[i]][[1]]
xk <- XKX[[i]][[2]]
# ans <- K[[i]]%*%XX[[i]]%*%xkx%*%xk
res1 <- bdwXw( XX[[i]], diag(K[[i]]), 'wX')
ans <- blockmult(blockmult(res1,xkx),xk, onmemory = T)
#..# ans <- blockmult(blockmult(bdwXw( XX[[i]], diag(K[[i]]), 'wX'),xkx),xk)
#..# ans <- mult_wX(XX[[i]], diag(K[[i]]))%*%xkx%*%xk
ans
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.