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
|
### ============================================================================
### colLogSumExps
###
### ----------------------------------------------------------------------------
### Non-exported methods
###
.DelayedMatrix_block_colLogSumExps <- function(lx, rows = NULL, cols = NULL,
na.rm = FALSE,
..., useNames = TRUE) {
# Check input type
stopifnot(is(lx, "DelayedMatrix"))
DelayedArray:::.get_ans_type(lx, must.be.numeric = TRUE)
# Subset
lx <- ..subset(lx, rows, cols)
# Compute result
val <- colblock_APPLY(x = lx,
FUN = colLogSumExps,
na.rm = na.rm,
...,
useNames = useNames)
if (length(val) == 0L) {
return(numeric(ncol(lx)))
}
unlist(val, recursive = FALSE, use.names = useNames)
}
### ----------------------------------------------------------------------------
### Exported methods
###
# ------------------------------------------------------------------------------
# General method
#
#' @inherit MatrixGenerics::colLogSumExps
#' @importMethodsFrom DelayedArray seed
#' @rdname colLogSumExps
#' @template common_params
#' @template lx
#' @export
#' @author Peter Hickey
#' @examples
#' x <- DelayedArray(matrix(runif(10), ncol = 2))
#' colLogSumExps(log(x))
setMethod("colLogSumExps", "DelayedMatrix",
function(lx, rows = NULL, cols = NULL, na.rm = FALSE,
force_block_processing = FALSE, ..., useNames = TRUE) {
.smart_seed_dispatcher(lx, generic = MatrixGenerics::colLogSumExps,
blockfun = .DelayedMatrix_block_colLogSumExps,
force_block_processing = force_block_processing,
rows = rows,
cols = cols,
na.rm = na.rm,
...,
useNames = useNames)
}
)
|
