require(EBSpatCGAL)
del3 <- Delaunay(3)
insert(del3,matrix(runif(3000),ncol=3))
cat("edgesMoins\n")
print(del3$rcpp()$conflicted_and_boundary_edges( c(.5,.5,.5))->edgesMoins)
insert(del3,c(.5,.5,.5))
del3$rcpp()$incident_edges( del3$point.index) -> edgesPlus
cat("edgesPlus\n")
print(edgesPlus)
cat("Check same boundary before and after insertion!\n")
unique(rbind(edgesMoins$boundary_edges[,1:3],edgesMoins$boundary_edges[,4:6])) -> moins
unique(rbind(edgesPlus[,c(1:3)],edgesPlus[,c(4:6)])) -> plus
plus <- plus[apply(plus,1,function(v) all(v!=c(.5,.5,.5))),]
moins[order(moins[,1]),] -> moins
plus[order(plus[,1]),] -> plus
print(moins);print(plus)
print(all(plus - moins == 0))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.