bm_repeat: Mapper for repeating a mapper

In inlabru-org/inlabru: Bayesian Latent Gaussian Modelling using INLA and Extensions

Mapper for repeating a mapper

Description

Defines a repeated-space mapper that sums the contributions for each copy. The ibm_n() method returns ibm_n(mapper) * n_rep, and ibm_values() returns seq_len(ibm_n(mapper)).

Usage

bm_repeat(mapper, n_rep, interleaved = FALSE)

bru_mapper_repeat(...)

## S3 method for class 'bm_repeat'
ibm_n(mapper, ...)

## S3 method for class 'bm_repeat'
ibm_n_output(mapper, ...)

## S3 method for class 'bm_repeat'
ibm_values(mapper, ...)

## S3 method for class 'bm_repeat'
ibm_jacobian(
  mapper,
  input,
  state = NULL,
  inla_f = FALSE,
  multi = FALSE,
  ...,
  sub_lin = NULL
)

## S3 method for class 'bm_repeat'
ibm_eval(mapper, input, state, multi = FALSE, ..., sub_lin = NULL)

## S3 method for class 'bm_repeat'
ibm_linear(mapper, input, state, ...)

## S3 method for class 'bm_repeat'
ibm_invalid_output(mapper, input, state, ...)

Arguments

mapper	The mapper to be repeated.
n_rep	The number of times to repeat the mapper. If a vector, the non-interleaved repeats are combined into a single repeat mapping, and combined with interleaved repeats via a bm_sum() of mappers.
interleaved	logical; if TRUE, the repeated mapping columns are interleaved; ⁠(x1[1], x2[1], ..., x1[2], x2[2], ...)⁠. If FALSE (default), the repeated mapping columns are contiguous, ⁠(x1[1], x1[2], ..., x2[1], x2[2], ...)⁠, and the Jacobian is a cbind() of the Jacobians of the repeated mappers. If n_rep is a vector, interleaved should either be a single logical, or a vector of the same length. Each element applies to the corresponding n_rep repetition specification.
...	Arguments passed on to other methods
input	Data input for the mapper.
state	A vector of latent state values for the mapping, of length ibm_n(mapper, inla_f = FALSE)
inla_f	logical; when TRUE for ibm_n() and ibm_values(), the result must be compatible with the INLA::f(...) and corresponding INLA::inla.stack(...) constructions. For ⁠ibm_{eval,jacobian,linear}⁠, the input interpretation may be different. Implementations do not normally need to do anything different, except for mappers of the type needed for hidden multicomponent models such as "bym2", which can be handled by bm_collect.
multi	logical; If TRUE (or positive), recurse one level into sub-mappers
sub_lin	Internal, optional pre-computed sub-mapper information

Value

A bm_repeat or bm_sum object, or the original input mapper.

Methods (by generic)

• ibm_jacobian(bm_repeat): The input should take the format of the repeated submapper.

• ibm_invalid_output(bm_repeat): Passes on the input to the corresponding method.

See Also

bru_mapper, bru_mapper_generics

Other mappers: bm_aggregate(), bm_collect(), bm_const(), bm_factor(), bm_fmesher(), bm_harmonics(), bm_index(), bm_linear(), bm_logsumexp(), bm_marginal(), bm_matrix(), bm_mesh_B(), bm_multi(), bm_pipe(), bm_scale(), bm_shift(), bm_sum(), bm_taylor(), bru_get_mapper(), bru_mapper(), bru_mapper.fm_mesh_1d(), bru_mapper.fm_mesh_2d(), bru_mapper_generics

Examples

(m0 <- bm_index(3))
(m <- bm_repeat(m0, 5))
ibm_n(m)
ibm_values(m)
ibm_jacobian(m, 1:3)
ibm_eval(m, 1:3, seq_len(ibm_n(m)))

# Interleaving and grouping
(m <- bm_repeat(m0, c(2, 1, 2), c(TRUE, FALSE, FALSE)))
ibm_n(m)
ibm_values(m)
ibm_jacobian(m, 1:3)
ibm_eval(m, 1:3, seq_len(ibm_n(m)))

inlabru-org/inlabru documentation built on July 17, 2025, 2:11 a.m.