Nothing
################################################################################
context("PROD_VEC")
set.seed(SEED)
################################################################################
# Simulating some data
N <- 73
M <- 43
x <- matrix(rnorm(N * M, mean = 100, sd = 5), N)
################################################################################
for (t in TEST.TYPES) {
X <- `if`(t == "raw", asFBMcode(x), big_copy(x, type = t))
test_that("equality with %*%", {
replicate(10, {
n <- sample(c(sample(N, size = 1), sample(0:1, size = 1)), size = 1)
m <- sample(c(sample(M, size = 1), sample(0:1, size = 1)), size = 1)
ind.row <- sample(N, size = n)
ind.col <- sample(M, size = m)
y.col <- rnorm(m)
expect_equal(big_prodVec(X, y.col, ind.row, ind.col, ncores = test_cores()),
drop(X[ind.row, ind.col, drop = FALSE] %*% y.col))
y.row <- rnorm(n)
expect_equal(big_cprodVec(X, y.row, ind.row, ind.col, ncores = test_cores()),
drop(crossprod(X[ind.row, ind.col, drop = FALSE], y.row)))
center <- rnorm(m); scale <- runif(m)
expect_equal(big_prodVec(X, y.col, ind.row, ind.col, center = center,
ncores = test_cores()),
drop(scale(X[ind.row, ind.col, drop = FALSE],
center = center, scale = FALSE) %*% y.col))
expect_equal(big_prodVec(X, y.col, ind.row, ind.col, ncores = test_cores(),
center = center, scale = scale),
drop(scale(X[ind.row, ind.col, drop = FALSE],
center = center, scale = scale) %*% y.col))
expect_equal(big_cprodVec(X, y.row, ind.row, ind.col, center = center,
ncores = test_cores()),
drop(crossprod(scale(X[ind.row, ind.col, drop = FALSE],
center = center, scale = FALSE), y.row)))
expect_equal(big_cprodVec(X, y.row, ind.row, ind.col, ncores = test_cores(),
center = center, scale = scale),
drop(crossprod(scale(X[ind.row, ind.col, drop = FALSE],
center = center, scale = scale), y.row)))
})
})
test_that("Incompatiblity between dimensions", {
ind.row <- sample(N, size = 21)
ind.col <- sample(M, size = 11)
y.col <- rnorm(21)
expect_error(big_prodVec(X, y.col, ind.row, ind.col), GET_ERROR_DIM())
y.row <- rnorm(11)
expect_error(big_cprodVec(X, y.row, ind.row, ind.col), GET_ERROR_DIM())
})
}
################################################################################
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.