Nothing
library(gpuR)
context("Ordering Methods")
# set seed
set.seed(123)
ORDER <- 4
# Base R objects
A <- matrix(rnorm(ORDER^2), nrow=ORDER)
idx <- sample(seq.int(ORDER), ORDER)
C <- A[idx,]
test_that("vclMatrix permute",
{
has_gpu_skip()
fgpuX <- vclMatrix(A, type="float")
gpuC <- permute(fgpuX, order = idx)
expect_is(gpuC, "fvclMatrix")
expect_equal(gpuC[], C, tolerance=1e-06,
info="float row permutations not equivalent")
has_double_skip()
fgpuX <- vclMatrix(A, type="double")
gpuC <- permute(fgpuX, order = idx)
expect_is(gpuC, "dvclMatrix")
expect_equal(gpuC[], C, tolerance=.Machine$double.eps^0.5,
info="double row permutations not equivalent")
})
test_that("vclBlockMatrix permute",
{
has_gpu_skip()
fgpuX <- vclMatrix(A, type="float")
blockX <- block(fgpuX, 1L, 4L, 1L, 3L)
gpuC <- permute(blockX, order = idx)
expect_is(gpuC, "fvclMatrixBlock")
expect_equal(gpuC[], C[,1:3], tolerance=1e-06,
info="float block row permutations not equivalent")
has_double_skip()
fgpuX <- vclMatrix(A, type="double")
blockX <- block(fgpuX, 1L, nrow(A), 1L, 3L)
gpuC <- permute(blockX, order = idx)
expect_is(gpuC, "dvclMatrixBlock")
expect_equal(gpuC[], C[,1:3], tolerance=.Machine$double.eps^0.5,
info="double block row permutations not equivalent")
})
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.