Running examples in ‘RCDT-Ex.R’ failed
The error most likely occurred in:

> ### Name: plotDelaunay
> ### Title: Plot 2D Delaunay triangulation
> ### Aliases: plotDelaunay
> 
> ### ** Examples
> 
> library(RCDT)
> # random points in a square ####
> square <- rbind(
+   c(-1, 1), c(1, 1), c(1, -1), c(-1, -1)
+ )
> library(uniformly)
> set.seed(314)
> pts_in_square <- runif_in_cube(10L, d = 2L)
> pts <- rbind(square, pts_in_square)
> d <- delaunay(pts)
> opar <- par(mar = c(0, 0, 0, 0))
> plotDelaunay(
+   d, type = "n", xlab = NA, ylab = NA, axes = FALSE, asp = 1, 
+   fillcolor = "random", luminosity = "dark", lwd_borders = 3
+ )
> par(opar)
> 
> # a constrained Delaunay triangulation: pentagram ####
> # vertices
> R <- sqrt((5-sqrt(5))/10)     # outer circumradius
> r <- sqrt((25-11*sqrt(5))/10) # circumradius of the inner pentagon
> k <- pi/180 # factor to convert degrees to radians
> X <- R * vapply(0L:4L, function(i) cos(k * (90+72*i)), numeric(1L))
> Y <- R * vapply(0L:4L, function(i) sin(k * (90+72*i)), numeric(1L))
> x <- r * vapply(0L:4L, function(i) cos(k * (126+72*i)), numeric(1L))
> y <- r * vapply(0L:4L, function(i) sin(k * (126+72*i)), numeric(1L))
> vertices <- rbind(
+   c(X[1L], Y[1L]),
+   c(x[1L], y[1L]),
+   c(X[2L], Y[2L]),
+   c(x[2L], y[2L]),
+   c(X[3L], Y[3L]),
+   c(x[3L], y[3L]),
+   c(X[4L], Y[4L]),
+   c(x[4L], y[4L]),
+   c(X[5L], Y[5L]),
+   c(x[5L], y[5L])
+ )
> # constraint edge indices (= boundary)
> edges <- cbind(1L:10L, c(2L:10L, 1L))
> # constrained Delaunay triangulation
> del <- delaunay(vertices, edges)
> # plot
> opar <- par(mar = c(0, 0, 0, 0))
> plotDelaunay(
+   del, type = "n", asp = 1, fillcolor = "distinct", lwd_borders = 3,
+   xlab = NA, ylab = NA, axes = FALSE
+ )
> par(opar)
> # interactive plot with 'rgl'
> mesh <- del[["mesh"]]
> library(rgl)
> open3d(windowRect = c(100, 100, 612, 612))
null 
   1 
> shade3d(mesh, color = "red", specular = "orangered")
> wire3d(mesh, color = "black", lwd = 3, specular = "black")
> # plot only the border edges - we could find them in `del[["edges"]]` 
>   # but we use the 'rgl' function `getBoundary3d` instead
> open3d(windowRect = c(100, 100, 612, 612))
null 
   2 
> shade3d(mesh, color = "darkred", specular = "firebrick")
> shade3d(getBoundary3d(mesh), lwd = 3)
> 
> # an example where `fillcolor` is a vector of colors ####
> n <- 50L # number of sides of the outer polygon
> angles1 <- head(seq(0, 2*pi, length.out = n + 1L), -1L)
> outer_points <- cbind(cos(angles1), sin(angles1))
> m <- 5L # number of sides of the inner polygon
> angles2 <- head(seq(0, 2*pi, length.out = m + 1L), -1L)
> phi <- (1+sqrt(5))/2 # the ratio  2-phi  will yield a perfect pentagram
> inner_points <- (2-phi) * cbind(cos(angles2), sin(angles2))
> points <- rbind(outer_points, inner_points)
> # constraint edges
> indices <- 1L:n
> edges_outer <- cbind(indices, c(indices[-1L], indices[1L]))
> indices <- n + 1L:m
> edges_inner <- cbind(indices, c(indices[-1L], indices[1L]))
> edges <- rbind(edges_outer, edges_inner)
> # constrained Delaunay triangulation
> del <- delaunay(points, edges) 
> # there are 55 triangles:
> del[["mesh"]]
 mesh3d object with 55 vertices, 55 triangles.
> # we make a cyclic palette of colors:
> colors <- viridisLite::turbo(28)
> colors <- c(colors, rev(colors[-1L]))
> # plot
> opar <- par(mar = c(0, 0, 0, 0))
> plotDelaunay(
+   del, type = "n", asp = 1, lwd_borders = 3, col_borders = "black", 
+   fillcolor = colors, col_edges = "black", lwd_edges = 1.5, 
+   axes = FALSE, xlab = NA, ylab = NA
+ )
Error in FUN(newX[, i], ...) : unused argument (simplify = TRUE)
Calls: plotDelaunay -> apply
Execution halted