scripts/temp-distance-functions/expand_matrix.R

In jtipton25/sgMRA: Multi-Resolution Gaussian Process Approximations with Stochastic Gradient Descent

expand_matrix <- function(A) {
    m <- nrow(A)
    n <- ncol(A)

    if (m > n) {
        stop("A must have more columns than rows")
    }
    # B <- Matrix(0, nrow = m, ncol = m, sparse = TRUE)
    # C <- Matrix(0, nrow = n, ncol = n, sparse = TRUE)
    # cbind(rbind(B, t(A)), rbind(A, C))
    rbind(A, Matrix(0, nrow = n-m, ncol = n, sparse = TRUE))
}