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
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
|
### ============================================================================
### colProds
###
### ----------------------------------------------------------------------------
### Non-exported methods
###
.DelayedMatrix_block_colProds <- function(x, rows = NULL, cols = NULL,
na.rm = FALSE,
method = c("direct", "expSumLog"),
..., useNames = TRUE) {
# Check input
method <- match.arg(method)
stopifnot(is(x, "DelayedMatrix"))
DelayedArray:::.get_ans_type(x, must.be.numeric = FALSE)
# Subset
x <- ..subset(x, rows, cols)
# Compute result
val <- colblock_APPLY(x = x,
FUN = colProds,
na.rm = na.rm,
method = method,
...,
useNames = useNames)
if (length(val) == 0L) {
return(numeric(ncol(x)))
}
unlist(val, recursive = FALSE, use.names = useNames)
}
### ----------------------------------------------------------------------------
### Exported methods
###
# ------------------------------------------------------------------------------
# General method
#
#' @inherit MatrixGenerics::colProds
#' @importMethodsFrom DelayedArray seed
#' @rdname colProds
#' @template common_params
#' @template lowercase_x
#' @author Peter Hickey
#' @export
setMethod("colProds", "DelayedMatrix",
function(x, rows = NULL, cols = NULL, na.rm = FALSE,
method = c("direct", "expSumLog"),
force_block_processing = FALSE, ..., useNames = TRUE) {
method <- match.arg(method)
.smart_seed_dispatcher(x, generic = MatrixGenerics::colProds,
blockfun = .DelayedMatrix_block_colProds,
force_block_processing = force_block_processing,
rows = rows,
cols = cols,
na.rm = na.rm,
#method = method, # Wait for fix on SMS's side.
...,
useNames = useNames)
}
)
# ------------------------------------------------------------------------------
# Seed-aware methods
#
#' @importMethodsFrom IRanges Views viewApply
#' @rdname colProds
#' @export
#' @template example_dm_matrix
#' @template example_dm_HDF5
#' @examples
#'
#' colProds(dm_matrix)
setMethod("colProds", "SolidRleArraySeed",
function(x, rows = NULL, cols = NULL, na.rm = FALSE,
method = c("direct", "expSumLog"), ..., useNames = TRUE) {
method <- match.arg(method)
if (method != "direct") {
stop("Only the 'direct' method is currently supported for ",
"DelayedMatrix with '", class(x), "' seed.")
}
message2(class(x), get_verbose())
irl <- get_Nindex_as_IRangesList(Nindex = list(rows, cols),
dim = dim(x))
views <- Views(subject = x@rle, start = unlist(irl))
val <- viewApply(X = views, FUN = prod, na.rm = na.rm)
if (length(irl) == 0) {
return(val)
}
n <- length(irl[[1]])
if (n == 1) {
if (useNames) {
nms <- colnames(x)
if (!is.null(cols)) {
nms <- nms[cols]
}
names(val) <- nms
}
return(val)
}
IDX <- rep(seq_along(irl), each = n)
val <- unlist(
lapply(X = split(val, IDX), FUN = prod, na.rm = na.rm))
if (useNames) {
nms <- colnames(x)
if (!is.null(cols)) {
nms <- nms[cols]
}
names(val) <- nms
}
val
}
)
|
