#' @useDynLib celda _rowSumByGroup
.rowSumByGroup <- function(x, group, L) {
    group <- factor(group, levels = seq(L))
    res <- .Call("_rowSumByGroup", x, group)
    return(res)
}

#' @useDynLib celda _rowSumByGroupChange
.rowSumByGroupChange <- function(x, px, group, pgroup, L) {
    group <- factor(group, levels = seq(L))
    pgroup <- factor(pgroup, levels = seq(L))
    res <- .Call("_rowSumByGroupChange", x, px, group, pgroup)
    return(res)
}

#' @useDynLib celda _colSumByGroup
.colSumByGroup <- function(x, group, K) {
    group <- factor(group, levels = seq(K))
    res <- .Call("_colSumByGroup", x, group)
    return(res)
}

#' @useDynLib celda _colSumByGroupChange
.colSumByGroupChange <- function(x, px, group, pgroup, K) {
    group <- factor(group, levels = seq(K))
    pgroup <- factor(pgroup, levels = seq(K))
    res <- .Call("_colSumByGroupChange", x, px, group, pgroup)
    return(res)
}


#' @useDynLib celda _rowSumByGroup_numeric
.rowSumByGroupNumeric <- function(x, group, L) {
    group <- factor(group, levels = seq(L))
    res <- .Call("_rowSumByGroup_numeric", x, group)
    return(res)
}

#' @useDynLib celda _colSumByGroup_numeric
.colSumByGroupNumeric <- function(x, group, K) {
    group <- factor(group, levels = seq(K))
    res <- .Call("_colSumByGroup_numeric", x, group)
    return(res)
}

#' @useDynLib celda _perplexityG
.perplexityGLogPx <- function(x, phi, psi, group, L) {
    group <- factor(group, levels = seq(L))
    res <- .Call("_perplexityG", x, phi, psi, group)
    return(res)
}