1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
|
### =============================================================================
### rowSums2
###
### ----------------------------------------------------------------------------
### Non-exported methods
###
.DelayedMatrix_block_rowSums2 <- function(x, rows = NULL, cols = NULL,
na.rm = FALSE, ..., useNames = TRUE) {
# Check input
stopifnot(is(x, "DelayedMatrix"))
DelayedArray:::.get_ans_type(x, must.be.numeric = FALSE)
# Subset
x <- ..subset(x, rows, cols)
# Compute result
val <- rowblock_APPLY(x = x,
FUN = rowSums2,
na.rm = na.rm,
...,
useNames = useNames)
if (length(val) == 0L) {
return(numeric(nrow(x)))
}
unlist(val, recursive = FALSE, use.names = useNames)
}
### ----------------------------------------------------------------------------
### Exported methods
###
# ------------------------------------------------------------------------------
# General method
#
#' @inherit MatrixGenerics::rowSums2
#' @importMethodsFrom DelayedArray seed
#' @rdname colSums2
#' @export
#' @examples
#'
#' # NOTE: Temporarily use verbose output to demonstrate which method is
#' # which method is being used
#' options(DelayedMatrixStats.verbose = TRUE)
#' # By default, this uses a seed-aware method for a DelayedMatrix with a
#' # 'SolidRleArraySeed' seed
#' rowSums2(dm_Matrix)
#' # Alternatively, can use the block-processing strategy
#' rowSums2(dm_Matrix, force_block_processing = TRUE)
#' options(DelayedMatrixStats.verbose = FALSE)
setMethod("rowSums2", "DelayedMatrix",
function(x, rows = NULL, cols = NULL, na.rm = FALSE,
force_block_processing = FALSE, ..., useNames = TRUE) {
.smart_seed_dispatcher(x, generic = MatrixGenerics::rowSums2,
blockfun = .DelayedMatrix_block_rowSums2,
force_block_processing = force_block_processing,
rows = rows,
cols = cols,
na.rm = na.rm,
...,
useNames = useNames)
}
)
# ------------------------------------------------------------------------------
# Seed-aware methods
#
#' @importMethodsFrom Matrix rowSums
#' @rdname colSums2
#' @export
setMethod("rowSums2", "Matrix",
function(x, rows = NULL, cols = NULL, na.rm = FALSE,
..., useNames = TRUE) {
message2(class(x), get_verbose())
x <- ..subset(x, rows, cols)
val <- rowSums(x = x, na.rm = na.rm)
if (!useNames) {
val <- unname(val)
}
val
}
)
|
