inst/tests/test_R_sim2_mark.R

In rcqls/EBSpatCGAL: EBSpat using CGAL

require(EBSpatCGAL)

# init del2
del2m <- Delaunay()

insert(del2m,x=runif(100,-350,350),y=runif(100,-350,350),m=runifDisc(100,1,2))


# init gibbs 
gdm <- SimGibbs(del2m ~ 2 + Del2(th[1]*(l<=20)+th[2]*(20<l & l<=80)*abs(v[[1]]$m-v[[2]]$m),th=c(2,4))|m ~ unifDisc(1,2)) 

run(gdm)

# init scene => TODO: atomatically created later
(scm <- Scene(del=del2m)) %<<% window2d(c(-350,350),c(-350,350),xlab="",ylab="",main="what a beautiful plot!") 
scm %<<% points(del,col=m)  %<<% lines(del,type="vor") 
scm %<<% lines(del)

plot(scm)

# delete some part
delete(del2m,inside=c(-300,-300,300,300))

plot(scm)

run(gdm)

plot(scm)