bm_repeat_indexing {inlabru}R Documentation

Re-indexing matrix for repeated mappers

Description

Helper methods for regular and interleaved state vector indexing, meant for use in bm_repeat mappers.

Usage

bm_repeat_indexing(n_map, n_rep, interleaved = FALSE)

bm_repeat_indexing_matrix(n_map, n_rep)

Arguments

n_map

The block mapper size

n_rep

The number of times the block is repeated

interleaved

logical; if TRUE, the state vector indexing is interleaved. Default is FALSE, for blockwise indexing.

Value

bm_repeat_indexing: Returns a list with a vector of offsets, offsets, and a vector of relative index values, index; block k indexing is given by offsets[k] + index.

bm_repeat_indexing_matrix: A sparseMatrix object.

Functions

Examples

(idx <- bm_repeat_indexing(3, 2, FALSE))
(idx <- bm_repeat_indexing(3, 2, TRUE))
(A <- bm_repeat_indexing_matrix(3, 2))
(x_interleaved <- 1:6)
(x_blockwise <- as.vector(A %*% x_interleaved))
(x_recovered <- as.vector(Matrix::t(A) %*% x_blockwise))


[Package inlabru version 2.13.0 Index]