<%@include file="includes/setup.md.rsp"%>

<%@string colname="colCumprods"%>
<%@string rowname="rowCumprods"%>
<%@string fcnname="colRowCumprods_subset"%>
<%@meta title="${colname}() and ${rowname}() benchmarks  on subsetted computation"%>
<%@meta author="Dongcan Jiang"%>
<%@meta date="2015-06-06"%>

<%@include file="${header}"%>


# <%@meta name="title"%>

This report benchmark the performance of <%=colname%>() and <%=rowname%>() on subsetted computation.


<% for (mode in c("integer", "double")) { %>

## Data type "<%=mode%>"

### Data
```r
<%=withCapture({
<%@include file="R/random-matrices.R"%>
data <- rmatrices(mode = mode, range = c(-1, 1))
})%>
```

### Results

<% for (dataLabel in names(data)) { %>
<% mprintf("%s: %s\n", mode, dataLabel)  %>
#### <%=dataLabel%> <%=mode%> matrix

```r
<%=withCapture({
X <- data[[.dataLabel.]]
rows <- sample.int(nrow(X), size = nrow(X)*0.7)
cols <- sample.int(ncol(X), size = ncol(X)*0.7)
X_S <- X[rows, cols]
gc()

colStats <- microbenchmark(
 "colCumprods_X_S"           = colCumprods(X_S),
 "colCumprods(X, rows, cols)"  = colCumprods(X, rows = rows, cols = cols),
 "colCumprods(X[rows, cols])" = colCumprods(X[rows, cols]),
 unit = "ms"
)

X <- t(X)
X_S <- t(X_S)
gc()

rowStats <- microbenchmark(
 "rowCumprods_X_S"           = rowCumprods(X_S),
 "rowCumprods(X, cols, rows)"  = rowCumprods(X, rows = cols, cols = rows),
 "rowCumprods(X[cols, rows])" = rowCumprods(X[cols, rows]),
 unit = "ms"
)
})%>
```

<% crBenchmarkResults(colStats, rowStats, tags=c(mode, dataLabel)) %>

<% } # for (dataLabel ...) %>

<% } # for (mode ...) %>


<%@include file="${footer}"%>


<%---------------------------------------------------------------------------
HISTORY:
2015-06-06
o Created.
---------------------------------------------------------------------------%>
