# \donttest{
# ipoints needs INLA
if (bru_safe_inla()) {
# Create integration points in dimension 'myDim' and 'myDiscreteDim'
ips1 <- ipoints(rbind(c(0, 3), c(3, 8)), 17, name = "myDim")
ips2 <- ipoints(domain = c(1, 2, 4), name = "myDiscreteDim")
# Calculate the cross product
ips <- cprod(ips1, ips2)
# Plot the integration points
plot(ips$myDim, ips$myDiscreteDim, cex = 10 * ips$weight)
}
# }
Run the code above in your browser using DataLab