secrfunc (version 1.0.0)

makegkPolygoncpp: Pre-compute Detection Hazard

Description

Hazards and probabilities are computed for each parameter combination, detector and mask point. Used by secr internal functions generalsecrloglik() and secr_makegk().

Usage

makegkPolygoncpp(detectfn, dim, convex, grain, ncores, gsbval, cumk, traps, mask)

Value

List with components

H

unbounded integrated hazard from radial function

gk

probability vector for elements of cc x detector x mask array

hk

hazard vector for elements of cc x detector x mask array

Arguments

detectfn

integer detector code (see secr)

dim

integer; 1 = transect, 2 = polygon

convex

logical; if TRUE all detector polygons are convex

grain

integer grain argument of RcppParallel::parallelFor

ncores

integer number of threads

gsbval

matrix of real parameter values, one row per combination

cumk

integer vector of start indices of vertices for each detector

traps

2-column matrix; coordinates of detector polygon vertices

mask

2-column matrix; coordinates of mask points

Examples

detectfn <- 14   # hazard half-normal
dim <- 2
convex <- TRUE
grain <- 1
ncores <- 2
# for a single parameter combination
gsbval <- matrix(c(0.2,1), nrow = 1)
cumk <- c(0,5)   
# centred unit polygon 
traps <- matrix(c(0,0,1,1,0,0,1,1,0,0), ncol = 2) - 0.5
mask <- as.matrix(expand.grid(x=-2:2, y=-2:2))

makegkPolygoncpp(detectfn, dim, convex, grain, ncores, gsbval, cumk, traps, mask)